Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,613,278 A→T 9.7% intergenic (‑27/+37) glsB ← / ← sad glutaminase 2/NAD(+)‑dependent succinate semialdehyde dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,613,2780AT 9.7% 380.1 / 23.7 131intergenic (‑27/+37)glsB/sadglutaminase 2/NAD(+)‑dependent succinate semialdehyde dehydrogenase
Reads supporting (aligned to +/- strand):  ref base A (58/60);  new base T (3/10);  total (61/70)
Fisher's exact test for biased strand distribution p-value = 8.61e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.87e-01

Only 100 of 138 total aligned reads displayed.

CAATCGGGAACCGTCTACTGTAGCCAGCGCCGGAATATAATCCGCGACTTTACCCTGACCAATGAGCGGCCGCACTTGCCGCAAGATGTTTTCTAAAATTGCATTATCCATGGCGACTGCCACTTTCTACTCCTGGACCGCAGGTCTGAAAAGACCTGCGAGTATATCAGAGCTGAATATGTCGCGTCAGATCCGGTCTTTCCACACCGTCTGGATATTACAGAATTCGTGTAAGCCGAAATGGGAAAGCTCACGACCAAAGCCACTCTTTTTCACGCCACCAAAGGCCACTCGCGCGT  >  NC_000913/1613129‑1613427
                                                                                                                                                     |                                                                                                                                                     
cAATCGGGAACCGTCTACTGTAGCCAGCGCCGGAATATAATCCGCGACTTTACCCTGACCAATGAGCGGCCGCACTTGCCGCAAGATGTTTTCTAAAATTGCATTATCCATGGCGACTGCCACTTTCTACTCCTGGACCGCAGGTCTGaa                                                                                                                                                       >  1:2352839/1‑150 (MQ=255)
  aTCGGGAACCGTCTACTGTAGCCAGCGCCGGAATATAATCCGCGACTTTACCCTGACCAATGAGCGGCCGCACTTGCCGCAAGATGTTTTCTAAAATTGCATTATCCATGGCGACTGCCACTTTCTACTCCTGGACCGCAGGTCTGaaaa                                                                                                                                                     <  2:364032/150‑1 (MQ=255)
  aTCGGGAACCGTCTACTGTAGCCAGCGCCGGAATATAATCCGCGACTTTACCCTGACCAATGAGCGGCCGCACTTGCCGCAAGATGTTTTCTAAAATTGCATTATCCATGGCGACTGCCACTTTCTACTCCTGGACCGCAGGTCTGaaaa                                                                                                                                                     >  1:645719/1‑150 (MQ=255)
   tCGGGAACCGTCTACTGTAGCCAGCGCCGGAATATAATCCGCGACTTTACCCTGACCAATGAGCGGCCGCACTTGCCGCAAGATGTTTTCTAAAATTGCATTATCCATGGCGACTGCCACTTTCTACTCCTGGACCGCAGGTCTGAAAAg                                                                                                                                                    <  2:2277323/150‑1 (MQ=255)
   tCGGGAACCGTCTACTGTAGCCAGCGCCGGAATATAATCCGCGACTTTACCCTGACCAATGAGCGGCCGCACTTGCCGCAAGATGTTTTCTAAAATTGCATTATCCATGGCGACTGCCACTTTCTACTCCTGGACCGCAGGTCTGAAAAg                                                                                                                                                    <  2:883579/150‑1 (MQ=255)
     gggAACCGTCTACTGTAGCCAGCGCCGGAATATAATCCGCGACTTTACCCTGACCAATGAGCGGCCGCACTTGCCGCAAGATGTTTTCTAAAATTGCATTATCCATGGCGACTGCCACTTTCTACTCCTGGACCGCAGGTCTGGAAAGAc                                                                                                                                                  >  1:1104665/1‑150 (MQ=255)
      ggAACCGTCTACTGTAGCCAGCGCCGGAATATAATCCGCGACTTTACCCTGACCAATGAGCGGCCGCACTTGCCGCAAGATGTTTTCTAAAATTGCATTATCCATGGCGACTGCCACTTTCTACTCCTGGACCGCAGGTCTGAAAAGAcc                                                                                                                                                 >  2:17848/1‑150 (MQ=255)
      ggAACCGTCTACTGTAGCCAGCGCCGGAATATAATCCGCGACTTTACCCTGACCAATGAGCGGCCGCACTTGCCGCAAGATGTTTTCTAAAATTGCATTATCCATGGCGACTGCCACTTTCTACTCCTGGACCGCAGGTCTGAAAAGAcc                                                                                                                                                 <  2:539534/150‑1 (MQ=255)
       gAACCGTCTACTGTAGCCAGCGCCGGAATATAATCCGCGACTTTACCCTGACCAATGAGCGGCCGCACTTGCCGCAAGATGTTTTCTAAAATTGCATTATCCATGGCGACTGCCACTTTCTACTCCTGGACCGCAGGTCTGAAAAGACCt                                                                                                                                                <  2:2371562/150‑1 (MQ=255)
        aaCCGTCTACTGTAGCCAGCGCCGGAATATAATCCGCGACTTTACCCTGACCAATGAGCGGCCGCACTTGCCGCAAGATGTTTTCTAAAATTGCATTATCCATGGCGACTGCCACTTTCTACTCCTGGACCGCAGGTCTGAAAAGACCTg                                                                                                                                               <  1:1119569/150‑1 (MQ=255)
         aCCGTCTACTGTAGCCAGCGCCGGAATATAATCCGCGACTTTACCCTGACCAATGAGCGGCCGCACTTGCCGCAAGATGTTTTCTAAAATTGCATTATCCATGGCGACTGCCACTTTCTACTCCTGGACCGCAGGTCTTTTCAGACCTGc                                                                                                                                              <  1:1369422/150‑1 (MQ=255)
          ccGTCTACTGTAGCCAGCGCCGGAATATAATCCGCGACTTTACCCTGACCAATGAGCGGCCGCACTTGCCGCAAGATGTTTTCTAAAATTGCATTATCCATGGCGACTGCCACTTTCTACTCCTGGACCGCAGGTCTGAAAAGACCTGCg                                                                                                                                             <  1:228645/150‑1 (MQ=255)
           cGTCTACTGTAGCCAGCGCCGGAATATAATCCGCGACTTTACCCTGACCAATGAGCGGCCGCACTTGCCGCAAGATGTTTTCTAAAATTGCATTATCCATGGCGACTGCCACTTTCTACTCCTGGACCGCAGGTCTTTTCAGACCTGCgg                                                                                                                                            <  1:2139585/150‑2 (MQ=255)
             tCTACTGTAGCCAGCGCCGGAATATAATCCGCGACTTTACCCTGACCAATGAGCGGCCGCACTTGCCGCAAGATGTTTTCTAAAATTGCATTATCCATGGCGACTGCCACTTTCTACTCCTGGACCGCAGGTCTTTTCAGACCTGCggtc                                                                                                                                          <  1:466263/150‑4 (MQ=255)
                aCTGTAGCCAGCGCCGGAATATAATCCGCGACTTTACCCTGACCAATGAGCGGCCGCACTTGCCGCAAGATGTTTTCTAAAATTGCATTATCCATGGCGACTGCCACTTTCTACTCCTGGACCGCAGGTCTTTTCAGACCTGCggtccag                                                                                                                                       <  2:486125/150‑7 (MQ=255)
                aCTGTAGCCAGCGCCGGAATATAATCCGCGACTTTACCCTGACCAATGAGCGGCCGCACTTGCCGCAAGATGTTTTCTAAAATTGCATTATCCATGGCGACTGCCACTTTCTACTCCTGGACCGCAGGTCTGAAAAGACCTGCGAGtata                                                                                                                                       >  2:641605/1‑150 (MQ=255)
                 cTGTAGCCAGCGCCGGAATATAATCCGCGACTTTACCCTGACCAATGAGCGGCCGCACTTGCCGCAAGATGTTTTCTAAAATTGCATTATCCATGGCGACTGCCACTTTCTACTCCTGGACCGCAGGTCTGAAAAGACCTGCGAGTatat                                                                                                                                      >  1:1734179/1‑150 (MQ=255)
                   gTAGCCAGCGCCGGAATATAATCCGCGACTTTACCCTGACCAATGAGCGGCCGCACTTGCCGCAAGATGTTTTCTAAAATTGCATTATCCATGGCGACTGCCACTTTCTACTCCTGGACCGCAGGTCTGAAAAGACCTGCGAGTATATCa                                                                                                                                    >  2:1451464/1‑150 (MQ=255)
                    tAGCCAGCGCCGGAATATAATCCGCGACTTTACCCTGACCAATGAGCGGCCGCACTTGCCGCAAGATGTTTTCTAAAATTGCATTATCCATGGCGACTGCCACTTTCTACTCCTGGACCGCAGGTCTGAAAAGACCTGCGAGTATATCag                                                                                                                                   <  1:929058/150‑1 (MQ=255)
                     aGCCAGCGCCGGAATATAATCCGCGACTTTACCCTGACCAATGAGCGGCCGCACTTGCCGCAAGATGTTTTCTAAAATTGAATTATCAATGGCGACTGCCACTTTCTACTCCTGGACCGCAGGTCTGAAAAGACCTGCGAGTATATCaga                                                                                                                                  >  2:1548671/1‑150 (MQ=255)
                      gCCAGCGCCGGAATATAATCCGCGACTTTACCCTGACCAATGAGCGGCCGCACTTGCCGCAAGATGTTTTCTAAAATTGCATTATCCATGGCGACTGCCACTTTCTACTCCTGGACCGCAGGTCTGAAAAGACCTGCGAGTATATCagag                                                                                                                                 <  2:1955602/150‑1 (MQ=255)
                      gCCAGCGCCGGAATATAATCCGCGACTTTACCCTGACCAATGAGCGGCCGCACTTGCCGCAAGATGTTTTCTAAAATTGCATTATCCATGGCGACTGCCACTTTCTACTCCTGGACCGCAGGTCTGAAAAGACCTGCGAGTATATCagag                                                                                                                                 <  1:638531/150‑1 (MQ=255)
                        cAGCGCCGGAATATAATCCGCGACTTTACCCTGACCAATGAGCGGCCGCACTTGCCGCAAGATGTTTTCTAAAATTGCATTATCCATGGCGACTGCCACTTTCTACTCCTGGACCGCAGGTCTGAAAAGACCTGCGAGTATATCAGAGCt                                                                                                                               <  1:1451290/150‑1 (MQ=255)
                             ccGGAATATAATCCGCGACTTTACCCTGACCAATGAGCGGCCGCACTTGCCGCAAGATGTTTTCTAAAATTGCATTATCCATGGCGACTGCCACTTTCTACTCCTGGACCGCAGGTCTGAAAAGACCTGCGAGTATATCAGAGCTGAata                                                                                                                          >  1:988138/1‑150 (MQ=255)
                             ccGGAATATAATCCGCGACTTTACCCTGACCAATGAGCGGCCGCACTTGCCGCAAGATGTTTTCTAAAATTGCATTATCCATGGCGACTGCCACTTTCTACTCCTGGACCGCAGGTCTGAAAAGACCTGCGAGTATATCAGAGCTGAata                                                                                                                          <  1:2266165/150‑1 (MQ=255)
                                gAATATAATCCGCGACTTTACCCTGACCAATGAGCGGCCGCACTTGCCGCAAGATGTTTTCTAAAATTGCATTATCCATGGCGACTGCCACTTTCTACTCCTGGACCGCAGGTCTGAAAAGACCTGCGAGTATATCAGAGCTGAATATGt                                                                                                                       <  2:812723/150‑1 (MQ=255)
                                 aaTATAATCCGCGACTTTACCCTGACCAATGAGCGGCCGCACTTGCCGCAAGATGTTTTCTAAAATTGCATTATCCATGGCGACTGCCACTTTCTACTCCTGGACCGCAGGTCTGAAAAGACCTGCGAGTATATCAGAGCTGAATATGTc                                                                                                                      >  1:2056314/1‑150 (MQ=255)
                                         ccGCGACTTTACCCTGACCAATGAGCGGCCGCACTTGCCGCAAGATGTTTTCTAAAATTGCATTATCCATGGCGACTGCCACTTTCTACTCCTGGACCGCAGGTCTGAAAAGACCTGCGAGTATATCAGAGCTGAATATGTCGCGTCAGa                                                                                                              >  2:709683/1‑150 (MQ=255)
                                           gcgACTTTACCCTGACCAATGAGCGGCCGCACTTGCCGCAAGATGTTTTCTAAAATTGCATTATCCATGGCGACTGCCACTTTCTACTCCTGGACCGCAGGTCTGAAAAGACCTGCGAGTATATCAGAGCTGAATATGTCGCGTCAGATc                                                                                                            <  1:342891/150‑1 (MQ=255)
                                                tttACCCTGACCAATGAGCGGCCGCACTTGCCGCAAGATGTTTTCTAAAATTGCATTATCCATGGCGACTGCCACTTTCTACTCCTGGACCGCAGGTCTGAAAAGACCTGCGAGTATATCAGAGCTGAATATGTCGCGTCAGATCCGGTc                                                                                                       <  1:509223/150‑1 (MQ=255)
                                                tttACCCTGACCAATGAGCGGCCGCACTTGCCGCAAGATGTTTTCTAAAATTGCATTATCCATGGCGACTGCCACTTTCTACTCCTGGACCGCAGGTCTGAAAAGACCTGCGAGTATATCAGAGCTGAATATGTCGCGTCAGATCCGGTc                                                                                                       <  2:2056314/150‑1 (MQ=255)
                                                tttACCCTGACCAATGAGCGGCCGCACTTGCCGCAAGATGTTTTCTAAAATTGCATTATCCATGGCGACTGCCACTTTCTACTCCTGGACCGCAGGTCTGAAAAGACCTGCGAGTATATCAGAGCTGAATATGTCGCGTCAGATCCGGTc                                                                                                       <  2:1718268/150‑1 (MQ=255)
                                                 ttACCCTGACCAATGAGCGGCCGCACTTGCCGCAAGATGTTTTCTAAAATTGCATTATCCATGGCGACTGCCACTTTCTACTCCTGGACCGCAGGTCTGAAAAGACCTGCGAGTATATCAGAGCTGAATATGTCGCGTCAGATCCGGTCt                                                                                                      >  1:2031363/1‑150 (MQ=255)
                                                   aCCCTGACCAATGAGCGGCCGCACTTGCCGCAAGATGTTTTCTAAAATTGCATTATCCATGGCGACTGCCACTTTCTACTCCTGGACCGCAGGTCTGAAAAGACCTGCGAGTATATCAGAGCTGAATATGTCGCGTCAGATCCGGTCttt                                                                                                    >  1:2307891/1‑150 (MQ=255)
                                                    cccTGACCAATGAGCGGCCGCACTTGCCGCAAGATGTTTTCTAAAATTGCATTATCCATGGCGACTGCCACTTTCTACTCCTGGACCGCAGGTCTGAAAAGACCTGCGAGTATATCAGCGCTGAATATGTCGCGTCAGATCCGGTCTTTc                                                                                                   >  2:121741/1‑150 (MQ=255)
                                                     ccTGACCAATGAGCGGCCGCACTTGCCGCAAGATGTTTTCTAAAATTGCATTATCCATGGCGACTGCCACTTTCTACTCCTGGACCGCAGGTCTGAAAAGACCTGCGAGTATATCAGAGCTGAATATGTCGCGTCAGATCCGGTCTTTcc                                                                                                  <  2:1348944/150‑1 (MQ=255)
                                                      cTGACCAATGAGCGGCCGCACTTGCCGCAAGATGTTTTCTAAAATTGCATTATCCATGGCGACTGCCACTTTCTACTCCTGGACCGCAGGTCTGAAAAGACCTGCGAGTATATCAGAGCTGAATATGTCGCGTCAGATCCGGTCTTTCca                                                                                                 >  1:2203505/1‑150 (MQ=255)
                                                         aCCAATGAGCGGCCGCACTTGCCGCAAGATGTTTTCTAAAATTGCATTATCCATGGCGACTGCCACTTTCTACTCCTGGACCGCAGGTCTGAAAAGACCTGCGAGTATATCAGAGCTGAATATGTCGCGTCAGATCCGGTCTTTCCacac                                                                                              <  2:1456148/150‑1 (MQ=255)
                                                            aaTGAGCGGCCGCACTTGCCGCAAGATGTTTTCTAAAATTGCATTATCCATGGCGACTGCCACTTTCTACTCCTGGACCGCAGGTCTGAAAAGACCTGCGAGTATATCAGAGCTGAATATGTCGCGTCAGATCCGGTCTTTCCACACCGt                                                                                           >  1:1943914/1‑150 (MQ=255)
                                                                       gCACTTGCCGCAAGATGTTTTCTAAAATTGCATTATCCATGGCGACTGCCACTTTCTACTCCTGGACCGCAGGTCTGAAAAGACCTGCGAGTATATCAGAGCTGAATATGTCGCGTCAGATCCGGTCTTTCCACACCGTCTGGATATTAc                                                                                <  1:413818/150‑1 (MQ=255)
                                                                       gCACTTGCCGCAAGATGTTTTCTAAAATTGCATTATCCATGGCGACTGCCACTTTCTACTCCTGGACCGCAGGTCTGAAAAGACCTGCGAGTATATCAGAGCTGAATATGTCGCGTCAGATCCGGTCTTTCCACACCGTCTGGATATTAc                                                                                <  1:398471/150‑1 (MQ=255)
                                                                         aCTTGCCGCAAGATGTTTTCTAAAATTGCATTATCCATGGCGACTGCCACTTTCTACTCCTGGACCGCAGGTCTGAAAAGACCTGCGAGTATATCAGAGCTGAATATGTCGCGTCAGATCCGGTCTTTCCACACCGTCTGGATATTACAg                                                                              <  2:1992073/150‑1 (MQ=255)
                                                                          cTTGCCGCAAGATGTTTTCTAAAATTGCATTATCCATGGCGACTGCCACTTTCTACTCCTGGACCGCAGGTCTGAAAAGACCTGCGAGTATATCAGAGCTGAATATGTCGCGTCAGATCCGGTCTTTCCACACCGTCTGGATATTACAGa                                                                             >  1:883020/1‑150 (MQ=255)
                                                                           ttGCCGCAAGATGTTTTCTAAAATTGCATTATCCATGGCGACTGCCACTTTCTACTCCTGGACCGCAGGTCTGAAAAGACCTGCGAGTATATCAGAGCTGAATATGTCGCGTCAGATCCGGTCTTTCCACACCGTCTGGATATTACAGaa                                                                            <  1:1883168/150‑1 (MQ=255)
                                                                            tGCCGCAAGATGTTTTCTAAAATTGCATTATCCATGGCGACTGCCACTTTCTACTCCTGGACCGCAGGTCTGAAAAGACCTGCGAGTATATCAGAGCTGAATATGTCGCGTCAGATCCGGTCTTTCCACACCGTCTGGATATTACAGAAt                                                                           <  2:244589/150‑1 (MQ=255)
                                                                              ccGCAAGATGTTTTCTAAAATTGCATTATCCATGGCGACTGCCACTTTCTACTCCTGGACCGCAGGTCTGAAAAGACCTGCGAGTATATCAGAGCTGAATATGTCGCGTCAGATCCGGTCTTTCCACACCGTCTGGATATTACAGAATTc                                                                         <  2:2307891/150‑1 (MQ=255)
                                                                              ccGCAAGATGTTTTCTAAAATTGCATTATCCATGGCGACTGCCACTTTCTACTCCTGGACCGCAGGTCTGAAAAGACCTGCGAGTATATCAGAGCTGAATATGTCGCGTCAGATCCGGTCTTTCCACACCGTCTGGATATTACAGAATTc                                                                         >  1:2231782/1‑150 (MQ=255)
                                                                                 cAATATGTTTTCTAAAATTGCATTATCCATGGAGACTGCCACTTTCTACTCCTGGACCGCTGGTCTGAAAAGACATGCGAGTATATCTGAGCTGAATATGTCGCGTCAGATCCGGTCTTTCCTCAACTTGTTGATATTACATAATTCgtg                                                                      <  2:1546017/150‑1 (MQ=255)
                                                                                 cAAGATGTTTTCTAAAATTGCATTATCCATGGCGACTGCCACTTTCTACTCCTGGACCGCAGGTCTGAAAAGACCTGCGAGTATATCAGAGCTGAATATGTCGCGTCAGATCCGGTCTTTCCACACCGTCTGGATATTACAGAATTCgtg                                                                      <  1:1780993/150‑1 (MQ=255)
                                                                                   aGATGTTTTCTAAAATTGCATTATCCATGGCGACTGCCACTTTCTACTCCTGGACCGCAGGTCTGAAAAGACCTGCGAGTATATCAGAGCTGAATATGTCGCGTCAGATCCGGTCTTTCCACACCGTCTGGATATTACAGAATTCGTGTa                                                                    >  2:1364336/1‑150 (MQ=255)
                                                                                       gTTTTCTAAAATTGCATTATCCATGGCGACTGCCACTTTCTACTCCTGGACCGCAGGTCTGAAAAGACCTGCGAGTATATCAGAGCTGAATATGTCGCGTCAGATCCGGTCTTTCCACACCGTCTGGATATTACAGAATTCGTGTAAGcc                                                                >  2:1666997/1‑150 (MQ=255)
                                                                                        ttttCTAAAATTGCATTAT‑‑ATGGCGACTGCCACTTTCTACTCCTGGACCGCAGGTCTGAAAAGACCTGCGAGTATATCAGAGCTGAATATGTCGCGTCAGATCCGGTCTTTCCACACCGTCTGGATATTACAGAATTCGTGTAAGCCGaa                                                             >  2:319455/1‑150 (MQ=255)
                                                                                             aaaaaTTGCATTATCCATGGCGACTGCCACTTTCTACTCCTGGACCGCAGGTCTGAAAAGACCTGCGAGTATATCAGAGCTGAATATGTCGCGTCAGATCCGGTCTTTCCACACCGTCTCGATATTACAGAATTCGTGTAAGCCGAAATg                                                          >  1:670726/2‑150 (MQ=255)
                                                                                              aaaaTTGCATTATCCATGGCGACTGCCACTTTCTACTCCTGGACCGCAGGTCTGAAAAGACCTGCGAGTATATCAGAGCTGAATATGTCGCGTCAGATCCGGTCTTTCCACACCGTCTGGATATTACAGAATTCGTGTAAGCCGAAAttg                                                         >  2:2103207/1‑148 (MQ=255)
                                                                                              aaaaTTGCATTATCCATGGCGACTGCCACTTTCTACTCCTGGACCGCAGGTCTGAAAAGACCTGCGAGTATATCAGAGCTGAATATGTCGCGTCAGATCCGGTCTTTCCACACCGTCTGGATATTACAGAATTCGTGTAAGCCGAAATgg                                                         <  2:1618740/150‑1 (MQ=255)
                                                                                              aaaaTTGCATTATCCATGGCGACTGCCACTTTCTACTCCTGGACCGCAGGTCTGAAAAGACCTGCGAGTATATCAGAGCTGAATATGTCGCGTCAGATCCGGTCTTTCCACACCGTCTGGATATTACAGAATTCGTGTAAGCCGAAATgg                                                         >  1:1629936/1‑150 (MQ=255)
                                                                                               aaaTTGCATTATCCATGGCGACTGCCACTTTCTACTCCTGGACCGCAGGTCTGAAAAGACCTGCGAGTATATCAGAGCTGAATATGTCGCGTCAGATCCGGTCTTTCCACACCGTCTGGATATTACAGAATTCGTGTAAGCCGAAATggg                                                        <  2:343422/150‑1 (MQ=255)
                                                                                                 aTTGCATTATCCATGGCGACTGCCACTTTCTACTCCTGGACCGCAGGTCTGAAAAGACCTGCGAGTATATCAGAGCTGAATATGTCGCGTCAGATCCGGTCTTTCCACACCGTCTGGATATTACAGAATTCGTGTAAGCCGAAATGGGaa                                                      <  1:1518718/150‑1 (MQ=255)
                                                                                                  ttGCATTATCCATGGCGACTGCCACTTTCTACTCCTGGACCGCAGGTCTGAAAAGACCTGCGAGTATATCAGAGCTGAATATGTCGCGTCAGATCCGGTCTTTCCACACCGTCTGGATATTACAGAATTCGTGTAAGCCGAAATGGGaaa                                                     <  2:597862/150‑1 (MQ=255)
                                                                                                   tGCATTATCCATGGCGACTGCCACTTTCTACTCCTGGACCGCAGGTCTGAAAAGACCTGCGAGTATATCAGAGCTGAATATGTCGCGTCAGATCCGGTCTTTCCACACCGTCTGGATATTACAGAATTCGTGTAAGCCGAAATGGGAAAg                                                    >  2:1425745/1‑150 (MQ=255)
                                                                                                    gCATTATCCATGGCGACTGCCACTTTCTACTCCTGGACCGCAGGTCTGAAAAGACCTGCGAGTATATCAGAGCTGAATATGTCGCGTCAGATCCGGTCTTTCCACACCGTCTGGATATTACAGAATTCGTGTAAGCCGAAATGGGAAAGc                                                   <  1:2368672/150‑1 (MQ=255)
                                                                                                       ttATCCATGGCGACTGCCACTTTCTACTCCTGGACCGCAGGTCTGAAAAGACCTGCGAGTATATCAGAGCTGAATATGTCGCGTCAGATCCGGTCTTTCCACACCGTCTGGATATTACAGAATTCGTGTAAGCCGAAATGGGAAAGCTCa                                                >  2:2307709/1‑150 (MQ=255)
                                                                                                        tATCCATGGCGACTGCCACTTTCTACTCCTGGACCGCAGGTCTGAAAAGACCTGCGAGTATATCAGAGCTGAATATGTCGCGTCAGATCCGGTCTTTCCACACCGTCTGGATATTACAGAATTCGTGTAAGCCGAAATGGGAAAGCTCAc                                               <  2:645719/150‑1 (MQ=255)
                                                                                                         aTCCATGGCGACTGCCACTTTCTACTCCTGGACCGCAGGTCTGAAAAGACCTGCGAGTATATCAGAGCTGAATATGTCGCGTCAGATCCGGTCTTTCCACACCGTCTGGATATTACAGAGTTCGTGTAAGCCGAAATGGGAAAGCTCACg                                              >  1:466216/1‑150 (MQ=255)
                                                                                                            cATGGCGACTGCCACTTTCTACTCCTGGACCGCAGGTCTGAAAAGACCTGCGAGTATATCAGAGCTGAATATGTCGCGTCAGATCCGGTCTTTCCACACCGTCTGGATATTACAGAATTCGTGTAAGCCGAAATGGGAAAGCTCACGAcc                                           <  2:742296/150‑1 (MQ=255)
                                                                                                            cATGGCGACTGCCACTTTCTACTCCTGGACCGCAGGTCTGAAAAGACCTGCGAGTATATCAGAGCTGAATATGTCGCGTCAGATCCGGTCTTTCCACACCGTCTGGATATTACAGAATTCGTGTAAGCCGAAATGGGAAAGCTCACGAcc                                           <  2:2231782/150‑1 (MQ=255)
                                                                                                               ggCGACTGCCACTTTCTACTCCTGGACCGCAGGTCTGAAAAGACCTGCGAGTATATCAGAGCTGAATATGTCGCGTCAGATCCGGTCTTTCCACACCGTCTGGATATTACAGAATTCGTGTAAGCCGAAATGGGAAAGCTCACGACCaaa                                        >  1:2165785/1‑150 (MQ=255)
                                                                                                                      gCCACTTTCTACTCCTGGACCGCAGGTCTGAAAAGACCTGCGAGTATATCAGAGCTGAATATGTCGCGTCAGATCCGGTCTTTCCACACCGTCTGGATATTACAGAATTCGTGTAAGCCGAAATGGGAAAGCTCACGACCAAAGCCActc                                 >  1:582374/1‑150 (MQ=255)
                                                                                                                      gCCACTTTCTACTCCTGGACCGCAGGTCTGAAAAGACCTGCGAGTATATCAGAGCTGAATATGTCGCGTCAGATCCGGTCTTTCCACACCGTCTGGATATTACAGAATTCGTGTAAGCCGAAATGGGAAAGCTCACGACCAAAGCCActc                                 >  1:2481184/1‑150 (MQ=255)
                                                                                                                       ccACTTTCTACTCCTGGACCGCAGGTCTGAAAAGACCTGCGAGTATATCAGAGCTGAATATGTCGCGTCAGATCCGGTCTTTCCACACCGTCTGGATATTACAGAATTCGTGTAAGCCGAAATGGGAAAGCTCACGCCCAAAGCCActct                                >  1:2255771/1‑150 (MQ=255)
                                                                                                                       ccACTTTCTACTCCTGGACCGCAGGTCTGAAAAGACCTGCGAGTATATCAGAGCTGAATATGTCGCGTCAGATCCGGTCTTTCCACACCGTCTGGATATTACAGAATTCGTGTAAGCCGAAATGGGAAAGCTCACGACCAAAGCCActct                                >  2:2231351/1‑150 (MQ=255)
                                                                                                                          cTTTCTACTCCTGGACCGCAGGTCTGAAAAGACCTGCGAGTATATCAGAGCTGAATATGTCGCGTCAGATCCGGTCTTTCCACACCGTCTGGATATTACAGAATTCGTGTAAGCCGAAATGGGAAAGCTCACGACCAAAGCCACTCtttt                             <  2:470468/150‑1 (MQ=255)
                                                                                                                           tttCTACTCCTGGACCGCAGGTCTGAAAAGACCTGCGAGTATATCAGAGCTGAATATGTCGCGTCAGATCCGGTCTTTCCACACCGTCTGGATATTACAGAATTCGTGTAAGCCGAAATGGGAAAGCTCACGACCAAAGCCACTCttttt                            >  2:954860/1‑150 (MQ=255)
                                                                                                                           tttCTACTCCTGGACCGCAGGTCTGAAAAGACCTGCGAGTATATCAGAGCTGAATATGTCGCGTCAGATCCGGTCTTTCCACACCGTCTGGATATTACAGAATTCGTGTAAGCCGAAATGGGAAAGCTCACGACCAAAGCCACTCttttt                            <  1:709683/150‑1 (MQ=255)
                                                                                                                            ttCTACTCCTGGACCGCAGGTCTGAAAAGACCTGCGAGTATATCAGAGCTGAATATGTCGCGTCAGATCCGGTCTTTCCACACCGTCTGGATATTACAGAATTCGTGTAAGCCGAAATGGGAAAGCTCACGACCAAAGCCACTCTTTTTc                           <  2:2193256/150‑1 (MQ=255)
                                                                                                                            ttCTACTCCTGGACCGCAGGTCTGAAAAGACCTGCGAGTATATCAGAGCTGAATATGTCGCGTCAGATCCGGTCTTTCCACACCGTCTGGATATTACAGAATTCGTGTAAGCCGAAATGGGAAAGCTCACGACCAAAGCCACTCTTTTTc                           <  2:1922907/150‑1 (MQ=255)
                                                                                                                                 cTCCTGGACCGCAGGTCTGAAAAGACCTGCGAGTATATCAGAGCTGAATATGTCGCGTCAGATCCGGTCTTTCCACACCGTCTGGATATTACAGAATTCGTGTAAGCCGAAATGGGAAAGCTCACGACCAAAGCCACTCTTTTTCACGcc                      <  1:1096927/150‑1 (MQ=255)
                                                                                                                                 cTCCTGGACCGCAGGTCTGAAAAGACCTGCGAGTATATCAGAGCTGAATATGTCGCGTCAGATCCGGTCTTTCCACACCGTCTGGATATTACAGAATTCGTGTAAGCCGAAATGGGAAAGCTCACGACCAAAGCCACTCTTTTTCACGcc                      <  1:759280/150‑1 (MQ=255)
                                                                                                                                   ccTGGACCGCAGGTCTGAAAAGACCTGCGAGTATATCAGAGCTGAATATGTCGCGTCAGATCCGGTCTTTCCACACCGTCTGGATATTACAGAATTCGTGTAAGCCGAAATGGGAAAGCTCACGACCAAAGCCACTCTTTTTCACGccac                    <  2:2203505/150‑1 (MQ=255)
                                                                                                                                    tatactcGCAGGTCTTTTCAGACCTGCGAGTATATCAGAGCTGAATATGTCGCGTCAGATCCGGTCTTTCCACACCGTCTGGATATTACAGAATTCGTGTAAGCCGAAATGGGAAAGCTCACGACCAAAGCCACTCTTTTTCACGccacc                   >  2:953355/7‑150 (MQ=255)
                                                                                                                                    cTGGACCGCAGGTCTGAAAAGACCTGCGAGTATATCAGAGCTGAATATGTCGCGTCAGATCCGGTCTTTCCACACCGTCTGGATATTACAGAATTCGTGTAAGCCGAAATGGGAAAGCTCACGACCAAAGCCACTCTTTTTCACGccacc                   >  1:752640/1‑150 (MQ=255)
                                                                                                                                        ctcGCAGGTCTTTTCAGACCTGCGAGTATATCAGAGCTGAATATGTCGCGTCAGATCCGGTCTTTCCACACCGTCTGGATATTACAGAATTCGTGTAAGCCGAAATGGGAAAGATCACGACCAAAGCCACTCTTTTTCACGCCACCAAAg               >  2:2123759/3‑150 (MQ=255)
                                                                                                                                        aCCGCAGGTCTGAAAAGACCTGCGAGTATATCAGAGCTGAATATGTCGCGTCAGATCCGGTCTTTCCACACCGTCTGGATATTACAGAATTCGTGTAAGCCGAAATGGGAAAGCTCACGACCAAAGCCACTCTTTTTCACGCCACCAAAg               >  1:588897/1‑150 (MQ=255)
                                                                                                                                         ccGCAGGTCTGAAAAGACCTGCGAGTATATCAGAGCTGAATATGTCGCGTCAGATCCGGTCTTTCCACACCGTCTGGATATTACAGAATTCGTGTAAGCCGAAATGGGAAAGCTCACGACCAAAGCCACTCTTTTTCACGCCACCAAAgg              >  1:767048/1‑150 (MQ=255)
                                                                                                                                         ccGCAGGTCTGAAAAGACCTGCGAGTATATCAGAGCTGAATATGTCGCGTCAGATCCGGTCTTTCCACACCGTCTGGATATTACAGAATTCGTGTAAGCCGAAATGGGAAAGATCACGACCAAAGCCACTCTTTTTCACGCCACCAAAgg              >  1:1785574/1‑150 (MQ=255)
                                                                                                                                          cGCAGGTCTTTTCAGACCTGCGAGTATATCAGAGCTGAATATGTCGCGTCAGATCCGGTCTTTCCACACCGTCTGGATATTACAGAATTCGTGTAAGCCGAAATGGGAAAGCTCACGACCAAAGCCACTCTTTTTCACGCCACCAAAGGc             >  1:965961/1‑150 (MQ=255)
                                                                                                                                          cGCAGGTCTGAAAAGACCTGCGAGTATATCAGAGCTGAATATGTCGCGTCAGATCCGGTCTTTCCACACCGTCTGGATATTACAGAATTCATGTAAGCCGAAATGGGAAAGCTCACGACCAAAGCCACTCTTTTTCACGCCACCAAAGGc             >  2:718210/1‑150 (MQ=255)
                                                                                                                                           gCAGGTCTGAAAAGACCTGCGAGTATATCAGAGCTGAATATGTCGCGTCAGATCCGGTCTTTCCACACCGTCTGGATATTACAGAATTCGTGTAAGCCGAAATGGGAAAGCTCACGACCAAAGCCACTCTTTTTCACGCCACCAAAGGcc            <  2:1194696/150‑1 (MQ=255)
                                                                                                                                           gCAGGTCTGAAAAGACCTGCGAGTATATCAGAGCTGAATATGTCGCGTCAGATCCGGTCTTTCCACACCGTCTGGATATTACAGAATTCGTGTAAGCCGAAATGGGAAAGCTCACGACCAAAGCCACTCTTTTTCACGCCACCAAAGGcc            >  1:404707/1‑150 (MQ=255)
                                                                                                                                              ggTCTGAAAAGACCTGCGAGTATATCAGAGCTGAATATGTCGCGTCAGATCCGGTCTTTCCACACCGTCTGGATATTACAGAATTCGTGTAAGCCGAAATGGGAAAGCTCACGACCAAAGCCACTCTTTTTCACGCCACCAAAGGCCACt         <  1:2425260/150‑1 (MQ=255)
                                                                                                                                              ggTCTCAAAAGACCTGCAACTATATCAGAGCTGAATATCTCGCGTCAGAACCCGTCTATCCACACCATCTGGATATTACAGAATTCGTGTAAGCCGAAATGGGAAAGCTCACGACCAAAGCCACTCTTTTTCACGCCACCAAAGGCCACt         >  2:530495/1‑150 (MQ=255)
                                                                                                                                               gTCTGAAAAGACCTGCGAGTATATCAGAGCTGAATATGTCGCGTCAGATCCTGTCTTTCCACACCGTCTGGATATTACAGAATTCGTGTAAGCCGAAATGGGAAAGCTCACGACCAAAGCCACTCTTTTTCACGCCACCAAAGGCCACTc        >  2:544454/1‑150 (MQ=255)
                                                                                                                                                tCTGAAAAGACCTGCGAGTATATCAGAGCTGAATATGTCGCGTCAGATCCGGTCTTTCCACACCGTCTGGATATTACAGAATTCGTGTAAGCCGAAATGGGAAAGCTCACGACCAAAGCCACTCTTTTTCACGCCACCAAAGGCCACTcg       <  2:2407188/150‑1 (MQ=255)
                                                                                                                                                tCTGAAAAGACCTGCGAGTATATCAGAGCTGAATATGTCGCGTCAGATCCGGTCTTTCCACACCGTCTGGATATTACAGAATTCGTGTAAGCCGAAATGGGAAAGCTCACGACCAAAGCCACTCTTTTTCACGCCACCAAAGGCCACTcg       >  2:1592411/1‑150 (MQ=255)
                                                                                                                                                 cTGAAAAGACCTGCGAGTATATCAGAGCTGAATATGTCGCGTCAGATCCGGTCTTTCCACACCGTCTGGATATTACAGAATTCGTGTAAGCCGAAATgg                                                         <  2:1642407/99‑1 (MQ=255)
                                                                                                                                                 cTGAAAAGACCTGCGAGTATATCAGAGCTGAATATGTCGCGTCAGATCCGGTCTTTCCACACCGTCTGGATATTACAGAATTCGTGTAAGCCGAAATgg                                                         >  1:1642407/1‑99 (MQ=255)
                                                                                                                                                  tGAAAAGACCTGCGAGTATATCAGAGCTGAATATGTCGCGTCAGATCCGGTCTTTCCACACCGTCTGGATATTACAGAATTCGTGTAAGCCGAAATGGGAAAGCTCACGACCAAAGCCACTCTTTTTCACGCCACCAAAGGCCACTcgcg     <  1:1862711/150‑1 (MQ=255)
                                                                                                                                                   gAAAAGACCTGCGAGTATATCAGAGCTGAATATGTCGCGTCAGATCCGGTCTTTCCACACCGTCTGGATATTACAGAATTCGTGTAAGCCGAAATGGGAAAGCTCACGACCAAAGCCACTCTTTTTCACGCCACCAAAGGCCACTcgcgc    >  1:880301/1‑150 (MQ=255)
                                                                                                                                                    aaaaGACCTGCGAGTATATCAGAGCTGAATATGTCGCGTCAGATCCGGTCTTTCCACACCGTCTGGATATTACAGAATTCGTGTAAGCCGAAATGGGAAAGCTCACGACCAAAGCCACTCTTTTTCACGCCACCAAAGGCCACTcgcgcg   <  1:794077/150‑1 (MQ=255)
                                                                                                                                                    aaaaGACCTGCGAGTATATCAGAGCTGAATATGTCGCGTCAGATCCGGTCTTTCCACACCGTCTGGATATTACAGAATTCGTGTAAGCCGAAATGGGAAAGCTCACGACCAAAGCCACTCTTTTTCACGCCACCAAAGGCCACTcgcgcg   >  2:1164947/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CAATCGGGAACCGTCTACTGTAGCCAGCGCCGGAATATAATCCGCGACTTTACCCTGACCAATGAGCGGCCGCACTTGCCGCAAGATGTTTTCTAAAATTGCATTATCCATGGCGACTGCCACTTTCTACTCCTGGACCGCAGGTCTGAAAAGACCTGCGAGTATATCAGAGCTGAATATGTCGCGTCAGATCCGGTCTTTCCACACCGTCTGGATATTACAGAATTCGTGTAAGCCGAAATGGGAAAGCTCACGACCAAAGCCACTCTTTTTCACGCCACCAAAGGCCACTCGCGCGT  >  NC_000913/1613129‑1613427

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