Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,408,285 T→A 8.8% V79D (GTC→GAC)  zntB → Zn(2(+)):H(+) symporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,408,2850TA 8.8% 221.2 / 11.1 80V79D (GTC→GAC) zntBZn(2(+)):H(+) symporter
Reads supporting (aligned to +/- strand):  ref base T (31/42);  new base A (4/3);  total (35/45)
Fisher's exact test for biased strand distribution p-value = 6.93e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.59e-01

GGAAAATACAGATGTGATTGATGAAGCGCATCCCTGTTGGCTCCACCTTAATTATGTACACCATGATAGCGCCCAATGGCTGGCGACAACACCGCTGCTTCCCAATAACGTACGTGATGCGCTGGCGGGCGAGAGCACGCGTCCCCGAGTCAGCCGTCTCGGTGAAGGCACGCTGATTACATTGCGCTGTATAAACGGCAGCACCGATGAACGCCCCGATCAACTGGTCGCCATGCGTGTATATATGGACGGGCGGTTAATTGTTTCGACCCGACAACGCAAAGTGCTGGCGCTG  >  NC_000913/1408136‑1408430
                                                                                                                                                     |                                                                                                                                                 
ggAAAATACAGATGTGATTGATGAAGCGCATCCCTGTTGGCTCCACCTTAATTATGTACACCATGATAGCGCCCAATGGCTGGCGACAACACCGCTGCTTCCCAATAACGTACGTGATGCGCTGGCGGGCGAGAGCACGCGTCCCCGAGt                                                                                                                                                   <  1:1285386/150‑1 (MQ=255)
ggAAAATACAGATGTGATTGATGAAGCGCATCCCTGTTGGCTCCACCTTAATTATGTACACCATGATAGCGCCCAATGGCTGGCGACAACACCGCTGCTTCCCAATAACGTACGTGATGCGCTGGCGGGCGAGAGCACGCGTCCCCGAGt                                                                                                                                                   <  1:1328793/150‑1 (MQ=255)
ggAAAATACAGATGTCATTGATGAAGCGCATCCCTGTTGGCTCCACCTTAATTATGTACACCATGATAGCGCCCAATGGCTGGCGACAACACCGCTGCTTCCCAATAACGTACGTGATGCGCTGGCGGGCGAGAGCACGCGTCCCCGAGt                                                                                                                                                   >  2:956031/1‑150 (MQ=255)
 gAAAATACAGATGTGATTGATGAAGCGCATCCCTGTTGGCTCCACCTTAATTATGTACACCATGATAGCGCCCAATGGCTGGCGACAACACCGCTGCTTCCCAATAACGTACGTGATGCGCTGGCGGGCGAGAGCACGCGTCCCCGAGTc                                                                                                                                                  <  1:302570/150‑1 (MQ=255)
   aaaTACAGATGTGATTGATGAAGCGCATCCCTGTTGGCTCCACCTTAATTATGTACACCATGATAGCGCCCAATGGCTGGCGACAACACCGCTGCTTCCCAATAACGTACGTGATGCGCTGGCGGGCGAGAGCACGCGTCCCCGAGTCAg                                                                                                                                                <  2:847030/150‑1 (MQ=255)
      tACAGATGTGATTGATGAAGCGCATCCCTGTTGGCTCCACCTTAATTATGTACACCATGATAGCGCCCAATGGCTGGCGACAACACCGCTGCTTCCCAATAACGTACGTGATGCGCTGGCGGGCGAGAGCACGCGTCCCCGAGTCAGCCg                                                                                                                                             >  1:1519742/1‑150 (MQ=255)
         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                                                                                                                        >  2:1006366/1‑150 (MQ=255)
                           gcATCCCTGTTGGCTCCACCTTAATTATGTACACCATGATAGCGCCCAATGGCTGGCGACAACACCGCTGCTTCCCAATAACGTACGTGATGCGCTGGCGGGCGAGAGCACGCGTCCCCGAGTCAGCCGTCTCGGTGAAGGCACGCTGAt                                                                                                                        >  1:733343/1‑150 (MQ=255)
                           gcATCCCTGTTGGCTCCACCTTAATTATGTACACCATGATAGCGCCCAATGGCTGGCGACAACACCGCTGCTTCCCAATAACGTACGTGATGCGCTGGCGGGCGAGAGCACGCGTCCCCGAGTCAGCCGTCTCGGTGAAGGCACGCTGAt                                                                                                                        <  1:555001/150‑1 (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:1509827/150‑1 (MQ=255)
                                                 aaTTATGTACACCATGATAGCGCCCAATGGCTGGCGACAACACCGCTGCTTCCCAATAACGTACGTGATGCGCTGGCGGGCGAGAGCACGCGTCCCCGAGTCAGCCGTCTCGGTGAAGGCACGCTGATTACATTGCGCTGTATAAACGgc                                                                                                  <  1:832676/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         >  1:874920/7‑150 (MQ=255)
                                                                                                                                          ccttcaccGAGACGGCTGACTCGGTGAAGGCACGCTGATTACATTGCGCTGTATAAACGGCAGCACCGATGAACGCCCCGATCAACTGGTCGCCATGCGTGTATATATGGACGGGCGGTTAATTGTTTCGACCCGACAACGCAAAGTGCt         >  2:256257/7‑150 (MQ=255)
                                                                                                                                          ccttcaccGAGACGGCTGACTCGGTGAAGGCACGCTGATTACATTGCGCTGTATAAACGGCAGCACCGATGAACGCCCCGATCAACTGGTCGCCATGCGTGTATATATGGACGGGCGGTTAATTGTTTCGACCCGACAACGCAAAGTGCt         >  2:1362657/7‑150 (MQ=255)
                                                                                                                                            ttcaccGAGACGGCTGACTCGGTGAAGGCACGCTGATTACATTGCGCTGTATAAACGGCAGCACCGATGAACGCCCCGATCAACTGGTCGCCATGCGTGTATATATGGACGGGCGGTTAATTGTTTCGACCCGACAACGCAAAGTGCTgg       <  1:256257/146‑1 (MQ=255)
                                                                                                                                            ttcaccGAGACGGCTGACTCGGTGAAGGCACGCTGATTACATTGCGCTGTATAAACGGCAGCACCGATGAACGCCCCGATCAACTGGTCGCCATGCGTGTATATATGGACGGGCGGTTAATTGTTTCGACCCGACAACGCAAAGTGCTgg       <  1:1362657/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)
                                                                                                                                                     |                                                                                                                                                 
GGAAAATACAGATGTGATTGATGAAGCGCATCCCTGTTGGCTCCACCTTAATTATGTACACCATGATAGCGCCCAATGGCTGGCGACAACACCGCTGCTTCCCAATAACGTACGTGATGCGCTGGCGGGCGAGAGCACGCGTCCCCGAGTCAGCCGTCTCGGTGAAGGCACGCTGATTACATTGCGCTGTATAAACGGCAGCACCGATGAACGCCCCGATCAACTGGTCGCCATGCGTGTATATATGGACGGGCGGTTAATTGTTTCGACCCGACAACGCAAAGTGCTGGCGCTG  >  NC_000913/1408136‑1408430

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