Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,470,951 A→C 8.2% intergenic (‑38/+35) ridA ← / ← pyrI enamine/imine deaminase, reaction intermediate detoxification/aspartate carbamoyltransferase, regulatory subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,470,9510AC 8.2% 177.5 / 6.1 61intergenic (‑38/+35)ridA/pyrIenamine/imine deaminase, reaction intermediate detoxification/aspartate carbamoyltransferase, regulatory subunit
Reads supporting (aligned to +/- strand):  ref base A (29/26);  new base C (3/2);  total (32/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.06e-01

TACGGTTGCGAAGTCGTTCAGATCTTTTACAAACACGGTAGTTTTAACGATGTCGCCCACTTTCAGGCCAGCGGCTTCGACGATCGCTTTTACGTTATCCAGCGACTGACGTGCCTGTGCAGCGACGTCTGCCGGTACTTCGCCCGTTTTCGGATTTACCGGGATCTGACCGGAGGTGATGATCATATTGCCCAGATCAACGCCCTGTACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTGGCGCGTTTTCGCACGGCAAAGCTGGATGAAACCGGTTCGGCATGGCTGATACAGTTGCTGTTCGGGCAGACCAGCACATTGTCGATGCGCTCCGGCAGACTTGGGCGCGATTTACCCACCACTTCATAGTTGTCGATACGGTTAACCGTGGCTTGCGGCGCATACAATGCCAGTTGATCTACTTGATCTTCACTC  >  NC_000913/4470653‑4471250
                                                                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                                                                           
tACGGTTGCGAAGTCGTTCAGATCTTTTACAAACACGGTAGTTTTAACGATGTCGCCCACTTTCAGGCCAGCGGCTTCGACGATCGCTTTTACGTTATCCAGCGACTGACGTGCCTGTGCAGCGACGTCTGCCGGTACTTCGCCCGTTTTCGGATTTACCGGGATCTGACCGGAGGTGATGATCATATTGCCCAGATCAACGCCCTGTACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCAtt                                                                                                                                                                                                                                                                                                           <  2:231396/301‑1 (MQ=255)
  cGGTTGCGAATTCGTTCAGATCTTTTACAAACACGGTAGTTTTAACGATGTCGCCCACTTTCAGGCCAGCGGCTTCGACGATCGCTTTTACGTTATCCAGCGACTGACGTGCCTGTGCAGCGACGTCTGCCGGTACTTCGCCCGTTTTCGGATTTACCGGGATCTGACCGGAGGTGATGATCATATTGCCCAGATCAACGCCCTGTACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATtat                                                                                                                                                                                                                                                                                                         <  2:558118/301‑1 (MQ=255)
        cGAAGTCGTTCAGATCTTTTACAAACACGGTAGTTTTAACGATGTCGCCCACTTTCAGGCCAGCGGCTTCGACGATCGCTTTTACGTTATCCAGCGACTGACGTGCCTGTGCAGCGACGTCTGCCGGTACTTCGCCCGTTTTCGGATTTACCGGGATCTGACCGGAGGTGATGATCATATTGCCCAGATCAACGCCCTGTACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGACTGGCTCagatcggaaga                                                                                                                                                                                                                                                                                                   >  2:49517/1‑291 (MQ=255)
          aaGTCGTTCAGATCTTTTACAAACACGGTAGTTTTAACGATGTCGCCCACTTTCAGGCCAGCGGCTTCGACGATCGCTTTTACGTTATCCAGCGACTGACGTGCCTGTGCAGCGACGTCTGCCGGTACTTCGCCCGTTTTCGGATTTACCGGGATCTGACCGGAGGTGATGATCATATTGCCCAGATCAACGCCCTGTACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGcc                                                                                                                                                                                                                                                                                                 <  2:84658/301‑1 (MQ=255)
                ttccctacacgtcgctcttccgatcttttAACGATGTCGCCCACTTTCAGGCCAGCGGCTTCGACGATCGCTTTTACGTTATCCAGCGACTGACGTGCCTGTGCAGCGACGTCTGCCGGTACTTCGCCCGTTTTCGGATTTACCGGGATCTGACCGGAGGTGATGATCATATTGCCCAGATCAACGCCCTGTACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCCCTATAATGAGCCAGACtt                                                                                                                                                                                                                                                                                           <  1:453100/276‑1 (MQ=255)
                   aGATCTTTTACAAACACGGTAGTTTTAACGATGTCGCCCACTTTCAGGCCAGCGGCTTCGACGATCGCTTTTACGTTATCCAGCGACTGACGTGCCTGTGCAGCGACGTCTGCCGGTACTTCGCCCGTTTTCGGATTTACCGGGATCTGACCGGAGGTGATGATCATATTGCCCAGATCAACGCCCTGTACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGAATTtta                                                                                                                                                                                                                                                                                        >  2:424358/1‑301 (MQ=255)
                                        gTTTTAACGATGTCGCCCACTTTCAGGCCAGCGGCTTCGACGATCGCTTTTACGTTATCCAGCGACTGACGTGCCTGTGCAGCGACGTCTGCCGGTACTTCGCCCGTTTTCGGATTTACCGGGATCTGACCGGAGGTGATGATCATATTGCCCAGATCAACGCCCTGTACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTgg                                                                                                                                                                                                                                                                   >  2:383782/1‑301 (MQ=255)
                                          tttAACGATGTCGCCCACTTTCAGGCCAGCGGCTTCGACGATCGCTTTTACGTTATCCAGCGACTGACGTGCCTGTGCAGCGACGTCTGCCGGTACTTCGCCCGTTTTCGGATTTACCGGGATCTGACCGGAGGTGATGATCATATTGCCCAGATCAACGCCCTGTACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCCCTATAATGAGCCAGACttagatcggaagagcacacgtctgaacc                                                                                                                                                                                                                                                                 >  2:453664/1‑275 (MQ=255)
                                               gtaTGTCGCCCACTGTCAGTCCTGCGGTTTCGACGTTGGCTTTTACGTTGTCCAGCGATTGACGTGCCTGTGCAGCGACGTTTGCCGGTACTTCGCCCGTTTTCGGATTTACCGGGATCTGACCGGAGGTGATGATCATATTGCCCAGATCAACGCCCTGTACGTAAGTACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGcac                                                                                                                                                                                                                                                            <  1:268261/299‑1 (MQ=255)
                                                    tCGCCCACTTTCAGGCCAGCGGCTTCGACGATCGCTTTTACGTTATCCAGCGACTGACGTGCCTGTGCAGCGACGTCTGCCGGTACTTCGCCCGTTTTCGGATTTACCGGGATCTGACCGGAGGTGATGATCATATTGCCCAGATCAACGCCCTGTACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCacaacaa                                                                                                                                                                                                                                                       >  2:560476/1‑300 (MQ=255)
                                                     cGCCCACTTTCAGGCCAGCGGCTTCGACGATCGCTTTTACGTTATCCAGCGACTGACGTGCCTGTGCAGCGACGTCTGCCGGTACTTCGCCCGTTTTCGGATTTACCGGGATCTGACCGGAGGTGATGATCATATTGCCCAGATCAACGCCCTGTACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACAtt                                                                                                                                                                                                                                                      >  2:119304/1‑301 (MQ=255)
                                                      gCCCACTTTCAGGCCAGCGGCTTCGACGATCGCTTTTACGTTATCCAGCGACTGACGTGCCTGTGCAGCGACGTCTGCCGGTACTTCGCCCGTTTTCGGATTTACCGGGATCTGACCGGAGGTGATGATCATATTGCCCAGATCAACGCCCTGTACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTa                                                                                                                                                                                                                                                     >  2:628013/1‑301 (MQ=255)
                                                         cACTTTTAGTCCAGCGGCTTCGACGATCGCTTTTACGTTATCCAGCGACTGACGTGCCTGTGCAGCGACGTCTGCCGGTACTTCGCCCGTTTTCGGATTTACCGGGATCTGACCGGAGGTGATGATCATATTGCCCAGATCAACGCCCTGTACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATgg                                                                                                                                                                                                                                                  <  2:296007/301‑1 (MQ=255)
                                                                      gCGGCTTCGACGATCGCTTTTACGTTATCCAGCGACTGACGTGCCTGTGCAGCGACGTCTGCCGGTACTTCGCCCGTTTTCGGATTTACCGGGATCTGACCGGAGGTGATGATCATATTGCCCAGATCAACGCCCTGTACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAACCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCAACACATTATGGAAAAAATCttttt                                                                                                                                                                                                                                     >  2:378243/1‑301 (MQ=255)
                                                                         tcTTCGACGATCGCTTTTACGTTATCCAGCGACTGACGTGCCTGTGCAGCGACGTCTGCCGGTACTTCGCCCGTTTTCGGATTTACCGGGATCTGACCGGAGGTGATGATCATATTGCCCAGATCAACGCCCTGTACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTcac                                                                                                                                                                                                                                  <  2:190636/300‑1 (MQ=255)
                                                                                cgaTCGCTTTTACGTTATCCAGCGACTGACGTGCCTGTGCAGCGACGTCTGCCGGTACTTCGCCCGTTTTCGGATTTACCGGGATCTGACCGGAGGTGATGATCATATTGCCCAGATCAACGCCCTGTACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACCCTTTTTCACAtttttt                                                                                                                                                                                                                           >  1:5659/1‑301 (MQ=255)
                                                                                 taTCGCTTTTACGTTATCCAGCTACTGACGTGCCTGTGCAGCGACGTCTGCCTGTACTTCGCCCGTTTTCGGATTTACCGGGATCTGACCGGAGGTGATGATCATATTACCCAGATCAACGCCCTGTACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTg                                                                                                                                                                                                                          <  2:62907/300‑1 (MQ=255)
                                                                                       ttttTCGTTGTCCAGCGACTGACGTGCCTGGGCGGCGACGTCTGCCGGTACTTCGCCCGTGTTCGGATTTACCGGGATCTGACCGGAGGTGATGATCATATTGCCCAGATCAACGCCCTGTACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTg                                                                                                                                                                                                                    <  1:478049/301‑1 (MQ=255)
                                                                                                  ccAGCGACTGACGTGCCTGTGCAGCGACGTCTGCCGGTACTTCGCCCGTTTTCGGATTTACCGGGATCTGACCGGAGGTGATGATCATATTGCCCAGATCAACGCCCTGTACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATc                                                                                                                                                                                                         <  2:5659/301‑1 (MQ=255)
                                                                                                           tcgtgtgctcttccgatctgTCTGCCGGTACTTCGCCCGTTTTCGGATTTACCGGGATCTGACCGGAGGTGATGATCATATTGCCCAGATCAACGCCCTGTACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTGGcgcg                                                                                                                                                                                                <  2:391867/282‑1 (MQ=255)
                                                                                                              gTGCTTGTGCAGCGTCGTCTGCCGGTACTTCTCCCGTTTTCGGATTTACCTGGATCTGACCGGAGGTGATGATCATATTGCCCAGATCAACGCCCTGTACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTGGCGCGttt                                                                                                                                                                                             <  1:423795/301‑1 (MQ=255)
                                                                                                                gCCTGTGCAGCGACGTCTGCCGGTACTTCGCCCGTTTTCGGATTTACCGGGATCTGACCGGAGGTGATGATCATATTGCCCAGATCAACGCCCTGTACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTGGCGAGttttt                                                                                                                                                                                           >  2:483065/1‑300 (MQ=255)
                                                                                                                 tctttccctacacgactctcttccgatctcccGTTGTCGGATTTACCGGGATCTGACCGGAGGTGATGATCATATTGCCCAGATCAACGCCCTGTACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTGGCGCGTTTTCg                                                                                                                                                                                          <  1:304624/272‑1 (MQ=255)
                                                                                                                       cAGCGACGTCTGCCGGTACTTCGCCCGTTTTCGGATTTACCGGGATCTGACCGGAGGTGATGATCATATTGCCCAGATCAACGCCCTGTACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTGGCGCGTTTTCGCACGGc                                                                                                                                                                                    >  2:179934/1‑301 (MQ=255)
                                                                                                                              gTCTGCCGGTACTTCGCCCGTTTTCGGATTTACCGGGATCTGTCCGGAGGTGATGATCATATTGCCCAGATCAACGCCCTGTACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGACAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCCATGGcgcgagatcggaagagcgtcgtg                                                                                                                                                                             >  1:391353/1‑282 (MQ=255)
                                                                                                                               tCTGCCGGTACTTCGCCCGTTTTCGGATTTACCGGGATCTGACCGGAGGTGATGATCATATTGCCCAGATCAACGCCCTGTACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTGGCGCGTTTTCGCACGGCAAAGCTgg                                                                                                                                                                            >  2:211162/1‑301 (MQ=255)
                                                                                                                                   ccGGTATTTCTGTCGTTTTCGGATTTACCGGGTTCTGACCGGAGGTGATGATCATATTGCCCAGATCAACGCCCTGTACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCGCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTGGCGCGTTTTCGCACGGCAAAGCTGGATGa                                                                                                                                                                        <  1:482464/301‑1 (MQ=255)
                                                                                                                                      gTACTTTGCCCGTTTTCGGATTTACCGGGATCTGACCGGAGGTGATGATCATATTGCCCAGATCAACGCCCTGTACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTGGCGCGTTTTCGCACGGCAAAGCTGGATGAAAc                                                                                                                                                                     <  1:119180/301‑1 (MQ=255)
                                                                                                                                             gCCCGTTTTCGGATTTACCGGGATCTGACCGGAGGTGATGATCATATTGCCCAGATCAACGCCCTGTACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTGGCGCGTTTTCGCACGGCAAAGCTGGATGAAACCGGTTCg                                                                                                                                                              >  2:488715/1‑301 (MQ=255)
                                                                                                                                              cccGTTTTCGGATTTACCGGGATCTGACCGGAGGTGATGATCATATTGCCCAGATCAACGCCCTGTACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTGGCGCGTTTTCgagatcggaagagcacacgtctgaactcca                                                                                                                                                             >  2:305028/1‑272 (MQ=255)
                                                                                                                                              cccGTTTTCGGATTTACCGGGATCTGACCGGAGGTGATGATCATATTGCCCAGATCAACGCCCTGTACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTGGCGCGTTTTCGCACGGCAAAGCTGGATGAAACCGGTTCgg                                                                                                                                                             >  1:85623/1‑301 (MQ=255)
                                                                                                                                                                 ttaTCTGATCGGAGGTGATTATCATATTTCCCAGTTCAACGCCCTGTACGTAAGTACCGATAGCTGCCGGTGCATTTTCCGTCGCGATTGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTGGCGCGTTTTCGCACGGCAAAGCTGGATGAAACCGGTTCGGCATGGCTGATACAGTTGCt                                                                                                                                          <  2:123112/299‑1 (MQ=255)
                                                                                                                                                                     cttaCCGGTGGTGTTGATCATATTGCCCAGATTAACGCCCTGTACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGGCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATAGCATTGGCGCGTTTTCGCACGGCAAAGCTGGATGAAACCGGTTCGGCATGGCTGATACAGTTGCTGTTc                                                                                                                                      <  1:511203/298‑1 (MQ=255)
                                                                                                                                                                          ggggggTGATTATCATATTGCCCAGATCAACGCCCTGTACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTGTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTGGCGCGTTTTCGCACGGCAAAGCTGGATGAAACCGGTTCGGCATGGCTGATACAGTTGCTGTTCGGGCa                                                                                                                                 <  1:383286/297‑1 (MQ=255)
                                                                                                                                                                           ggAGGTGATGATCATATTGCCCAGATCAACGCCCTGTACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGCTAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTGGCGCGTTTTCGCACGGCAAAGCTGGATGAAACCGGTTCGGCATGGCTGATACAGTTACTGTTCGGGCAg                                                                                                                                >  1:539587/1‑301 (MQ=255)
                                                                                                                                                                                   tgatCATATTGCCCAGATCAACGCCCTGTACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTGGCGCGTTTTCGCACGGCAAAGCTGGATGAAACCGGTTCGGCATGGCTGATACAGTTGCTGTTCGGGAAGACCAgaaa                                                                                                                        >  2:396916/1‑298 (MQ=255)
                                                                                                                                                                                     ttCATATTGCCCAGATCAACGCCCTGTACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTGGCGCGTTTTCGCACGGCAAAGCTGGATGAAACCGGTTCGGCATGGCTGATACAGTTGCTGTTCGGGCAGACCAGCACAt                                                                                                                      <  2:85747/300‑1 (MQ=255)
                                                                                                                                                                                     atCATATTGCCCAGATCAACGCCCTGTACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTGGCGCGTTTTCGCACGGCAAAGCTGGATGAAACCGGTTCGGCATGGCTGATAAAGTTGCTGTTCGGGCAGACCAGcacaa                                                                                                                      >  2:237278/1‑300 (MQ=255)
                                                                                                                                                                                        atatTGCCCAGATCAACGCCCTGTACGTAAGGACCGATAGCTGCCGGTGCATTTACCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTGGCGCGTTTTCGCACGGCAAAGCTGGATGAAACCGGTTCGGCATGGCTGATACAGTTGCTGTTCGGGCAGACCAGCACCTTGt                                                                                                                   >  1:90814/1‑301 (MQ=255)
                                                                                                                                                                                         tatTGCCCAGATCAACGCCCTGTACGTAAGGACCGATAGCTGCCGGTGCATTTACCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTGGCGCGTTTTCGCACGGCAAAGCTGGATGAAACCGGTTCGGCATGGCTGATACAGTTGCTGTTCGGGCAGACCAGCACATTGTc                                                                                                                  <  2:90938/301‑1 (MQ=255)
                                                                                                                                                                                                     taaCGCCCTTTACGTAAGGACCGATAGCTGCCGTTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGTGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTGGCGCGTTTTCGCACGGCAAAGCTGGATGAAACCGGTTCGGCATGGCTGATACAGTTGCTGTTCGGGCAGACCAGCACATTGTCGATGCGCTCCgg                                                                                                      <  2:273215/300‑1 (MQ=255)
                                                                                                                                                                                                     cAACGCCCTGTACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTGGCGCGTTTTCGCACGGCAAAGCTGGATGAAACCGGTTCGGCATGGCTGATACAGTTGCTGTTCGGGCAGACCAGCACATTGTCGATGCGCTCCgg                                                                                                      >  2:175869/1‑301 (MQ=255)
                                                                                                                                                                                                     aaaCGCCCTGTACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTGGCGCGTTTTCGCACGGCAAAGCTGGATGAAACCGGTTCGGCATGGCTGATACAGTTGCTGTTCGGGCAGACCAGCACATTGTCGATGCGCTCCgg                                                                                                      <  1:210919/300‑1 (MQ=255)
                                                                                                                                                                                                                                     ggTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTGGCGCGTTTTCGCACGGCAAAGCTGGATGAAACCGGTTCGGCATGGCTGATACAGTTGCTGTTCGGGCAGACCAGCACATTGTCGATGCGCTCCGGCAGACTTGGGCGCGATTTACCCACCACTTCaa                                                                      >  2:533621/1‑300 (MQ=255)
                                                                                                                                                                                                                                         cATTTTCCGTCGCGATAGTTTTGCTCATGATTTCCCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTGGCGCGTTTTCGCACGGCAAAGCTGGATGAAACCGGTTCGGCATGGCTGATACAGTTGCTGTTCGGGCAGACCAGCACATTGTCGATGCGCTCCGGCAGACTTGGGCGCGATTTACCCACCACTTCATAGtt                                                                  >  2:299076/1‑301 (MQ=255)
                                                                                                                                                                                                                                            tttCCGTCTCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTGGCGCGTTTTCGCACGGCAAAGCTGGATGAAACCGGTTCGGCATGGCTGATACAGTTGCTGTTCGGGCAGACCAGCACATTGTCGATGCGCTCCGGCAGACTTGGGCGCGATTTACCCACCACTTCATAGTTGTc                                                               <  1:532948/301‑1 (MQ=255)
                                                                                                                                                                                                                                                   cgcgATAGTTTTGCTCAAGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTGGCGCGTTTTCGCACGGCAAAGCTGGATGAAACCGGTTCGGCATGGCTGATACAGTTGCTGTTCGGGCAGACCAGCACATTGTCGATGCGCTCCGGCAAACTTGGGCGCGATTTACCCAACACTTCATAGTTGTCGATACgg                                                        >  1:457338/1‑301 (MQ=255)
                                                                                                                                                                                                                                                     ggttAGTTTTTCTCATGATTTCTCCGTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTAGTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTGGCGCGTTTTCGCACGGCAAAGCTGGATGAAACCGGTTCGGCATGGCTGATACAGTTGCTGTTCGGGCAGACCAGCACATTGTCGATGCGCTCCGGCAGACTTGGGCGCGATTTACCCACCACTTCATAGTTGTCGATACGGtt                                                      <  1:559762/298‑1 (MQ=255)
                                                                                                                                                                                                                                                                cTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACATTATTTGCATTTGAGCGCGATATCATTGGCGCGTTTTCGCACGGCAAAGCTGGATGAAACCGGTTCGGCATGGCTGATACAGTTGCTGTTCGGGCAGACCAGCACATTGTCGATGCGCTCCGGCAGACTTGGGCGCGATTTACCCACCACTTCATAGTTGTCGATACGGTTAACCGTGGCtt                                           >  2:275039/1‑301 (MQ=255)
                                                                                                                                                                                                                                                                cTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTGGCGCGTTTTCGCACGGCAAAGCTGGATGAAACCGGTTCGGCATGGCTGATACAGTTGCTGTTCGGGCAGACCAGCACATTGTCGATGCGCTCCGGCAGACTTGGGCGCGATTTACCCACCACTTCATAGTTGTCGATACGGTTAACCGTGGCtt                                           >  2:569977/1‑301 (MQ=255)
                                                                                                                                                                                                                                                                  cATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTGGCGCGTTTTCGCACGGCAAAGCTGGATGAAACCGGTTCGGCATGGCTGATACAGTTGCTGTTCGGGCAGACCAGCACATTGTCGATGCGCTCCGGCAGACTTGGGCGCGATTTACCCACCACTTCATAGTTATCGATACGGTTAACCGTGGCTTgc                                         >  1:378759/1‑301 (MQ=255)
                                                                                                                                                                                                                                                                                 ccgtttgttaataaaaGTTTGGCTTCCTATATTGTGCCAGATTTTTATTACCAACCGCAATTAATTGGACAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTGGCGCGTTTTCGCACGGCAAAGCTGGATGAAACCGGTTCGGCATGGCTGATACAGTTGCTGTTCGGGCAGACCAGCACATTGTCGATGCGCTCCGGCAGACTTGGGCGCGATTTACCCACCACTTCATAGTTGTCGATACGGTTAACCGTGGCTTGCGGCGCATACAATGcc                          <  1:612746/289‑1 (MQ=255)
                                                                                                                                                                                                                                                                                      aGCGGTAAAAAGTCTGGTTCATTATAGGGGTCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACAGTATGTGAAAACTCTTTTTCACATTATTTGCATTTGAGCGCGATATCATTGGCGCGTTTTCGCACGGCAAAGCTGGATGAAACCGGTTCGGCATGGCTGATACAGTTGCTGTTCGGGCAGACCAGCACATTGTCGATGCGCTCCGGCAGACTTGGGCGCGATTTACCCACCACTTCATAGTTGTCGATACGGTTAACCGTGGCTTGCGGCGCATACAATGCCAGTTg                     <  1:274736/301‑1 (MQ=255)
                                                                                                                                                                                                                                                                                        gtaataaaaGTCTGGCTCCCTATAATGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTGGCGCGTTTTCGCACGGCAAAGCTGGATGAAACCGGTTCGGCATGGCTGATACAGTTGCTGTTCGGGCAGACCAGCACATTGTCGATGCGCTCCGGCAGACTTGGGCGCGATTTACCCACCACTTCATAGTTGTCGATACGGTTAACCGTGGCTTGCGGCGCATACAATGCCAgatcgg                   >  2:613604/6‑296 (MQ=255)
                                                                                                                                                                                                                                                                                             aaaaGTCTGGCTCTTTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTGGCGCGTTTTCGCACGGCAAAGCTGGATGAAACCGGTTCGGCATGGCTGATACAGTTGCTGTTCGGGCAGACCAGCACATTGTCGATGCGCTCCGGCAGACTTGGGCGCGATTTACCCACCACTTCATAGTTGTCGATACGGTTAACCGTGGCTTGGGGCGCATACAATGCCAGTTGATCTACt              <  2:308313/301‑1 (MQ=255)
                                                                                                                                                                                                                                                                                              aaaGTCTGGCTCCCTATAATGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTGGCGCGTTTTCGCACGGCAAAGCTGGATGAAACCGGTTCGGCATGGCTGATACAGTTGCTGTTCGGGCAGACCAGCACATTGTCGATGCGCTCCGGCAGACTTGGGCGCGATTTACCCACCACTTCATAGTTGTCGATACGGTTAACCGTAGCTTGCGGCCCATACAATACCAGttcagaacatt             >  1:360243/1‑292 (MQ=255)
                                                                                                                                                                                                                                                                                              aaaGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTGGCGCGTTTTCGCACGGCAAAGCTGGATGAAACCGGTTCGGCATGGCTGATACAGTTGCTGTTCGGGCAGACCCGCACATTGTCGATGCGCTCCGGCCGCCTTGGGCACGATTTACCCACCACTTCACAGGTGTCGATACGGTTAAACGTGGCTTGCGGCCCATACCATCCCCGTTGATCTCCtt             >  1:117316/1‑301 (MQ=255)
                                                                                                                                                                                                                                                                                               aaGTTTGGGTCATTATAGGGGGCCAGACTTTTATTACCAATCGCAATTAATTGGCCAGCACCACATTATGGGAATACTCTGTTTCACAGTAGTTGCATTTGAGCGCGATATCATTGGCGCGTTTTCGCACGGCAAAGCTGGATGAAACCGGTTCGGCATGGCGGATACAGTTGCTGTTCGGGCAGACCAGCACATTGTCGATGCGCTCCGGCAGACTTGGGCGCGATTTACCCACCACTTCATAGTTGTCGATACGGTTAACCGTGGCTTGCGGCGCATACAATGCCAGTTGATCTACTTg            <  1:569263/301‑1 (MQ=255)
                                                                                                                                                                                                                                                                                                     ggttCATTTAAGGGAGTCAGACTTTTATTACCAACAGAAATTAATTGGCCAGAACCACATTATGGGAAAACGCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTGGCGCGTTTTCGCACGGCAAAGCTGGATGAAACCGGTTCGGCATGGCTGATACAGTTGCTGTGCGGGCAGACCAGCACATTGTCGATGCGCTCCGGCAGACTTGGGCGCGATTTACCCACCACTTCATAGTTGTCGATACGGTTAACCGTGGCTTGCGGCGCATACAATGCCAGTTGATCTACTTGATCTTc      <  1:175672/298‑1 (MQ=255)
                                                                                                                                                                                                                                                                                                     ggCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTGGCGCGTTTTCGCACGGCAAAGCTGGATGAAACCGGTTCGGCATGGCTGATACAGTTGCTGTTCGGGCAGACCAGCACATTGTCGATGCGCTCCGGCAGACTTGGGCGCGATTTACCCACCACTTCATAGTTGTCGATACGGTTAACCGTGGCTTGCGGCGCATACAATGCCAGTTGATATACTTGATCtta      >  2:590269/1‑300 (MQ=255)
                                                                                                                                                                                                                                                                                                      gCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTGGCGCGTTTTCGCACGGCAAAGCTGGATGAAACCGGTTCGGCATGGCTGATACAGTTGCTGTTCGGGCAGACCAGCACATTGTCGATGCGCTCCGGCAGACTTGGGCGCGATTTACCCACCACTTCATAGTTGTCGATACGGTTAACCGTGGCTTGCGGCACATACAATGCCAGTTAAACTACTTAATCTTct     >  1:446223/1‑300 (MQ=255)
                                                                                                                                                                                                                                                                                                         cATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTGGCGCGTTTTCGCACGGCAAAGCTGGATGAAACCGGTTCGGCATGGCTGATACAGTTGCTGTTCGGGCAGACCAGCACATTGTCGATGCGCTCCGGCAGACTTGGGCGCGATTTACCCACCACTTCATAGTTGTCGATACGGTTAACCGTGGCTTGCGGCGCATACAATGCCAGTTGATCTACTTGATCTTCACTc  >  2:383164/1‑301 (MQ=255)
                                                                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                                                                           
TACGGTTGCGAAGTCGTTCAGATCTTTTACAAACACGGTAGTTTTAACGATGTCGCCCACTTTCAGGCCAGCGGCTTCGACGATCGCTTTTACGTTATCCAGCGACTGACGTGCCTGTGCAGCGACGTCTGCCGGTACTTCGCCCGTTTTCGGATTTACCGGGATCTGACCGGAGGTGATGATCATATTGCCCAGATCAACGCCCTGTACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTGGCGCGTTTTCGCACGGCAAAGCTGGATGAAACCGGTTCGGCATGGCTGATACAGTTGCTGTTCGGGCAGACCAGCACATTGTCGATGCGCTCCGGCAGACTTGGGCGCGATTTACCCACCACTTCATAGTTGTCGATACGGTTAACCGTGGCTTGCGGCGCATACAATGCCAGTTGATCTACTTGATCTTCACTC  >  NC_000913/4470653‑4471250

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 9 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑