Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,002,259 G→A 5.7% M278I (ATG→ATA)  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,2590GA 5.7% 557.1 / 13.8 176M278I (ATG→ATA) corANi(2(+))/Co(2(+))/Mg(2(+)) transporter
Reads supporting (aligned to +/- strand):  ref base G (84/81);  new base A (4/6);  total (88/88)
Fisher's exact test for biased strand distribution p-value = 5.35e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.68e-01

Only 100 of 178 total aligned reads displayed.

CATAACGAATCCCTGTTCCAGAAGGTGAACTTCCTGATGCAGGCGGCAATGGGTTTTATCAACATCGAGCAGAACCGCATCATCAAAATCTTCTCGGTGGTATCCGTGGTATTCCTGCCGCCGACGCTCGTTGCTTCCAGCTATGGCATGAACTTTGAGTTTATGCCAGAACTGAAGTGGAGCTTCGGCTACCCTGGCGCGATTATCTTTATGATCCTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAGTCTGGCTAAG  >  NC_000913/4002110‑4002407
                                                                                                                                                     |                                                                                                                                                    
  tAACGAATCCCTGTTCCAGAAGGTGAACTTCCTGATGCAGGCGGCAATGGGTTTTATCAACATCGAGCAGAACCGCATCATCAAAATCTTCTCGGTGGTATCCGTGGTATTCCTGCCGCCGACGCTCGTTGCTTCCAGCTATGGCATGaa                                                                                                                                                    >  2:1259664/1‑150 (MQ=255)
    aCGAATCCCTGTTCCAGAAGGTGAACTTCCTGATGCAGGCGGCAATGGGTTTTATCAACATCGAGCAGAACCGCATCATCAAAATCTTCTCGGTGGTATCCGTGGTATTCCTGCCGCCGACGCTCGTTGCTTCCAGCTATGGCATGAACt                                                                                                                                                  <  2:27977/150‑1 (MQ=255)
      gAATCCCTGTTCCAGAAGGTGAACTTCCTGATGCAGGCGGCAATGGTTTTTATCAACATCGAGCAGAACCGCATCATCAAAATCTTCTCGGTGGTATCCGTGGTATTCCTGCCGCCGACGCTCGTTGCTTCCAGCTATGGCATGAACttt                                                                                                                                                <  1:928694/150‑1 (MQ=255)
        aTCCCTGTTCCAGAAGGTGAACTTCCTGATGCAGGCGGCAATGGGTTTTATCAACATCGAGCAGAACCGCATCATCAAAATCTTCTCGGTGGTATCCGTGGTATTCCTGCCGCCGACGCTCGTTGCTTCCAGCTATGGCATGAACTTTGa                                                                                                                                              <  2:2480188/150‑1 (MQ=255)
         tCCCTGTTCCAGAAGGTGAACTTCCTGATGCAGGCGGCAATGGGTTTTATCAACATCGAGCAGAACCGCATCATCAAAATCTTCTCGGTGGTATCCGTGGTATTCCTGCCGCCGACGCTCGTTGCTTCCAGCTATGGCATGAACTTTGAg                                                                                                                                             <  1:1030299/150‑1 (MQ=255)
         tCCCTGTTCCAGAAGGTGAACTTCCTGATGCAGGCGGCAATGGGTTTTATCAACATCGAGCAGAACCGCATCATCAAAATCTTCTCGGTGGTATCCGTGGTATTCCTGCCGCCGACGCTCGTTGCTTCCAGCTATGGCATGAACTTTGAg                                                                                                                                             <  1:1738168/150‑1 (MQ=255)
         tCCCTGTTCCAGAAGGTGAACTTCCTGATGCAGGCGGCAATGGGTTTTATCAACATCGAGCAGAACCGCATCATCAAAATCTTCTCGGTGGTATCCGTGGTATTCCTGCCGCCGACGCTCGTTGCTTCCAGCTATGGCATGAACTTTGAg                                                                                                                                             <  1:129410/150‑1 (MQ=255)
             tGTTCCAGAAGGTGAACTTCCTGATGCAGGCGGCAATGGGTTTTATCAACATCGAGCAGAACCGCATCATCAAAATCTTCTCGGTGGTATCCGTGGTATTCCTGCCGCCGACGCTCGTTGCTTCCAGCTATGGCATGAACTTTGAGt                                                                                                                                            >  1:460237/1‑147 (MQ=255)
                tCCAGAAGGTGAACTTCCTGATGCAGGCGGCAATGGGTTTTATCAACATCGAGCAGAACCGCATCATCAAAATCTTCTCGGTGGTATCCGTGGTATTCCTGCCGCCGACGCTCGTTGCTTCCAGCTATGGCATGAACTTTGAGTTTATGc                                                                                                                                      >  1:1297457/1‑150 (MQ=255)
                 ccAGATGGTGAACTTCCTGTTGCAGGCGGAAATGGGTTTTATCAACATCGAGCAGAACCGCATCATCAAAATCTTCTCGGTGGTATCCGTTGTATTCCTGCCGCCGACGCTCGTTGCTTCCAGCTATGGCATGAACTTTGAGTTTATGcc                                                                                                                                     <  2:1618777/150‑1 (MQ=255)
                  cAGAAGGTGAACTTCCTGATGCAGGCGGCAATGGGTTTTATCAACATCGAGCAGAACCGCATCATCAAAATCTTCTCGGTGGTATCCGTGGTATTCCTGCCGCCGACGCTCGTTGCTTCCAGCTATGGCATAAACTCAAAGTTCATGCCa                                                                                                                                    <  1:1034038/150‑1 (MQ=255)
                    gAAGGTGAACTTCCTGATGCAGGCGGCAATGGGTTTTATCAACATCGAGCAGAACCGCATCATCAAAATCTTCTCGGTGGTATCCGTGGTATTCCTGCCGCCGACGCTCGTTGCTTCCAGCTATGGCATAAACTCAAAGTTCATGCCata                                                                                                                                  <  1:1460356/150‑3 (MQ=255)
                      aGGTGAACTTCCTGATGCAGGCGGCAATGGGTTTTATCAACATCGAGCAGAACCGCATCATCAAAATCTTCTCGGTGGTATCCGTGGTATTCCTGCCGCCGACGCTCGTTGCTTCCAGCTATGGCATGAACTTTGAGTTTATGCCAGAAc                                                                                                                                >  2:2231229/1‑150 (MQ=255)
                       ggTGAACTTCCTGATGCAGGCGGCAATGGGTTTTATCAACATCGAGCAGAACCGCATCATCAAAATCTTCTCGGTGGTATCCGTGGTATTCCTGCCGCCGACGCTCGTTGCTTCCAGCTATGGCATAAACTCAAAGTTCATGCCatagct                                                                                                                               <  1:2228314/150‑6 (MQ=255)
                        gTGAACTTCCTGATGCAGGCGGCAATGGGTTATATCAACATCGAGCAGAACCGCATCAGCAAAATCTTCTCGGGGGTATCCGTGGTATTCCTGCCGCCGACGCTCGTTGCTTCCAGCTATGGCATAAACTCAAAGTTCATGCCATAGCTg                                                                                                                              <  2:1403935/150‑1 (MQ=255)
                         ggAACTTCCTGATGCAGGCGGCAATGGGTTTTATCAACATCGAGCAGAACCGCATCATCAAAATCTTCTCGGTGGTATCCGTTTTTTTCCTGCCGCCGACGCTCGTTGCTTCCAGCTATGGCATAAACTCATAGTTCATGCCATAGCTgg                                                                                                                             <  2:1005977/149‑2 (MQ=255)
                            aCTTCCTGATGCAGGCGGCAATGGGTTTTATCAACATCGAGCAGAACCGCATCATCAAAATCTTCTCGGTGGTATCCGTGGTATTCCTGCCGCCGACGCTCGTTGCTTCCAGCTATGGCATGAACTTTGAGTTTATGCCAGAACTGAAGt                                                                                                                          <  2:2024121/150‑1 (MQ=255)
                             cTTCCTGATGCAGGCGGCAATGGGTTTTATCAACATCGAGCAGAACCGCATCATCAAAATCTTCTCGGTGGTATCCGTGGTATTCCTGCCGCCGACGCTCGTTGCTTCCAGCTATGGCATGAACTTTGAGTTTATGCCAGAACTGAAGTg                                                                                                                         >  1:1236997/1‑150 (MQ=255)
                             cTTCCTGATGCAGGCGGCAATGGGTTTTATCAACATCGAGCAGAACCGCATCATCAAAATCTTCTCGGTGGTATCCGTGGTATTCCTGCCGCCGACGCTCGTTGCTTCCAGCTATGGCATGAACTTTGAGTTTATGCCAGAACTGAAGTg                                                                                                                         >  2:846591/1‑150 (MQ=255)
                                 cTGATGCAGGCGGCAATGGGTTTTATCAACATCGAGCAGAACCGCATCATCAAAATCTTCTCGGTGGTATCCGTGGTATTCCTGCCGCCGACGCTCGTTGCTTCCAGCTATGGCATGAACTTTGAGTTTATGCCAGAACTGAAGTGGAGc                                                                                                                     <  1:1864604/150‑1 (MQ=255)
                                  tGATGGAGGCGGCAATGGGTTTTATCAACATCGAGCAGAACCGCATCATCAAAATCTTCTCGGTGGTATCCGTGGTATTCCTGCCGCCGACGCTCGTTGCTTCCAGCTATGGCATGAACTTTGAGTTTATGCCAGAACTGAAGTGGAGCt                                                                                                                    <  2:229957/150‑1 (MQ=255)
                                  tGATGCAGGCGGCAATGGGTTTTATCAACATCGAGCAGAACCGCATCATCAAAATCTTCTCGGTGGTATCCGCGGTATTCCTGCCGCCGACGCTCGTTGCTTCCAGCTATGGCATGAACTTTGAGTTTATGCCAGAACTGAAGTGGAGCt                                                                                                                    >  1:2309248/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)
                                       cAGGCGGCAATGGGTTTTATCAACATCGAGCAGAACCGCATCATCAAAATCTTCTCGGTGGTATCCGTGGTATTCCTGCCGCCGACGCTCGTTGCTTCCAGCTATGGCATGAACTTTGAGTTTATGCCAGAACTGAAGTGGAGCTTCGGc                                                                                                               <  1:1109711/150‑1 (MQ=255)
                                         ggcggcAATGGGTTTTATCAACATCGAGCAGAACCGCATCATCAAAATCTTCTCGGTGGTATCCGTGGTATTCCTGCCGCCGACGCTCGTTGCTTCCAGCTATGGCATGAACTTTGAGTTTATGCCAGAACTGAAGTGGAGCTTCGGCTa                                                                                                             >  2:1288205/1‑150 (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)
                                                  gggTTTTATCAACATCGAGCAGAACCGCATCATCAAAATCTTCTCGGTGGTATCCGTGGTATTCCTGCCGCCGACGCTCGTTGCTTCCAGCTATGGCATGAACTTTGAGTTTATGCCAGAACTGAAGTGGAGCTTCGGCTACCCTGgcgc                                                                                                    >  1:1134213/1‑150 (MQ=255)
                                                  gggTTTTATCAACATCGAGCAGAACCGCATCATCAAAATCTTCTCGGTGGTATCCGTGGTATTCCTGCCGCCGACGCTCGTTGCTTCCAGCTATGGCATGAACTTTGAGTTTATGCCAGAACTGAAGTGGAGCTTCGGCTACCCTGgcgc                                                                                                    <  1:1733538/150‑1 (MQ=255)
                                                   ggTTTTATCAACATCGAGCAGAACCGCATCATCAAAATCTTCTCGGTGGTATCCGTGGTATTCCTGCCGCCGACGCTCGTTGCTTCCAGCTATGGCATGAACTTTGAGTTTATGCCAGAACTGAAGTGGAGCTTCGGCTACCCTGGcgcg                                                                                                   >  2:182693/1‑150 (MQ=255)
                                                     ttttATCAACATCGAGCAGAACCGCATCATCAAAATCTTCTCGGTGGTATCCGGGGTATTCCTGCCGCCGACGCTCGTTGCTTCCAGCTATGGCATGAACTTTGAGTTTATGCCAGAACTGAAGTGGAGCTTCGGCTACCCTGGCGCGAt                                                                                                 <  2:1500590/150‑1 (MQ=255)
                                                          tCAACATCGAGCAGAACCGCATCATCAAAATCTTCTCGGTGGTATCCGTGGTATTCCTGCCGCCGACGCTCGTTGCTTCCAGCTATGGCATGAACTTTGAGTTTATGCCAGAACTGAAGTGGAGCTTCGGCTACCCTGGCGCGATTATCt                                                                                            >  1:1022606/1‑150 (MQ=255)
                                                           cAACATCGAGCAGAACCGCATCATCAAAATCTTCTCGGTGGTATCCGTGGTATTCCTGCCGCCGACGCTCGTTGCTTCCAGCTATGGCATGAACTTTGAGTTTATGCCAGAACTGAAGTGGAGCTTCGGCTACCCTGGCGCGATTATCtt                                                                                           <  1:2388186/150‑1 (MQ=255)
                                                            aaCATCGAGCAGAACCGCATCATCAAAATCTTCTCGGTGGTATCCGTGGTATTCCTGCCGCCGACGCTCGTTGCTTCCAGCTATGGCATGAACTTTGAGTTTATGCCAGAACTGAAGTGGAGCTTCGGCTACCCTGGCGCGATTATCttt                                                                                          >  1:1047322/1‑150 (MQ=255)
                                                            aaCATCGAGCAGAACCGCATCATCAAAATCTTCTCGGTGGTATCCGTGGTATTCCTGCCGCCGACGCTCGTTGCTTCCAGCTATGGCATGAACTTTGAGTTTATGCCAGAACTGAAGTGGAGCTTCGGCTACCCTGGCGCGATTATCttt                                                                                          >  1:905823/1‑150 (MQ=255)
                                                             aCATCGAGCAGAACCGCATCATCAAAATCTTCTCGGTGGTATCCGTGGTATTCCTGCCGCCGACGCTCGTTGCTTCCAGCTATGGCATGAACTTTGAGTTTATGCCAGAACTGAAGTGGAGCTTCGGCTACCCTGGCGCGATTATCTTTa                                                                                         >  1:1300427/1‑150 (MQ=255)
                                                               aTCGAGCAGAACCGCATCATCAAAATCTTCTCGGTGGTATCCGTGGTATTCCTGCCGCCGACGCTCGTTGCTTCCAGCTATGGCATGAACTTTGAGTTTATGCCAGAACTGAAGTGGAGCTTCGGCTACCCTGGCGCGATTATCTTTATg                                                                                       <  2:1115196/150‑1 (MQ=255)
                                                                 cGAGCAGAACCGCATCATCAAAATCTTCTCGGTGGTATCCGTGGTATTCCTGCCGCCGACGCTCGTTGCTTCCAGCTATGGCATGAACTTTGAGTTTATGCCAGAACTGAAGTGGAGCTTCGGCTACCCTGGCGCGATTATCTTTATGAt                                                                                     >  2:984759/1‑150 (MQ=255)
                                                                  gAGCAGAACCGCATCATCAAAATCTTCTCGGTGGTATCCGTGGTATTCCTGCCGCCGACGCTCGTTGCTTCCAGCTATGGCATGAACTTTGAGTTTATGCCAGAACTGAAGTGGAGCTTCGGCTACCCTGGCGCGATTATCTTTATGATc                                                                                    >  2:495295/1‑150 (MQ=255)
                                                                   aGCAGAACCGCATCATCAAAATCTTCTCGGTGGTATCCGTGGTATTCCTGCCGCCGACGCTCGTTGCTTCCAGCTATGGCATGAACTTTGAGTTTATGCCAGAACTGAAGTGGAGCTTCGGCTACCCTGGCGCGATTATCTTTATGATcc                                                                                   <  1:1284419/150‑1 (MQ=255)
                                                                    gCAGAACCGCATCATCAAAATCTTCTCGGTGGTATCCGTGGTATTCCTGCCGCCGACGCTCGTTGCTTCCAGCTATGGCATGAACTTTGAGTTTATGCCAGAACTGAAGTGGAGCTTCGGCTACCCTGGCGCGATTATCTTTATGATCCt                                                                                  >  1:536269/1‑150 (MQ=255)
                                                                      aGAACCGCATCATCAAAATCTTCTCGGTGGTATCCGTGGTATTCCTGCCGCCGACGCTCGTTGCTTCCAGCTATGGCATGAACTTTGAG‑TTATGCCAGAACTGAAGTGGAGCTTCGGCTACCCTGGCGCGATTATCTTTATGATCCTcgc                                                                               >  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:278321/1‑150 (MQ=255)
                                                                             catcatCAAAATCTTCTCGGTGGTATCCGTGGTATTCCTGCCGCCGACGCTCGTTGCTTCCAGCTATGGCATGAACTTTGAGTTTATGCCAGAACTGAAGTGGAGCTTCGGCTACCCTGGCGCGATTATCTTTATGATCCTCGCGGGCCt                                                                         <  1:370040/150‑1 (MQ=255)
                                                                              atcatcAAAATCTTCTCGGTGGTATCCGTGGTATTCCTGCCGCCGACGCTCGTTGCTTCCAGCTATGGCATGAACTTTGAGTTTATGCCAGAACTGAAGTGGAGCTTCGGCTACCCTGGCGCGATTATCTTTATGATCCTCGCGGGCCTg                                                                        <  2:103968/150‑1 (MQ=255)
                                                                                  tcaAAATCTTCTCGGTGGTATCCGTGGTATTCCTGCCGCCGACGCTCGTTGCTTCCAGCTATGGCATGAACTTTGAGTTTATGCCAGAACTGAAGTGGAGCTTCGGCTACCCTGGCGCGATTATCTTTATGATCCTCGCGGGCCTGGCAc                                                                    <  1:1956685/150‑1 (MQ=255)
                                                                                   caAAATCTTCTCGGTGGTATCCGTGGTATTCCTGCCGCCGACGCTCGTTGCTTCCAGCTATGGCATGAACTTTGAGTTTATGCCAGAACTGAAGTGGAGCTTCGGCTACCCTGGCGCGATTATCTTTATGATCCTCGCGGGCCTGGCAcc                                                                   >  2:45364/1‑150 (MQ=255)
                                                                                    aaaaTCTTCTCGGTGGTATCCGTGGTATTCCTGCCGCCGACGCTCGTTGCTTCCAGCTATGGCATGAACTTTGAGTTTATGCCAGAACTGAAGTGGAGCTTCGGCTACCCTGGCGCGATTATCTTTATGATCCTCGCGGGCCTGGCACCg                                                                  <  1:1657554/150‑1 (MQ=255)
                                                                                    aaaaTCTTCTCGGTGGTATCCGTGGTATTCCTGCCGCCGACGCTCGTTGCTTCCAGCTATGGCATGAACTTTGAGTTTATGCCAGAACTGAAGTGGAGCTTCGGCTACCCTGGCGCGATTATCTTTATGATCCTCGCGGGCCTGGCACCg                                                                  >  2:948298/1‑150 (MQ=255)
                                                                                      aaTCTTCTCGTTTGTATCCGGGGTATTCCTGCCGCCGACGCTAGTTGCTTCCAGCTATGGCATGAACTTTGAGTTTTTGCCAGAACTGAAGTGGAGCTTCGGCTACCCTGGCGCGATTATCTTTATGATCCTCGCGGGCCTGGCACCGTa                                                                <  2:1859261/150‑1 (MQ=255)
                                                                                      aaTCTTCTCGGTGGTATCCGCGGTATTCCTGCCGCCGACGCTCGTTGCTTCCAGCTATGGCATGAACTTTGAGTTTATGCCAGAACTGAAGTGGAGCTTCGGCTACCCTGGCGCGATTATCTTTATGATCCTCGCGGGCCTGGCACCGTa                                                                <  2:2309248/150‑1 (MQ=255)
                                                                                         cttctCGGTGGTATCCGTGGTATTCCTGCCGCCGACGCTCGTTGCTTCCAGCTATGGCATGAACTTTGAGTTTATGCCAGAACTGAAGTGGAGCTTCGGCTACCCTGGCGCGATTATCTTTATGATCCTCGCGGGCCTGGCACCGTATCt                                                             <  2:895297/150‑1 (MQ=255)
                                                                                         cttctCGGTGGTATCCGTGGTATTCCTGCCGCCGACGCTCGTTGCTTCCAGCTATGGCATGAACTTTGAGTTTATGCCAGAACTGAAGTGGAGCTTCGGCTACCCTGGCGCGATTATCTTTATGATCCTCGCGGGCCTGGCACCGTATCt                                                             >  1:7233/1‑150 (MQ=255)
                                                                                          ttgtcGGTGGTATCCGTGGTATTCCTGCCGCCGACGCTCGTTGCTTCCAGCTATGGCATGAACTTTGAGTTTATGCCAGAACTGAAGTGGAGCTTCGGCTACCCTGGCGCGATTATCTTTATGATCCTCGCGGGCCTGGCACCGTATCTg                                                            <  2:689093/147‑1 (MQ=255)
                                                                                            ctcGGTGGTATCCGTGGTATTCCTGCCGCCGACGCTCGTTGCTTCCAGCTATGGCATGAACTTTGAGTTTATGCCAGAACTGAAGTGGAGCTTCGGCTACCCTGGCGCGATTATCTTTATGATCCTCGCGGGCCTGGCACCGTATCTGTa                                                          <  1:121132/150‑1 (MQ=255)
                                                                                              cGGTGGTATCCGTGGTATTCCTGCCGCCGACGCTCGTTGCTTCCAGCTATGGCATGAACTTTGAGTTTATGCCAGAACTGAAGTGGAGCTTCGGCTACCCTGGCGCGATTATCTTTATGATCCTCGCGGGCCTGGCACCGTATCTGTACt                                                        <  2:2449274/150‑1 (MQ=255)
                                                                                               ggtggtATCCGTGGTATTCCTGCCGCCGACGCTCGTTGCTTCCAGCTATGGCATGAACTTTGAGTTTATGCCAGAACTGAAGTGGAGCTTCGGCTACCCTGGCGCGATTATCTTTATGATCCTCGCGGGCCTGGCACCGTATCTGTACtt                                                       >  1:1531423/1‑150 (MQ=255)
                                                                                                    tATCCGTGGTATTCCTGCCGCCGACGCTCGTTGCTTCCAGCTATGGCATGAACTTTGAGTTTATGCCAGAACTGAAGTGGAGCTTCGGCTACCCTGGCGCGATTATCTTTATGATCCTCGCGGGCCTGGCACCGTATCTGTACTTTAAGc                                                  <  2:819867/150‑1 (MQ=255)
                                                                                                     aTCCGTGGTATTCCTGCCGCCGACGCTCGTTGCTTCCAGCTATGGCATGAACTTTGAGTTTATGCCAGAACTGAAGTGGAGCTTCGGCTACCCTGGCGCGATTATCTTTATGATCCTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCg                                                 <  1:779050/150‑1 (MQ=255)
                                                                                                      tCCGTGGTATTCCTGCCGCCGACGCTCGTTGCTTCCAGCTATGGCATGAACTTTGAGTTTATGCCAGAACTGAAGTGGAGCTTCGGCTACCCTGGCGCGATTATCTTTATGATCCTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCgg                                                >  2:2271736/1‑150 (MQ=255)
                                                                                                       ccGTGGTATTCCTGCCGCCGACGCTCGTTGCTTCCAGCTATGGCATGAACTTTGAGTTTATGCCAGAACTGAAGTGGAGCTTCGGCTACCCTGGCGCGATTATCTTTATGATCCTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGga                                               >  2:414819/1‑150 (MQ=255)
                                                                                                       ccGTGGTATTCCTGCCGCCGACGCTCGTTGCTTCCAGCTATGGCATGAACTTTGAGTTTATGCCAGAACTGAAGTGGAGCTTCGGCTACCCTGGCGCGATTATCTTTATGATCCTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGga                                               >  2:1878060/1‑150 (MQ=255)
                                                                                                         gTGGTATTCCTGCCGCCGACGCTCGTTGCTTCCAGCTATGGCATGAACTTTGAGTTTATGCCAGAACTGAAGTGGAGCTTCGGCTACCCTGGCGCGATTATCTTTATGATCCTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGgaag                                             >  1:1082738/1‑150 (MQ=255)
                                                                                                          tGGTATTCCTGCCGACGACGCTCCTTGCTTCCAGCTATGGCATGAACTTTTAGTTTATGCCAGAACTGAAGTGGAGCTTCGGCTACCCTGGCGCGATTATCTTTATGATCCTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGgaaga                                            <  2:1602985/150‑1 (MQ=255)
                                                                                                                 ccTGCCGCCGACGCTCGTTGCTTCCAGCTATGGCATGAACTTTGAGTTTATGCCAGAACTGAAGTGGAGCTTCGGCTACCCTGGCGCGATTATCTTTATGATCCTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGtt                                     <  2:541175/150‑1 (MQ=255)
                                                                                                                 ccTGCCGCCGACGCTCGTTGCTTCCAGCTATGGCATGAACTTTGAGTTTATGCCAGAACTGAAGTGGAGCTTCGGCTACCCTGGCGCGATTATCTTTATGATCCTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGtt                                     <  1:1492074/150‑1 (MQ=255)
                                                                                                                 ccTGCCGCCGACGCTCGTTGCTTCCAGCTATGGCATGAACTTTGAGTTTATGAAAGAACTGAAGTTGAGCTTCGTCTACCCTGGCGCGATTATCTTTATGATCCTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGtt                                     >  2:1910286/1‑150 (MQ=255)
                                                                                                                  cTGCCGCCGACGCTCGTTGCTTCCAGCTATGGCATGAACTTTGAGTTTATGCCAGAACTGAAGTGGAGCTTCGGCTACCCTGGCGCGATTATCTTTATGATCCTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGTTg                                    <  2:1467248/150‑1 (MQ=255)
                                                                                                                      cgccgACGCTCGTTGCTTCCAGCTATGGCATGAACTTTGAGTTTATGCCAGAACTGAAGTGGAGCTTCGGCTACCCTGGCGCGATTATCTTTATGATCCTCGCGGGCCTGGCACCGTATCt                                                             <  2:349787/121‑1 (MQ=255)
                                                                                                                      cgccgACGCTCGTTGCTTCCAGCTATGGCATGAACTTTGAGTTTATGCCAGAACTGAAGTGGAGCTTCGGCTACCCTGGCGCGATTATCTTTATGATCCTCGCGGGCCTGGCACCATATCTGTACTTTAAGCGGAAGAACTGGTTGTaaa                                <  2:1231122/150‑1 (MQ=255)
                                                                                                                            cGCTCGTTGCTTCCAGCTATGGCATGAACTTTGAGTTTATGCCAGAACTGAAGTGGAGCTTCGGCTACCCTGGCGCGATTATCTTTATGATCCTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACga                          >  2:2236546/1‑150 (MQ=255)
                                                                                                                             gCTCGTTGCTTCCAGCTATGGCATGAACTTTGAGTTTATGCCAGAACTGAAGTGGAGCTTCGGCTACCCTGGCGCGATTATCTTTATGATCCTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACgag                         <  2:1166045/150‑1 (MQ=255)
                                                                                                                              cTCGTTGCTTCCAGCTATGGCATGAACTTTGAGTTTATGCCAGAACTGAAGTGGAGCTTCGGCTACCCTGGCGCGATTATCTTTATGATCCTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACgaga                        >  1:1221223/1‑150 (MQ=255)
                                                                                                                                     cTTCCAGCTATGGCATGAACTTTGAGTTTATGCCAGAACTGAAGTGGAGCTTCGGCTACCCTGGCGCGATTATCTTTATGATCCTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTgg                 <  1:2029412/150‑1 (MQ=255)
                                                                                                                                     cTTCCAGCTATGGCATGAACTTTGAGTTTATGCCAGAACTGAAGTGGAGCTTCGGCTACCCTGGCGCGATTATCTTTATGATCCTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTgg                 >  2:723935/1‑150 (MQ=255)
                                                                                                                                     cTTCCAGCTATGGCATGAACTTTGAGTTTATGCCAGAACTGAAGTGGAGCTTCGGCTACCCTGGCGCGATTATCTTTATGATCCTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTgg                 >  1:2183069/1‑150 (MQ=255)
                                                                                                                                      ttCCAGCTATGGCATGAACTTTGAGTTTATGCCAGAACTGAAGTGGAGCTTCGGCTACCCTGGCGCGATTATCTTTATGATCCTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGc                >  1:1045662/1‑150 (MQ=255)
                                                                                                                                       tCCAGCTATGGCATGAACTTTGAGTTTATGCCAGAACTGAAGTGGAGCTTCGGCTACCCTGGCGCGATTATCTTTATGATCCTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCt               >  2:169174/1‑150 (MQ=255)
                                                                                                                                       tCCAGCTATGGCATGAACTTTGAGTTTATGCCAGAACTGAAGTGGAGCTTCGGCTACCCTGGCGCGATTATCTTTATGATCCTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCt               >  1:1920901/1‑150 (MQ=255)
                                                                                                                                       tCCAGCTATGGCATGAACTTTGAGTTTATGCCAGAACTGAAGTGGAGCTTCGGCTACCCTGGCGCGATTATCTTTATGATCCTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCt               >  2:636770/1‑150 (MQ=255)
                                                                                                                                       tCCAGCTATGGCATGAACTTTGAGTTTATGCCAGAACTGAAGTGGAGCTTCGGCTACCCTGGCGCGATTATCTTTATGATCCTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCt               >  1:926739/1‑150 (MQ=255)
                                                                                                                                        ccAGCTATGGCCTCACCTTTGAGTTTATGCCAGAACTGAAGTCGAGCTTCTGCTACCCTCGCGCGATTATCTTTATGATCCTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCtt              <  1:424999/150‑1 (MQ=255)
                                                                                                                                          agttctGGCATAAACTCAAAGTTCATGCCAGAACTGAAGTGGAGCTTCGGCTACCCTGGCGCGATTATCTTTATGATCCTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAg            >  1:2091037/6‑150 (MQ=255)
                                                                                                                                           gttctGGCATAAACTCAAAGTTCAAGACAGAACTTAAATGCACATTCGGCTACCCTCCCGAGATTAACATAATGATCCACGCGGGCCTGGCACCGTATCTGAACAAAAAGACGAAGCAAAGGTTGTAAAAAAATAGAGCTGTGGCTTact           >  2:1401374/5‑148 (MQ=255)
                                                                                                                                            cTATGGCATGAACTTTGAGTTTATGCCAGAACTGAAGTGGAGCTTCGGCTACCCTGGCGCGATTATCTTTATGATCCTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAGTc          <  2:2019909/150‑1 (MQ=255)
                                                                                                                                            cTATGGCATGAACTTTGAGTTTATGCCAGAACTGAAGTGGAGCTTCGGCTACCCTGGCGCGATTATCTTTATGATCCTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAGTc          >  2:795655/1‑150 (MQ=255)
                                                                                                                                              aTGGCATGAACTTTGAGTTTATGCCAGAACTGAAGTGGAGCTTCGGCTACCCTGGCGCGATTATCTTTATGATCCTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAGTCTg        <  1:2438181/150‑1 (MQ=255)
                                                                                                                                               tGGCATAAACTCAAAGTTCATGCCAGAACTGAAGTGGAGCTTCGGCTACCCTGGCGCGATTATCTTTATGATCCTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAGTCTgg       >  2:319994/1‑150 (MQ=255)
                                                                                                                                                ggCATGAACTTTGAGTTTATGCCATAACTGAAGTGGAGCTTCGGCTACGCTGGCGCGATTATCTTTATGATCCTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAGTCTGGc      <  1:1048017/150‑1 (MQ=255)
                                                                                                                                                ggCATGAACTTTGAGTTTATGCCAGAACTGAAGTGGAGCTTCTGCTACCCTGGCGAGATTATCTTTATGATCCTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAGTCTGGc      <  1:1819303/150‑1 (MQ=255)
                                                                                                                                                ggCATGAACTTTGAGTTTATGCCAGAACTGAAGTGGAGCTTCGGCTACCCTGGCGCGATTATCTTTATGATCCTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAGTCTggg      >  1:2433783/1‑149 (MQ=255)
                                                                                                                                                  cATGAACTTTGAGTTTATGCCAGAACTGAAGTGGAGCTTCGGCTACCCTGGCGCGATTATCTTTATGATCCTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAGTCTGGCTa    >  2:2165352/1‑150 (MQ=255)
                                                                                                                                                    tGAACTTTGAGTTTATGCCAGAACTGAAGTGGAGCTTCGGCTACCCTGGCGCGATTATCTTTATGATCCTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAGTCTGGCTAAg  <  1:1363791/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
CATAACGAATCCCTGTTCCAGAAGGTGAACTTCCTGATGCAGGCGGCAATGGGTTTTATCAACATCGAGCAGAACCGCATCATCAAAATCTTCTCGGTGGTATCCGTGGTATTCCTGCCGCCGACGCTCGTTGCTTCCAGCTATGGCATGAACTTTGAGTTTATGCCAGAACTGAAGTGGAGCTTCGGCTACCCTGGCGCGATTATCTTTATGATCCTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAGTCTGGCTAAG  >  NC_000913/4002110‑4002407

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: ‑