Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,408,537 A→T 6.1% E163V (GAG→GTG)  zntB → Zn(2(+)):H(+) symporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,408,5370AT 6.1% 281.0 / 5.7 95E163V (GAG→GTG) zntBZn(2(+)):H(+) symporter
Reads supporting (aligned to +/- strand):  ref base A (47/42);  new base T (4/2);  total (51/44)
Fisher's exact test for biased strand distribution p-value = 6.83e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.86e-01

GGTTAATTGTTTCGACCCGACAACGCAAAGTGCTGGCGCTGGACGATGTGGTGAGCGATCTGGAAGAGGGCACGGGTCCGACCGATTGCGGGGGATGGCTGGTGGATGTGTGCGATGCGTTGACCGATCATTCCAGTGAATTTATCGAGCAGCTGCACGATAAAATTATCGACCTTGAAGATAATCTCCTTGATCAGCAAATTCCACCGCGTGGATTCCTGGCTCTGCTGCGCAAACAATTAATTGTGATGCGTCGCTATATGGCACCGCAACGTGATGTTTATGCTCGTCTTGCC  >  NC_000913/1408390‑1408685
                                                                                                                                                   |                                                                                                                                                    
ggTTAATTGTTTCGACCCGACAACGCAAAGTGCTGGCGCTGGACGATGTGGTGAGCGATCTGGAAGAGGGCACGGGTCCGACCGATTGCGGGGGATGGCTGGTGGATGTGTGCGATGCGTTGACCGATCATTCCAGTGAATTTATCGagc                                                                                                                                                    >  1:384314/1‑150 (MQ=255)
   tAATTGTTTCGCCCCGACAACGCAAAGTGCTGGCGCTGGACGATGTGGTGAGCGATCTGGAAGAGGGCACGGGTCCGACCGATTGCGGGGGATGGCTGGTGGATGTGTGCGATGCGTTGACCGATCATTCCAGTGAATTTATCGagcagc                                                                                                                                                 >  2:510194/1‑150 (MQ=255)
    aaTTGTTTCGACCCGACAACGCAAAGTGCTGGCGCTGGACGATGTGGTGAGCGATCTGGAAGAGGGCACGGGTCCGACCGATTGCGGGGGATGGCTGGTGGATGTGTGCGATGCGTTGACCGATCATTCCAGTGAATTTATCGAGCAGCt                                                                                                                                                >  1:79901/1‑150 (MQ=255)
    aaTTGTTTCGACCCGACAACGCAAAGTGCTGGCGCTGGACGATGTGGTGAGCGATCTGGAAGAGGGCACGGGTCCGACCGATTGCGGGGGATGGCTGGTGGATGTGTGCGATGCGTTCACCGATCATTCCAGTGAATTTATCGAGCAGCt                                                                                                                                                >  1:1336774/1‑150 (MQ=255)
     aTTGTTTCGACCCGACAACGCAAAGTGCTGGCGCTGGACGATGTGGTGAGCGATCTGGAAGAGGGCACGGGTCCGACCGATTGCGGGGGATGGCTGGTGGATGTGTGCGATGCGTTGACCGATCATTCCAGTGAATTTATCGAGCAGCTg                                                                                                                                               >  1:233016/1‑150 (MQ=255)
       tGTTTCGACCCGACAACGCAAAGTGCTGGCGCTGGACGATGTGGTGAGCGATCTGGAAGAGGGCACGGGTCCGACCGATTGCGGGGGATGGCTGGTGGATGTGTGCGATGCGTTGACCGATCATTCCAGTGAATTTATCGAGCAGCTGCa                                                                                                                                             <  1:158335/150‑1 (MQ=255)
        gTTTCGACCCGACAACGCAAAGTGCTGGCGCTGGACGATGTGGTGAGCGATCTGGAAGAGGGCACGGGTCCGACCGATTGCGGGGGATGGCTGGTGGATGTGTGCGATGCGTTGACCGATCATTCCAGTGAATTTATCGAGCAGCTGCAc                                                                                                                                            <  2:233016/150‑1 (MQ=255)
        gTTTCGACCCGACAACGCAAAGTGCTGGCGCTGGACGATGTGGTGAGCGATCTGGAAGAGGGCACGGGTCCGACCGATTGCGGGGGATGGCTGGTGGATGTGTGCGATGCGTTGACCGATCATTCCAGTGAATTTATCGAGCAGCTGCAc                                                                                                                                            >  1:188469/1‑150 (MQ=255)
        gTTTCGACCCGACAACGCAAAGTGCTGGCGCTGGACGATGTGGTGAGCGAGCTGGAAGAGGGCACGGGTCCGACCGATTGCGGGGGATGGCTGGTTTATGTGTGCGATGCGTTGACCGATCATTCCAGTGAATTTATCGAGCAGCTGCAc                                                                                                                                            <  2:181489/150‑1 (MQ=255)
         tttCGACCCGACAAGGCAAAGTGCTGGCGCTGGACGATGTGGTGAGCGATCTGGAAGAGGGCACGGGTCCGACCGATTGCGGGGGATGGCTGGTGGATGTGTGCGATGCGTTGACCGATCATTCCAGTGAATTTATCGAGCAGCTGCACg                                                                                                                                           >  2:1241642/1‑150 (MQ=255)
           tCGACCCGACAACGCAAAGTGCTGGCGCTGGACGATGTGGTGAGCGATCTGGAAGAGGGCACGGGTCCGACCGATTGCGGGGGATGGCTGGTGGATGTGTGCGATGCGTTGACCGATCATTCCAGTGAATTTATCGAGCAGCTGCACGAt                                                                                                                                         >  2:118915/1‑150 (MQ=255)
            cGACCCGACAACGCAAAGTGCTGGCGCTGGACGATGTGGTGAGCGATCTGGAAGAGGGCACGGGTCCGACCGATTGCGGGGGATGGCTGGTGGATGTGTGCGATGCGTTGACCGATCATTCCAGTGAATTTATCGAGCAGCTGCACGATa                                                                                                                                        <  2:874920/150‑1 (MQ=255)
              aCCCGACAACGCAAAGTGCTGGCGCTGGACGATGTGGTGAGCGATCTGGAAGAGGGCACGGGTCCGACCGATTGCGGGGGATTGCTGGTGGATGTGTGCGATGCGTTGACCGATCATTCCAGTGAATTTATCGAGCAGCTGCACGATaaa                                                                                                                                      <  1:1432083/150‑1 (MQ=255)
               cccGACAACTCAACGGTCTGGCGCTGGACGATGTTTTGAGCGATCTTGAACAGGGCAGGGTTCCGACCGATTGCGGGGGATGGCTGGTGGATGTGTGCGATGCGTTGACCGATCATTCCAGTGAATTTATCGTGCAGCTGCTCGATaaat                                                                                                                                     <  1:55361/150‑2 (MQ=255)
                ccGACAACGCAAAGTGCTGGCGCTGGACGATGTGGTGAGCGATCTGGAAGAGGGCACGGGTCCGACCGATTGCGGGGGATGGCTGGTGGATGTGTGCGATGCGTTGACCGATCATTCCAGTGAATTTATCGAGCAGCTGCACGATAAAAt                                                                                                                                    <  1:623605/150‑1 (MQ=255)
                 cGACAACGCAAAGTGCTGGCGCTGGACGATGTGGTGAGCGATCTGGAAGAGGGCACGGGTCCGACCGATTGCGGGGGATGGCTGGTGGATGTGTGCGATGCGTTGACCGATCATTCCAGTGAATTTATCGAGCAGCTGCACGATAAAAtt                                                                                                                                   >  1:839937/1‑150 (MQ=255)
                  gACAACGCAAAGTGCTGGCGCTGGACGATGTGGTGAGCGATCTGGAAGAGGGCACGGGTCCGACCGATTGCGGGGGATGGCTGGTGGATGTGTGCGATGCGTTGACCGATCATTCCAGTGAATTTATCGTGCAGCTGCTCGATaaattca                                                                                                                                  <  2:1358048/150‑5 (MQ=255)
                  gACAACGCAAAGTGCTGGCGCTGGACGATGTGGTGAGCGATCTGGAAGAGGGCACGGGTCCGACCGATTGCGGGGGATGGCTGGTGGATGTGTGCGATGCGTTGACCGATCATTCCAGTGAATTTATCGAGCAGCTGCACGATAAAATTa                                                                                                                                  >  2:357511/1‑150 (MQ=255)
                    cAACGCAAAGTGCTGGCGCTGGACGATGTGGTGAGCGATCTGGAAGAGGGCACGGGTCCGACCGATTGCGGGGGATGGCTGGTGGATGTGTGCGATGCGTTGACCGATCATTCCAGTGAATTTATCGAGCAGCTGCACGATAAAATTATc                                                                                                                                >  2:1329979/1‑150 (MQ=255)
                    cAACGCAAAGTGCTGGCGCTGGACGATGTGGTGAGCGATCTGGAAGAGGGCACGGGTCCGACCGATTGCGGGGGATGGCTGGTGGATGTGTGCGATGCGTTGACCGATCATTCCAGTGAATTTATCGAGCAGCTGCACGATAAAATTATc                                                                                                                                >  2:1403521/1‑150 (MQ=255)
                     aaCGCAAAGTGCTGGCGCTGGACGATGTGGTGAGCGATCTGGAAGAGGGCACGGGTCCGACCGATTGCGGGGGATGGCTGGTGGATGTGTGCGATGCGTTGACCGATCATTCCAGTGAATTTATCGAGCAGCTGCACGATAAAATTATCg                                                                                                                               <  2:842960/150‑1 (MQ=255)
                     aaCGCAAAGTGCTGGCGCTGGACGATGTGGTGAGCGATCTGGAAGAGGGCACGGGTCCGACCGATTGCGGGGGATGGCTGGTGGATGTGTGCGATGCGTTGACCGATCATTCCAGTGAATTTATCGAGCAGCTGCACGATAAAATTATCg                                                                                                                               <  2:402451/150‑1 (MQ=255)
                       cGCAAAGTGCTGGCGCTGGACGATGTGGTGAGCGATCTGGAAGAGGGCACGGGTCCGACCGATTGCGGGGGATGGCTGGTGGATGTGTGCGATGCGTTGACCGATCATTCCAGTGAATTTATCGAGCAGCTGCACGATAAAATTATCGAc                                                                                                                             <  1:400840/150‑1 (MQ=255)
                         cAAAGTGCTGGCGCTGGACGATGTGGTGAGCGATCTGGAAGAGGGCACGGGTCCGACCGATTGCGGGGGATGGCTGGTGGATGTGTGCGATGCGTTGACCGATCATTCCAGTGAATTTATCGAGCAGCTGCACGATAAAATTATCGACCt                                                                                                                           >  1:782317/1‑150 (MQ=255)
                           aaGTGCTGGCGCTGGACGATGTGGTGAGCGATCTGGAAGAGGGCACGGGTCCGACCGATTGCGGGGGATGGCTGGTGGATGTGTGCGATGCGTTGACCGATCATTCCAGTGAATTTATCGAGCAGCTGCACGATAAAATTATCGACCTTg                                                                                                                         >  1:959758/1‑150 (MQ=255)
                            aGTGCTGGCGCTGGACGATGTGGTGAGCGATCTGGAAGAGGGCACGGGTCCGACCGATTGCGGGGGATGGCTGGTGGATGTGTGCGATGCGTTGACCGATCATTCCAGTGAATTTATCGAGCAGCTGCACGATAAAATTATCGACCTTGa                                                                                                                        <  2:384314/150‑1 (MQ=255)
                             gTGCTGGCGCTGGACGATGTGGTGAGCGATCTGGAAGAGGGCACGGGTCCGACCGATTGCGGGGGATGGCTGGTGGATGTGTGCGATGCGTTGACCGATCATTCCAGTGAATTTATCGAGCAGCTGCACGATAAAATTATCGACCTTGaa                                                                                                                       >  1:962821/1‑150 (MQ=255)
                                    cgcTGGACGATGTGGTGAGCGATCTGGAAGAGGGCACGGGTCCGACCGATTGCGGGGGATGGCTGGTGGATGTGTGCGATGCGTTGACCGATCATTCCAGTGAATTTATCGAGCAGCTGCACGATAAAATTATCGACCTTGAAGATAAtc                                                                                                                <  1:975610/150‑1 (MQ=255)
                                      cTGGACGATGTGGTGAGCGATCTGGAAGAGGGCACGGGTCCGACCGATTGCGGGGGATGGCTGGTGGATGTGTGCGATGCGTTGACCGATCATTCCAGTGAATTTATCGAGCAGCTGCACGATAAAATTATCGACCTTGAAGATAAtctc                                                                                                              <  2:1307846/150‑1 (MQ=255)
                                        ggACGATGTGGTGAGCGATCTGGAAGAGGGCACGGGTCCGACCGATTGCGGGGGATGGCTGGTGGATGTGTGCGATGCGTTGACCGATCATTCCAGTGAATTTATCGAGCAGCTGCACGATAAAATTATCGACCTTGAAGATAATCTCCt                                                                                                            >  2:744570/1‑150 (MQ=255)
                                        ggACGATGTGGTGAGCGATCTGGAAGAGGGCACGGGTCCGACCGATTGCGGGGGATGGCTGGTGGATGTGTGCGATGCGTTGACCGATCATTCCAGTGAATTTATCGAGCAGCTGCACGATAAAATTATCGACCTTGAAGATAATCTCCt                                                                                                            >  2:750159/1‑150 (MQ=255)
                                        ggACGATGTGGTGAGCGATCTGGAAGAGGGCACGGGTCCGACCGATTGCGGGGGATGGCTGGTGGATGTGTGCGATGCGTTGACCGATCATTCCAGTGAATTTATCGAGCAGCTGCACGATAAAATTATCGACCTTGAAGATAATCTCCt                                                                                                            >  2:1212197/1‑150 (MQ=255)
                                         gACGATGTGGTGAGCGATCTGGAAGAGGGCACGGGTCCGACCGATTGCGGGGGATGGCTGGTGGATGTGTGCGATGCGTTGACCGATCATTCCAGTGAATTTATCGAGCAGCTGCACGATAAAATTATCGACCTTGAAGATAATCTGCtt                                                                                                           <  2:289616/150‑1 (MQ=255)
                                         gACGATGTGGTGAGCGATCTGGAAGAGGGCACGGGTCCGACCGATTGCGGGGGATGGCTGGTGGATGTGTGCGATGCGTTGACCGATCATTCCAGTGAATTTATCGAGCAGCTGCACGATAAAATTATCGACCTTGAAGATAATCTCCtt                                                                                                           <  2:1002071/150‑1 (MQ=255)
                                         gACGATGTGGTGAGCGATCTGGAAGAGGGCACGGGTCCGACCGATTGCGGAGGATGGCTGGTGGATGTGTGCGATGCGTTGACCGATCATTCCAGTGAATTTATCGAGCAGCTGCACGATAAAATTATCGACCTTGAAGATAATCTCCtt                                                                                                           <  1:485224/150‑1 (MQ=255)
                                          aCGATGTGGTGAGCGATCTGGAAGAGGGCACGGGTCCGACCGATTGCGGGGGATGGCTGGTGGATGTGTGCGATGCGTTGACCGATCATTCCAGTGAATTTATCGAGCAGCTGCACGATAAAATTATCGACCTTGAAGATAATCTCCTTg                                                                                                          <  2:959758/150‑1 (MQ=255)
                                          aCGATGTGGTGAGCGATCTGGAAGAGGGCACGGGTCCGACCGATTGCGGGGGATGGCTGGTGGATGTGTGCGATGCGTTGACCGATCATTCCAGTGAATTTATCGAGCAGCTGCACGATAAAATTATCGACCTTGAAGATAATCTCCTTg                                                                                                          <  2:1056343/150‑1 (MQ=255)
                                            gATGTGGTGAGCGATCTGGAAGAGGGCACGGGTCCGACCGATTGCGGGGGATGGCTGGTGGATGTGTGCGATGCGTTGACCGATCATTCCAGTGAATTTATCGAGCAGCTGCACGATAAAATTATCGACCTTGAAGATAATCTCCTTGAt                                                                                                        <  1:690565/150‑1 (MQ=255)
                                             aTGTGGTGAGCGATCTGGAAGAGGGCACGGGTCCGACCGATTGCGGGGGATGGCTGGTGGATGTGTGCGATGCGTTGACCGATCATTCCAGTGAATTTATCGAGCAGCTGCACGATAAAATTATCGACCTTGAAGATAATCTCCTTGATc                                                                                                       <  2:525612/150‑1 (MQ=255)
                                             aTGTGGTGAGCGATCTGGAAGAGGGCACGGGTCCGACCGATTGCGGGGGATGGCTGGTGGATGTGTGCGATGCGTTGACCGATCATTCCAGTGAATTTATCGAGCAGCTGCACGATAAAATTATCGACCTTGAAGATAATCTCCTTGATc                                                                                                       <  2:1399994/150‑1 (MQ=255)
                                               gtggtgAGCGATCTGGAAGAGGGCACGGGTCCGACCGATTGCGGGGGATGGCTGGTGGATGTGTGCGATGCGTTGACCGATCATTCCAGTGAATTTATCGAGCAGCTGCACGATAAAATTATCGACCTTGAAGATAATCTCCTTGATCAg                                                                                                     <  1:76508/150‑1 (MQ=255)
                                                 ggtgAGCGATCTGGAAGAGGGCACGGGTCCGACCGATTGCGGGGGATGGCTGGTGGATGTGTGCGATGTGTTGACCGATCATTCCAGTGAATTTATCGAGCAGCTGCACGATAAAATTATCGACCTTGAAGATAATCTCCTTGATAAGCa                                                                                                   >  2:1178020/1‑150 (MQ=255)
                                                 ggtgAGCGATCTGGAAGAGGGCACGGGTCCGACCGATTGCGGGGGATGGCTGGTGGATGTGTGCGATGCGTTGACCGATCATTCCAGTGAATTTATCGAGCAGCTGCACGATAAAATTATCGACCTTGAAGATAATCTCCTTGATCAGCa                                                                                                   >  1:722259/1‑150 (MQ=255)
                                                 ggtgAGCGATCTGGAAGAGGGCACGGGTCCGACCGATTGCGGGGGATGGCTGGTGGATGTGTGCGATGCGTTGACCGATCATTCCAGTGAATTTATCGAGCAGCTGCACGATAAAATTATCGACCTTGAAGATAATCTCCTTGATCAGCa                                                                                                   >  2:1494853/1‑150 (MQ=255)
                                                  gtgAGCGATCTGGAAGAGGGCACGGGTCCGACCGATTGCGGGGGATGGCTGGTGGATGTGTGCGATGCGTTGACCGATCATTCCAGTGAATTTATCGAGCAGCTGCACGATAAAATTATCGACCTTGAAGATAATCTCCTTGATCAGCaa                                                                                                  >  1:1415473/1‑150 (MQ=255)
                                                    gAGCGATCTGGAAGAGGGCACGGGTCCGACCGATTGCGGGGGATGGCTGGTGGATGTGTGCGATGCGTTGACCGCTCATTCCAGTGAATTTATCGAGCAGCTGCACGATAAAATTATCGCCCTTGAAGATAATCTCCTTGATCAGCAAAt                                                                                                >  2:945035/1‑150 (MQ=255)
                                                     aGCGATCTGGAAGAGGGCACGGGTCCGACCGATTGCGGGGGATGGCTGGTGGATGTGTGCGATGCGTTTACCGATCATTCCAGTGAATTTATCGAGCAGCTGCACGATAAAATTATCGACCTTGAAGATAATCTCCTTGATCAGCAAAtt                                                                                               <  2:1336774/150‑1 (MQ=255)
                                                     aGCGATCTGGAAGAGGGCACGGGTCCGACCGATTGCGGGGGATGGCTGGTGGATGTGTGCGATGCGTTGACCGATCATTCCAGTGAATTTATCGAGCAGCTGCACGATAAAATTATCGACCTTGAAGATAATCTCCTTGATCAGCAAAtt                                                                                               >  1:594601/1‑150 (MQ=255)
                                                       cGATCTGGAAGAGGGCACGGGTCCGACCGATTGCGGGGGATGGCTGGTGGATGTGTGCGATGCGTTGACCGATCATTCCAGTGAATTTATCGAGCAGCTGCACGATAAAATTATCGACCTTGAAGATAATCTCCTTGATCAGCAAATTcc                                                                                             <  2:129885/150‑1 (MQ=255)
                                                       cGATCTGGAAGAGGGCACGGGTCCGACCGATTGCGGGGGATGGCTGGTGGATGTGTGCGATGCGTTGACCGATCATTCCAGTGAATTTATCGAGCAGCTGCACGATAAAATTATCGACCTTGAAGATAATCTCCTTGATCAGCAAATTcc                                                                                             <  2:962821/150‑1 (MQ=255)
                                                         aTCTGGAAGAGGGCACGGGTCCGACCGATTGCGGGGGATGGCTGGTGGATGTGTGCGATGCGTTGACCGATCATTCCAGTGAATTTATCGAGCAGCTGCACGATAAAATTATCGACCTTGAAGATAATCTCCTTGATCAGCAAATTCCAc                                                                                           >  1:873212/1‑150 (MQ=255)
                                                          tCTGGAAGAGGGCACGGGTCCGACCTCTTGCGGTGGCTGGCTGGTGGATGTGTGCGATGCGTTGACCGATCATTCCAGGGAATTTATCGAGCAGCTGCACGATAAAATTATCGCCCTTGTAGATAATCTCCTTGATCAGCAAATTCCAcc                                                                                          <  2:1559975/150‑1 (MQ=255)
                                                           cTGGAAGAGGGCACTGGTCCGACCGATTGCGGGGGATGGCTGGTGGATGTGTGCGATGCGTTGACCGATCATTCCAGTGAATTTATCGAGCAGCTGCACGATAAAATTATCGACCTTGAAGATAATCTCCTTGATCAGCAAATTCCACcg                                                                                         >  2:432762/1‑150 (MQ=255)
                                                           cTGGAAGAGGGCACGGGTCCGACCGATTGCGGGGGATGGCTGGTGGATGTGTGCGATGCGTTGACCGATCATTCCAGTGAATTTATCGAGCAGCTGCACGATAAAATTATCGACCTTGAAGATAATCTCCTTGATCAGCAAATTCCACcg                                                                                         >  2:1588329/1‑150 (MQ=255)
                                                             ggAAGAGGGCACGGGTCCGACCGATTGCGGGGGATGGCTGGTGGATGTGTGCGATGCGTTGACCGATCATTCCAGTGAATTTATCGAGCAGCTGCACGATAAAATTATCGACCTTGAAGATAATCTCCTTGATCAGCAAATTCCACcgcg                                                                                       <  2:594601/150‑1 (MQ=255)
                                                               aaGAGGGCACGGGTCCGACCGATTGCGGGGGATGGCTGGTGGATGTGTGCGATGCGTTGACCGATCATTCCAGTGAATTTATCGAGCAGCTGCACGATAAAATTATCGACCTTGAAGATAATCTCCTTGATCAGCAAATTCCACCGCGTg                                                                                     <  1:432762/150‑1 (MQ=255)
                                                                     gaaCGGGTCCGACCGATAGCGGGGGATGGCTGGTGGATGTGTGCGATGCGTTGACCGATCATTCCAGTGAATATATCGAGCAGCTGCACGATAAAATTATCGAGCTTGAAGATAATCTCCTTGATCAGCAAATTCCACCGCGAGGATTcc                                                                               >  2:1400107/3‑150 (MQ=255)
                                                                      cACGGGTCCGACCGATTGCGGGGGATGGCTGGTGGATGTGTGCGATGCGTTGACCGATCATTCCAGTGAATTTATCGAGCAGCTGCACGATAAAATTATCGACCTTGAAGATAATCTCCTTGATCAGCAAATTCCACCGCGTGGATTCCt                                                                              <  1:385687/150‑1 (MQ=255)
                                                                         gggTCCGACCGATTGCGGGGGATGGCTGGTGGATGTGTGCGATGCGTTGACCGATCATTCCAGTGAATTTATCGAGCAGCTGCACGATAAAATTATCGACCTTGAAGATAATCTCCTTGATCAGCAAATTCCACCGCGTGGATTCCTGGc                                                                           <  1:186470/150‑1 (MQ=255)
                                                                                    aTTGCGGGGGATGGCTGGTGGATGTGTGCGATGCGTTGACCGATCATTCCAGTGAATTTATCGAGCAGCTGCACGATAAAATTATCGACCTTGAAGATAATCTCCTTGATCAGCAAATTCCACCGCGTGGATTCCTGGCTCTGCTGCGCa                                                                >  1:1051355/1‑150 (MQ=255)
                                                                                    aTTGCGGGGGATGGCTGGTGGATGTGTGCGATGCGTTGACCGATCATTCCAGTGAATTTATCGAGCAGCTGCACGATAAAATTATCGACCTTGAAGATAATCTCCTTGATCAGCAAATTCCACCGCGTGGATTCCTGGCTCTGCTGCGCa                                                                <  1:1494853/150‑1 (MQ=255)
                                                                                      tGCGGGGGATGGCTGGTGGATGTGTGCGATGCGTTGACCGATCATTCCAGTGAATTTATCGAGCAGCTGCACGATAAAATTATCGACCTTGAAGATAATCTCCTTGATCAGCAAATTCCACCGCGTGGATTCCTGGCTCTGATGCGCaaa                                                              >  2:1029272/1‑150 (MQ=255)
                                                                                       gCGGGGGATGGCTGGTGGATGTGTGCGATGCGTTGACCGATCATTCCAGTGAATTTATCGAGCAGCTGCACGATAAAATTATCGACCTTGAAGATAATCTCCTTGATCAGCAAATTCCACCGCGTGGATTCCTGGCTCTGCTGCGCAAAc                                                             <  1:1366274/150‑1 (MQ=255)
                                                                                           gggATGGCTGGTGGATGTGTGCGATGCGTTGACCGATCATTCCAGTGAATTTATCGAGCAGCTGCACGATAAAATTATCGACCTTGAAGATAATCTCCTTGATCAGCAAATTCCACCGCGTGGATTCCTGGCTCTGCTGCGCAAACaatt                                                         <  2:188469/150‑1 (MQ=255)
                                                                                               tGGCTGGTGGATGTGTGCGATGCGTTGACCGATCATTCCAGTGAATTTATCGAGCAGCTGCACGATAAAATTATCGACCTTGAAGATAATCTCCTTGATCAGCAAATTCCACCGCGTGGATTCCTGGCTCTGCTGCGCAAACaattaatt                                                     >  2:1158007/1‑150 (MQ=255)
                                                                                                ggCTGGTGGATGTGTGCGATGCGTTGACCGATCATTCCAGTGAATTTATCGAGCAGCTGCACGATAAAATTATCGACCTTGAAGATAATCTCCTTGATCAGCAAATTCCACCGCGTGGATTCCTGGTTCTGCTGCGCAAACAATTAATtg                                                    >  2:1317622/1‑150 (MQ=255)
                                                                                                 gCTGGTGGATGTGTGCGATGCGTTGACCGATCATTCCAGTGAATTTATCGAGCAGCTGCACGATAAAATTATCGACCTTGAAGATAATCTCCTTGATCAGCAAATTCCACCGCGTGGATTCCTGGCTCTGCTGCGCAAACAATTAATtgt                                                   <  2:722259/150‑1 (MQ=255)
                                                                                                 gCTGGTGGATGTGTGCGATGCGTTGACCGATCATTCCAGTGAATTTATCGAGCAGCTGCACGATAAAATTATCGACCTTGAAGATAATCTCCTTGATCAGCAAATTCCACCGCGTGGATTCCTGGCTCTGCTGCGCAAACAATTAATtgt                                                   <  1:118915/150‑1 (MQ=255)
                                                                                                  cTGGGGGATGTGTGCGATGCGTTGACCGATCATTCCAGTGAATTTATCGAGCAGCTGCACGATAAAATTATCGACCTTGAAGATAATCTCCTTGATCAGCAAATTCCACCGCGTGGATTCCTGGCTCTGCTGCGCAAACAATTAATtgtg                                                  <  2:942503/150‑1 (MQ=255)
                                                                                                   tggtggATGTGTGCGATGCGTTGACCGATCATTCCAGTGAATTTATCGAGCAGCTGCACGATAAAATTATCGACCTTGAAGATAATCTCCTTGATCAGCAAATTCCACCGCGTGGATTCCTGGCTCTGCTGCGCAAACAATTAATTGTGa                                                 >  1:1306170/1‑150 (MQ=255)
                                                                                                    ggtggATGTGTGCGATGCATTGACCGATCATTCCAGTGAATTTATCGAGCAGCTGCACGATAAAATTATCGACCTTGAAGATAATCTCCTTGATCAGCAAATTCCACCGCGTGGATTCCTGGCTCTGCTGCGCAAACAATTAATTGTGAt                                                >  2:1067683/1‑150 (MQ=255)
                                                                                                          tgtgtgCGATGCGTTGACCGATCATTCCAGTGAATTTATCGAGCAGCTGCACGATAAAATTATCGACCTTGAAGATAATCTCCTTGATCAGCAAATTCCACCGCGTGGATTCCTGGCTCTGCTGCGCAAACAATTAATTGTGATGCGTCg                                          >  2:1509776/1‑150 (MQ=255)
                                                                                                           gtgtgCGATGCGTTGACCGATCATTCCAGTGAATTTATCGAGCAGCTGCACGATAAAATTATCGACCTTGAAGATAATCTCCTTGATCAGCAAATTCCACCGCGTGGATTCCTGGCTCTGCTGCGCAAACAATTAATTGTGATGCGTCGc                                         >  1:558125/1‑150 (MQ=255)
                                                                                                            tgtgCGATGCGTTGACCGATCATTCCAGTGAATTTATCGAGCAGCTGCACGATAAAATTATCGACCTTGAAGATAATCTCCTTGATCAGCAAATTCCACCGCGTGGATTCCTGGCTCTGCTGCGCAAACAATTAATTGTGATGCGTCGCt                                        >  1:1108043/1‑150 (MQ=255)
                                                                                                             gtgCGATGCGTTGACCGATCATTCCAGTGAATTTATCGAGCAGCTGCACGATAAAATTATCGACCTTGAAGATAATCTCCTTGATCAGCAAATTCCACCGCGTGGATTCCTGGCTCTGCTGCGCAAACAATTAATTGTGATGCGTCGCta                                       >  1:477019/1‑150 (MQ=255)
                                                                                                               gCGATGCGTTGACCGATCATTCCAGTGAATTTATCGAGCAGCTGCACGATAAAATTATCGACCTTGAAGATAATCTCCTTGATCAGCAAATTCCACCGCGTGGATTCCTGGCTCTGCTGCGCAAACAATTAATTGTGATGc                                              >  2:1389257/1‑141 (MQ=255)
                                                                                                               gCGATGCGTTGACCGATCATTCCAGTGAATTTATCGAGCAGCTGCACGATAAAATTATCGACCTTGAAGATAATCTCCTTGATCAGCAAATTCCACCGCGTGGATTCCTGGCTCTGCTGCGCAAACAATTAATTGTGATGc                                              <  1:1389257/141‑1 (MQ=255)
                                                                                                               gCGATGCGTTGACCGATCATTCCAGTGAATTTATCGAGCAGCTGCACGATAAAATTATCGACCTTGAAGATAATCTCCTTGATCAGCAAATTCCACCGCGTGGATTCCTGGCTCTGCTGCGCAAACAATTAATTGTGATGCGTCGCtata                                     >  1:545226/1‑150 (MQ=255)
                                                                                                               gCGATGCGTTGACCGATCATTCCAGTGAATTTATCGAGCAGCTGCACGATAAAATTATCGACCTTGAAGATAATCTCCTTGATCAGCAAATTCCACCGCGTGGATTCCTGGCTCTGCTGCGCAAACAATTAATTGTGATGCGTCGCtata                                     >  1:276444/1‑150 (MQ=255)
                                                                                                                    gCGTTGACCGATCATTCCAGTGAATTTATCGAGCAGCTGCACGATAAAATTATCGACCTTGAAGATAATCTCCTTGATCAGCAAATTCCACCGCGTTGATTCCTGGCTCTGCTGCGCAAACAATTAATTGTGATGCGTCGCTATATGGCa                                <  2:805988/150‑1 (MQ=255)
                                                                                                                     cGTTGACCGATCATTCCAGTGAATTTATCGAGCAGCTGCACGATAAAATTATCGACCTTGAAGATAATCTCCTTGATCAGCAAATTCCACCGCGTGGATTCCTGGCTCTGCTGCGCAAACAATTAATTGTGATGCGTCGCTATATGGCAc                               >  1:821318/1‑150 (MQ=255)
                                                                                                                        tGACCGATCATTCCAGTGAATTTATCGAGCAGCTGCACGATAAAATTATCGACCTTGAAGATAATCTCCTTGATCAGCAAATTCCACCGCGTTGATTCCTGGCTCTGCTGCGCAAACAATTAATTGTGATGCGTCGCTATATGGCACCGc                            <  2:839937/150‑1 (MQ=255)
                                                                                                                         gACCGATCATTCCAGTGAATTTATCGAGCAGCTGCACGATAAAATTATCGACCTTGAAGATAATCTCCTTGATCAGCAAATTCCACCGCGTGGATTCCTGGCTCTGCTGCACAAACAATTAATTGTGATGCGTCGCTATATGGCACCGCa                           >  1:1321053/1‑150 (MQ=255)
                                                                                                                         gACCGATCATTCCAGTGAATTTATCGAGCAGCTGCACGATAAAATTATCGACCTTGAAGATAATCTCCTTGATCAGCAAATTCCACCGCGTGGATTCCTGGCTCTGCTGAGCAAACAATTAATTGTGATGCGTCGCTATATGGCACCGCa                           <  2:873212/150‑1 (MQ=255)
                                                                                                                            cGATCATTCCAGTGAATTTATCGAGCAGCTGCACGATAAAATTATCGACCTTGAAGATAATCTCCTTGATCAGCAAATTCCACCGCGTGGATTCCTGGCTCTGCTGCGCAAACAATTAATTGTGATGCGTCGCTATATGGCACCGCAACg                        <  2:1306170/150‑1 (MQ=255)
                                                                                                                             gATCATTCCAGTGAATTTATCGAGCAGCTGCACGATAAAATTATCGACCTTGAAGATAATCTCCTTGATCAGCAAATTCCACCGCGTGGATTCCTGGCTCTGCTGCGCAAACAATTAATTGTGATGCGTAGCTATATGGCACCGCAACGt                       >  2:589095/1‑150 (MQ=255)
                                                                                                                                    ccAGTGAATTTATCGAGCAGCTGCACGATAAAATTATCGACCTTGAAGATAATCTCCTTGATCAGCAAATTCCACCGCGTGGATTCCTGGCTCTGCTGCGCAAACAATTAATTGTGATGCGTCGCTATATGGCACCGCAACGTGATGttt                >  1:36114/1‑150 (MQ=255)
                                                                                                                                     cAGTGAATTTATCGAGCAGCTGCACGATAAAATTATCGACCTTGAAGATAATCTCCTTGATCAGCAAATTCCACCGCGTGGATTCCTGGCTCTGCTGCGCAAACAATTAATTGTGATGCGTCGCTATATGGCACCGCAACGTGATGTTTa               <  1:357511/150‑1 (MQ=255)
                                                                                                                                      gataattttATCGTGCAGCTGCTCGATAAAATTATCGACCTTGAAGATAATCTCCTTGATCAGCAAATTCCACCGCGTGGATTCCTGGCTCTGCTGCGCAAACAATTAATTGTGATGCGTCGCTATATGGCACCGCAACGTGATGTTTAt              >  2:395062/7‑150 (MQ=255)
                                                                                                                                       gTGAATTTATCGAGCAGCTGCACGATAAAATTATCGACCTTGAAGATAATCTCCTTGATCAGCAAATTCCACCGCGTGGATTCCTGGCTCTGCTGCGCAAACAATTAATTGTGATGCGTCGCTATATGGCACCGCAACGTGATGTTTATg             <  1:1212197/150‑1 (MQ=255)
                                                                                                                                       gTGAATTTATCGAGCAGCTGCACGATAAAATTATCGACCTTGAAGATAATCTCCTTGATCAGCAAATTCCACCGCGTGGATTCCTGGCTCTGCTGCGCAAACAATTAATTGTGATGCGTCGCTATATGGCACCGCAACGTGATGTTTATg             >  1:1113908/1‑150 (MQ=255)
                                                                                                                                       ataattttATCGTGCAGCTGCTCGATAAAATTATCGACCTTGAAGATAATCTCCTTGATCAGCAAATTCCACCGCGTGGATTCCTGGCTCTGCTGCGCAAACAATTAATTGTGATGCGTCGCTATATGGCACCGCAACGTGATGTTTATg             >  1:412199/6‑150 (MQ=255)
                                                                                                                                       ataattttATCGTGCAGCTGCTCGATAAAATTATCGACCTTGAAGATAATCTCCTTGATCAGCAAATTCCACCGCGTGGATTCCTGGCTCTGCTGCGCAAACAATTAATTGTGATGCGTCGCTATATGGCACCGCAACGTGATGTTTATg             >  1:984515/6‑150 (MQ=255)
                                                                                                                                          aaTTTATCGAGCAGCTGCACGATAAAATTATCGACCTTGAAGATAATCTCCTTGATCAGCAAATTCCACCGCGTGGATTCCTGGCTCTGCTGCGCAAACAATTAATTGTGATGCGTCGCTATATGGCACCGCAACGTGATGTTTATGCTc          >  2:866913/1‑150 (MQ=255)
                                                                                                                                              tATCGAGCAGCTGCACGATAAAATTATCGACCTTGAAGATAATCTCCTTGATCAGCAAATTCCACCGCGTGGATTCCTGGCTCTGCTGCGCAAACAATTAATTGTGATGCGTCGCTATATGGCACCGCAACGTGATGTTTATGCTCGTCt      >  2:380991/1‑150 (MQ=255)
                                                                                                                                               aTCGTGCAGCTGCTCGATAAAATTATCGACCTTGAAGATAATCTCCTTGATCAGCAAATTCCACCGCGTGGATTCCTGGCTCTGCTGCGCAAACAATTAATTGTGATGCGTCGCTATATGGCACCGCAACGTGATGTTTATGCTCGTCtt     >  1:775422/1‑150 (MQ=255)
                                                                                                                                                  gAGCAGCTGCACGATAAAATTATCGACCTTGAAGATAATCTCCTTGATCAGCAAATTCCACCGCGTGGATTCCTGGCTCTGCTGCGCAAACAATTAATTGTGATGAGTCGCTATATGGCACCGCAACGTGATGTTTATGCTCGTCTTGcc  >  2:667131/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
GGTTAATTGTTTCGACCCGACAACGCAAAGTGCTGGCGCTGGACGATGTGGTGAGCGATCTGGAAGAGGGCACGGGTCCGACCGATTGCGGGGGATGGCTGGTGGATGTGTGCGATGCGTTGACCGATCATTCCAGTGAATTTATCGAGCAGCTGCACGATAAAATTATCGACCTTGAAGATAATCTCCTTGATCAGCAAATTCCACCGCGTGGATTCCTGGCTCTGCTGCGCAAACAATTAATTGTGATGCGTCGCTATATGGCACCGCAACGTGATGTTTATGCTCGTCTTGCC  >  NC_000913/1408390‑1408685

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