Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,613,277 A→T 9.4% intergenic (‑26/+38) 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,2770AT 9.4% 387.8 / 21.9 131intergenic (‑26/+38)glsB/sadglutaminase 2/NAD(+)‑dependent succinate semialdehyde dehydrogenase
Reads supporting (aligned to +/- strand):  ref base A (57/61);  new base T (3/9);  total (61/70)
Fisher's exact test for biased strand distribution p-value = 1.42e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.01e-01

Only 100 of 137 total aligned reads displayed.

CAATCGGGAACCGTCTACTGTAGCCAGCGCCGGAATATAATCCGCGACTTTACCCTGACCAATGAGCGGCCGCACTTGCCGCAAGATGTTTTCTAAAATTGCATTATCCATGGCGACTGCCACTTTCTACTCCTGGACCGCAGGTCTGAAAAGACCTGCGAGTATATCAGAGCTGAATATGTCGCGTCAGATCCGGTCTTTCCACACCGTCTGGATATTACAGAATTCGTGTAAGCCGAAATGGGAAAGCTCACGACCAAAGCCACTCTTTTTCACGCCACCAAAGGCCACTCGCGCG  >  NC_000913/1613129‑1613426
                                                                                                                                                    |                                                                                                                                                     
  aTCGGGAACCGTCTACTGTAGCCAGCGCCGGAATATAATCCGCGACTTTACCCTGACCAATGAGCGGCCGCACTTGCCGCAAGATGTTTTCTAAAATTGCATTATCCATGGCGACTGCCACTTTCTACTCCTGGACCGCAGGTCTGaaaa                                                                                                                                                    >  1:645719/1‑150 (MQ=255)
   tCGGGAACCGTCTACTGTAGCCAGCGCCGGAATATAATCCGCGACTTTACCCTGACCAATGAGCGGCCGCACTTGCCGCAAGATGTTTTCTAAAATTGCATTATCCATGGCGACTGCCACTTTCTACTCCTGGACCGCAGGTCTGAAAAg                                                                                                                                                   <  2:883579/150‑1 (MQ=255)
    cGGGAACCGTCTACACTAGCCAGCCCCGGAATATAATCCGCGACTTTACCCTGACCAATGAGCGGCCGCACTTGCCGCAAGATGTTTTCTAAAATTGCATTATCCATGGCGACTGCCACTTTCTACTCCTGGACCGCAGGTCTGAAAAGa                                                                                                                                                  >  2:398471/1‑150 (MQ=255)
     gggAACCGTCTACTGTAGCCAGCGCCGGAATATAATCCGCGACTTTACCCTGACCAATGAGCGGCCGCACTTGCCGCAAGATGTTTTCTAAAATTGCATTATCCATGGCGACTGCCACTTTCTACTCCTGGACCGCAGGTCTGGAAAGAc                                                                                                                                                 >  1:1104665/1‑150 (MQ=255)
      ggAACCGTCTACTGTAGCCAGCGCCGGAATATAATCCGCGACTTTACCCTGACCAATGAGCGGCCGCACTTGCCGCAAGATGTTTTCTAAAATTGCATTATCCATGGCGACTGCCACTTTCTACTCCTGGACCGCAGGTCTGAAAAGAcc                                                                                                                                                <  2:539534/150‑1 (MQ=255)
      ggAACCGTCTACTGTAGCCAGCGCCGGAATATAATCCGCGACTTTACCCTGACCAATGAGCGGCCGCACTTGCCGCAAGATGTTTTCTAAAATTGCATTATCCATGGCGACTGCCACTTTCTACTCCTGGACCGCAGGTCTGAAAAGAcc                                                                                                                                                >  2:17848/1‑150 (MQ=255)
       gAACCGTCTACTGTAGCCAGCGCCGGAATATAATCCGCGACTTTACCCTGACCAATGAGCGGCCGCACTTGCCGCAAGATGTTTTCTAAAATTGCATTATCCATGGCGACTGCCACTTTCTACTCCTGGACCGCAGGTCTGAAAAGACCt                                                                                                                                               <  2:2371562/150‑1 (MQ=255)
        aaCCGTCTACTGTAGCCAGCGCCGGAATATAATCCGCGACTTTACCCTGACCAATGAGCGGCCGCACTTGCCGCAAGATGTTTTCTAAAATTGCATTATCCATGGCGACTGCCACTTTCTACTCCTGGACCGCAGGTCTGAAAAGACCTg                                                                                                                                              <  1:1119569/150‑1 (MQ=255)
         aCCGTCTACTGTAGCCAGCGCCGGAATATAATCCGCGACTTTACCCTGACCAATGAGCGGCCGCACTTGCCGCAAGATGTTTTCTAAAATTGCATTATCCATGGCGACTGCCACTTTCTACTCCTGGACCGCAGGTCTTTTCAGACCTGc                                                                                                                                             <  1:1369422/150‑1 (MQ=255)
         aCCGTCTACTGTAGCCAGCGCCGGAATATAATCCGCGACTTTACCCTGACCAATGAGCGGCCGCACTTGCCGCAAGATGTTTTCTAAAATTGCATTATCCATGGCGACTGCCACTTTCTACTCCTGGACCGCAGGTCTTTTCAGACCTGc                                                                                                                                             <  2:1506275/150‑1 (MQ=255)
         aCCGTCTACTGTAGCCAGCGCCGGAATATAATCCGCGACTTTACCCTGACCAATGAGCGGCCGCACTTGCCGCAAGATGTTTTCTAAAATTGCATTATCCATGGCGACTGCCACTTTCTACTCCTGGACCGCAGGTCTTTTCAGACCTGc                                                                                                                                             <  2:1447183/150‑1 (MQ=255)
          ccGTCTACTGTAGCCAGCGCCGGAATATAATCCGCGACTTTACCCTGACCAATGAGCGGCCGCACTTGCCGCAAGATGTTTTCTAAAATTGCATTATCCATGGCGACTGCCACTTTCTACTCCTGGACCGCAGGTCTGAAAAGACCTGCg                                                                                                                                            <  1:228645/150‑1 (MQ=255)
          ccGTCTACTGTAGCCAGCGCCGGAATATAATCCGCGACTTTACCCTGACCAATGAGCGGCCGCACTTGCCGCAAGATGTTTTCTAAAATTGCATTATCCATGGCGACTGCCACTTTCTACTCCTGGACCGCAGGTCTGAAAAGACCTGCg                                                                                                                                            <  1:2061027/150‑1 (MQ=255)
           cGTCTACTGTAGCCAGCGCCGGAATATAATCCGCGACTTTACCCTGACCAATGAGCGGCCGCACTTGCCGCAAGATGTTTTCTAAAATTGCATTATCCATGGCGACTGCCACTTTCTACTCCTGGACCGCAGGTCTTTTCAGACCTGCgg                                                                                                                                           <  1:567751/150‑2 (MQ=255)
             tCTACTGTAGCCAGCGCCGGAATATAATCCGCGACTTTACCCTGACCAATGAGCGGCCGCACTTGCCGCAAGATGTTTTCTAAAATTGCATTATCCATGGCGACTGCCACTTTCTACTCCTGGACCGCAGGTCTTTTCAGACCTGCggtc                                                                                                                                         <  1:466263/150‑4 (MQ=255)
               tACTGTAGCCAGCGCCGGAATATAATCCGCGACTTTACCCTGACCAATGAGCGGCCGCACTTGCCGCAAGATGTTTTCTAAAATTGCATTATCCATGGCGACTGCCACTTTCTACTCCTGGACCGCAGGTCTGAAAAGACCTGCGAGtat                                                                                                                                       >  2:1053722/1‑150 (MQ=255)
                aCTGTAGCCAGCGCCGGAATATAATCCGCGACTTTACCCTGACCAATGAGCGGCCGCACTTGCCGCAAGATGTTTTCTAAAATTGCATTATCCATGGCGACTGCCACTTTCTACTCCTGGACCGCAGGTCTTTTCAGACCTGCggtccag                                                                                                                                      <  2:486125/150‑7 (MQ=255)
                aCTGTAGCCAGCGCCGGAATATAATCCGCGACTTTACCCTGACCAATGAGCGGCCGCACTTGCCGCAAGATGTTTTCTAAAATTGCATTATCCATGGCGACTGCCACTTTCTACTCCTGGACCGCAGGTCTTTTCAGACCTGCggtccag                                                                                                                                      <  1:2309245/150‑7 (MQ=255)
                aCTGTAGCCAGCGCCGGAATATAATCCGCGACTTTACCCTGACCAATGAGCGGCCGCACTTGCCGCAAGATGTTTTCTAAAATTGCATTATCCATGGCGACTGCCACTTTCTACTCCTGGACCGCAGGTCTGAAAAGACCTGCGAGtata                                                                                                                                      >  2:641605/1‑150 (MQ=255)
                    tAGCCAGCGCCGGAATATAATCCGCGACTTTACCCTGACCAATGAGCGGCCGCACTTGCCGCAAGATGTTTTCTAAAATTGCATTATCCATGGCGACTGCCACTTTCTACTCCTGGACCGCAGGTCTGAAAAGACCTGCGAGTATATCag                                                                                                                                  <  1:929058/150‑1 (MQ=255)
                     aGCCAGCGCCGGAATATAATCCGCGACTTTACCCTGACCAATGAGCGGCCGCACTTGCCGCAAGATGTTTTCTAAAATTGCATTATCCATGGCGACTGCCACTTTCTACTCCTGGACCGCAGGTCTGAAAAGACCTGCGAGTATATCaga                                                                                                                                 <  2:1617548/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)
                                          cgcgACTTTACCCTGACCAATGAGCGGCCGCACTTGCCGCAAGATGTTTTCTAAAATTGCATTATCCATGGCGACTGCCACTTTCTACTCCTGGACCGCAGGTCTGAAAAGACCTGCGAGTATATCAGAGCTGAATATGTCGCGTCAGAt                                                                                                            >  1:343422/1‑150 (MQ=255)
                                           gcgACTTTACCCTGACCAATGAGCGGCCGCACTTGCCGCAAGATGTTTTCTAAAATTGCATTATCCATGGCGACTGCCACTTTCTACTCCTGGACCGCAGGTCTGAAAAGACCTGCGAGTATATCAGAGCTGAATATGTCGCGTCAGATc                                                                                                           <  1:342891/150‑1 (MQ=255)
                                            cgACTTTACCCTGACCAATGAGCGGCCGCACTTGCCGCAAGATGTTTTCTAAAATTGCATTATCCATGGCGACTGCCACTTTCTACTCCTGGACCGCAGGTCTGAAAAGACCTGCGAGTATATCAGAGCTGAATATGTCGCGTCAGATcc                                                                                                          <  2:2285283/150‑1 (MQ=255)
                                            cgACTTTACCCTGACCAATGAGCGGCCGCAATTGCCGCAAGATGTTTTCTAAAATTGCATTATCCATGGCGACTGCCACTTTCTACTCCTGGACCGCAGGTCTGAAAAGACCTGCGAGTATATCAGAGCTGAATATGTCGCGTCAGATcc                                                                                                          <  1:1730165/150‑1 (MQ=255)
                                                tttACCCTGACCAATGAGCGGCCGCACTTGCCGCAAGATGTTTTCTAAAATTGCATTATCCATGGCGACTGCCACTTTCTACTCCTGGACCGCAGGTCTGAAAAGACCTGCGAGTATATCAGAGCTGAATATGTCGCGTCAGATCCGGTc                                                                                                      <  2:2056314/150‑1 (MQ=255)
                                                tttACCCTGACCAATGAGCGGCCGCACTTGCCGCAAGATGTTTTCTAAAATTGCATTATCCATGGCGACTGCCACTTTCTACTCCTGGACCGCAGGTCTGAAAAGACCTGCGAGTATATCAGAGCTGAATATGTCGCGTCAGATCCGGTc                                                                                                      <  1:509223/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)
                                                               gAGCGGCCGCACTTGCCGCAAGATGTTTTCTAAAATTGCATTATCCATGGCGACTGCCACTTTCTACTCCTGGACCGCAGGTCTGAAAAGACCTGCGAGTATATCAGAGCTGAATATGTCGCGTCAGATCCGGTCTTTCCACACCGTCTg                                                                                       >  1:593098/1‑150 (MQ=255)
                                                                       gCACTTGCCGCAAGATGTTTTCTAAAATTGCATTATCCATGGCGACTGCCACTTTCTACTCCTGGACCGCAGGTCTGAAAAGACCTGCGAGTATATCAGAGCTGAATATGTCGCGTCAGATCCGGTCTTTCCACACCGTCTGGATATTAc                                                                               <  1:413818/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)
                                                                             gCCGCAAGATGTTTTCTAAAATTGCATTATCCATGGCGACTGCCACTTTCTACTCCTGGACCGCAGGTCTGAAAAGACCTGCGAGTATATCAGAGCTGAATATGTCGCGTCAGATCCGGTCTTTCCACACCGTCTGGATATTACAGAAtt                                                                         >  1:702857/1‑150 (MQ=255)
                                                                              ccGCAAGATGTTTTCTAAAATTGCATTATCCATGGCGACTGCCACTTTCTACTCCTGGACCGCAGGTCTGAAAAGACCTGCGAGTATATCAGAGCTGAATATGTCGCGTCAGATCCGGTCTTTCCACACCGTCTGGATATTACAGAATTc                                                                        >  1:2231782/1‑150 (MQ=255)
                                                                              ccGCAAGATGTTTTCTAAAATTGCATTATCCATGGCGACTGCCACTTTCTACTCCTGGACCGCAGGTCTGAAAAGACCTGCGAGTATATCAGAGCTGAATATGTCGCGTCAGATCCGGTCTTTCCACACCGTCTGGATATTACAGAATTc                                                                        <  2:2307891/150‑1 (MQ=255)
                                                                                 cAATATGTTTTCTAAAATTGCATTATCCATGGAGACTGCCACTTTCTACTCCTGGACCGCTGGTCTGAAAAGACATGCGAGTATATCTGAGCTGAATATGTCGCGTCAGATCCGGTCTTTCCTCAACTTGTTGATATTACATAATTCgtg                                                                     <  2:1546017/150‑1 (MQ=255)
                                                                                 cAAGATGTTTTCTAAAATTGCATTATCCATGGCGACTGCCACTTTCTACTCCTGGACCGCAGGTCTGAAAAGACCTGCGAGTATATCAGAGCTGAATATGTCGCGTCAGATCCGGTCTTTCCACACCGTCTGGATATTACAGAATTCgtg                                                                     <  1:1780993/150‑1 (MQ=255)
                                                                                   tgATGTTTTCTAAAATTGCATTATCCATGGCGACTGCCACTTTGTACTCCTGGAGCGCAGGTCTGAAAAGACCTGCGAGTATATCAGAGCTGAATATGTCGCGTCAGATCCGGTCTTTCCACACCGTCTGGATATTACAGAATTCGTGTa                                                                   <  2:1701189/149‑1 (MQ=255)
                                                                                       gTTTTCTAAAATTGCATTATCCATGGCGACTGCCACTTTCTACTCCTGGACCGCAGGTCTGAAAAGACCTGCGAGTATATCAGAGCTGAATATGTCGCGTCAGATCCGGTCTTTCCACACCGTCTGGATATTACAGAATTCGTGTAAGcc                                                               >  2:1666997/1‑150 (MQ=255)
                                                                                        ttttCTAAAATTGCATTAT‑‑ATGGCGACTGCCACTTTCTACTCCTGGACCGCAGGTCTGAAAAGACCTGCGAGTATATCAGAGCTGAATATGTCGCGTCAGATCCGGTCTTTCCACACCGTCTGGATATTACAGAATTCGTGTAAGCCGaa                                                            >  2:319455/1‑150 (MQ=255)
                                                                                             aaaaaTTGCATTATCCATGGCGACTGCCACTTTCTACTCCTGGACCGCAGGTCTGAAAAGACCTGCGAGTATATCAGAGCTGAATATGTCGCGTCAGATCCGGTCTTTCCACACCGTCTCGATATTACAGAATTCGTGTAAGCCGAAATg                                                         >  1:670726/2‑150 (MQ=255)
                                                                                              aaaaTTGCATTATCCATGGCGACTGCCACTTTCTACTCCTGGACCGCAGGTCTGAAAAGACCTGCGAGTATATCAGAGCTGAATATGTCGCGTCAGATCCGGTCTTTCCACACCGTCTGGATATTACAGAATTCGTGTAAGCCGAAATgg                                                        <  2:1618740/150‑1 (MQ=255)
                                                                                              aaaaTTGCATTATCCATGGCGACTGCCACTTTCTACTCCTGGACCGCAGGTCTGAAAAGACCTGCGAGTATATCAGAGCTGAATATGTCGCGTCAGATCCGGTCTTTCCACACCGTCTGGATATTACAGAATTCGTGTAAGCCGAAATgg                                                        >  1:1629936/1‑150 (MQ=255)
                                                                                               aaaTTGCATTATCCATGGCGACTGCCACTTTCTACTCCTGGACCGCAGGTCTGAAAAGACCTGCGAGTATATCAGAGCTGAATATGTCGCGTCAGATCCGGTCTTTCCACACCGTCTGGATATTACAGAATTCGTGTAAGCCGAAATggg                                                       <  1:763666/150‑1 (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)
                                                                                                            cATGGCGACTGCCACTTTCTACTCCTGGACCGCAGGTCTGAAAAGACCTGCGAGTATATCAGAGCTGAATATGTCGCGTCAGATCCGGTCTTTCCACACCGTCTGGATATTACAGAATTCGTGTAAGCCGAAATGGGAAAGCTCACGAcc                                          <  2:742296/150‑1 (MQ=255)
                                                                                                            cATGGCGACTGCCACTTTCTACTCCTGGACCGCAGGTCTGAAAAGACCTGCGAGTATATCAGAGCTGAATATGTCGCGTCAGATCCGGTCTTTCCACACCGTCTGGATATTACAGAATTCGTGTAAGCCGAAATGGGAAAGCTCACGAcc                                          <  2:2231782/150‑1 (MQ=255)
                                                                                                             aTGGCGACTGCCACTTTCTACTCCTGGACCGCAGGTCTGAAAAGACCTGCGAGTATATCAGAGCTGAATATGTCGCGTCAGATCCGGTCTTTCCACACCGTCTGGATATTACAGAATTCGTGTAAGCCGAAATGGGAAAGCTCACGACCa                                         >  1:814849/1‑150 (MQ=255)
                                                                                                               ggCGACTGCCACTTTCTACTCCTGGACCGCAGGTCTGAAAAGACCTGCGAGTATATCAGAGCTGAATATGTCGCGTCAGATCCGGTCTTTCCACACCGTCTGGATATTACAGAATTCGTGTAAGCCGAAATGGGAAAGCTCACGACCaaa                                       >  1:2165785/1‑150 (MQ=255)
                                                                                                                  gACTGCCACTTTCTACTCCTGGACCGCAGGTCTTAAAAGACCTGCGAGTATATCAGAGCTTAATATGTCGCGTCAGATCCGGTCTTTGTACTCCGTCTGGTTATTACAGTATTGTTGTAAGCCTAAATGTTTTAGCTCACGTCCTAAGcc                                    <  2:2203586/150‑1 (MQ=255)
                                                                                                                       ccACTTTCTACTCCTGGACCGCAGGTCTGAAAAGACCTGCGAGTATATCAGAGCTGAATATGTCGCGTCAGATCCGGTCTTTCCACACCGTCTGGATATTACAGAATTCGTGTAAGCCGAAATGGGAAAGCTCACGCCCAAAGCCActct                               >  1:2255771/1‑150 (MQ=255)
                                                                                                                       ccACTTTCTACTCCTGGACCGCAGGTCTGAAAAGACCTGCGAGTATATCAGAGCTGAATATGTCGCGTCAGATCCGGTCTTTCCACACCGTCTGGATATTACAGAATTCGTGTAAGCCGAAATGGGAAAGCTCACGACCAAAGCCActct                               >  2:2231351/1‑150 (MQ=255)
                                                                                                                       ccACTTTCTACTCCTGGACCGCAGGTCTGAAAAGACCTGCGAGTATATCAGAGCTGAATATGTCGCGTCAGATCCGGTCTTTCCACACCGTCTGGATATTACAGAATTCGTGTAAGCCGAAATGGGAAAGCTCACGACCAAAGCCActct                               >  2:12339/1‑150 (MQ=255)
                                                                                                                          cTTTCTACTCCTGGACCGCAGGTCTGAAAAGACCTGCGAGTATATCAGAGCTGAATATGTCGCGTCAGATCCGGTCTTTCCACACCGTCTGGATATTACAGAATTCGTGTAAGCCGAAATGGGAAAGCTCACGACCAAAGCCACTCtttt                            <  2:470468/150‑1 (MQ=255)
                                                                                                                           tttCTACTCCTGGACCGCAGGTCTGAAAAGACCTGCGAGTATATCAGAGCTGAATATGTCGCGTCAGATCCGGTCTTTCCACACCGTCTGGATATTACAGAATTCGTGTAAGCCGAAATGGGAAAGCTCACGACCAAAGCCACTCttttt                           >  2:954860/1‑150 (MQ=255)
                                                                                                                            ttCTACTCCTGGACCGCAGGTCTGAAAAGACCTGCGAGTATATCAGAGCTGAATATGTCGCGTCAGATCCGGTCTTTCCACACCGTCTGGATATTACAGAATTCGTGTAAGCCGAAATGGGAAAGCTCACGACCAAAGCCACTCTTTTTc                          <  2:2193256/150‑1 (MQ=255)
                                                                                                                            ttCTACTCCTGGACCGCAGGTCTGAAAAGACCTGCGAGTATATCAGAGCTGAATATGTCGCGTCAGATCCGGTCTTTCCACACCGTCTGGATATTACAGAATTCGTGTAAGCCGAAATGGGAAAGCTCACGACCAAAGCCACTCTTTTTc                          <  2:1922907/150‑1 (MQ=255)
                                                                                                                              cTACTCCTGGACCGCAGGTCTGAAAAGACCTGCGAGTATATCAGAGCTGAATATGTCGCGTCAGCTCCGGTCTTTCCCCACCGTCTGGATATTACAGAAATCGTGTAAGCCGAAATGGGAAAGCTCACGCCCAAAGCCACTCGTTTTCAc                        >  1:2474398/1‑150 (MQ=255)
                                                                                                                                 cTCCTGGACCGCAGGTCTGAAAAGACCTGCGAGTATATCAGAGCTGAATATGTCGCGTCAGATCCGGTCTTTCCACACCGTCTGGATATTACAGAATTCGTGTAAGCCGAAATGGGAAAGCTCACGACCAAAGCCACTCTTTTTCACGcc                     <  1:759280/150‑1 (MQ=255)
                                                                                                                                 cTCCTGGACCGCAGGTCTGAAAAGACCTGCGAGTATATCAGAGCTGAATATGTCGCGTCAGATCCGGTCTTTCCACACCGTCTGGATATTACAGAATTCGTGTAAGCCGAAATGGGAAAGCTCACGACCAAAGCCACTCTTTTTCACGcc                     <  1:1096927/150‑1 (MQ=255)
                                                                                                                                   ccTGGACCGCAGGTCTGAAAAGACCTGCGAGTATATCAGAGCTGAATATGTCGCGTCAGATCCGGTCTTTCCACACCGTCTGGATATTACAGAATTCGTGTAAGCCGAAATGGGAAAGCTCACGACCAAAGCCACTCTTTTTCACGccac                   >  2:2011429/1‑150 (MQ=255)
                                                                                                                                   ccTGGACCGCAGGTCTGAAAAGACCTGCGAGTATATCAGAGCTGAATATGTCGCGTCAGATCCGGTCTTTCCACACCGTCTGGATATTACAGAATTCGTGTAAGCCGAAATGGGAAAGCTCACGACCAAAGCCACTCTTTTTCACGccac                   >  2:1881361/1‑150 (MQ=255)
                                                                                                                                   ccTGGACCGCAGGTCTGAAAAGACCTGCGAGTATATCAGAGCTGAATATGTCGCGTCAGATCCGGTCTTTCCACACCGTCTGGATATTACAGAATTCGTGTAAGCCGAAATGGGAAAGCTCACGACCAAAGCCACTCTTTTTCACGccac                   <  2:702857/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)
                                                                                                                                         ccGCAGGTCTGAAAAGACCTGCGAGTATATCAGAGCTGAATATGTCGCGTCAGATCCGGTCTTTCCACACCGTCTGGATATTACAGAATTCGTGTAAGCCGAAATGGGAAAGATCACGACCAAAGCCACTCTTTTTCACGCCACCAAAgg             >  1:1785574/1‑150 (MQ=255)
                                                                                                                                          cGCAGGTCTGAAAAGACCTGCGAGTATATCAGAGCTGAATATGTCGCGTCAGATCCGGTCTTTCCACACCGTCTGGATATTACAGAATTCGTGTAAGCCGAAATGGGAAAGCTCACGACCAAAGCCACTCTTTTTCACGCCACCAAAGGc            >  1:439517/1‑150 (MQ=255)
                                                                                                                                          cGCAGGTCTGAAAAGACCTGCGAGTATATCAGAGCTGAATATGTCGCGTCAGATCCGGTCTTTCCACACCGTCTGGATATTACAGAATTCATGTAAGCCGAAATGGGAAAGCTCACGACCAAAGCCACTCTTTTTCACGCCACCAAAGGc            >  2:718210/1‑150 (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:1592411/1‑150 (MQ=255)
                                                                                                                                                tCTGAAAAGACCTGCGAGTATATCAGAGCTGAATATGTCGCGTCAGATCCGGTCTTTCCACACCGTCTGGATATTACAGAATTCGTGTAAGCCGAAATGGGAAAGCTCACGACCAAAGCCACTCTTTTTCACGCCACCAAAGGCCACTcg      <  2:2407188/150‑1 (MQ=255)
                                                                                                                                                 cTGAAAAGACCTGCGAGTATATCAGAGCTGAATATGTCGCGTCAGATCCGGTCTTTCCACACCGTCTGGATATTACAGAATTCGTGTAAGCCGAAATgg                                                        >  1:1642407/1‑99 (MQ=255)
                                                                                                                                                  tGAAAAGACCTGCGAGTATATCAGAGCTGAATATGTCGCGTCAGATCCGGTCTTTCCACACCGTCTGGATATTACAGAATTCGTGTAAGCCGAAATGGGAAAGCTCACGACCAAAGCCACTCTTTTTCACGCCACCAAAGGCCACTcgcg    <  1:1862711/150‑1 (MQ=255)
                                                                                                                                                   gAAAAGACCTGCGAGTATATCAGAGCTGAATATGTCGCGTCAGATCCGGTCTTTCCACACCGTCTGGATATTACAGAATTCGTGTAAGCCGAAATGGGAAAGCTCACGACCAAAGCCACTCTTTTTCACGCCACCAAAGGCCACTcgcgc   >  1:942304/1‑150 (MQ=255)
                                                                                                                                                    aaaaGACCTGCGAGTATATCAGAGCTGAATATGTCGCGTCAGATCCGGTCTTTCCACACCGTCTGGATATTACAGAATTCGTGTAAGCCGAAATGGGAAAGCTCACGACCAAAGCCACTCTTTTTCACGCCACCAAAGGCCACTcgcgcg  <  1:794077/150‑1 (MQ=255)
                                                                                                                                                    aaaaGACCTGCGAGTATATCAGAGCTGAATATGTCGCGTCAGATCCGGTCTTTCCACACCGTCTGGATATTACAGAATTCGTGTAAGCCGAAATGGGAAAGCTCACGACCAAAGCCACTCTTTTTCACGCCACCAAAGGCCACTcgcgcg  >  2:1164947/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
CAATCGGGAACCGTCTACTGTAGCCAGCGCCGGAATATAATCCGCGACTTTACCCTGACCAATGAGCGGCCGCACTTGCCGCAAGATGTTTTCTAAAATTGCATTATCCATGGCGACTGCCACTTTCTACTCCTGGACCGCAGGTCTGAAAAGACCTGCGAGTATATCAGAGCTGAATATGTCGCGTCAGATCCGGTCTTTCCACACCGTCTGGATATTACAGAATTCGTGTAAGCCGAAATGGGAAAGCTCACGACCAAAGCCACTCTTTTTCACGCCACCAAAGGCCACTCGCGCG  >  NC_000913/1613129‑1613426

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: