Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,613,279 A→C 6.9% intergenic (‑28/+36) glsB ← / ← sad glutaminase 2/succinate semialdehyde dehydrogenase, NAD(P)+‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,613,2790AC 6.9% 447.7 / 18.5 130intergenic (‑28/+36)glsB/sadglutaminase 2/succinate semialdehyde dehydrogenase, NAD(P)+‑dependent
Reads supporting (aligned to +/- strand):  ref base A (59/62);  new base C (4/5);  total (63/67)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.82e-01

Only 100 of 137 total aligned reads displayed.

AATCGGGAACCGTCTACTGTAGCCAGCGCCGGAATATAATCCGCGACTTTACCCTGACCAATGAGCGGCCGCACTTGCCGCAAGATGTTTTCTAAAATTGCATTATCCATGGCGACTGCCACTTTCTACTCCTGGACCGCAGGTCTGAAAAGACCTGCGAGTATATCAGAGCTGAATATGTCGCGTCAGATCCGGTCTTTCCACACCGTCTGGATATTACAGAATTCGTGTAAGCCGAAATGGGAAAGCTCACGACCAAAGCCACTCTTTTTCACGCCACCAAAGGCCACTCG  >  NC_000913/1613130‑1613422
                                                                                                                                                     |                                                                                                                                               
aaTCGGGAACCGTCTACTGTAGCCAGCGCCGGAATATAATCCGCGACTTTACCCTGACCAATGAGCGGCCGCACTTGCCGCAAGATGTTTTCTAAAATTGCATTATCCATGGCGACTGCCACTTTCTACTCCTGGACCGCAGGTCTGaaa                                                                                                                                                 >  1:825328/1‑150 (MQ=255)
aaTCGGGAACCGTCTACTGTAGCCAGCGCCGGAATATAATCCGCGACTTTACCCTGACCAATGAGCGGCCGCACTTGCCGCAAGATGTTTTCTAAAATTGCATTATCCATGGCGACTGCCACTTTCTACTCCTGGACCGCAGGTCTGaaa                                                                                                                                                 >  1:117331/1‑150 (MQ=255)
aaTCGGGAACCGTCTACTGTAGCCAGCGCCGGAATATAATCCGCGACTTTACCCTGACCAATGAGCGGCCGCACTTGCCGCAAGATGTTTTCTAAAATTGCATTATCCATGGCGACTGCCACTTTCTACTCCTGGACCGCAGGTCTGaaa                                                                                                                                                 >  2:1573372/1‑150 (MQ=255)
aaTCGGGAACCGTCTACTGTAGCCAGCGCCGGAATATAATCCGCGACTTTACCCTGACCAATGAGCGGCCGCACTTGCCGCAAGATGTTTTCTAAAATTGCATTATCCATGGCGACTGCCACTTTCTACTCCTGGACCGCAGGTCTGaaa                                                                                                                                                 >  2:1573375/1‑150 (MQ=255)
 aTCGGGAACCGTCTACTGTAGCCAGCGCCGGAATATAATCCGCGACTTTACCCTGACCAATGAGCGGCCGCACTTGCCGCAAGATGTTTTCTAAAATTGCATTATCCATGGCGACTGCCACTTTCTACTCCTGGACCGCAGGTCTGaaaa                                                                                                                                                >  1:417902/1‑150 (MQ=255)
 aTCGGGAACCGTCTACTGTAGCCAGCGCCGGAATATAATCCGCGACTTTACCCTGACCAATGAGCGGCCGCACTTGCCGCAAGATGTTTTCTAAAATTGCATTATCCATGGCGACTGCCACTTTCTACTCCTGGACCGCAGGTCTGaaaa                                                                                                                                                >  2:1973422/1‑150 (MQ=255)
    gggAACCGTCTACTGTAGCCAGCGCCGGAATATAATCCGCGACTTTACCCTGACCAATGAGCGGCCGCACTTGCCGCAAGATGTTTTCTAAAATTGCATTATCCATGGCGACTGCCACTTTCTACTCCTGGACCGCAGGTCTGAAAAGAc                                                                                                                                             >  2:2451800/1‑150 (MQ=255)
       aaCCGTCTACTGTAGCCAGCGCCGGAATATAATCCGCGACTTTACCCTGACCAATGAGCGGCCGCACTTGCCGCAAGATGTTTTCTAAAATTGCATTATCCATGGCGACTGCCACTTTCTACTCCTGGACCGCAGGTCTGAAAAGACCTg                                                                                                                                          >  1:1893727/1‑150 (MQ=255)
       aaCCGTCTACTGTAGCCAGCGCCGGAATATAATCCGCGACTTTACCCTGACCAATGAGCGGCCGCACTTGCCGCAAGATGTTTTCTAAAATTGCATTATCCATGGCGACTGCCACTTTCTACTCCTGGACCGCAGGTCTGAAAAGACCTg                                                                                                                                          >  1:398033/1‑150 (MQ=255)
         ccGTCTACTGTAGCCAGCGCCGGAATATAATCCGCTAATTTACCCTGACCAATGAGCGGCCGCACTTGCCGCAAGATTTTTTCTAAATTTGCATTATCCATGGCGACTGCCACTTTCTACTCCTGGACCGCAGGTCTGAAAAGACCTGCg                                                                                                                                        <  1:1653514/150‑1 (MQ=255)
          cGTCTACTGTAGCCAGCGCCTGAATATAAGCCGCGACTTTACCCTTACAAATGAGCGGCCGCACTTGCCGCAATATGTTTTCTAAAATTGCATTATCCATGGCGACTGCCACTTTCTACTCCTGGACCGCAGGTCTTTTCAGACCTGCgg                                                                                                                                       <  1:1043316/150‑2 (MQ=255)
          cGTCTACTGTAGCCAGCGCCGTAATATAATCCGCGACTTTACCCTGACCAATGCGCGGCCGCACTTGCCGCAAGATGTTTTCTAAAATTGCATTATCCATGGCGACTGCCACTTTCTACTCCTGGACCGCAGGTCTTTTCAGACCTGCgg                                                                                                                                       <  1:1234358/150‑2 (MQ=255)
          cGTCTACTGTAGCCAGCGCCGGAATATAATCCGCGACTTTACCCTGACCAATGAGCGGCCGCACTTGCCGCAAGATGTTTTCTAAAATTGCATTATCCATGGCGACTGCCACTTTCTACTCCTGGACCGCAGGTCTTTTCAGACCTGCgg                                                                                                                                       <  1:958579/150‑2 (MQ=255)
          cGTCTACTGTAGCCAGCGCCGGAATATAATCCGCGACTTTACCCTGACCAATGAGCGGCCGCACTTGCCGCAAGATGTTTTCTAAAATTGCATTATCCATGGCGACTGCCACTTTCTACTCCTGGACCGCAGGTCTTTTCAGACCTGCgg                                                                                                                                       <  1:1043302/150‑2 (MQ=255)
          cGTCTACTGTAGCCAGCGCCGGAATATAATCCGCGACTTTACCCTGACCAATGAGCGGCCGCACTTGCCGCAAGATGTTTTCTAAAATTGCATTATCCATGGCGACTGCCACTTTCTACTCCTGGACCGCAGGTCTGAAAAGACCTGCGa                                                                                                                                       <  2:684788/150‑1 (MQ=255)
          cGTCTACTGTAGCCAGCGCCGGAATATAATCCGCGACTTTACCCTGACCAATGAGCGGCCGCACTTGCCGCAAGATGTTTTCTAAAATTGCATTATCCATGGCGACTGCCACTTTCTACTCCTGGACCGCAGGTCTGAAAAGACCTGCGa                                                                                                                                       >  1:430880/1‑150 (MQ=255)
          cGTCTACTGTAGCCAGCGCCGGAATATAATCCGCGACTTTACCCTGACCAATGAGCGGCCGCACTTGCCGCAAGATGTTTTCTAAAATTGCATTATCCATGGCGACTGCCACTTTCTACTCCTGGACCGCAGGTCTGAAAAGACCTGCGa                                                                                                                                       >  1:706205/1‑150 (MQ=255)
            tCTACTGTAGCCAGCGCCGGAATATAATCCGCGACTTTACCCTGACCAATGAGCGGCCGCACTTGCCGCAAGATGTTTTCTAAAATTGCATTATCCATGGCGACTGCCACTTTCTACTCCTGGACCGCAGGTCTGAAAAGACCTGCGAGt                                                                                                                                     <  1:1808664/150‑1 (MQ=255)
                 tGTAGCCAGCGCCGGAATATAATCCGCGACTTTACCCTGACCAATGAGCGGCCGCACTTGCCGCAAGATGTTTTCTAAAATTGCATTATCCATGGCGACTGCCACTTTCTACTCCTGGACCGCAGGTCTGAAAAGACCTGCGAGTATATc                                                                                                                                >  1:158676/1‑150 (MQ=255)
                   tAGCCAGCGCCGGAATATAATCCGCGACTTTACCCTGACCAATGAGCGGCCGCACTTGCCGCAAGATGTTTTCTAAAATTGCATTATCCATGGCGACTGCCACTTTCTACTCCTGGACCGCAGGTCTGAAAAGACCTGCGAGTATATCag                                                                                                                              >  2:2363203/1‑150 (MQ=255)
                   tAGCCAGCGCCGGAATATAATCCGCGACTTTACCCTGACCAATGAGCGGCCGCACTTGCCGCAAGATGTTTTCTAAAATAGCATTATCAATGGCGACTGCCACTTTCTACTCCTGGACCGCAGGTCTGAAAAGACCTGCGAGTATATCag                                                                                                                              >  2:2363187/1‑150 (MQ=255)
                    aGCCAGCGCCGGAATATAATCCGCGACTTTACCCTGACCAATGAGCGGCCGCACTTGCCGCAAGATGTTTTCTAAAATTGCATTATCCATGGCGACTGCCACTTTCTACTCCTGGACCGCAGGTCTGAAAAGACCTGCGAGTATATCaga                                                                                                                             <  1:2736108/150‑1 (MQ=255)
                    aGCCAGCGCCGGAATATAATCCGCGACTTTACCCTGACCAATGAGCGGCCGCACTTGCCGCAAGATGTTTTCTAAAATTGCATTATCCATGGCGACTGCCACTTTCTACTCCTGGACCGCAGGTCTGAAAAGACCTGCGAGTATATCaga                                                                                                                             <  1:2664001/150‑1 (MQ=255)
                      ccAGCGCCGTAATATAATCCGCGACTTTACCCTGACCAATGAGCGGCCGCACTTGCCGCAAGATGTTTTCTAAAATTGCATTATCCATGGCGACTGCCACTTTCTACTCCTGGACCGCAGGTCTGAAAAGACCTGCGAGTATATCAGAGc                                                                                                                           <  1:1452445/150‑1 (MQ=255)
                      ccAGCGCCGGAATATAATCCGCGACTTTACCCTGACCAATGAGCGGCCGCACTTGCCGCAAGATGTTTTCTAAAATTGCATTATCCATGGCGACTGCCACTTTCTACTCCTGGACCGCAGGTCTGAAAAGACCTGCGAGTATATCAGAGc                                                                                                                           <  1:383796/150‑1 (MQ=255)
                       cAGCGCCGGAATATAATCCGCGACTTTACCCTGACCAATGAGCGGCCGCACTTGCCGCAAGATGTTTTCTAAAATTGCATTATCCATGGCGACTGCCACTTTCTACTCCTGGACCGCAGGTCTGAAAAGACCTGCGAGTATATCAGAGCt                                                                                                                          >  2:957070/1‑150 (MQ=255)
                                 atatAATCCGCGACTTTACCCTGACCAATGAGCGGCCGCACTTGCCGCAAGATGTTTTCTAAAATTGCATTATCCATGGCGACTGCCACTTTCTACTCCTGGACCGCAGGTCTGAAAAGACCTGCGAGTATATCAGAGCTGAATATGTcg                                                                                                                <  2:1730310/150‑1 (MQ=255)
                                      aTCCGCGACTTTACCCTGACCACTGAGCGGCCGCACTTGCCGCAAGATGTTTTCTAAAATTGCATTATCCATGGCGACTGCCACTTTCTACTCCTGGACCGCAGGTCTGAAAAGACCTGCGAGTATATCAGAGCTGAATATGTCGCGTCa                                                                                                           <  2:786568/150‑1 (MQ=255)
                                       tCCGCGACTTTACCCTGACCAATGAGCGGCCGCACTTGCCGCAAGATGTTTTCTAAAATTGCATTATCCATGGCGACTGCCACTTTCTACTCCTGGACCGCAGGTCTGAAAAGACCTGCGAGTATATCAGAGCTGAATATGTCGCGTCAg                                                                                                          >  2:1353761/1‑150 (MQ=255)
                                       tCCGCGACTTTACCCTGACCAATGAGCGGCCGCACTTGCCGCAAGATGTTTTCTAAAATTGCATTATCCATGGCGACTGCCACTTTCTACTCCTGGACCGCAGGTCTGAAAAGACCTGCGAGTATATCAGAGCTGAATATGTCGCGTCAg                                                                                                          >  2:836844/1‑150 (MQ=255)
                                       tCCGCGACTTTACCCTGACCAATGAGCGGCCGCACTTGCCGCAAGATGTTTTCTAAAATTGCATTATCCATGGCGACTGCCACTTTCTACTCCTGGACCGCAGGTCTGAAAAGACCTGCGAGTATATCAGAGCTGAATATGTCGCGTCAg                                                                                                          >  2:1591386/1‑150 (MQ=255)
                                          gcgACTTTACCCTGACCAATGAGCGGCCGCACTTGCCGCAAGATGTTTTCTAAAATTGCATTATCCATGGCGACTGCCACTTTCTACTCCTGGACCGCAGGTCTGAAAAGACCTGCGAGTATATCAGAGCTGAATATGTCGCGTCAGATc                                                                                                       >  1:2088806/1‑150 (MQ=255)
                                          gcgACTTTACCCTGACCAATGAGCGGCCGCACTTGCCGCAAGATGTTTTCTAAAATTGCATTATCCATGGCGACTGCCACTTTCTACTCCTGGACCGCAGGTCTGAAAAGACCTGCGAGTATATCAGAGCTGAATATGTCGCGTCAGATc                                                                                                       >  2:848296/1‑150 (MQ=255)
                                          gcgACTTTACCCTGACCAATGAGCGGCCGCACTTGCCGCAAGATGTTTTCTAAAATTGCATTATCCATGGCGACTGCCACTTTCTACTCCTGGACCGCAGGTCTGAAAAGACCTGCGAGTATATCAGAGCTGAATATGTCGCGTCAGATc                                                                                                       >  2:554313/1‑150 (MQ=255)
                                               tttACCCTGACCAATGAGCGGCCGCACTTGCCGCAAGATGTTTTCTAAAATTGCATTATCCATGGCGACTGCCACTTTCTACTCCTGGACCGCAGGTCTGAAAAGACCTGCGAGTATATCAGAGCTGAATATGTCGCGTCAGATCCGGTc                                                                                                  >  2:1697568/1‑150 (MQ=255)
                                               tttACCCTGACCAATGAGCGGCCGCACTTGCCGCAAGATGTTTTCTAAAATTGCATTATCCATGGCGACTGCCACTTTCTACTCCTGGACCGCAGGTCTGAAAAGACCTGCGAGTATATCAGAGCTGAATATGTCGCCTCAGATCCGGTc                                                                                                  >  2:2098373/1‑150 (MQ=255)
                                                 tACCCTGACCAATGAGCGGCCGCACTTGCCGCAAGATGTTTTCTAAAATTGCATTATCCATGGCGACTGCCACTTTCTACTCCTGGACCGCAGGTCTGAAAAGACCTGCGAGTATATCAGAGCTGAATATGTCGCGTCAGATCCGGTCtt                                                                                                <  1:302331/150‑1 (MQ=255)
                                                      tGACCAATGAGCGGCCGCACTTGCCGCAAGATGTTTTCTAAAATTGCATTATCCATGGCGACTGCCACTTTCTACTCCTGGACCGCAGGTCTGAAAAGACCTGCGAGTCTATCAGAGCTGAATATGTCGCGTCAGATCCGGTCTTTCcac                                                                                           >  2:1426264/1‑150 (MQ=255)
                                                      tGACCAATGAGCGGCCGCACTTGCCGCAAGATGTTTTCTAAAATTGCATTATCCATGGCGACTGCCACTTTCTACTCCTGGACCGCAGGTCTGAAAAGACCTGCGAGTATATCAGAGCTGAATATGTCGCGTCAGATCCGGTCTTTCcac                                                                                           >  1:548552/1‑150 (MQ=255)
                                                      tGACCAATGAGCGGCCGCACTTGCCGCAAGATGTTTTCTAAAATTGCATTATCCATGGCGACTGCCACTTTCTACTCCTGGACCGCAGGTCTGAAAAGACCTGCGAGTATATCAGAGCTGAATATGTCGCGTCAGATCCGGTCTTTCcac                                                                                           >  1:1616402/1‑150 (MQ=255)
                                                      tGACCAATGAGCGGCCGCACTTGCCGCAAGATGTTTTCTAAAATTGCATTATCCATGGCGACTGCCACTTTCTACTCCTGGACCGCAGGTCTGAAAAGACCTGCGAGTATATCAGAGCTGAATATGTCGCGTCAGATCCGGTCTTTCcac                                                                                           >  1:1115399/1‑150 (MQ=255)
                                                      tGACCAATGAGCGGCCGCACTTGCCGCAAGATGTTTTCTAAAATTGCATTATCCATGGCGACTGCCACTTTCTACTCCTGGACCGCAGGTCTGAAAAGACCTGCGAGTATATCAGAGCTGAATATGTCGCGTCAGATCCGGTCTTTCcac                                                                                           >  2:50206/1‑150 (MQ=255)
                                                      tGACCAATGAGCGGCCGCACTTGCCGCAAGATGTTTTCTAAAATTGCATTATCCATGGCGACTGCCACTTTCTACTCCTGGACCGCAGGTCTGAAAAGACCTGCGAGTATATCAGAGCTGAATATGTCGCGTCAGATCCGGTCTGTCcac                                                                                           >  2:2232396/1‑150 (MQ=255)
                                                      tGACCAATGAGCGGCCGCACTTGCCGCAAGATGTTGTCTAAAATTGCATTATCCATGGCGACTGCCACTTTCTACTCCTGGACCGCAGGTCTGAAAAGACCTGCGAGTATATCAGAGCTGAATATGTCGCGTCAGATCCGGTCTTTCcac                                                                                           >  2:2232364/1‑150 (MQ=255)
                                                                gCGGCCGCACTTGCCGCAAGATGTTTTCTAAAATTGCATTATCCATGGCGACTGCCACTTTCTACTCCTGGACCGCAGGTCTGAAAAGACCTGCGAGTATATCAGAGCTGAATATGTCGCGTCAGATCCGGTCTTTCCACACCGTCTGGa                                                                                 >  2:1992029/1‑150 (MQ=255)
                                                                gCGGCCGCACTTGCCGCAAGATGTTTTCTAAAATTGCATTATCCATGGCGACTGCCACTTTCTACTCCTGGACCGCAGGTCTGAAAAGACCTGCGAGTATATCAGAGCTGAATATGTCGCGTCAGATCCGGTCTTTCCACACCGTCTGGa                                                                                 >  2:203072/1‑150 (MQ=255)
                                                                 cGGCCGCACTTGCCGCAAGATGTTTTCTAAAATTGCATTATCCATGGCGACTGCCACTTTCTACTCCTGGACCGCAGGTCTGAAAAGACCTGCGAGTATATCAGAGCTGAATATGTCGCGTCAGATCCGGTCTTTCCACACCGTCTGGat                                                                                >  1:2334120/1‑150 (MQ=255)
                                                                    ccGCACTTGCCGCAAGATGTTTTCTAAAATTGCATTATCCATGGCGACTGCCACTTTCTACTCCTGGACCGCAGGTCTGAAAAGACCTGCGAGTATATCAGAGCTGAATATGTCGCGTCAGATCCGGTCTTTCCACACCGTCTGGATAtt                                                                             <  1:1143204/150‑1 (MQ=255)
                                                                    ccGCACTTGCCGCAAGATGTTTTCTAAAATTGCATTATCCATGGCGACTGCCACTTTCTACTCCTGGACCGCAGGTCTGAAAAGACCTGCGAGTATATCAGAGCTGAATATGTCGCGTCAGATCCGGTCTTTCCACACCGTCTGGATAtt                                                                             <  1:1928385/150‑1 (MQ=255)
                                                                       cACTTGCCGCAAGTTTTTTTCTAAAATTGCATTATCGATGGCGACTGCCACTTTCTACTCCTGGAGCGCAGGTCTGAAAAGACCTGCGAGTATATCAGAGCTGTATATGTCGCGTCAGATCCGGTCTTTCCACACCGTCTGGATATTACa                                                                          <  1:2425110/150‑1 (MQ=255)
                                                                       cACTTGCCGCAAGCTGTTTTCTAACATTGCATGATCCCTGGCGACTGCCACTTTCTCCTCCTGGACCGCAGGTCTGAAAAGACCTGCGAGGATATCAGAGCTGAATATGTCGCGTCAGCTCCGGTCTTTCCACACCGTCTGGATATTACa                                                                          <  2:898513/150‑1 (MQ=255)
                                                                       cACTTGCCGCAAGATGTTTTCTAAAATTGCATTATCCATGGTGACTGCCACTTTCTACTCCTGGACCGCAGGTCTGAAAAGACCTGCGAGTATATCAGAGCTGAATATGTCGCGTCAGATCCGGTCTTTCCACACCGTCTGGATATTACa                                                                          <  2:2634600/150‑1 (MQ=255)
                                                                       cACTTGCCGCAAGATGTTTTCTAAAATTGCATTATCCATGGCGACTGCCACTTTCTACTCCTGGACCGCAGGTCTGAAAAGACCTGCGAGTATATCAGAGCTGAATATGTCGCGTCAGATCCGGTCTTTCCACACCGTCTGGATATTACa                                                                          <  1:95105/150‑1 (MQ=255)
                                                                       cACTTGCCGCAAGATGTTTTCTAAAATTGCATTATCCATGGCGACTGCCACTTTCTACTCCTGGACCGCAGGTCTGAAAAGACCTGCGAGTATATCAGAGCTGAATATGTCGCGTCAGATCCGGTCTTTCCACACCGTCTGGATATTACa                                                                          <  2:1838941/150‑1 (MQ=255)
                                                                       cACTTGCCGCAAGATGTTTTCTAAAATTGCATTATCCATGGCGACTGCCACTTTCTACTCCTGGACCGCAGGTCTGAAAAGACCTGCGAGTATATCAGAGCTGAATATGTCGCGTCAGATCCGGTCTTTCCACACCGTCTGGATATTACa                                                                          <  2:1782314/150‑1 (MQ=255)
                                                                       cACTTGCCGCAAGATGGTTTCTAAAAGTGCATTATCCCTGGCGACTGCCACTTTCTACTCCTGGACCGCAGGTCTGAAAAGACCTGCGAGTATATCAGAGCTGAATATGTCGCGTCAGATCCGGTCTTTCCACACCGTCTGGATATTACa                                                                          <  2:1813211/150‑1 (MQ=255)
                                                                       cACTTGCCACAAGATGTTTTCTAAAATTGCATTATCCATGGCGACTGCCACTTTCTACTCCTGGACCGCAGGTCTGAAAAGACCTGCGAGTATATCAGAGCTGAATATGTCGCGTCAGATCCGGTCTTTCCACACCGTCTGGATATTACa                                                                          <  2:2115116/150‑1 (MQ=255)
                                                                                   gATGTTTTCTAAAATTGCATTATCCATGGCGACTGCCACTTTCTACTCCTGGACCGCAGGTCTGAAAAGACCTGCGAGTATATCAGAGCTGAATATGTCGCGTCAGATCCGGTCTTTCCACACCGTCTGGATATTACAGAATTCGTGTaa                                                              <  2:1744143/150‑1 (MQ=255)
                                                                                     tGTTTTCTAAAATTGCATTATCCATGGCGACTGCCACTTTCTACTCCTGGACCGCAGGTCTGAAAAGACCTGCGAGTATATCAGAGCTGAATATGTCGCGTCAGATCCGGTCTTTCCACACCGTCTGGATATTACAGAATTCGTGTAAGc                                                            <  1:1050063/150‑1 (MQ=255)
                                                                                     tGTTTTCTAAAATTGCATTATCCATGGCGACTGCCACTTTCTACTCCTGGACCGCAGGTCTGAAAAGACCTGCGAGTATATCAGAGCTGAATATGTCGCGTCAGATCCGGTCTTTCCACACCGTCTGGATATTACAGAATTCGTGTAAGc                                                            <  1:335568/150‑1 (MQ=255)
                                                                                       ttttCTAAAATTGCATTATCCATGGCGACTGCCACTTTCTACTCCTGGACCGCAGGTCTGAAAAGACCTGCGAGTATATCAGAGCTGAATATGTCGCGTCAGATCCGGTCTTTCCACACCGTCTGGATATTACAGAATTCGTGTAAGCCg                                                          >  1:484557/1‑150 (MQ=255)
                                                                                               aaTTGCATTATCCATGGCGACTGCCACTTTCTACTCCTGGACCGCAGGTCTGAAAAGACCTGCGAGTATATCAGAGCTGAATATGTCGCGTCAGATCCGGTCTTTCCACACCGTCTGGATATTACAGAATTCGTGTAAGCCGAAATGGGa                                                  <  1:2710910/150‑1 (MQ=255)
                                                                                                 ttGCATTATCCATGGCGACTGCCACTTTCTACTCCTGGACCGCAGGTCTGAAAAGACCTGCGAGTATATCAGAGCTGAATATGTCGCGTCAGATCCGGTCTTTCCACACCGTCTGGATATTACAGAATTCGTGTAAGCCGAAATGGGaaa                                                <  1:594441/150‑1 (MQ=255)
                                                                                                 ttGCATTATCCATGGCGACTGCCACTTTCTACTCCTGGACCGCAGGTCTGAAAAGACCTGCGAGTATATCAGAGCTGAATATGTCGCGTCAGATCCGGTCTTTCCACACCGTCTGGATATTACAGAATTCGTGTAAGCCGAAATGGGaaa                                                <  2:1375759/150‑1 (MQ=255)
                                                                                                 gtGCATTATCCATGGCGACTGCCACTTTCTACTCCTGGACCGCAGGTCTGAAAAGACCTGCGAGTATATCAGAGCTGAATATGTCGCGTCAGATCCGGTCTTTCCACACCGTCTGGATATTACAGAATTCGTGTAAGCCGAAATGGGaaa                                                <  2:2280197/149‑1 (MQ=255)
                                                                                                     aTTATCCATGGCGACTGCCACTTTCTACTCCTGGACCGCAGGTCTGAAAAGACCTGCGAGTATATCAGAGCTGAATATGTCGCGTCAGATCCGGTCTTTCCACACCGTCTGGATATTACAGAATTCGTGTAAGCCGAAATGGGAAAGCTc                                            >  2:28921/1‑150 (MQ=255)
                                                                                                     aTTATCCATGGCGACTGCCACTTTCTACTCCTGGACCGCAGGTCTGAAAAGACCTGCGAGTATATCAGAGCTGAATATGTCGCGTCAGATCCGGTCTTTCCACACCGTCTGGATATTACAGAATTCGTGTAAGCCGAAATGGGAAAGCTc                                            >  2:29121/1‑150 (MQ=255)
                                                                                                         tCCATGGCGACTGCCACTTTCTACTCCTGGACCGCAGGTCTGAAAAGACCTGCGAGTATATCAGAGCTGAATATGTCGCGTCAGATCCGGTCTTTCCACACCGTCTGGATATTACAGAATTCGTGTAAGCCGAAATGGGAAAGCTCACGa                                        <  1:1886657/150‑1 (MQ=255)
                                                                                                                  aCTGCCACTTTCTACTCCTGGACCGCAGGTCTGAAAAGACCTGCGAGTATATCAGAGCTGAATATGTCGCGTCAGATCCGGTCTTTCCACACCGTCTGGATATTACAGAATTCGTGTAAGCCGAAATGGGAAAGCTCACGACCAAAGCCa                               <  2:2627292/150‑1 (MQ=255)
                                                                                                                  aCTGCCACTTTCTACTCCTGGACCGCAGGTCTGAAAAGACCTGCGAGTATATCAGAGCTGAATATGTCGCGTCAGATCCGGTCTTTCCACACCGTCTGGATATTACAGAATTCGTGTAAGCCGAAATGGGAAAGCTCACGACCAAAGCCa                               <  2:2070427/150‑1 (MQ=255)
                                                                                                                  aCTGCCACTTTCTACTCCTGGACCGCAGGTCTGAAAAGACCTGCGAGTATATCAGAGCTGAATATGTCGCGTCAGATCCGGTCTTTCCACACCGTCTGGATATTACAGAATTCGTGTAAGCCGAAATGGGAAAGCTCACGACCAAAGCCa                               <  1:1882622/150‑1 (MQ=255)
                                                                                                                       cACTTTCTACTCCTGGACCGCAGGTCTGAAAAGACCTGCGAGTATATCAGAGCTGAATATGTCGCGTCAGATCCGGTCTTTCCACACCGTCTGGATATTACAGAATTCGTGTAAGCCGAAATGGGAAAGCTCACGACCAAAGCCACTCtt                          >  1:1880204/1‑150 (MQ=255)
                                                                                                                       cACTTTCTACTCCTGGACCGCAGGTCTGAAAAGACCTGCGAGTATATCAGAGCTGAATATGTCGCGTCAGATCCGGTCTTTCCACACCGTCTGGATATTACAGAATTCGTGTAAGCCGAAATGGGAAAGCTCACGACCAAAGCCACTCtt                          >  2:898527/1‑150 (MQ=255)
                                                                                                                          tttCTACTCCTGGACCGCAGGTCTGAAAAGACCTGCTTGTATATCAGAGCTTAATATGTCGCGTCAGATCCGGTCTTTCCACACCGTCTGGATATTACATATTTCGTGTAAGCCGAAATGGGAAAGCTCACGACCAAAGCCACTCttttt                       <  1:604706/150‑1 (MQ=255)
                                                                                                                          tttCTACTCCTGGACCGCAGGTCTGAAAAGACCTGCGAGTATATCAGAGCTGAATATGTCGCGTCAGATCCGGTCTTTCCACACCGTCTGGATATTACAGAATTCGTGTAAGCCGAAATGGGAAAGCTCACGACCAAAGCCACTCttttt                       <  2:936704/150‑1 (MQ=255)
                                                                                                                              tACTCCTGGACCGCAGGTCTGAAAAGACCTGCGAGTATATCAGAGCTGAATATGTCGCGTCAGATCCGGTCTTTCCACACCGTCTGGATATTACAGAATTCGTGTAAGCCGAAATGGGAAAGCTCACGACCAAAGCCACTCTTTTTCACg                   <  2:2176678/150‑1 (MQ=255)
                                                                                                                              tACTCCTGGACCGCAGGTCTGAAAAGACCTGCGAGTATATCAGAGCTGAATATGTCGCGTCAGATCCGGTCTTTCCACACCGTCTGGATATTACAGAATTCGTGTAAGCCGAAATGGGAAAGCTCACGACCAAAGCCACTCTTTTTCACg                   <  2:484557/150‑1 (MQ=255)
                                                                                                                              tACTCCTGGACCGCAGGTCTGAAAAGACCTGCGAGTATATCAGAGCTGAATATGTCGCGTCAGATCCGGTCTTTCCACACCGTCTGGATATTACAGAATTCGTGTAAGCCGAAATGGGAAAGCTCACGACCAAAGCCACTCTTTTTCACg                   <  2:1868319/150‑1 (MQ=255)
                                                                                                                               aCTCCTGGACCGCAGGTCTGAAAAGACCTGCGAGTATATCAGAGCTGAATATGTCGCGTCAGATCCGGTCTTTCCACACCGTCTGGATATTACAGAATTCGTGTAAGCCGAAATGGGAAAGCTCACGACCAAAGCCACTCTTTTTCACGc                  >  2:111783/1‑150 (MQ=255)
                                                                                                                               aCTCCTGGACCGCAGGTCTGAAAAGACCTGCGAGTATATCAGAGCTGAATATGTCGCGTCAGATCCGGTCTTTCCACACCGTCTGGATATTACAGAATTCGTGTAAGCCGAAATGGGAAAGCTCACGACCAAAGCCACTCTTTTTCACGc                  >  2:2182505/1‑150 (MQ=255)
                                                                                                                                   cTGGACCGCAGGTCTGAAAAGACCTGCGAGTATATCAGAGCTGAATATGTCGCGTCAGATCCGGTCTTTCCACACCGTCTGGATATTACAGCATTCGTGTAAGCCGAAATGGGAAAGCTCACGACCAAAGCCACTCTTTTTCACGccacc              <  2:1780866/150‑1 (MQ=255)
                                                                                                                                        tcGCAGGTCTTTTCAGACCTGCGAGTATATCAGAGCTGAATATGTCGCGTCAGATCCGGTCTTTCCACACCGTCTGGATATTACAGAATTCGTGTAAGCCGAAATGGGAAAGCTCACGACCAAAGCCACTCTTTTTCACGCCACCAAAgg         >  2:2325430/2‑150 (MQ=255)
                                                                                                                                        tcGCAGGTCTTTTCAGACCTGCGAGTATATCAGAGCTGAATATGTCGCGTCAGATCCGGTCTTTCCACACCGTCTGGATATTACAGAATTCGTGTAAGCCGAAATGGGAAAGCTCACGACCAAAGCCACTCTTTTTCACGCCACCAAAgg         >  2:1020780/2‑150 (MQ=255)
                                                                                                                                        tcGCAGGTCTTTTCAGACCTGCGAGTATATCAGAGCTGAATATGTCGCGTCAGATCCGGTCTTTCCACACCGTCTGGATATTACAGAATTCGTGTAAGCCGAAATGGGAAAGCTCACGACCAAAGCCACTCTTTTTCACGCCACCAAAgg         >  2:1427511/2‑150 (MQ=255)
                                                                                                                                         cGCAGGTCTGAAAAGCCCTGGGATTATATCAGAGCTGATTATGTCGCGTCAGATCCGGTCTTTCCACACCGTCTGTATATTACAGATTTCGTGTAAGCCGAAATGTTAAAGCTCACGACCAAAGCCACTCTTTTTCACGCCACCAAAGGc        <  1:836844/150‑1 (MQ=255)
                                                                                                                                         cGCAGGTCTGAAAAGACCTGCGAGTATATCAGAGCTGAATATGTCGCGTCAGATCCGGTCTTTCCACACCGTCTGGATATTACAGAATTCGTGTAAGCCGAAATGGGAAAGCTCACGACCAAAGCCACTCTTTTTCACGCCACCAAAGGc        <  1:1794617/150‑1 (MQ=255)
                                                                                                                                         cGCAGGTCTGAAAAGACCTGCGAGTATATCAGAGCTGAATATGTCGCGTCAGATCCGGTCTTTCCACACCGTCTGGATATTACAGAATTCGTGTAAGCCGAAATGGGAAAGCTCACGACCAAAGCCACTCTTTTTCACGCCACCAAAGGc        >  1:1409469/1‑150 (MQ=255)
                                                                                                                                         cGCAGGTCTGAAAAGACCTGCGAGTATATCAGAGCTGAATATGTCGCGTCAGATCCGGTCTTTCCACACCGTCTGGATATTACAGAATTCGTGTAAGCCGAAATGGGAAAGCTCACGACCAAAGCCACTCTTTTTCACGCCACCAAAGGc        <  1:1353761/150‑1 (MQ=255)
                                                                                                                                           cAGGTCTGAAAAGACCTGCGAGTATATCAGAGCTGAATATGTCGCGTCAGATCCGGTCTTTCCACACCGTCTGGATATTACAGAATTCGTGTAAGCCGAAATGGGAAAGCTCACGACCAAAGCCACTCTTTTTCACGCCACCAAAGGCCa      <  2:2483649/150‑1 (MQ=255)
                                                                                                                                           cAGGTCTGAAAAGACCTGCGAGTATATCAGAGCTGAATATGTCGCGTCAGATCCGGTCTTTCCACACCGTCTGGATATTACAGAATTCGTGTAAGCCGAAATGGGAAAGCTCACGACCAAAGCCACTCTTTTTCACGCCACCAAAGGCCa      <  2:2682465/150‑1 (MQ=255)
                                                                                                                                           cAGGTCTGAAAAGACCTGCGAGTATATCAGAGCTGAATATGTCGCGTCAGATCCGGTCTTTCCACACCGTCTGGATATTACAGAATTCGTGTAAGCCGAAATGGGAAAGCTCACGACCAAAGCCACTCTTTTTCACGCCACCAAAGGCCa      <  2:581115/150‑1 (MQ=255)
                                                                                                                                            aGGTCTGAAAAGACCTGCGAGTATATCAGAGCTGAATATGTCGCGTCAGATCCGGTCTTTCCACACCGTCTGGATATTACAGAATTCGTGTAAGCCGAAATGGGAAAGCTCACGACCAAAGCCACTCTTTTTCACGCCACCAAAg          <  1:164224/145‑1 (MQ=255)
                                                                                                                                              gTCTGAAAAGACCTGCGAGTATATCAGAGCTGAATATGTCGCGTCAGATCCGGTCTTTCCACACCGTCTGGATATTACAGAATTCGTGTAAGCCGAAATGGGAAAGCTCACGACCAAAGCCACTCTTTTTCACGCCACCAAAGGCCACTc   >  1:1065202/1‑150 (MQ=255)
                                                                                                                                              gTCTGAAAAGACCTGCGAGTATATCAGAGCTGAATATGTCGCGTCAGATCCGGTCTTTCCACACCGTCTGGATATTACAGAATTCGTGTAAGCCGAAATGGGAAAGCTCACGACCAAAGCCACTCTTTTTCACGCCACCAAAGGCCACTc   >  2:738512/1‑150 (MQ=255)
                                                                                                                                              gTCTGAAAAGACCTGCGAGTATATCAGAGCTGAATATGTCGCGTCAGATCCGGTCTTTCCACACCGTCTGGATATTACAGAATTCGTGTAAGCCGAAATGGGAAAGCTCACGACCAAAGCCACTCTTTTTCACGCCACCAAAGGCCACTc   >  1:1065186/1‑150 (MQ=255)
                                                                                                                                              gTCTGAAAAGACCTGCGAGTATATCAGAGCTGAATATGTCGCGTCAGATCCGGTCTTTCCACACCGTCTGGATATTACAGAATTCGTGTAAGCCGAAATGGGAAAGCTCACGACCAAAGCCACTCTTTTTCACGCCACCAAAGGCCACTc   >  2:2477615/1‑150 (MQ=255)
                                                                                                                                              gTCTGAAAAGACCTGCGAGTATATCAGAGCTGAATATGTCGCGTCAGATCCGGTCTTTCCACACCGTCTGGATATTACAGAATTCGTGTAAGCCGAAATGGGAAAGCTCACGACCAAAGCCACTCTTTTTCACGCCACCAAAGGCCACTc   >  2:193782/1‑150 (MQ=255)
                                                                                                                                              gTCTGAAAAGACCTGCGAGTATATCAGAGCTGAATATGTCGCGTCAGATCCGGTCTTTCCACACCGTCTGGATATTACAGAATTCGTGTAAGCCGAAATGGGAAAGCTCACGACCAAAGCCACTCTTTTTCACGCCACCAAAGGCCACTc   >  2:193754/1‑150 (MQ=255)
                                                                                                                                               tCTGAAAAGACCTGCGAGTATATCAGAGCTGAATATGTCGCGTCAGATGCGGTCTTTCCACACCGTCTGGATATTACAGAATTCGTGTAAGCCGAAATGGGAAAGCTCACGACCAAAGCCACTCTTTTTCACGCCACCAAAGGCCACTcg  >  1:1335710/1‑150 (MQ=255)
                                                                                                                                               tCTGAAAAGACCTGCGAGTATATCAGAGCTGAATATGTCGCGTCAGATCCGGTCTTTCCACACCGTCTGGATATTACAGAATTCGTGTAAGCCGAAATGGGAAAGCTCACGACCAAAGCCACTCTTTTTCACGCCACCAAAGGCCACTcg  >  1:2234675/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                               
AATCGGGAACCGTCTACTGTAGCCAGCGCCGGAATATAATCCGCGACTTTACCCTGACCAATGAGCGGCCGCACTTGCCGCAAGATGTTTTCTAAAATTGCATTATCCATGGCGACTGCCACTTTCTACTCCTGGACCGCAGGTCTGAAAAGACCTGCGAGTATATCAGAGCTGAATATGTCGCGTCAGATCCGGTCTTTCCACACCGTCTGGATATTACAGAATTCGTGTAAGCCGAAATGGGAAAGCTCACGACCAAAGCCACTCTTTTTCACGCCACCAAAGGCCACTCG  >  NC_000913/1613130‑1613422

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: