Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,002,265 T→A 5.4% F280L (TTT→TTA) ‡ corA → Ni(2(+))/Co(2(+))/Mg(2(+)) transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,002,2650TA 5.4% 564.8 / 13.6 175F280L (TTT→TTA) ‡corANi(2(+))/Co(2(+))/Mg(2(+)) transporter
Reads supporting (aligned to +/- strand):  ref base T (83/82);  new base A (4/6);  total (87/88)
Fisher's exact test for biased strand distribution p-value = 7.47e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.81e-01

Only 100 of 186 total aligned reads displayed.

GAATCCCTGTTCCAGAAGGTGAACTTCCTGATGCAGGCGGCAATGGGTTTTATCAACATCGAGCAGAACCGCATCATCAAAATCTTCTCGGTGGTATCCGTGGTATTCCTGCCGCCGACGCTCGTTGCTTCCAGCTATGGCATGAACTTTGAGTTTATGCCAGAACTGAAGTGGAGCTTCGGCTACCCTGGCGCGATTATCTTTATGATCCTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAGTCTGGCTAAGCCACCT  >  NC_000913/4002116‑4002413
                                                                                                                                                     |                                                                                                                                                    
gAATCCCTGTTCCAGAAGGTGAACTTCCTGATGCAGGCGGCAATGGTTTTTATCAACATCGAGCAGAACCGCATCATCAAAATCTTCTCGGTGGTATCCGTGGTATTCCTGCCGCCGACGCTCGTTGCTTCCAGCTATGGCATGAACttt                                                                                                                                                      <  1:928694/150‑1 (MQ=255)
 aaTCCCTGTTCCAGAAGGTGAACTTCCTGATGCAGGCGGCAATGGGTTTTATCAACATCGAGCAGAACCGCATCATCAAAATCTTCTCGGTGGTATCCGTGGTATTCCTGCCGCCGACGCTCGTTGCTTCCAGCTATGGCATGAACTTTg                                                                                                                                                     <  2:2234737/150‑1 (MQ=255)
   tCCCTGTTCCAGAAGGTGAACTTCCTGATGCAGGCGGCAATGGGTTTTATCAACATCGAGCAGAACCGCATCATCAAAATCTTCTCGGTGGTATCCGTGGTATTCCTGCCGCCGACGCTCGTTGCTTCCAGCTATGGCATGAACTTTGAg                                                                                                                                                   <  1:1030299/150‑1 (MQ=255)
   tCCCTGTTCCAGAAGGTGAACTTCCTGATGCAGGCGGCAATGGGTTTTATCAACATCGAGCAGAACCGCATCATCAAAATCTTCTCGGTGGTATCCGTGGTATTCCTGCCGCCGACGCTCGTTGCTTCCAGCTATGGCATGAACTTTGAg                                                                                                                                                   <  2:1182551/150‑1 (MQ=255)
   tCCCTGTTCCAGAAGGTGAACTTCCTGATGCAGGCGGCAATGGGTTTTATCAACATCGAGCAGAACCGCATCATCAAAATCTTCTCGGTGGTATCCGTGGTATTCCTGCCGCCGACGCTCGTTGCTTCCAGCTATGGCATGAACTTTGAg                                                                                                                                                   <  1:1738168/150‑1 (MQ=255)
   tCCCTGTTCCAGAAGGTGAACTTCCTGATGCAGGCGGCAATGGGTTTTATCAACATCGAGCAGAACCGCATCATCAAAATCTTCTCGGTGGTATCCGTGGTATTCCTGCCGCCGACGCTCGTTGCTTCCAGCTATGGCATGAACTTTGAg                                                                                                                                                   <  1:129410/150‑1 (MQ=255)
    cccTGTTCCAGAAGGTTAACTTCCTGATGCAGGCGGCAATGGGTTTTATCAACATCGAGCAGAACCGCATCATCAAAATCTTCTCGGTGGTATCCGTGGTATTCCTGCCGCCGACGCTCGTTGCTTCCAGCTATGGCATGAACTTTGAGt                                                                                                                                                  <  2:1289244/150‑1 (MQ=255)
      cTGTTCCAGAAGGTGAACTTCCTGATGCAGGCGGCAATGGGTTTTATCAACATCGAGCAGAACCGCATCATCAAAATCTTCTCGGTGGTATCCGTGGTATTCCTGCCGCCGACGCTCGTTGCTTCCAGCTATGGCATGAACTTTGAGttt                                                                                                                                                >  2:1888789/1‑150 (MQ=255)
       tGTTCCAGAAGGTGAACTTCCTGATGCAGGCGGCAATGGGTTTTATCAACATCGAGCAGAACCGCATCATCAAAATCTTCTCGGTGGTATCCGTGGTATTCCTGCCGCCGACGCTCGTTGCTTCCAGCTATGGCATGAACTTTGAGt                                                                                                                                                  <  2:460237/147‑1 (MQ=255)
       tGTTCCAGAAGGTGAACTTCCTGATGCAGGCGGCAATGGGTTTTATCAACATCGAGCAGAACCGCATCATCAAAATCTTCTCGGTGGTATCCGTGGTATTCCTGCCGCCGACGCTCGTTGCTTCCAGCTATGGCATGAACTTTGAGt                                                                                                                                                  >  1:460237/1‑147 (MQ=255)
         ttCCAGAAGGTGAACTTCCTGATGCAGGCGGCAATGGGTTTTATCAACATCGAGCAGAACCGCATCATCAAAATCTTCTCGGTGGTATCCGTGGTATTCCTGCCGCCGACGCTCGTTGCTTCCAGCTATGGCATGAACTTTGAGTTTATg                                                                                                                                             <  2:534142/150‑1 (MQ=255)
           ccAGAAGGTGAACTTCCTGATGCAGGCGGCAATGGGTTTTATCAACATCGAGCAGAACCGCATCATCAAAATCTTCTCGGTGGTATCCGTGGTATTCCTGCCGCCGACGCTCGTTGCTTCCAGCTATGGCATGAACTTTGAGTTTATGcc                                                                                                                                           >  2:1109711/1‑150 (MQ=255)
            cAGAAGGTGAACTTCCTGATGCAGGCGGCAATGGGTTTTATCAACATCGAGCAGAACCGCATCATCAAAATCTTCTCGGTGGTATCCGTGGTATTCCTGCCGCCGACGCTCGTTGCTTCCAGCTATGGCATAAACTCAAAGTTCATGCCa                                                                                                                                          <  1:1034038/150‑1 (MQ=255)
                aGGTGAACTTCCTGATGCAGGCGGCAATGGGTTTTATCAACATCGAGCAGAACCGCATCATCAAAATCTTCTCGGTGGTATCCGTGGTATTCCTGCCGCCGACGCTCGTTGCTTCCAGCTATGGCATGAACTTTGAGTTTATGCCAGAAc                                                                                                                                      >  2:2231229/1‑150 (MQ=255)
                 ggTGAACTTCCTGATGCAGGCGGCAATGGGTTTTATCAACATCGAGCAGAACCGCATCATCAAAATCTTCTCGGTGGTATCCGTGGTATTCCTGCCGCCGACGCTCGTTGCTTCCAGCTATGGCATGAACTTTGAGTTTATGCCAGAACt                                                                                                                                     >  1:279423/1‑150 (MQ=255)
                  gTGAACTTCCTGATGCAGGCGGCAATGGGTTATATCAACATCGAGCAGAACCGCATCAGCAAAATCTTCTCGGGGGTATCCGTGGTATTCCTGCCGCCGACGCTCGTTGCTTCCAGCTATGGCATAAACTCAAAGTTCATGCCATAGCTg                                                                                                                                    <  2:1403935/150‑1 (MQ=255)
                   tGAACTTCCTGATGCAGGCGGCAATGGGTTTTATCAACATCGAGCAGAACCGCATCATCAAAATCTTCTCGGTGGTATCCGTGGTATTCCTGCCGCCGACGCTCGTTGCTTCCAGCTATGGCATGAACTTTGAGTTTATGCCAGAACTGa                                                                                                                                   >  1:257144/1‑150 (MQ=255)
                   ggAACTTCCTGATGCAGGCGGCAATGGGTTTTATCAACATCGAGCAGAACCGCATCATCAAAATCTTCTCGGTGGTATCCGTTTTTTTCCTGCCGCCGACGCTCGTTGCTTCCAGCTATGGCATAAACTCATAGTTCATGCCATAGCTgg                                                                                                                                   <  2:1005977/149‑2 (MQ=255)
                      aCTTCCTGATGCAGGCGGCAATGGGTTTTATCAACATCGAGCAGAACCGCATCATCAAAATCTTCTCGGTGGTATCCGTGGTATTCCTGCCGCCGACGCTCGTTGCTTCCAGCTATGGCATGAACTTTGAGTTTATGCCAGAACTGAAGt                                                                                                                                >  1:2305020/1‑150 (MQ=255)
                      aCTTCCTGATGCAGGCGGCAATGGGTTTTATCAACATCGAGCAGAACCGCATCATCAAAATCTTCTCGGTGGTATCCGTGGTATTCCTGCCGCCGACGCTCGTTGCTTCCAGCTATGGCATGAACTTTGAGTTTATGCCAGAACTGAAGt                                                                                                                                <  1:111616/150‑1 (MQ=255)
                      aCTTCCTGATGCAGGCGGCAATGGGTTTTATCAACATCGAGCAGAACCGCATCATCAAAATCTTCTCGGTGGTATCCGTGGTATTCCTGCCGCCGACGCTCGTTGCTTCCAGCTATGGCATGAACTTTGAGTTTATGCCAGAACTGAAGt                                                                                                                                <  2:2024121/150‑1 (MQ=255)
                       cTTCCTGATGCAGGCGGCAATGGGTTTTATCAACATCGAGCAGAACCGCATCATCAAAATCTTCTCGGTGGTATCCGTGGTATTCCTGCCGCCGACGCTCGTTGCTTCCAGCTATGGCATGAACTTTGAGTTTATGCCAGAACTGAAGTg                                                                                                                               >  1:1236997/1‑150 (MQ=255)
                           cTGATGCAGGCGGCAATGGGTTTTATCAACATCGAGCAGAACCGCATCATCAAAATCTTCTCGGTGGTATCCGTGGTATTCCTGCCGCCGACGCTCGTTGCTTCCAGCTATGGCATGAACTTTGAGTTTATGCCAGAACTGAAGTGGAGc                                                                                                                           <  1:1864604/150‑1 (MQ=255)
                            tGATGGAGGCGGCAATGGGTTTTATCAACATCGAGCAGAACCGCATCATCAAAATCTTCTCGGTGGTATCCGTGGTATTCCTGCCGCCGACGCTCGTTGCTTCCAGCTATGGCATGAACTTTGAGTTTATGCCAGAACTGAAGTGGAGCt                                                                                                                          <  2:229957/150‑1 (MQ=255)
                             gATGCAGGCGGCAATGGGTTTTATCAACATCGAGCAGAACCGCATCATCAAAATCTTCTCGGTGGTATCCGTGGTATTCCTCCCGCCGACGCTCGTTGCTTCCAGCTATGGCATGAACTTTGAGTTTATGCCAGAACTGAAGTGGAGCtt                                                                                                                         >  2:2010044/1‑150 (MQ=255)
                               tGCAGGCGGCAATGGGTTTTATCAACATCGAGCAGAACCGCATCATCAAAATCTTCTCGGTGGTATCCGTGGTATTCCTGCCGCCGACGCTCGTTGCTTCCAGCTATGGCATGAACTTTGAGTTTATGCCAGAACTGAAGTGGAGCTTCg                                                                                                                       >  1:1889853/1‑150 (MQ=255)
                                gCAGGCGGCAATGGGTTTTATCAACATCGAGCAGAACCGCATCATCAAAATCTTCTCGGTGGTATCCGTGGTATTCCTGCCGCCGACGCTCGTTGCTTCCAGCTATGGCATGAACTTTGAGTTTATGCCAGAACTGAAGTGGAGCTTCgg                                                                                                                      <  2:1784599/150‑1 (MQ=255)
                                 cAGGCGGCAATGGGTTTTATCAACATCGAGCAGAACCGCATCATCAAAATCTTCTCGGTGGTATCCGTGGTATTCCTGCCGCCGACGCTCGTTGCTTCCAGCTATGGCATGAACTTTGAGTTTATGCCAGAACTGAAGTGGAGCTTCGGc                                                                                                                     <  2:2099544/150‑1 (MQ=255)
                                      ggcAATGGGTTTTATCAACATCGAGCAGAACCGCATCATCAAAATCTTCTCGGTGGTATCCGTGGTATTCCTGCCGCCGACGCTCGTTGCTTCCAGCTATGGCATGAACTTTGAGTTTATGCCAGAACTGAAGTGGAGCTTCGGCTAccc                                                                                                                <  1:1851674/150‑1 (MQ=255)
                                      ggcAATGGGTTTTATCAACATCGAGCAGAACCGCATCATCAAAATCTTCTCGGTGGTATCCGTGGTATTCCTGCCGCCGACGCTCGTTGCTTCCAGCTATGGCATGAACTTTGAGTTTATGCCAGAACTGAAGTGGAGCTTCGGCTAccc                                                                                                                >  2:1737076/1‑150 (MQ=255)
                                         aaTGGGTTTTATCAACATCGAGCAGAACCGCATCATCAAAATCTTCTCGGTGGTATCCGTGGTATTCCTGCCGCCGACGCTCGTTGCTTCCAGCTATGGCATGAACTTTGAGTTTATGCCAGAACTGAAGTGGAGCTTCGGCTACCCTgg                                                                                                             <  1:735569/150‑1 (MQ=255)
                                           tGGGTTTTATCAACATCGAGCAGAACCGCATCATCAAAATCTTCTCGGTGGTATCCGTGGTATTCCTGCCGCCGACGCTCGTTGCTTCCAGCTATGGCATGAACTTTGAGTTTATGCCAGAACTGAAGTGGAGCTTCGGCTACCCTGgcg                                                                                                           <  1:846591/150‑1 (MQ=255)
                                             ggTTTTATCAACATCGAGCAGAACCGCATCATCAAAATCTTCTCGGTGGTATCCGTGGTATTCCTGCCGCCGACGCTCGTTGCTTCCAGCTATGGCATGAACTTTGAGTTTATGCCAGAACTGAAGTGGAGCTTCGGCTACCCTGGcgcg                                                                                                         >  2:182693/1‑150 (MQ=255)
                                               ttttATCAACATCGAGCAGAACCGCATCATCAAAATCTTCTCGGTGGTATCCGGGGTATTCCTGCCGCCGACGCTCGTTGCTTCCAGCTATGGCATGAACTTTGAGTTTATGCCAGAACTGAAGTGGAGCTTCGGCTACCCTGGCGCGAt                                                                                                       <  2:1500590/150‑1 (MQ=255)
                                                   aTCAACATCGAGCAGAACCGCATCATCAAAATCTTCTCGGTGGTATCCGTGGTATTCCTGCCGCCGACGCTCGTTGCTTCCAGCTATGGCATGAACTTTGAGTTTATGCCAGAACTGAAGTGGAGCTTCGGCTACCCTGGCGCGATTATc                                                                                                   <  1:1204417/150‑1 (MQ=255)
                                                    tCAACATCGAGCAGAACCGCATCATCAAAATCTTCTCGGTGGTATCCGTGGTATTCCTGCCGCCGACGCTCGTTGCTTCCAGCTATGGCATGAACTTTGAGTTTATGCCAGAACTGAAGTGGAGCTTCGGCTACCCTGGCGCGATTATCt                                                                                                  >  1:1941713/1‑150 (MQ=255)
                                                    tCAACATCGAGCAGAACCGCATCATCAAAATCTTCTCGGTGGTATCCGTGGTATTCCTGCCGCCGACGCTCGTTGCTTCCAGCTATGGCATGAACTTTGAGTTTATGCCAGAACTGAAGTGGAGCTTCGGCTACCCTGGCGCGATTATCt                                                                                                  >  1:1022606/1‑150 (MQ=255)
                                                     cAACATCGAGCAGAACCGCATCATCAAAAGCTTCTCGGTGGTATCCGTGGTATTCCTGCCGCCGACGCTCGTTGCTTCCAGCTATGGCATGAACTTTGAGTTTATGCCAGAACTGAAGTGGAGCTTCGGCTACCCTGGCGCGATTATCtt                                                                                                 <  2:1866964/150‑1 (MQ=255)
                                                      aaCATCGAGCAGAACCGCATCATCAAAATCTTCTCGGTGGTATCCGTGGTATTCCTGCCGCCGACGCTCGTTGCTTCCAGCTATGGCATGAACTTTGAGTTTATGCCAGAACTGAAGTGGAGCTTCGGCTACCCTGGCGCGATTATCttt                                                                                                >  1:905823/1‑150 (MQ=255)
                                                      aaCATCGAGCAGAACCGCATCATCAAAATCTTCTCGGTGGTATCCGTGGTATTCCTGCCGCCGACGCTCGTTGCTTCCAGCTATGGCATGAACTTTGAGTTTATGCCAGAACTGAAGTGGAGCTTCGGCTACCCTGGCGCGATTATCttt                                                                                                >  1:1047322/1‑150 (MQ=255)
                                                       aCATCGAGCAGAACCGCATCATCAAAATCTTCTCGGTGGTATCCGTGGTATTCCTGCCGCCGACGCTCGTTGCTTCCAGCTATGGCATGAACTTTGAGTTTATGCCAGAACTGAAGTGGAGCTTCGGCTACCCTGGCGCGATTATCTTTa                                                                                               >  1:1300427/1‑150 (MQ=255)
                                                           cGAGCAGAACCGCATCATCAAAATCTTCTCGGTGGTATCCGTGGTATTCCTGCCGCCGACGCTCGTTGCTTCCAGCTATGGCATGAACTTTGAGTTTATGCCAGAACTGAAGTGGAGCTTCGGCTACCCTGGCGCGATTATCTTTATGAt                                                                                           >  2:984759/1‑150 (MQ=255)
                                                             aGCAGAACCGCATCATCAAAATCTTCTCGGTGGTATCCGTGGTATTCCTGCCGCCGACGCTCGTTGCTTCCAGCTATGGCATGAACTTTGAGTTTATGCCAGAACTGAAGTGGAGCTTCGGCTACCCTGGCGCGATTATCTTTATGATcc                                                                                         <  1:1284419/150‑1 (MQ=255)
                                                             aGCAGAACCGCATCATCAAAATCTTCTCGGTGGTATCCGTGGTATTCCTGCCGCCCACGCTCGTTGCTTCCAGCTATGGCATGAACTTTGAGTTTATGCCAGAACTGAAGTGGAGCTTCGGCTACCCTGGCGGGATTATCTTTATGATcc                                                                                         >  2:809933/1‑150 (MQ=255)
                                                                aGAACCGCATCATCAAAATCTTCTCGGTGGTATCCGTGGTATTCCTGCCGCCGACGCTCGTTGCTTCCAGCTATGGCATGAACTTTGAGTTATGCCAGAACTGAAGTGGAGCTTCGGCTACCCTGGCGCGATTATCTTTATGATCCTcgc                                                                                     >  1:456286/1‑150 (MQ=255)
                                                                aGAACCGAATCATCAAAATCTTCTCGGTGGTATCCGTGGTATTCCTGCCGCCGACGCTCGTTGCTTCCAGCTATGGCATGAACTTTGAGTTTATGCCAGAACTGAAGTGGAGCTTCGGCTACCCTGGCGCGATTATCTTTATGATCCTcg                                                                                      <  2:411143/150‑1 (MQ=255)
                                                                   aCCGCATCATCAAAATCTTCTCGGTGGTATCCGTGGTATTCCTGCCGCCGACGCTCGTTGCTTCCAGCTATGGCATGAACTTTGAGTTTATGCCAGAACTGAAGTGGAGCTTCGGCTACCCTGGCGCGATTATCTTTATGATCCTCGCgg                                                                                   <  1:338938/150‑1 (MQ=255)
                                                                    ccGCATCATCAAAATCTTCTCGGTGGTATCCGTGGTATTCCTGCCGCCGACGCTCGTTGCTTCCAGCTATGGCATGAACTTTGAGTTTATGCCAGAACTGAAGTGGAGCTTCGGCTACCCTGGCGCGATTATCTTTATGATCCTCGCggg                                                                                  <  1:1613455/150‑1 (MQ=255)
                                                                      gCATCATCAAAATCTTCTCGGTGGTATCCGTGGTATTCCTGCCGCCGACGCTCGTTGCTTCCAGCTATGGCATGAACTTTGAGTTTATGCCAGAACTGAAGTGGAGCTTCGGCTACCCTGGCGCGATTATCTTTATGATCCTCGCGGGcc                                                                                <  1:2206090/150‑1 (MQ=255)
                                                                       catcatCAAAATCTTCTCGGTGGTATCCGTGGTATTCCTGCCGCCGACGCTCGTTGCTTCCAGCTATGGCATGAACTTTGAGTTTATGCCAGAACTGAAGTGGAGCTTCGGCTACCCTGGCGCGATTATCTTTATGATCCTCGCGGGCCt                                                                               <  1:370040/150‑1 (MQ=255)
                                                                            tcaAAATCTTCTCGGTGGTATCCGTGGTATTCCTGCCGCCGACGCTCGTTGCTTCCAGCTATGGCATGAACTTTGAGTTTATGCCAGAACTGAAGTGGAGCTTCGGCTACCCTGGCGCGATTATCTTTATGATCCTCGCGGGCCTGGCAc                                                                          <  1:1956685/150‑1 (MQ=255)
                                                                              aaaaTCTTCTCGGTGGTATCCGTGGTATTCCTGCCGCCGACGCTCGTTGCTTCCAGCTATGGCATGAACTTTGAGTTTATGCCAGAACTGAAGTGGAGCTTCGGCTACCCTGGCGCGATTATCTTTATGATCCTCGCGGGCCTGGCACCg                                                                        >  2:948298/1‑150 (MQ=255)
                                                                              aaaaTCTTCTCGGTGGTATCCGTGGTATTCCTGCCGCCGACGCTCGTTGCTTCCAGCTATGGCATGAACTTTGAGTTTATGCCAGAACTGAAGTGGAGCTTCGGCTACCCTGGCGCGATTATCTTTATGATCCTCGCGGGCCTGGCACCg                                                                        >  2:1278144/1‑150 (MQ=255)
                                                                              aaaaTCTTCTCGGTGGTATCCGTGGTATTCCTGCCGCCGACGCTCGTTGCTTCCAGCTATGGCATGAACTTTGAGTTTATGCCAGAACTGAAGTGGAGCTTCGGCTACCCTGGCGCGATTATCTTTATGATCCTCGCGGGCCTGGCACCg                                                                        <  1:1657554/150‑1 (MQ=255)
                                                                                aaTCTTCTCGTTTGTATCCGGGGTATTCCTGCCGCCGACGCTAGTTGCTTCCAGCTATGGCATGAACTTTGAGTTTTTGCCAGAACTGAAGTGGAGCTTCGGCTACCCTGGCGCGATTATCTTTATGATCCTCGCGGGCCTGGCACCGTa                                                                      <  2:1859261/150‑1 (MQ=255)
                                                                                aaTCTTCTCGGTGGTATCCGTGGTATTCCTGCCGCCGACACTCGTTGCTTCCAGCTATGGCATGAACTTTGAGTTTATGCCAGAACTGAAGTGGAGCTTCGGCTACCCTGGCGCGATTATCTTTATGATCCTCGCGGGCCTGGCACCGTa                                                                      <  2:326056/150‑1 (MQ=255)
                                                                                aaTCTTCTCGGTGGTATCCGCGGTATTCCTGCCGCCGACGCTCGTTGCTTCCAGCTATGGCATGAACTTTGAGTTTATGCCAGAACTGAAGTGGAGCTTCGGCTACCCTGGCGCGATTATCTTTATGATCCTCGCGGGCCTGGCACCGTa                                                                      <  2:2309248/150‑1 (MQ=255)
                                                                                  tcttctCGGTGGTCTCCGTGGGATTCCTGGCGCCGACGCTCGTTGCTTCCAGCTATGGCCTGAACTTTGAGTTTATGCCAGAACTGAAGTGGAGCTTCGGCTACCCTGGCGCGATTATCTTTATGATCCTCGCGGGCCTGGCACCATATc                                                                    <  2:689999/150‑1 (MQ=255)
                                                                                        cGGTGGTATCCGTGGTATTCCTGCCGCCGACGCTCGTTGCTTCCAGCTATGGCATGAACTTTGAGTTTATGCCAGAACTGAAGTGGAGCTTCGGCTACCCTGGCGCGATTATCTTTATGATCCTCGCGGGCCTGGCACCGTATCTGTACt                                                              <  2:2449274/150‑1 (MQ=255)
                                                                                         ggtggtATCCGTGGTATTCCTGCCGCCGACGCTCGTTGCTTCCAGCTATGGCATGAACTTTGAGTTTATGCCAGAACTGAAGTGGAGCTTCGGCTACCCTGGCGCGATTATCTTTATGATCCTCGCGGGCCTGGCACCGTATCTGTACtt                                                             >  1:1482060/1‑150 (MQ=255)
                                                                                         ggtggtATCCGTGGTATTCCTGCCGCCGACGCTCGTTGCTTCCAGCTATGGCATGAACTTTGAGTTTATGCCAGAACTGAAGTGGAGCTTCGGCTACCCTGGCGCGATTATCTTTATGATCCTCGCGGGCCTGGCACCGTATCTGTACtt                                                             >  1:1531423/1‑150 (MQ=255)
                                                                                               aTCCGTGGTATTCCTGCCGCCGACGCTCGTTGCTTCCAGCTATGGCATGAACTTTGAGTTTATGCCAGAACTGAAGTGGAGCTTCGGCTACCCTGGCGCGATTATCTTTATGATCCTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCg                                                       >  1:1696021/1‑150 (MQ=255)
                                                                                               aTCCGTGGTATTCCTGCCGCCGACGCTCGTTGCTTCCAGCTATGGCATGAACTTTGAGTTTATGCCAGAACTGAAGTGGAGCTTCGGCTACCCTGGCGCGATTATCTTTATGATCCTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCg                                                       <  1:779050/150‑1 (MQ=255)
                                                                                               aTCCGTGGTATTCCTGCCGCCGACGCTCGTTGCTTCCAGCTATGGCATGAACTTTGAGTTTATGCCAGAACTGAAGTGGAGCTTCGGCTACCCTGGCGCGATTATCTTTATGATCCTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCg                                                       >  1:1104404/1‑150 (MQ=255)
                                                                                                tCCGTGGTATTCCTGCCGCCGACGCTCGTTGCTTCCAGCTATGGCATGAACTTTGAGTTTATGCCAGAACTGAAGTGGAGCTTCGGCTACCCTGGCGCGATTATCTTTATGATCCTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCgg                                                      >  2:2271736/1‑150 (MQ=255)
                                                                                                 ccGTGGTATTCCTGCCGCCGACGCTCGTTGCTTCCAGCTATGGCATGAACTTTGAGTTTATGCCAGAACTGAAGTGGAGCTTCGGCTACCCTGGCGCGATTATCTTTATGATCCTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGga                                                     <  1:495295/150‑1 (MQ=255)
                                                                                                 ccGTGGTATTCCTGCCGCCGACGCTCGTTGCTTCCAGCTATGGCATGAACTTTGAGTTTATGCCAGAACTGAAGTGGAGCTTCGGCTACCCTGGCGCGATTATCTTTATGATCCTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGga                                                     >  2:437500/1‑150 (MQ=255)
                                                                                                   gTGGTATTCCTGCCGCCGACGCTCGTTGCTTCCAGCTATGGCATGAACTTTGAGTTTATGCCAGAACTGAAGTGGAGCTTCGGCTACCCTGGCGCGATTATCTTTATGATCCTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGgaag                                                   >  2:27904/1‑150 (MQ=255)
                                                                                                   gTGGTATTCCTACCGCCGACGCTCGTTGCTTCCAGCTATGGCATGAACTTTGAGTTTATGCCAGAACTGAAGTGGAGCTTCGGCTACCCTGGCGCGATTATCTTTATGATCCTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGgaag                                                   >  1:2357819/1‑150 (MQ=255)
                                                                                                    tGGTATTCCTGCCGACGACGCTCCTTGCTTCCAGCTATGGCATGAACTTTTAGTTTATGCCAGAACTGAAGTGGAGCTTCGGCTACCCTGGCGCGATTATCTTTATGATCCTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGgaaga                                                  <  2:1602985/150‑1 (MQ=255)
                                                                                                           ccTGCCGCCGACGCTCGTTGCTTCCAGCTATGGCATGAACTTTGAGTTTATGCCAGAACTGAAGTGGAGCTTCGGCTACCCTGGCGCGATTATCTTTATGATCCTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGtt                                           <  2:1017928/150‑1 (MQ=255)
                                                                                                            cTGCCGCCGACGCTCGTTGCTTCCAGCTATGGCATGAACTTTGAGTTTATGCCAGAACTGAAGTGGAGCTTCGGCTACCCTGGCGCGATTATCTTTATGATCCTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGTTg                                          <  2:1467248/150‑1 (MQ=255)
                                                                                                                cgccgACGCTCGTTGCTTCCAGCTATGGCATGAACTTTGAGTTTATGCCAGAACTGAAGTGGAGCTTCGGCTACCCTGGCGCGATTATCTTTATGATCCTCGCGGGCCTGGCACCGTATCt                                                                   <  2:349787/121‑1 (MQ=255)
                                                                                                                cgccgACGCTCGTTGCTTCCAGCTATGGCATGAACTTTGAGTTTATGCCAGAACTGAAGTGGAGCTTCGGCTACCCTGGCGCGATTATCTTTATGATCCTCGCGGGCCTGGCACCGTATCt                                                                   >  1:349787/1‑121 (MQ=255)
                                                                                                                cgccgACGCTCGTTGCTTCCAGCTATGGCATGAACTTTGAGTTTATGCCAGAACTGAAGTGGAGCTTCGGCTACCCTGGCGCGATTATCTTTATGATCCTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGTTGTaaa                                      >  1:2395324/1‑150 (MQ=255)
                                                                                                                cgccgACGCTCGTTGCTTCCAGCTATGGCATGAACTTTGAGTTTATGCCAGAACTGAAGTGGAGCTTCGGCTACCCTGGCGCGATTATCTTTATGATCCTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGTTGTaaa                                      <  1:920472/150‑1 (MQ=255)
                                                                                                                cgccgACGCTCGTTGCTTCCAGCTATGGCATGAACTTTGAGTTTATGCCAGAACTGAAGTGGAGCTTCGGCTACCCTGGCGCGATTATCTTTATGATCCTCGCGGGCCTGGCACCATATCTGTACTTTAAGCGGAAGAACTGGTTGTaaa                                      <  2:1231122/150‑1 (MQ=255)
                                                                                                                 gccgACGCTCGTTGCTTCCAGCTATGGCATGAACTTTGAGTTTATGCCAGAACTGAAGTGGAGCTTCGGCTACCCTGGCGCGATTATCTTTATGATCCTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGTTGTaaaa                                     >  2:880350/1‑150 (MQ=255)
                                                                                                                        cTCGTTGCTTCCAGCTATGGCATGAACTTTGAGTTTATGCCAGAACTGAAGTGGAGCTTCGGCTACCCTGGCGCGATTATCTTTATGATCCTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACgaga                              >  1:1221223/1‑150 (MQ=255)
                                                                                                                           gTTGCTTCCAGCTATGGCATGAACTTTGAGTTTATGCCAGAACTGAAGTGGAGCTTCGGCTACCCTGGCGCGATTATCTTTATGATCCTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCg                           <  1:1404454/150‑1 (MQ=255)
                                                                                                                              gCTTCCAGCTATGGCATGAACTTTGAGTTTATGCCAGAACTGAAGTGGAGCTTCGGCTACCCTGGCGCGATTATCTTTATGATCCTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTg                        >  1:1299864/1‑150 (MQ=255)
                                                                                                                               cTTCCAGCTATGGCATGAACTTTGAGTTTATGCCAGAACTGAAGTGGAGCTTCGGCTACCCTGGCGCGATTATCTTTATGATCCTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTgg                       <  1:2029412/150‑1 (MQ=255)
                                                                                                                               cTTCCAGCTATGGCATGAACTTTGAGTTTATGCCAGAACTGAAGTGGAGCTTCGGCTACCCTGGCGCGATTATCTTTATGATCCTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTgg                       >  2:723935/1‑150 (MQ=255)
                                                                                                                                ttCCAGCTATGGCATGAACTTTGAGTTTATGCCAGAACTGAAGTGGAGCTTCGGCTACCCTGGCGCGATTATCTTTATGATCCTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGc                      >  1:1045662/1‑150 (MQ=255)
                                                                                                                                 tCCAGCTATGGCATGAACTTTGAGTTTATGCCAGAACTGAAGTGGAGCTTCGGCTACCCTGGCGCGATTATCTTTATGATCCTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCt                     >  2:169174/1‑150 (MQ=255)
                                                                                                                                 tCCAGCTATGGCATGAACTTTGAGTTTATGCCAGAACTGAAGTGGAGCTTCGGCTACCCTGGCGCGATTATCTTTATGATCCTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCt                     >  2:636770/1‑150 (MQ=255)
                                                                                                                                   cAGCTATGGCATGAACTTTGAGTTTATGCCAGAACTGAAGTGGAGCTTCGGCTACCCTGGCGCGATTATCTTTATGATCCTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTa                   <  1:1055528/150‑1 (MQ=255)
                                                                                                                                    agttctGGCATAAACTCAAAGTTCATGCCAGAACTGAAGTGGAGCTTCGGCTACCCTGGCGCGATTATCTTTATGATCCTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAg                  >  1:2091037/6‑150 (MQ=255)
                                                                                                                                       tATGGCATGAACTTTGAGTTTATGCCAGAACTGAAGTGGAGCTTCGGCTACCCTGGCGCGATTATCTTTATGATCCTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAGTCt               >  1:590049/1‑150 (MQ=255)
                                                                                                                                          ggCATGAACTTTGAGTTTATGCCAGAACTGAAGTGGAGCTTCTGCTACCCTGGCGAGATTATCTTTATGATCCTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAGTCTGGc            <  1:1819303/150‑1 (MQ=255)
                                                                                                                                          ggCATGAACTTTGAGTTTATGCCAGAACTGAAGTGGAGCTTCGGCTACCCTGGCGCGATTATCTTTATGATCCTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAGTCTggg            >  1:2433783/1‑149 (MQ=255)
                                                                                                                                           gCATGAACTTTGAGTTTATGCCAGAACTGAAGTGGAGCTTCGGCTACCCTGGCGCGATTATCTTTATGATCCTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAGTCTGGCt           >  2:1633813/1‑150 (MQ=255)
                                                                                                                                           gCATGAACTTTGAGTTTATGCCAGAACTGAAGTGGAGCTTCGGCTACCCTGGCGCGATTATCTTTATGATCCTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAGTCTGGCt           <  1:45364/150‑1 (MQ=255)
                                                                                                                                            cATGAACTTTGAGTTTATGCCAGAACTGAAGTGGAGCTTCGGCTACCCTGGCGCGATTATCTTTATGATCCTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAGTCTGGCTa          >  2:2165352/1‑150 (MQ=255)
                                                                                                                                              tGAACTTTGAGTTTATGCCAGAACTGAAGTGGAGCTTCGGCTACCCTGGCGCGATTATCTTTATGATCCTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAGTCTGGCTAAg        <  1:1363791/150‑1 (MQ=255)
                                                                                                                                                  cTTTGAGTTTATGCCAGAACTGAAGTGGAGCTTCGGCTACCCTGGCGCGATTATCTTTATGATCCTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAGTCTGGCTAAGCCAc    >  2:905857/1‑150 (MQ=255)
                                                                                                                                                   tttGAGTTTATGCCAGAACTGAAGTGGAGCTTCGGCTACCCTGGCGCGATTATCTTTATGATCCTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAGTCTGGCTAAGCCAcc   <  2:1134213/150‑1 (MQ=255)
                                                                                                                                                    ttGAGTTTATGCCAGAACTGAAGTGGAGCTTCGGCTACCCTGGCGCGATTATCTTTATGATCCTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAGCCAGACTAAGCCAccg  <  2:1590632/150‑2 (MQ=255)
                                                                                                                                                    ttGAGTTTATGCCAGAACTGAAGTGGAGCTTCGGCTACCCTGGCGCGATTATCTTTATGATCCTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAGCCAGACTAAGCCAccg  <  1:1121080/150‑2 (MQ=255)
                                                                                                                                                    ttGAGTTTATGCCAGAACTGAAGTGGAGCTTCGGCTACCCTGGCGCGATTATCTTTATGATCCTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAGCCAGACTAAGCCAccg  <  2:821269/150‑2 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
GAATCCCTGTTCCAGAAGGTGAACTTCCTGATGCAGGCGGCAATGGGTTTTATCAACATCGAGCAGAACCGCATCATCAAAATCTTCTCGGTGGTATCCGTGGTATTCCTGCCGCCGACGCTCGTTGCTTCCAGCTATGGCATGAACTTTGAGTTTATGCCAGAACTGAAGTGGAGCTTCGGCTACCCTGGCGCGATTATCTTTATGATCCTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAGTCTGGCTAAGCCACCT  >  NC_000913/4002116‑4002413

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