Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,597,272 T→A 7.5% pseudogene (815/5422 nt) yneO ← AIDA‑I family autotransporter YneO

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,597,2720TA 7.5% 470.7 / 20.3 158pseudogene (815/5422 nt)yneOAIDA‑I family autotransporter YneO
Reads supporting (aligned to +/- strand):  ref base T (74/72);  new base A (4/8);  total (78/80)
Fisher's exact test for biased strand distribution p-value = 3.69e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.91e-01

Only 100 of 159 total aligned reads displayed.

GAACGCCCACCAATGAATATTCTTTGTTATTAATCAGTCCATTATTGGAAATATTGACTGTGCCATTACCGAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTTTATTAGCCACTCGCCACCCTTTTCAACGACAACCTGGCC  >  NC_000913/1597123‑1597420
                                                                                                                                                     |                                                                                                                                                    
     cccACCAAAGAATATTCTTAGTTATTAATCACTCCATTATTGGAAATATTGACTGTGCCAATACCGAAAAATCCATTATAGAGTCAGCGTACAGTGATCACAGAGTCTAGATCCAGCACAATCAGGGTTCCGATACCGCTGGCATTGCCa                                                                                                                                                 >  2:19250/1‑150 (MQ=255)
       caccaATGAATATTGTTTGTTATTAATCAGTCCATTATTGGAAATATTGACTGTGCCATTACCGAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCAcc                                                                                                                                               <  2:2444135/150‑1 (MQ=255)
       caccaATGAATATTCTTTGTTATTAATCAGTCCATTATTGGAAATATTGACTGTGCCATTACCGAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCAcc                                                                                                                                               >  2:1237770/1‑150 (MQ=255)
       caccaATGAATATTCTTTGTTATTAATCAGTCCATTATTGGAAATATTGACTGTGCCATTACCGAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCAcc                                                                                                                                               <  2:2435804/150‑1 (MQ=255)
       caccaATGAATATTCTTTGTTATTAATCAGTCCATTATTGGAAATATTGACTGTGCCATTACCGAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCAcc                                                                                                                                               <  2:1704781/150‑1 (MQ=255)
        accaATGAATATTCTTTGTTATTAATCAGTCCATTATTGGAAATATTGACTGTGCCATTACCGAAATATCCATTATCGAGTCTGCGTACCGTGATGACAGAGGCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCg                                                                                                                                              <  1:2036927/150‑1 (MQ=255)
        accaATGAATATTCTTTGTTATTAATCAGTCCATTATTGGAAATATTGACTGTGCCATTACCGAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTTCCACCg                                                                                                                                              >  2:1784356/1‑150 (MQ=255)
        accaATGAATATTCTTTGTTATTAATCAGTCCATTATTGGAAATATTGACTGTGCCATTACCGAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCg                                                                                                                                              <  2:499269/150‑1 (MQ=255)
         ccaATGAATATTCTTTGTTATTAATCAGTCCATTATTGGAAATATTGACTGTGCCATTACCGAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCAATGCCACCgg                                                                                                                                             <  2:1784829/150‑2 (MQ=255)
          caATGAATATTCTTTGTTATTAATCAGTCCATTATTGGAAATATTGACTGTGCCATTACCGAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGAt                                                                                                                                            <  1:2064841/150‑1 (MQ=255)
           aaTGAATATTCTTTGTTATTAATCAGTCCATTATTGGAAATATTGACTGTGCCATTACCGAAATATCCATAATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATa                                                                                                                                           >  2:1651957/1‑150 (MQ=255)
            aTGAATATTCTTTGTTATTAATCAGTCCATTATTGGAAATATTGACTGTGCCATTACCGAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCAATGCCACCGGTat                                                                                                                                          <  2:821242/150‑2 (MQ=255)
             tGAATATTCTTTGTTATTAATCAGTCCATTATTGGAAATATTGACTGTGCCATTACCGAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCAATGCCACCGGTatc                                                                                                                                         <  2:1378518/150‑3 (MQ=255)
              gAATATTCTTTGTTCTAAATCAGGCCTAAATTGGAAATATTGACTGTGCCATTACCGAATAATCCATTATAGAGTCTCCGTACCGTGATGACAGAGGCTTGATCCTGCACTTTCAGGGTTCCGATACCGGTGGCAATGCCACCGGTatcg                                                                                                                                        <  2:1997498/150‑4 (MQ=255)
              gAATATTCTTTGTTATTAATCAGTCCATTATTGGAAATATTGACTGTGCCATTACCGAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATc                                                                                                                                        >  2:2396063/1‑150 (MQ=255)
               aaTATTCTTTGTTATTAATCAGTCCATTATTGGAAATATTGACTGTGCCATTACCGAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCg                                                                                                                                       <  2:190046/150‑1 (MQ=255)
               aaTATTCTTTGTTATTAATCAGTCCATTATTGGAAATATTGACTGTGCCATTACCGAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCg                                                                                                                                       <  2:351965/150‑1 (MQ=255)
               aaTATTCTTTGTTATTAATCAGTCCATTATTGGAAATATTGACTGTGCCATTACCGAAATATCAATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCg                                                                                                                                       >  1:760151/1‑150 (MQ=255)
                 tatTCTTTGTTATTAATCAGTCCATTATTGGAAATATTGACTGTGCCATTACCGAAATATCCATTATAGAGTCTTCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTa                                                                                                                                     <  1:1143384/150‑1 (MQ=255)
                 tatTCTTTGTTATTAATCAGTCCATTATTGGAAATATTGACTGTGCCATTACCGAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTa                                                                                                                                     <  1:556081/150‑1 (MQ=255)
                 tatTCTTTGTTATTAATCAGTCCATTATTGGAAATATTGACTGTGCCATTACCGAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTa                                                                                                                                     <  1:1549938/150‑1 (MQ=255)
                    tCTTTGTTATTAATCAGTCCATTATTGGAAATATTGACTGTGCCATTACCGAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTAttt                                                                                                                                  >  1:1826896/1‑150 (MQ=255)
                          ttattaATCAGTCCATTATTGGAAATATTGACTGTGCCATTACCGAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGcc                                                                                                                            >  1:856453/1‑150 (MQ=255)
                          ttattaATCAGTCCATTATTGGAAATATTGACTGTGCCATTACCGAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGcc                                                                                                                            >  2:992785/1‑150 (MQ=255)
                           tattaATCAGTCCATTATTGGAAATATTGACTGTGCCATTACCGAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCg                                                                                                                           >  1:2157083/1‑150 (MQ=255)
                           tattaATCAGTCCATTATTGGAAATATTGACTGTGCCATTACCGAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCg                                                                                                                           <  1:1288820/150‑1 (MQ=255)
                            attaATCAGTCCATTATTGGAAATATTGACTGTGCCATTACCGAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGt                                                                                                                          >  2:2110996/1‑150 (MQ=255)
                            attaATCAGTCCATTATTGGAAATATTGACTGTGCCATTACCGAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGt                                                                                                                          >  2:71729/1‑150 (MQ=255)
                                 tCAGTCCATTATTGGAAATATTGACTGTGCCATTACCGAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTaacca                                                                                                                     >  1:757723/1‑150 (MQ=255)
                                       cATTATTGGAAATATTGACTGTGCCATTACCGAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAAccccc                                                                                                               <  2:645233/150‑1 (MQ=255)
                                       cATTATTGGAAATATTGACTGTGCCATTACCGAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAAccccc                                                                                                               >  1:869805/1‑150 (MQ=255)
                                       cATTATTGGAAATATTGACTGTGCCATTACCGAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCATCCGTAACCAAcccccc                                                                                                              >  2:2085423/1‑150 (MQ=255)
                                             tGGAAATATTGACTGTGCCATTACCGAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCACCCCCCCCTcgc                                                                                                         >  1:309824/1‑150 (MQ=255)
                                             tGGAAATATTGACTGTGCCATTACCGAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCATGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGAAACCAACCCCCACTcgc                                                                                                         >  2:1488227/1‑150 (MQ=255)
                                              ggAAATATTGACTGTGCCATTACCGAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTcgcg                                                                                                        <  2:2178152/150‑1 (MQ=255)
                                              ggAAATATTGACTGTGCCATTACCGAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTcgcg                                                                                                        <  2:601397/150‑1 (MQ=255)
                                              ggAAATATTGACTGTGCCATTACCGAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTcgcg                                                                                                        <  1:2258026/150‑1 (MQ=255)
                                              ggAAATATTGACTGTGCCATTACCGAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAAACCCCCCTcgcg                                                                                                        >  1:2472766/1‑150 (MQ=255)
                                               gAAATATTGACTGTGCCATTACCGAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGa                                                                                                       <  1:1647121/150‑1 (MQ=255)
                                                 aaTATTGACTGTGCCATTACCGAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAAt                                                                                                     <  1:2210409/150‑1 (MQ=255)
                                                    atTGACTGTGCCATTACCGAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGt                                                                                                  >  2:2213245/1‑150 (MQ=255)
                                                      tGACTGTGCCATTACCGAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCg                                                                                                <  2:562850/150‑1 (MQ=255)
                                                         cTGTGCCATTACCGAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGccg                                                                                             >  1:1840254/1‑149 (MQ=255)
                                                          tgtgCCATTACCGAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTc                                                                                            <  1:29560/150‑1 (MQ=255)
                                                           gtgCCATTACCGAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTcc                                                                                           >  1:27957/1‑150 (MQ=255)
                                                           gtgCCATTACCGAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTcc                                                                                           >  1:1009868/1‑150 (MQ=255)
                                                            tgCCATTACCGAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTccc                                                                                          >  2:1025244/1‑150 (MQ=255)
                                                             gCCATTACCGAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCTTTGTCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTcccc                                                                                         <  2:1162877/150‑1 (MQ=255)
                                                             gCCATTACCGAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTcccc                                                                                         <  2:3186/150‑1 (MQ=255)
                                                             gCCATTACCGAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGCCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTcccc                                                                                         <  1:796069/150‑1 (MQ=255)
                                                               cATTACCGAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAg                                                                                       >  2:1917106/1‑150 (MQ=255)
                                                                   aCCGAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTcc                                                                                   <  1:540578/150‑1 (MQ=255)
                                                                      gAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTg                                                                                <  1:636369/150‑1 (MQ=255)
                                                                      gAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCGTg                                                                                >  2:747480/1‑150 (MQ=255)
                                                                      gAAATATCCATTATAGAGTCTGCGTACCGTGAGGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTg                                                                                <  1:2045285/150‑1 (MQ=255)
                                                                       aaaTATCCATTATAGAGTCTGATTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAAGCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGa                                                                               <  2:463020/150‑1 (MQ=255)
                                                                             ccATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCa                                                                         >  2:2121428/1‑150 (MQ=255)
                                                                                 tataGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAAttt                                                                     >  1:1681146/1‑150 (MQ=255)
                                                                                 tataGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAAttt                                                                     >  1:1214066/1‑150 (MQ=255)
                                                                                 tataGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGCTGTCATTGACACCGATAATCGTATTTTCAGCCGTAACCAACCCCACCTCGCGAAAAGTCGCCTCCCCAGTTCCTTGCTTACCAAttt                                                                     >  2:644294/1‑150 (MQ=255)
                                                                                 tataGAGTCTGCGTACCGTGATGACAGAGTCTTGAACCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAAttt                                                                     <  1:380384/150‑1 (MQ=255)
                                                                                    agagTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGaa                                                                  <  1:1975326/150‑1 (MQ=255)
                                                                                      tgTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACTTTCAGGGTTCCTATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAAt                                                                <  2:730094/149‑1 (MQ=255)
                                                                                          tGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCaa                                                            >  2:330685/1‑150 (MQ=255)
                                                                                           gCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGCATAGTCGGCTCCCCAGTTCCTTGATTACCAATTTGAAATTGAAt                                                           >  1:2279158/1‑150 (MQ=255)
                                                                                           gCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAAt                                                           <  1:2121428/150‑1 (MQ=255)
                                                                                             gTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTg                                                         <  1:2396063/150‑1 (MQ=255)
                                                                                                ccGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCGCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAgg                                                      >  2:1043761/1‑150 (MQ=255)
                                                                                                             gTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCACCCCCCCCTCGCGAATAGTCGCGTCCCCAGTTCCTTGATTTCCAATTTGAAATTCAATTGAGGAAGCAtttttttt                                         >  1:2362475/1‑150 (MQ=255)
                                                                                                              tCTTGATCATGCACATTCAGGGTTCCGATACCGGTGGCTTTGCCACCGATAATCGTATTTTCAGCCGTAACAAACCCCCCCTCGCGAATAGTGGCGTCCCCAGTTCGTTGATAACAAATTTGAAATTCTATTGAGTAATCATtttttttt                                        >  2:1465147/1‑150 (MQ=255)
                                                                                                                ttGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCCACCCCCCCTCGCGACTAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATGTTttat                                      >  1:2222858/1‑150 (MQ=255)
                                                                                                                ttGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTttat                                      <  2:2011204/150‑1 (MQ=255)
                                                                                                                ttGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTttat                                      <  2:1991145/150‑1 (MQ=255)
                                                                                                                 tGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTttatt                                     >  2:1001569/1‑150 (MQ=255)
                                                                                                                 tGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTttatt                                     >  2:449047/1‑150 (MQ=255)
                                                                                                                   aTCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTTTATTAg                                   <  2:1052821/150‑1 (MQ=255)
                                                                                                                    tCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTTTATTAGc                                  >  1:155589/1‑150 (MQ=255)
                                                                                                                     ccTTCACATTCAGGGGTCCGATAACGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACAAACACCCCCTCGCTAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTGAATTGAGGAATCATTATTTTTTATTAgga                                 <  1:779238/150‑3 (MQ=255)
                                                                                                                     ccTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTTTATTAGcc                                 <  1:635072/150‑1 (MQ=255)
                                                                                                                     ccTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTTTATTAGcc                                 <  2:2157083/150‑1 (MQ=255)
                                                                                                                     ccTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTTTATTAGcc                                 >  2:1769261/1‑150 (MQ=255)
                                                                                                                          acaTTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTTTATTAGCCACTCg                            <  2:757723/150‑1 (MQ=255)
                                                                                                                           caTTCAGGGCTCCGATACCGGTGGCATTGCCACCGATACTCGTATTTTCAGCCGTAACCACCCCCCCCTCGCGAATAGTCGCCTCCCCCGTTCGTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTTTATTGGCCACTCGc                           >  2:1638885/1‑150 (MQ=255)
                                                                                                                            aTTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCCCCCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCGTGATTACCAATTTGAAATTCAATTGAGGACTCATTATTTTTTATTAGCCACTCGcc                          >  1:1625097/1‑150 (MQ=255)
                                                                                                                            aTTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTTTTTTTTATTAGCCACTCGcc                          >  2:1758648/1‑150 (MQ=255)
                                                                                                                            aTTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTTTATTAGCCACTCGcc                          <  1:2321448/150‑1 (MQ=255)
                                                                                                                               cAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTTTATTAGCCACTCGCCAcc                       <  1:92314/150‑1 (MQ=255)
                                                                                                                                 gggTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAAcccc                                                                                                                <  2:2405313/59‑1 (MQ=255)
                                                                                                                                 gggTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAAcccc                                                                                                                >  1:2405313/1‑59 (MQ=255)
                                                                                                                                 gggTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTTTATTAGCCACTCGCCACCCt                     >  2:458653/1‑150 (MQ=255)
                                                                                                                                  ggTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTTTATTAGCCACTCGCCACCCtt                    >  2:429027/1‑150 (MQ=255)
                                                                                                                                     tCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGAGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTTTATTAGcc                                 >  2:1108380/1‑134 (MQ=255)
                                                                                                                                     tCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGAGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTTTATTAGcc                                 <  1:1108380/134‑1 (MQ=255)
                                                                                                                                        gATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTTTATTAGCCACTCGCCACCCTTTTCAAc              <  1:1237770/150‑1 (MQ=255)
                                                                                                                                        gATACCGGTGGCATTGCCACCGATAATCGTAATTTCAGCCGTAAAAAACCAAAACTATCGAAAAGTAGCATCCCCAGATCGTTGAATAAAAATTTTAAATACAAATGAGGAAAAATTATTTTTTATTAGCCACTCGCCACCCTTTTAAAc              >  2:1274866/1‑150 (MQ=255)
                                                                                                                                         aTACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTTTATTAGCCACTCGCCACCCTTTTCAACg             <  1:103853/150‑1 (MQ=255)
                                                                                                                                               gTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTTTATTAGCCACTCGCCACCCTTTTCAACGACAAcc       <  1:1917106/150‑1 (MQ=255)
                                                                                                                                                  gCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTTTATTAGCCACTCGCCACCCTTTTCAACGACAACCTgg    >  2:1433658/1‑150 (MQ=255)
                                                                                                                                                  gCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTTTATTAGCCACTCGCCACCCTTTTCAACGACAACCTgg    >  1:644392/1‑150 (MQ=255)
                                                                                                                                                    aTTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCACCCTCGCGAATAGTCGCCTCCCCACTTCCTTGATTACCAAATTGAAATTAAATTGAGGAATCATTATTTTTTATTAGCCAATAGCCACCCTTTTCAACGACAACCAGGcc  >  2:2136427/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
GAACGCCCACCAATGAATATTCTTTGTTATTAATCAGTCCATTATTGGAAATATTGACTGTGCCATTACCGAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTTTATTAGCCACTCGCCACCCTTTTCAACGACAACCTGGCC  >  NC_000913/1597123‑1597420

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