Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,408,292 T→A 8.9% R81R (CGT→CGA) ‡ zntB → Zn(2(+)):H(+) symporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,408,2920TA 8.9% 217.0 / 10.9 79R81R (CGT→CGA) ‡zntBZn(2(+)):H(+) symporter
Reads supporting (aligned to +/- strand):  ref base T (31/41);  new base A (4/3);  total (35/44)
Fisher's exact test for biased strand distribution p-value = 6.93e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.04e-01

AGATGTGATTGATGAAGCGCATCCCTGTTGGCTCCACCTTAATTATGTACACCATGATAGCGCCCAATGGCTGGCGACAACACCGCTGCTTCCCAATAACGTACGTGATGCGCTGGCGGGCGAGAGCACGCGTCCCCGAGTCAGCCGTCTCGGTGAAGGCACGCTGATTACATTGCGCTGTATAAACGGCAGCACCGATGAACGCCCCGATCAACTGGTCGCCATGCGTGTATATATGGACGGGCGGTTAATTGTTTCGACCCGACAACGCAAAGTGCTGGCGCTGGACGATGTGG  >  NC_000913/1408145‑1408440
                                                                                                                                                   |                                                                                                                                                    
aGATGTGATTGATGAAGCGCATCCCTGTTGGCTCCACCTTAATTATGTACACCATGATAGCGCCCAATGGCTGGCGACAACACCGCTGCTTCCCAATAACGTACGTGATGCGCTGGCGGGCGAGAGCACGCGTCCCCGAGTCAGCCGtct                                                                                                                                                    <  1:563556/150‑1 (MQ=255)
  aTGTGATTGATGAAGCGCATCCCTGTTGGCTCCACCTTAATTATGTACACCATGATAGCGCCCAATGGCTGGCGACAACACCGCTGCTTCCCAATACCGTACGTGATGCGCTGGCGGGCGAGAGCCCCCGTCCCCGAGTCAGCCGTCTCg                                                                                                                                                  >  2:657404/1‑150 (MQ=255)
   tgtgATTGATGAAGCGCATCCCTGTTGGCTCCACCTTAATTATGTACACCATGATAGCGCCCAATGGCTGGCGACAACACCGCTGCTTCCCAATAACGTACGTGATGCGCTGGCGGGCGAGAGCACGCGTCCCCGAGTCAGCCGTCTCgg                                                                                                                                                 >  2:994427/1‑150 (MQ=255)
   tgtgATTGATGAAGCGCATCCCTGTTGGCTCCACCTTAATTATGTACACCATGATAGCGCCCAATGGCTGGCGACAACACCGCTGCTTCCCAATAACGTACGTGATGCGCTGGCGGGCGAGAGCACGCGTCCCCGAGTCAGCCGTCTCgg                                                                                                                                                 <  2:233675/150‑1 (MQ=255)
     tgattgatGAAGCGCATCCCTGTTGGCTCCACCTTAATTATGTACACCATGATAGCGCCCAATGGCTGGCGACAACACCGCTGCTTCCCAATAACGTACGTGATGCGCTGGCGGGCGAGAGCACGCGTCCCCGAGTCAGCCGTCTCGGTg                                                                                                                                               <  2:525811/150‑1 (MQ=255)
     tgattgatGAAGCGCATCCCTGTTGGCTCCACCTTAATTATGTACACCATGATAGCGCCCAATGGCTGGCGACAACACCGCTGCTTCCCAATAACGTACGTGATGCGCTGGCGGGCGAGAGCACGCGTCCCCGAGTCAGCCGTCTCGGTg                                                                                                                                               <  2:520343/150‑1 (MQ=255)
      gattgatGAAGCGCATCCCTGTTGGCTCCACCTTAATTATGTACACCATGATAGCGCCCAATGGCTGGCGACAACACCGCTGCTTCCCAATAACGTACGTGATGCGCTGGCGGGCGAGAGCACGCGTCCCCGAGTCAGCCGTCTCGGTGa                                                                                                                                              >  2:790597/1‑150 (MQ=255)
      gattgatGAAGCGCATCCCTGTTGGCTCCACCTTAATTATGTACACCATGATAGCGCCCAATGGCTGGCGACAACACCGCTGCTTCCCAATAACGTACGTGATGCGCTGGCGGGCGAGAGCACGCGTCCCCGAGACGGCTGACTCGGGGa                                                                                                                                              <  2:1345081/150‑1 (MQ=255)
         tgatgaAGCGCATCCCTGTTGGCTCCACCTTAATTATGTACACCATGATAGCGCCCAATGGCTGGCGACAACACCGCTGCTTCCCAATAACGTACGTGATGCGCTGGCGGGCGAGAGCACGCGTCCCCGAGTCAGCCGTCTCTGTGAAgg                                                                                                                                           <  1:1392344/150‑1 (MQ=255)
         tgatgaAGCGCATCCCTGTTGGCTCCACCTTAATTATGTACACCATGATAGCGCCCAATGGCTGGCGACAACACCGCTGCTTCCCAATAACGTACGTGATGCGCTGGCGGGCGAGAGCACGCGTCCCCGAGTCAGCCGTCTCGGTGAAgg                                                                                                                                           <  2:463565/150‑1 (MQ=255)
           atgaAGCGCATCCCTGTTGGCTCCACCTTAATTATGTACACCATGATAGCGCCCAATGGCTGGCGACAACACCGCTGCTTCCCAATAACGTACGTGATGCGCTGGCGGGCGAGAGCACGCGTCCCCGAGTCAGCCGTCTCGGTGAAGGCa                                                                                                                                         >  1:1559975/1‑150 (MQ=255)
            tgaAGCGCATCCCTGTTGGCTCCACCTTAATTATGTACACCATGATAGCGCCCAATGGCTGGCGACAACACCGCTGCTTCCCAATAACGTACGTGATGCGCTGGCGGGCGAGAGCACGCGTCCCCGAGTCAGCCGTCTCGGTGAAGGCAc                                                                                                                                        >  2:535703/1‑150 (MQ=255)
               aGCGCATCCCTGTTGGCTCCACCTTAATTATGTACACCATGATAGCGCCCAATGGCTGGCGACAACACCGCTGCTTCCCAATAACGTACGTGATGCGCTGGCGGGCGAGAGCACGCGTCCCCGAGTCAGCCGTCTCGGTGAAGGCACGCt                                                                                                                                     >  2:1009335/1‑150 (MQ=255)
               aGCGCATCCCTGTTGGCTCCACCTTAATTATGTACACCATGATAGCGCCCAATGGCTGGCGACAACACCGCTGCTTCCCAATAACGTACGTGATGCGCTGGCGGGCGAGAGCACGCGTCCCCGAGTCAGCCGTCTCGGTGAAGGCACGCt                                                                                                                                     <  2:1273472/150‑1 (MQ=255)
               aGCGCATCCCTGTTGGCTCCACCTTAATTATGTACACCATGATAGCGCCCAATGGCTGGCGACAACACCGCTGCTTCCCAATAACGTACGTGATGCGCTGGCGGGAGAGAGCACGCGTCCCCGAGTCAGCCGTCTCGGTGAAGGCACGCt                                                                                                                                     <  2:954651/150‑1 (MQ=255)
                gcgcATCCCTGTTGGCTCCACCTTAATTATGTACACCATGATAGCGCCCAATGGCTGGCGACAACACCGCTGCTTCCCAATAACGTACGTGATGCGCTGGCGGGCGAGAGCACGCGTCCCCGAGTCAGCCGTCTCGGTGAAGGCACGCTg                                                                                                                                    <  2:312305/150‑1 (MQ=255)
                gcgcATCCCTGTTGGCTCCACCTTAATTATGTACACCATGATAGCGCCCAATGGCTGGCGACAACACCGCTGCTTCCCAATAACGTACGTGATGCGCTGGCGGGCGAGAGCACGCGTCCCCGAGTCAGCCGTCTCGGTGAAGGCACGCTg                                                                                                                                    <  2:1466559/150‑1 (MQ=255)
                 cgcATCCCTGTTGGCTCCACCTTAATTATGTACACCATGATAGCGCCCAATGGCTGGCGACAACACCGCTGCTTCCCAATAACGTACGTGATGCGCTGGCGGGCGAGAGCACGCGTCCCCGAGTCAGCCGTCTCGGTGAAGGCACGCTGa                                                                                                                                   >  1:289616/1‑150 (MQ=255)
                  gcATCCCTGTTGGCTCCACCTTAATTATGTACACCATGATAGCGCCCAATGGCTGGCGACAACACCGCTGCTTCCCAATAACGTACGTGATGCGCTGGCGGGCGAGAGCACGCGTCCCCGAGTCAGCCGTCTCGGTGAAGGCACGCTGAt                                                                                                                                  <  1:555001/150‑1 (MQ=255)
                  gcATCCCTGTTGGCTCCACCTTAATTATGTACACCATGATAGCGCCCAATGGCTGGCGACAACACCGCTGCTTCCCAATAACGTACGTGATGCGCTGGCGGGCGAGAGCACGCGTCCCCGAGTCAGCCGTCTCGGTGAAGGCACGCTGAt                                                                                                                                  >  2:1006366/1‑150 (MQ=255)
                  gcATCCCTGTTGGCTCCACCTTAATTATGTACACCATGATAGCGCCCAATGGCTGGCGACAACACCGCTGCTTCCCAATAACGTACGTGATGCGCTGGCGGGCGAGAGCACGCGTCCCCGAGTCAGCCGTCTCGGTGAAGGCACGCTGAt                                                                                                                                  >  1:733343/1‑150 (MQ=255)
                    aTCCCTGTTGGCTCCACCTTAATTATGTACACCATGATAGCGCCCAATGGCTGGCGACAACACCGCTGCTTCCCAATAACGTACGTGATGCGCTGGCGGGCGAGAGCACGCGTCCCCGAGTCAGCCGTCTCGGTGAAGGCACGCTGATTa                                                                                                                                <  1:1383668/150‑1 (MQ=255)
                      cccTGTTGGCTCCACCTTAATTATGTACACCATGATAGCGCCCAATGGCTGGCGACAACACCGCTGCTTCCCAATAACGTACGTGATGCGCTGGCGGGCGAGAGCACGCGTCCCCGAGTCAGCCGTCTCGGTGAAGGCACGCTGATTACa                                                                                                                              >  2:1390001/1‑150 (MQ=255)
                              gCTCCACCTTAATTATGTACACCATGTTAGCGCCCAATGGCTGGCGACAACACCGCTGCTTCCCAATAACGTACGTGATGCGCTGGCGGGCGAGAGCACGCGTCCCCGAGTCAGCCGTCTCGGTGAAGGCACGCTGATTACATTGCGCTg                                                                                                                      <  2:234303/150‑1 (MQ=255)
                              gCTCCACCTTAATTATGTACACCATGATAGCGCCCAATGGCTGGCGACAACACCGCTGCTTCCCAATAACGTACGTGATGCGCTGGCGGGCGAGAGCACGCGTCCCCGAGTCAGCCGTCTCGGTGAAGGCACGCTGATTACATTGCGCTg                                                                                                                      <  1:1009335/150‑1 (MQ=255)
                                  cACCTTAATTATGTACACCATGATAGCGCCCAATGGCTGGCGACAACACCGCTGCTTCCCAATAACGTACGTGATGCGCTGGCGGGCGAGAGCACGCGTCCCCGAGTCAGCCGTCTCGGTGAAGGCACGCTGATTACATTGCGCTGtata                                                                                                                  <  2:443297/150‑1 (MQ=255)
                                   aCCTTAATTATGTACACCATGATAGCGCCCAATGGCTGGCGACAACACCGCTGCTTCCCAATAACGTACGTGATGCGCTGGCGGGCGAGAGCACGCGTCCCCGAGTCAGCCGTCTCGGTGAAGGCACGCTGATTACATTGCGCTGTATaa                                                                                                                 <  2:1521433/150‑1 (MQ=255)
                                     cTTAATTATGTACACCATGATAGCGCCCAATGGCTGGCGACAACACCGCTGCTTCCCAATAACGTACGTGATGCGCTGGCGGGCGAGAGCACGCGTCCCCGAGTCAGCCGTCTCGGTGAAGGCACGCTGATTACATTGCGCTGTATAAAc                                                                                                               >  1:1002071/1‑150 (MQ=255)
                                        aaTTATGTACACCATGATAGCGCCCAATGGCTGGCGACAACACCGCTGCTTCCCAATAACGTACGTGATGCGCTGGCGGGCGAGAGCACGCGTCCCCGAGTCAGCCGTCTCGGTGAAGGCACGCTGATTACATTGCGCTGTATAAACGgc                                                                                                            <  1:832676/150‑1 (MQ=255)
                                        aaTTATGTACACCATGATAGCGCCCAATGGCTGGCGACAACACCGCTGCTTCCCAATAACGTACGTGATGCGCTGGCGGGCGAGAGCACGCGTCCCCGAGTCAGCCGTCTCGGTGAAGGCACGCTGATTACATTGCGCTGTATAAACGgc                                                                                                            <  1:1509827/150‑1 (MQ=255)
                                         aTTATGTACACCATGATAGCGCCCAATGGCTGGCGACAACACCGCTGCTTCCCAATAACGTACGTGATGCGCTGGCGGGCGAGAGCACGCGTCCCCGAGTCAGCCGTCTCGGTGAAGGCACGCTGATTACATTGCGCTGTATAAACGgca                                                                                                           <  1:695413/150‑1 (MQ=255)
                                           tATGTACACCATGATAGCGCCCAATGGCTGGCGACAACACCGCTGCTTCCCAATAACGTACGTGATGCGCTGGCGGGCGAGAGCACGCGTCCCCGAGTCAGCCGTCTCGGTGAAGGCACGCTGATTACATTGCGCTGTATAAACGgcagc                                                                                                         >  2:405244/1‑150 (MQ=255)
                                                 cacCATGATAGCGCCCAATGGCTGGCGACAACACCGCTGCTTCCCAATAACGTACGTGATGCGCTGGCGGGCGAGAGCACGCGTCCCCGAGTCAGCCGTCTCGGTGAAGGCGCGCTGATTACATTGCGCTGTATAAACGGCAGCACCGAt                                                                                                   >  2:148785/1‑150 (MQ=255)
                                                 cacCATGATAGCGCCCAATGGCTGGCGACAACACCGCTGCTTCCCAATAACGTACGTGATGCGCTGGCGGGCGAGAGCACGCGTCCCCGAGTCAGCCGTCTCGGTGAAGGCACGCTGATTACATTGCGCTGTATAAACGGCAGCACCGAt                                                                                                   >  2:943127/1‑150 (MQ=255)
                                                 cacCATGATAGCGCCCAATGGCTGGCGACAACACCGCTGCTTCCCAATAACGTACGTGATGCGCTGGCGGGCGAGAGCACGCGTCCCCGAGTCAGCCGTCTCGGTGAAGGCACGCTGATTACATTGCGCTGTATAAACGGCAGCACCGAt                                                                                                   >  1:900532/1‑150 (MQ=255)
                                                 cacCATGATAGCGCCCAATGGCTGGCGACAACACCGCTGCTTCCCAATAACGTACGTGATGCGCTGGCGGGCGAGAGCACGCGTCCCCGAGTCAGCCGTCTCGGTGAAGGCACGCTGATTACATTGCGCTGTATAAACGGCAGCACCGAt                                                                                                   >  1:938406/1‑150 (MQ=255)
                                                   ccATGATAGCGCCCAATGGCTGGCGACAACACCGCTGCTTCCCAATAACGTACGTGATGCGCTGGCGGGCGAGAGCACGCGTCCCCGAGTCAGCCGTCTCGGTGAAGGCACGCTGATTACATTGCGCTGTATAAACGGCAGCACCGATGa                                                                                                 >  1:181489/1‑150 (MQ=255)
                                                    cATGATAGCGCCCAATGGCTGGCGACAACACCGCTGCTTCCCAATAACGTACGTGATGCGCTGGCGGGCGAGAGCACGCGTCCCCGAGTCAGCCGTCTCGGTGAAGGCACGCTGATTACATTGCGCTGTATAAACGGCAGCACCGATGaa                                                                                                >  2:1086939/1‑150 (MQ=255)
                                                      tGATAGCGCCCAATGGCTGGCGACAACACCGCTGCTTCCCAATAACGTACGTGATGCGCTGGCGGGCGAGAGCACGCGTCCCCGAGTCAGCCGTCTCGGTGAAGGCACGCTGATTACATTGCGCTGTATAAACGGCAGCACCGATGAACg                                                                                              <  1:507324/150‑1 (MQ=255)
                                                      tGATAGCGCACAATGGCTGGCTACAACACCGCTGCTTCCCAATCACGTACGTGATGCGCTGGCGGGCGAGAGCACGCGTCCACGAGTCAGCCGTCTCGGTGAAGGCACGCTGATTACATTGCGCTGTATAAACGGCAGCACCGATGAACg                                                                                              <  2:388040/150‑1 (MQ=255)
                                                          aGCGCCCAATGGCTGGCGACAAGACCGCTGCTTCCCAATAACGTACGTGATGCGCTGGCGGGCGAGAGCACGCGGCCCCGAGTCAGCCGTCTCGGTGAAGGCACGCTGATTACATTGCGCTGTATAAACGGCAGCACCGATGAACGcccc                                                                                          <  2:1108750/150‑1 (MQ=255)
                                                              cccAATGGCTGGCGACAACACCGCTGCTTCCCAATAACGTACGTGATGCGCTGGCGGGCGAGAGCACGCGTCCCCGAGTCAGCCGTCTCGGTGAAGGCACGCTGATTACATTGCGCTGTATAAACGGCAGCACCGATGAACGCCCCGATc                                                                                      >  1:805988/1‑150 (MQ=255)
                                                                 aaTGGCTGGCGACAACACCGCTGCTTCCCAATAACGTACGTGATGCGCTGGCGGGCGAGAGCACGCGTCCCCGAGTCAGCCGTCTCGGTGAAGGCACGCTGATTACATTGCGCTGTATAAACGGCAGCACCGATGAACGCCCCGATCAAc                                                                                   <  2:96738/150‑1 (MQ=255)
                                                                  aTGGCTGGCGACAACACCGCTGCTTCCCAATAACGTACGTGATGCGCTGGCGGGCGAGAGCACGCGTCCCCGAGTCAGCCGTCTCGGTGAAGGCACGCTGATTACATTGCGCTGTATAAACGGCAGCACCGATGAACGCCCCGATCAACt                                                                                  <  1:790597/150‑1 (MQ=255)
                                                                  aTGGCTGGCGACAACACCGCTGCTTCCCAATAACGTACGTGATGCGCTGGCGGGCGAGAGCACGCGTCCCCGAGTCAGCCGTCTCGGTGAAGGCACGCTGATTACATTGCGCTGTATAAACGGCAGCACCGATGAACGCCCCGATCAACt                                                                                  <  2:733343/150‑1 (MQ=255)
                                                                   tggctggcGACAACACCGCTGCTTCCCAATAACGTACGTGATGCGCTGGCGGGCGAGAGCACGCGTCCCCGAGTCAGCCGTCTCGGTGAAGGCACGCTGATTACATTGCGCTGTATAAACGGCAGCACCGATGAACGCCCCGATCAACTg                                                                                 >  1:61894/1‑150 (MQ=255)
                                                                           gACAACACAGCTGCTTCCCAATAACGTACGTGATGCGCTGGCGGGCGAGAGCACGCGTCCCCGAGTCAGCCGTCTCGGTGAAGGCACGCTGATTACATTGCGCTGTATAAACGGCAGCACCGATGAACGCCCCGATCAACTGGTCGCCAt                                                                         >  2:390111/1‑150 (MQ=255)
                                                                             caacaCCGCTGCTTCCCAATAACGTACGTGATGCGCTGGCGGGCGAGAGCACGCGTCCCCGAGTCAGCCGTCTCGGTGAAGGCACGCTGATTACATTGCGCTGTATAAACGGCAGCACCGATGAACGCCCCGATCAACTGGTCGCCATGc                                                                       >  1:815607/1‑150 (MQ=255)
                                                                              aacaCCGCTGCTTCCCAATAACGTACGTGATGCGCTGGCGGGCGAGAGCACGCGTCCCCGAGTCAGCCGTCTCGGTGAAGGCACGCTGATTACATTGCGCTGTATAAACGGCAGCACCGATGAACGCCCCGATCAACTGGTCGCCATGCg                                                                      <  2:61894/150‑1 (MQ=255)
                                                                                 acCGCTGCTTCCCAATAACGTACGTGATGCGCTGGCGGGCGAGAGCACGCGTCCCCGAGTCAGCCGTCTCGGTGAAGGCACGCTGATTACATTGCGCTGTATAAACGGCAGCACCGATGAACGCCCCGATCAACTGGTCGCCATGCgtgt                                                                   >  2:335583/1‑150 (MQ=255)
                                                                                   cGCTGCTTCCCAATAACGTACGTGATGCGCTGGCGGGCGAGAGCACGCGTCCCCGAGTCAGCCGTCTCGGTGAAGGCACGCTGATTACATTGCGCTGTATAAACGTCAGCACCGATGAACGCCCCGATCAACTGGTCGCCATGCGTGtat                                                                 <  2:1148343/150‑1 (MQ=255)
                                                                                   cGCTGCTTCCCAATAACGTACGTGATGCGCTGGCGGGCGAGAGCACGCGTCCCCGAGTCAGCCGTCTCGGTGAAGGCACGCTGATTACATTGCGCTGTATAAACGGCAGCACCGATGAACGCCCCGATCAACTGGTCGCCATGCGTGtat                                                                 <  1:535703/150‑1 (MQ=255)
                                                                                      tgctTCCCAATAACGTACGTTATGCGCTGGCGGGCGAGAGCACGCGTCCCCGAGTCAGCCGTCTCGGTGAAGGCACGCTGATTACATTGCGCTGTATAAACGGCAGCACCGATGAACGCCCCGATCAACTGGTCGCCATGCGTGtatata                                                              <  2:451404/150‑1 (MQ=255)
                                                                                       gctTCCCAATAACGTACGTGATGCGCTGGCGGGCGAGAGCACGCGTCCCCGAGTCAGCCGTCTCGGTGAAGGCACGCTGATTACATTGCGCTGTATAAACGGCAGCACCGATGAACGCCCCGATCAACTGGTCGCCATGCGTGTatatat                                                             >  2:735887/1‑150 (MQ=255)
                                                                                          tCCCAATAACGTACGTGATGCGCTGGCGGGCGAGAGCACGCGTCCCCGAGTCAGCCGTCTCGGTGAACGCACGCTGATTACATTGCGCTGTATAAACGGCAGCACCGATGAACGCCCCGATCAACTGGTGGCCATGCGTGTATATATGGa                                                          >  2:1432083/1‑150 (MQ=255)
                                                                                           ccaaaaaaCGTACGTGATGCGATGGCGGGCCAGAGCACGCGTCCCCGAGTCAGCCGTCTCGGTGAAGGCACGCTGAATACATTGCGCTCTATAAAAGGCAGCACCGATGAACGCCCCGATCAACTGGTCGCCATGCGTGTATATATGGAc                                                         >  2:610673/7‑150 (MQ=255)
                                                                                              aaTAACGTACGGGTTGCGCTGGCGGGCGAGAGCACGCGTCCCCGAGTCAGCCGTCTCGGTGAAGGCACGCTGATTACATTGCGCTGTATAAACGGCAGCACCGATGAACGCCCCGATCAACTGGTCGCCATGCGTGTATATATGGACggg                                                      <  2:238025/150‑1 (MQ=255)
                                                                                               aTAACGTACGTGATGCGCTGGCGGGCGAGAGCACGCGTCCCCGAGTCAGCCGTCTCGGTGAAGGCGCGCTGATTACATTGCGCTGTATAAACGGCAGCACCGATGAACGCCCCGATCAACTGGTCGCCATGCGTGTATATATGGACGGGc                                                     <  1:148785/150‑1 (MQ=255)
                                                                                                      acgtGATGCGCTGGCGGGCGAGAGCACGCGTCCCCGAGTCAGCCGTCTCGGTGAAGGCACGCTGATTACATTGCGCTGTATAAACGGCAGCACCGATGAACGCCCCGATCAACTGGTCGCCATGCGTGTATATATGGACGGGCGGTTAAt                                              >  1:87818/1‑150 (MQ=255)
                                                                                                       cgtGATGCGCTGGCGGGCGAGAGCACGCGTCCCCGAGTCAGCCGTCTCGGTGAAGGCACGCTGATTACATTGCGCTGTATAAACGGCAGCACCGATGAACGCCCCGATCAACTGGTCGCCATGCGTGTATATATGGACGGGCGGTTAAtt                                             <  1:994427/150‑1 (MQ=255)
                                                                                                         tGATGCGCTGGCGGGCGAGAGCACGCGTCCCCGAGTCAGCCGTCTCGGTGAAGGCACGCTGATTACATTGCGCTGTATAAACGGCAGCACCGATGAACGCCCCGATCAACTGGTCGCCATGCGTGTATATATGGACGGGCGGTTAATTGt                                           <  1:405244/150‑1 (MQ=255)
                                                                                                            tGCGCTGGCGGGCGAGAGCACGCGTCCCCGAGTCAGCCGTCTCGGTGAAGGCACGCTGATTACATTGCGCTGTATAAACGGCAGCACCGATGAACGCCCCGATCAACTGGTCGCCATGCGTGTATATATGGACGGGCGGTTAATTGTTTc                                        <  1:943127/150‑1 (MQ=255)
                                                                                                             gcgcTGGCGGGCGAGAGCACGCGTCCCCGAGTCAGCCGTCTCGGTGAAGGCACGCTGATTACATTGCGCTGTATAAACGGCAGCACCGATGAACGCCCCGATCAACTGGTCGCCATGCGTGTATATATGGACGGGCGGTTAATTGTTTCg                                       <  1:390111/150‑1 (MQ=255)
                                                                                                                  ggcgggcgAGAGCACGCGTCCCCGAGTCAGCCGTCTCGGTGAAGGCACGCTGATTACATTGCGCTGTATAAACGGCAGCACCGATGAACGCCCCGATCAACTGGTCGCCATGCGTGTATATATGGACGGGCGGTTAATTGTTTCGACCCg                                  >  1:694032/1‑150 (MQ=255)
                                                                                                                        cgAGAGCACGCGTCCCCGAGTCAGCCGTCTCGGTGAAGGCACGCTGATTACATTGCGCTGTATAAACGGCAGCACCGATGAACGCCCCGATCAACTGGTCGCCATGCGTGTATATATGGTCGGTCGGTTAATTGTTTCTACCCGACAACg                            <  2:1519742/150‑1 (MQ=255)
                                                                                                                             gCACGCGTCCCCGAGTCAGCCGTCTCGGTGAAGGCACGCTGATTACATTGCGCTGTATAAACGGCAGCACCGATGAACGCCCCGATCAACTGGTCGCCATGCGTGTATATATGGACGGGCGGTTAATTGTTTCGACCCGACAACGCAAAg                       >  1:324436/1‑150 (MQ=255)
                                                                                                                              cACGCGTCCCCGAGTCAGCCGTCTCGGTGAAGGCACGCTGATTACATTGCGCTGTATAAACGGCAGCACCGATGAACGCCCCGATCAACTGGTCGCCATGCGTGTATATATGGACGGGCGGTTAATTGTTTCGACCCGACAACGCAAAGt                      <  1:956031/150‑1 (MQ=255)
                                                                                                                                 gcgTCCCCGAGTCAGCCGTCTCGGTGAAGGCACGCTGATTACATTGCGCTGTATAAACGGCAGCACCGATGAACGCCCCGATCAACTGGTCGCCATGCGTGTATATATGGACGGGCGGTTAATTGTTTCGACCCGACAACGCAAAGTGCt                   <  1:187697/150‑1 (MQ=255)
                                                                                                                                 ccttcaccGAGACGGCTGACTCGGTGAAGGCACGCTGATTACATTGCGCTGTATAAACGGCAGCACCGATGAACGCCCCGATCAACTGGTCGCCATGCGTGTATATATGGACGGGCGGTTAATTGTTTCGACCCGACAACGCAAAGTGCt                   >  2:256257/7‑150 (MQ=255)
                                                                                                                                 ccttcaccGAGACGGCTGACTCGGTGAAGGCACGCTGATTACATTGCGCTGTATAAACGGCAGCACCGATGAACGCCCCGATCAACTGGTCGCCATGCGTGTATATATGGACGGGCGGTTAATTGTTTCGACCCGACAACGCAAAGTGCt                   >  1:874920/7‑150 (MQ=255)
                                                                                                                                 ccttcaccGAGACGGCTGACTCGGTGAAGGCACGCTGATTACATTGCGCTGTATAAACGGCAGCACCGATGAACGCCCCGATCAACTGGTCGCCATGCGTGTATATATGGACGGGCGGTTAATTGTTTCGACCCGACAACGCAAAGTGCt                   >  2:1362657/7‑150 (MQ=255)
                                                                                                                                   ttcaccGAGACGGCTGACTCGGTGAAGGCACGCTGATTACATTGCGCTGTATAAACGGCAGCACCGATGAACGCCCCGATCAACTGGTCGCCATGCGTGTATATATGGACGGGCGGTTAATTGTTTCGACCCGACAACGCAAAGTGCTgg                 <  1:1362657/146‑1 (MQ=255)
                                                                                                                                   ttcaccGAGACGGCTGACTCGGTGAAGGCACGCTGATTACATTGCGCTGTATAAACGGCAGCACCGATGAACGCCCCGATCAACTGGTCGCCATGCGTGTATATATGGACGGGCGGTTAATTGTTTCGACCCGACAACGCAAAGTGCTgg                 <  1:256257/146‑1 (MQ=255)
                                                                                                                                    tcaccGAGACGGCTGACTCGGTGAAGGCACGCTGATTACATTGCGCTGTATAAACGGCAGCACCGATGAACGCCCCGATCAACTGGTCGCCATGCGTGTATATATGGACGGGCGGTTAATTGTTTCGACCCGACAACGCAAAGTGCTGgc                >  2:690565/4‑150 (MQ=255)
                                                                                                                                    tCCCCGAGTCAGCCGTCTCGGTGAAGGCACGCTGATTACATTGCGCTGTATAAACGGCAGCACCGATGAACGCCCCGATCAACTGGTCGCCATGCGTGTATATATGGTCGGGCGGTTAATTGTTTCGACCCGACAACGCAAAGTGCTGgc                <  2:938406/150‑1 (MQ=255)
                                                                                                                                     ccccGAGTCAGCCGTCTCGGTGAAGGCACGCTGATTACATTGCGCTGTATAAACGGCAGCACCGATGAACGCCCCGATCAACTGGTCGCCATGCGTGTATATATGGACGGGCGGTTAATTGTTTCGACCCGACAACGCAAAGTGCTGgcg               >  2:544116/1‑150 (MQ=255)
                                                                                                                                        cGAGTCAGCCGTCTCGGTGAAGGCACGCTGATTACATTGCGCTGTATAAACGGCAGCACCGATGAACGCCCCGATCAACTGGTCGCCATGCGTGTATATATGGACGGGCGGTTAATTGTTTCGACCCGACAACGCAAAGTGCTGGCGCTg            <  1:735887/150‑1 (MQ=255)
                                                                                                                                             cAGCCGTCTCGGTGAAGGCACGCTGATTACATTGCGCTGTATAAACGGCAGCACCGATGAACGCCCCGATCAACTGGTCGCCATGCGTGTATATATGGACGGGCGGTTAATTGTTTCGACCCGACAACGCAAAGTGCTGGCGCTGGACGa       >  2:76508/1‑150 (MQ=255)
                                                                                                                                                ccGTCTCGGTGAAGGCACGCTGATTACATTGCGCTGTATAAACGGCAGCACCGATGAACGCCCCGATCAACTGGTCGCCATGCGTGTATATATGGACGGGCGGTTAATTGTTTCGACCCGACAACGCAAAGTGCTGGCGCTGGACGAtgt    <  1:232783/150‑1 (MQ=255)
                                                                                                                                                  gTCTCGGTGAAGGCACGCTGATTACATTGCGCTGTATAAACGGCAGCACCGATGAACGCCCCGATCAACTGGTCGCCATGCGTGTATATATGGACGGGCGGTTAATTGTTTCGACCCGACAACGCAAAGTGCTGGCGCTGGACGATgtgg  <  1:1086939/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
AGATGTGATTGATGAAGCGCATCCCTGTTGGCTCCACCTTAATTATGTACACCATGATAGCGCCCAATGGCTGGCGACAACACCGCTGCTTCCCAATAACGTACGTGATGCGCTGGCGGGCGAGAGCACGCGTCCCCGAGTCAGCCGTCTCGGTGAAGGCACGCTGATTACATTGCGCTGTATAAACGGCAGCACCGATGAACGCCCCGATCAACTGGTCGCCATGCGTGTATATATGGACGGGCGGTTAATTGTTTCGACCCGACAACGCAAAGTGCTGGCGCTGGACGATGTGG  >  NC_000913/1408145‑1408440

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