Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 963,764 G→A 6.6% intergenic (+96/‑64) rpsA → / → ihfB 30S ribosomal subunit protein S1/integration host factor (IHF), DNA‑binding protein, beta subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913963,7640GA 6.6% 412.4 / 12.1 121intergenic (+96/‑64)rpsA/ihfB30S ribosomal subunit protein S1/integration host factor (IHF), DNA‑binding protein, beta subunit
Reads supporting (aligned to +/- strand):  ref base G (65/48);  new base A (2/6);  total (67/54)
Fisher's exact test for biased strand distribution p-value = 1.37e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 123 total aligned reads displayed.

AACTTCTCCAACAACGCAATGGCTGAAGCTTTCAAAGCAGCTAAAGGCGAGTAATTCTCTGACTCTTCGGGATTTTTATTCCGAAGTTTGTTGAGTTTACTTGACAGATTGCAGGTTTCGTCCTGTAATCAAGCACTAAGGGCGGCTACGGCCGCCCTTAATCAATGCAGCAACAGCAGCCGCTTAATTTGCCTTTAAGGAACCGGAGGAATCATGACCAAGTCAGAATTGATAGAAAGACTTGCCACCCAGCAATCGCACATTCCCGCCAAGACGGTTGAAGATGCAGTAAAAGAGA  >  NC_000913/963615‑963912
                                                                                                                                                     |                                                                                                                                                    
 aCTTCTCCAACAACGCAATGGCTGAAGCTTTCAAAGCAGCTAAAGGCGAGTAATTCTCTGACTCTTCGGGATTTTTATTCCGAAGTTTGTTGAGTTTACTTGACAGATTGCAGGTTTCGTCCTGTAATCAAGCACTAAGGGCGGCTACgg                                                                                                                                                     <  2:1687960/150‑1 (MQ=255)
   ttCTCCAACAACGCAATGGCTGAAGCTTTCAAAGCAGCTAAAGGCGAGTAATTCTCTGACTCTTCGGGATTTTTATTCCGAAGTTTGTTGAGTTTACTTGACAGATTGCAGGTTTCGTCCTGTAATCAAGCACTAAGGGCGGCTACGgcc                                                                                                                                                   >  1:1664013/1‑150 (MQ=255)
     ctcCAACAACGCAATGGCTGAAGCTTTCAAAGCAGCTAAAGGCGAGTAATTCTCTGACTCTTCGGGATTTTTATTCCGAAGTTTGTTGAGTTTACTTGACAGATTGCAGGTTTCGTCCTGTAATCAAGCACTAAGGGCGGCTACGgccgc                                                                                                                                                 <  2:422204/150‑1 (MQ=255)
      tcCAACAACGCAATGGCTGAAGCTTTCAAAGCAGCTAAAGGCGAGTAATTCTCTGACTCTTCGGGATTTTTATTCCGAAGTTTGTTGAGTTTACTTGACAGATTGCAGGTTTCGTCCTGTAATCAAGCACTAAGGGCGGCTACGgccgcc                                                                                                                                                >  1:587975/1‑150 (MQ=255)
      tcCAACAACGCAATGGCTGAAGCTTTCAAAGCAGCTAAAGGCGAGTAATTCTCTGACTCTTCGGGATTTTTATTCCGAAGTTTGTTGAGTTTACTTGACAGATTGCAGGTTTCGTCCTGTAATCAAGCACTAAGGGCGGCTACGgccgcc                                                                                                                                                >  2:375998/1‑150 (MQ=255)
        caacaaCGCAATGGCTGAAGCTTTCAAAGCAGCTAAAGGCGAGTAATTCTCTGACTCTTCGGGATTTTTATTCCGAAGTTTGTTGAGTTTACTTGACAGATTGCAGGTTTCGTCCTGTAATCAAGCACTAAGGGCGGCTACGGCCGCCCt                                                                                                                                              >  1:1341506/1‑150 (MQ=255)
        caacaaCGCAATGGCTGAAGCTTTCAAAGCAGCTAAAGGCGAGTAATTCTCTGACTCTTCGGGATTTTTATTCCGAAGTTTGTTGAGTTTACTTGACAGATTGCAGGTTTCGTCCTGTAATCAAGCACTAAGGGCGGCTACGGCCGCCCt                                                                                                                                              >  1:1427392/1‑150 (MQ=255)
        caacaaCGCAATGGCTGAAGCTTTCAAAGCAGCTAAAGGCGAGTAATTCTCTGACTCTTCGGGATTTTTATTCCGAAGTTTGTTGAGTTTACTTGACAGATTGCAGGTTTCGTCCTGTAATCAAGCACTAAGGGCGGCCGTAGCCGCCCt                                                                                                                                              <  2:664766/150‑1 (MQ=255)
        caacaaCGCAATGGCTGAAGCTTTCAAAGCAGCTAAAGGCGAGTAATTCTCTGACTCTTCGGGATTTTTATTCCGAAGTTTGTTGAGTTTACTTGACAGATTGCAGGTTTCGTCCTGTAATCAAGCACTAAGGGCGGCCGTAGCCGCCCt                                                                                                                                              <  2:664829/150‑1 (MQ=255)
          acaacGCAATGGCTGAAGCTTTCAAAGCAGCTAAAGGCGAGTAATTCTCTGACTCTTCGGGATTTTTATTCCGAAGTTTGTTGAGTTTACTTGACAGATTGCAGGTTTCGTCCTGTAATCAAGCACTAAGGGCGGCCGTAGCCGCCCTTa                                                                                                                                            <  1:799738/150‑1 (MQ=255)
            aacGCAATGGCTGAAGCTTTCAAAGCAGCTAAAGGCGAGTAATTCTCTGACTCTTCGGGATTTTTATTCCGAAGTTTGTTGAGTTTACTTGACAGATTGCAGGTTTCGTCCTGTAATCAAGCACTAAGGGCGGCTACGGCCGCCCTTAAt                                                                                                                                          >  2:349919/1‑150 (MQ=255)
            aacGCAATGGCTGAAGCTTTCAAAGCAGCTAAAGGCGAGTAATTCTCTGACTCTTCGGGATTTTTATTCCGAAGTTTGTTGAGTTTACTTGACAGATTGCAGGTTTCGTCCTGTAATCAAGCACTAAGGGCGGCTACGGCCGCCCTTAAt                                                                                                                                          >  2:349914/1‑150 (MQ=255)
               gCACTGGCTGAAGCTTTCAAAGCAGCTACAGGCGAGTAATTCTCTGACTCTTCGGGATTTTTATTCCGAAGTTTGTTGAGTTTACTTGACAGATTGCAGGTTTCGTCCTGTAATCAAGCACTAAGGGCGGCTACGGCCGCCCTTAATCaa                                                                                                                                       <  2:165022/150‑1 (MQ=255)
               gCAATGGCTGAAGCTTTCAAAGCAGCTAAAGGCGAGTAATTCTCTGACTCTTCGGGATTTTTATTCCGAAGTTTGTTGAGTTTACTTGACAGATTGCAGGTTTCGTCCTGTAATCAAGCATTAAGGGCGGCCGTAGCCGCCCTTagtgct                                                                                                                                       <  2:908099/150‑6 (MQ=255)
               gCAATGGCTGAAGCTTTCAAAGCAGCTAAAGGCGAGTAATTCTCTGACTCTTCGGGATTTTTATTCCGAAGTTTGTTGAGTTTACTTGACAGATTGCAGGTTTCGTCCTGTAATCAAGCATTAAGGGCGGCCGTAGCCGCCCTTagtgct                                                                                                                                       <  2:908096/150‑6 (MQ=255)
               gCAATGGCTGAAGCTTTCAAAGCAGCTAAAGGCGAGTAATTCTCTGACTCTTCGGGATTTTTATTCCGAAGTTTGTTGAGTTTACTTGACAGATTGCAGGTTTCGTCCTGTAATCAAGCACTAAGGGCGGCTACGGCCGCCCTTAATCaa                                                                                                                                       <  2:936760/150‑1 (MQ=255)
               gCAATGGCTGAAGCTTTCAAAGCAGCTAAAGGCGAGTAATTCTCTGACTCTTCGGGATTTTTATTCCGAAGTTTGTTGAGTTTACTTGACAGATTGCAGGTTTCGTCCTGTAATCAAGCACTAAGGGCGGCTACGGCCGCCCTTAATCaa                                                                                                                                       <  1:1423375/150‑1 (MQ=255)
                   tGGCTGAAGCTTTCAAAGCAGCTAAAGGCGAGTAATTCTCTGACTCTTCGGGATTTTTATTCCGAAGTTTGTTGAGTTTACTTGACAGATTGCAGGTTTCGTCCTGTAATCAAGCACTAAGGGCGGCTACGGCCGCCCTTAATCAATgca                                                                                                                                   <  1:1694059/150‑1 (MQ=255)
                     gCTGAAGCTTTCAAAGCAGCTAAAGGCGAGTAATTCTCTGACTCTTCGGGATTTTTATTCCGAAGTTTGTTGAGTTTACTTGACAGATTGCAGGTTTCGTCCTGTAATCAAGCACTAAGGGCGGCTACGGCCGCCCTTAATCAATgcagc                                                                                                                                 >  2:385369/1‑150 (MQ=255)
                       tGAAGCTTTCAAAGCAGCTAAAGGCGAGTAATTCTCTGACTCTTCGGGATTTTTATTCCGAAGTTTGTTGAGTTTACTTGACAGATTGCAGGTTTCGTCCTGTAATCAAGCACTAAGGGCGGCTACGGCCGCCCTTAATCAATGCAGCaa                                                                                                                               <  1:24706/150‑1 (MQ=255)
                        gAAGCTTTCAAAGCAGCTAAAGGCGAGTAATTCTCTGACTCTTCGGGATTTTTATTCCGAAGTTTGTTGAGTTTACTTGACAGATTGCAGGTTTCGTCCTGTAATCAAGCACTAAGGGCGGCTACGGCCGCCCTTAATCAATGCAGCAAc                                                                                                                              <  2:1437155/150‑1 (MQ=255)
                        gAAGCTTTCAAAGCAGCTAAAGGCGAGTAATTCTCTGACTCTTCGGGATTTTTATTCCGAAGTTTGTTGAGTTTACTTGACAGATTGCAGGTTTCGTCCTGTAATCAAGCACTAAGGGCGGCTACGGCCGCCCTTAATCAATGCAGCAAc                                                                                                                              >  2:67191/1‑150 (MQ=255)
                         aaGCTTTCAAAGCAGCTAAAGGCGAGTAATTCTCTGACTCTTCGGGATTTTTATTCCGAAGTTTGTTGAGTTTACTTGACAGATTGCAGGTTTCGTCCTGTAATCAAGCACTAAGGGCGGCTACGGCCGCCCTTAATCAATGCAGCAAca                                                                                                                             <  1:500121/150‑1 (MQ=255)
                           gCTTTCAAAGCAGCTAAAGGCGAGTAATTCTCTGACTCTTCGGGATTTTTATTCCGAAGTTTGTTGAGTTTACTTGACAGATTGCAGGTTTCGTCCTGTAATCAAGCACTAAGGGCGGCTACGGCCGCCCTTAATCAATGCAGCAAcagc                                                                                                                           <  1:865239/150‑1 (MQ=255)
                           gCTTTCAAAGCAGCTAAAGGCGAGTAATTCTCTGACTCTTCGGGATTTTTATTCCGAAGTTTGTTGAGTTTACTTGACAGATTGCAGGTTTCGTCCTGTAATCAAGCACTAAGGGCGGCTACGGCCGCCCTTAATCAATGCAGCAAcagc                                                                                                                           <  1:865255/150‑1 (MQ=255)
                             tttCAAAGCAGCTAAAGGCGAGTAATTCTCTGACTCTTCGGGATTTTTATTCCGAAGTTTGTTGAGTTTACTTGACAGATTGCAGGTTTCGTCCTGTAATCAAGCACTAAGGGCGGCTACGGCCGCCCTTAATCAATGCAGCAAcagcag                                                                                                                         >  2:695787/1‑150 (MQ=255)
                             tttCAAAGCAGCTAAAGGCGAGTAATTCTCTGACTCTTCGGGATTTTTATTCCGAAGTTTGTCGAGTTTACTTGACAGATTGCAGGTTTCGTCCTGTAATCAAGCACTAAGGGCGGCTACGGCCGCCCTTAATCAATGCAGCAAcagcag                                                                                                                         >  1:1064984/1‑150 (MQ=255)
                              ttCAAAGCAGCTAAAGGCGAGTAATTCTCTGACTCTTCGGGATTTTTATTCCGAAGTTTGTTGAGTTTACTTGACAGATTGCAGGTTTCGTCCTGTAATCAAGCACTAAGGGCGGCTACGGCCGCCCTTAATCAATGCAGCAACagcagc                                                                                                                        >  2:1395470/1‑150 (MQ=255)
                                  aaGCAGCTAAAGGCGAGTAATTCTCTGACTCTTCGGGATTTTTATTCCGAAGTTTGTTGAGTTTACTTGACAGATTGCAGGTTTCGTCCTGTAATCAAGCACTAAGGGCGGCTACGGCCGCCCTTAATCAATGCAGCAACAGCAGCCGCt                                                                                                                    >  2:1645998/1‑150 (MQ=255)
                                    gcagcTAAAGGCGAGTAATTCTCTGACTCTTCGGGATTTTTATTCCGAAGTTTGTTGAGTTTACTTGACAGATTGCAGGTTTCGTCCTGTAATCAAGCACTAAGGGCGGCTACGGCCGCCCTTAATCAATGCAGCAAcagcag                                                                                                                         >  1:121045/1‑143 (MQ=255)
                                    gcagcTAAAGGCGAGTAATTCTCTGACTCTTCGGGATTTTTATTCCGAAGTTTGTTGAGTTTACTTGACAGATTGCAGGTTTCGTCCTGTAATCAAGCACTAAGGGCGGCTACGGCCGCCCTTAATCAATGCAGCAAcagcag                                                                                                                         >  1:121057/1‑143 (MQ=255)
                                    gcagcTAAAGGCGAGTAATTCTCTGACTCTTCGGGATTTTTATTCCGAAGTTTGTTGAGTTTACTTGACAGATTGCAGGTTTCGTCCTGTAATCAAGCACTAAGGGCGGCTACGGCCGCCCTTAATCAATGCAGCAAcagcag                                                                                                                         <  2:1577001/143‑1 (MQ=255)
                                    gcagcTAAAGGCGAGTAATTCTCTGACTCTTCGGGATTTTTATTCCGAAGTTTGTTGAGTTTACTTGACAGATTGCAGGTTTCGTCCTGTAATCAAGCACTAAGGGCGGCTACGGCCGCCCTTAATCAATGCAGCAAcagcag                                                                                                                         <  2:121057/143‑1 (MQ=255)
                                    gcagcTAAAGGCGAGTAATTCTCTGACTCTTCGGGATTTTTATTCCGAAGTTTGTTGAGTTTACTTGACAGATTGCAGGTTTCGTCCTGTAATCAAGCACTAAGGGCGGCTACGGCCGCCCTTAATCAATGCAGCAAcagcag                                                                                                                         >  1:1577001/1‑143 (MQ=255)
                                     cagcTAAAGGCGAGTAATTCTCTGACTCTTCGGGATTTTTATTCCGAAGTTTGTTGAGTTTACTTGACAGATTGCAGGTTTCGTCCTGTAATCAAGCACTAAGGGCGGCTACGGCCGCCCTTAATCAATGCAGCAACAGCAGCCGCTTaa                                                                                                                 <  1:1093405/150‑1 (MQ=255)
                                      agcTAAAGGCGAGTAATTCTCTGACTCTTCGGGATTTTTATTCCGAAGTTTGTTGAGTTTACTTGACAGATTGCAGGTTTCGTCCTGTAATCAAGCACTAAGGGCGGCTACGGCCGCCCTTAATCAATGCAGCAACAGCAGCCGCTTAAt                                                                                                                <  1:341295/150‑1 (MQ=255)
                                         tAAAGGCGAGTAATTCTCTGACTCTTCGGGATTTTTATTCCGAAGTTTGTTGAGTTTACTTGACAGATTGCAGGTTTCGTCCTGTAATCAAGCACTAAGGGCGGCTACGGCCGCCCTTAATCAATGCAGCAACAGCAGCCGCTTAATTTg                                                                                                             >  1:750966/1‑150 (MQ=255)
                                              gCGAGTAATTCTCTGACTCTTCGGGATTTTTATTCCGAAGTTTGTTGAGTTTACTTGACAGATTGCAGGTTTCGTCCTGTAATCAAGCACTAAGGGCGGCTACGGCCGCCCTTAATCAATGCAGCAACAGCAGCCGCTTAATTTGCCttt                                                                                                        >  1:1175693/1‑150 (MQ=255)
                                              gCGAGTAATTCTCTGACTCTTCGGGATTTTTATTCCGAAGTTTGTTGAGTTTACTTGACAGATTGCAGGTTTCGTCCTGTAATCAAGCACTAAGGGCGGCTACGGCCGCCCTTAATCAATGCAGCAACAGCAGCCGCTTAATTTGCCttt                                                                                                        >  1:1175695/1‑150 (MQ=255)
                                               cGAGTAATTCTCTGACTCTTCGGGATTTTTATTCCGAAGTTTGTTGAGTTTACTTGACAGATTGCAGGTTTCGTCCTGTAATCAAGCACTAAGGGCGGCTACGGCCGCCCTTAATCAATGCAGCAACAGCAGCCGCTTAATTTGCCTTTa                                                                                                       <  1:1352765/150‑1 (MQ=255)
                                                gAGTAATTCTCTGACTCTTCGGGATTTTTATTCCGAAGTTTGTTGAGTTTACTTGACAGATTGCAGGTTTCGTCCTGTAATCAAGCACTAAGGGCGGCTACGGCCGCCCTTAATCAATGCAGCAACAGCAGCCGCTTAATTTGCCTTTaa                                                                                                      <  1:1527473/150‑1 (MQ=255)
                                                      ttCTCTGACTCTTCGGGATTTTTATTCCGAAGTTTGGTGAGTTTACTTGACAGATTGCAGGTTTCGTCCTGTAATCAAGCACTAAGGGCGGCTACGGCCGCCCTTAATCAATGCTGCACCAGCAGCCGCTTAATTTGCCTTTAAGGAccc                                                                                                >  2:115389/1‑150 (MQ=255)
                                                        ctctGACTCTTCGGGATTTTTATTCCGAAGTTTGTTGAGTTTACTTGACAGATTGCAGGTTTCGTCCTGTAATCAAGCACTAAGGGCGGCTACGGCCGCCCTTAATCAATGCAGCAACAGCAGCCGCTTAATTTGCCTTTAAGGAACCgg                                                                                              <  2:1373458/150‑1 (MQ=255)
                                                         tctGACTCTTCGGGATTTTTATTCCGAAGTTTGTTGAGTTTACTTGACAGATTGCAGGTTTCGTCCTGTAATCAAGCACTAAGGGCGGCTACGGCCGCCCTTAATCAATGCAGCAACAGCAGCCGCTTAATTTGCCTTTAAGGAACCgga                                                                                             >  1:456416/1‑150 (MQ=255)
                                                         tctGACTCTTCGGGATTTTTATTCCGAAGTTTGTTGAGTTTACTTGACAGATTGCAGGTTTCGTCCTGTAATCAAGCACTAAGGGCAGCTACGGCCGCCCTTAATCAATGCAGCAACAGCAGCCGCTTAATTTGCCTTTAAGGAACCgga                                                                                             >  1:246265/1‑150 (MQ=255)
                                                         tctGACTCTTCGGGATTTTTATTCCGAAGTTTGTTGAGTTTACATGACAGATTGCACGTTTCGTCCTGTAATCTAGCACTAAGGGCGGCTACGGCCGCCCTTAATCAATGCAGCAACAGCAGCCGCTTAATTTGCCTTTAAGGAACCgga                                                                                             >  1:1303849/1‑150 (MQ=255)
                                                           tGACTCTTCGGGATTTTTATTCCGAAGTTTGTTGAGTTTACTTGACAGATTGCAGGTTTCGTCCTGTAATCAAGCACTAAGGGCGGCTACGGCCGCCCTTAATCAATGCAGCAACAGCAGCCGCTTAATTTGCCTTTAAGGAACCggagg                                                                                           >  1:370390/1‑150 (MQ=255)
                                                           tGACTCTTCGGGATTTTTATTCCGAAGTTTGTTGAGTTTACTTGACAGATTGCAGGTTTCGTCCTGTAATCAAGCACTAAGGGCGGCTACGGCCGCCCTTAATCAATGCAGCAACAGCAGCCGCTTAATTTGCCTTTAAGGAACCggagg                                                                                           >  1:1486971/1‑150 (MQ=255)
                                                                 ttCGGGATTTTTATTCCGAAGTTTGTTGAGTTTACTTGACAGATTGCAGGTTTCGTCCTGTAATCAAGCACTAAGGGCGGCTACGGCCGCCCTTAATCAATGCAGCAACAGCAGCCGCTTAATTTGCCTTTAAGGAACCGGAGGAATCAt                                                                                     >  2:1269467/1‑150 (MQ=255)
                                                                 ttCGGGATTTTTATTCCGAAGTTTGTTGAGTTTACTTGACAGATTGCAGGTTTCGTCCTGTAATCAAGCACTAAGGGCGGCTACGGCCGCCCTTAATCAATGCAGCAACAGCAGCCGCTTAATTTGCCTTTAAGGAACCGGAGGAATCAt                                                                                     <  2:1643769/150‑1 (MQ=255)
                                                                 ttCGGGATTTTTATTCCGAAGTTTGTTGAGTTTACTTGACAGATTGCAGGTTTCGTCCTGTAATCAAGCACTAAGGGCGGCTACGGCCGCCCTTAATCAATGCAGCAACAGCAGCCGCTTAATTTGCCTTTAAGGAACCGGAGGAATCAt                                                                                     >  2:717295/1‑150 (MQ=255)
                                                                 ttCGGGATTTTTATTCCGAAGTTTGTTGAGTTTACTTGACAGATTGCAGGTTTCGTCCTGTAATCAAGCACTAAGGGCGGCTACGGCCGCCCTTAATCAATGCAGCAACAGCAGCCGCTTAATTTGCCTTTAAGGAACCGGAGGAATCAt                                                                                     <  2:790256/150‑1 (MQ=255)
                                                                 ttCGGGATTTTTATTCCGAAGTTTGTTGAGTTTACTTGACAGATTGCAGGTTTCGTCCTGTAATCAAGCACTAAGGGCGGCTACGGCCGCCCTTAATCAATGCAGCAACAGCAGCCGCTTAATTTGCCTTTAAGGAACCGGAGGAATCAt                                                                                     >  2:1181651/1‑150 (MQ=255)
                                                                 ttCGGGATTTTTATTCCGAAGTTTGTTGAGTTTACTTGACAGATTGCAGGTTTCGTCCTGTAATCAAGCACTAAGGGCGGCTACGGCCGCCCTTAATCAATGCAGCAACAGCAGCCGCTTAATTTGCCTTTAAGGAACCGGAGGAATCAt                                                                                     >  2:1692676/1‑150 (MQ=255)
                                                                 ttCGGGATTTTTATTCCGAAGTTTGTTGAGTTTACTTGACAGATTGCAGGTTTCGTCCTGTAATCAAGCACTAAGGGCGGCTACGGCCGCCCTTAATCAATGCAGCAACAGCAGCCGCTTAATTTGCCTTTAAGGAACCGGAGGAATCAt                                                                                     <  2:1431182/150‑1 (MQ=255)
                                                                    gggATTTTTATTCCGAAGTTTGTTGAGTTTACTTGACAGATTGCAGGTTTCGTCCTGTAATCAAGCACTAAGGGCGGCTACGGCCGCCCTTAATCAATGCAGCAACAGCAGCCGCTTAATTTGCCTTTAAGGAACCGGAGGAATCATGAc                                                                                  >  2:1611332/1‑150 (MQ=255)
                                                                    gggATTTTTATTCCGAAGTTTGTTGAGTTTACTTGACAGATTGCAGGTTTCGTCCTGTAATCAAGCACTAAGGGCGGCTACGGCCGCCCTTAATCAATGCAGCAACAGCAGCCGCTTAATTTGCCTTTAAGGAACCGGAGGAATCATGAc                                                                                  >  2:1611280/1‑150 (MQ=255)
                                                                           ttATTCCGAAGTTTGTTGAGTTTACTTGACAGATTGCAGGTTTCGTCCTGTAATCAAGCACTAAGGGCGGCTACGGCCGCCCTTAATCAATGCAGCAACAGCAGCCGCTTAATTTGCCTTTAAGGAACCGGAGGAATCATGACCAAGTCa                                                                           >  1:1279995/1‑150 (MQ=255)
                                                                                ccGAAGTTTGTTGAGTTTACTTGACAGATTGCAGGTTTCGTCCTGTAATCAAGCACTAAGGGCGGCTACGGCCGCCCTTAATCAATGCAGCAACAGCAGCCGCTTAATTTGCCTTTAAGGAACCGGAGGAATCATGACCAAGTCAGAAtt                                                                      <  1:651321/150‑1 (MQ=255)
                                                                                 cGAAGTTTGTTGAGTTTACTTGACAGATTGCAGGTTTCGTCCTGTAATCAAGCACTAAGGGCGGCTACGGCCGCCCTTAATCAATGCAGCAACAGCAGCCGCTTAATTTGCCTTTAAGGAACCGGAGGAATCATGACCAAGTCAGAATTg                                                                     >  1:1583357/1‑150 (MQ=255)
                                                                                  gAAGTTTGTTGAGTTTACTTGACAGATTGCAGGTTTCGTCCTGTAATCAAGCACTAAGGGCGGCTACGGCCGCCCTTAATCAATGCAGCAACAGCAGCCGCTTAATTTGCCTTTAAGGAACCGGAGGAATCATGACCAAGTCAGAATTGa                                                                    >  1:1694093/1‑150 (MQ=255)
                                                                                  gAAGTTTGTTGAGTTTACTTGACAGATTGCAGGTTTCGTCCTGTAATCAAGCACTAAGGGCGGCTACGGCCGCCCTTAATCAATGCAGCAACAGCAGCCGCTTAATTTGCCTTTAAGGAACCGGAGGAATCATGACCAAGTCAGAATTGa                                                                    >  1:39875/1‑150 (MQ=255)
                                                                                  gAAGTTTGTTGAGTTTACTTGACAGATTGCAGGTTTCGTCCTGTAATCAAGCACTAAGGGCGGCTACGGCCGCCCTTAATCAATGCAGCAACAGCAGCCGCTTAATTTGCCTTTAAGGAACCGGAGGAATCATGACCAAGTCAGAATTGa                                                                    >  2:123121/1‑150 (MQ=255)
                                                                                  gAAGTTTGTTGAGTTTACTTGACAGATTGCAGGTTTCGTCCTGTAATCAAGCACTAAGGGCGGCTACGGCCGCCCTTAATCAATGCAGCAACAGCAGCCGCTTAATTTGCCTTTAAGGAACCGGAGGAATCATGACCAAGTCAGAATTGa                                                                    >  1:1705082/1‑150 (MQ=255)
                                                                                       ttgttgAGTTTACTTGACAGATTGCAGGTTTCGTCCTGTAATCAAGCACTAAGGGCGGCTACGGCCGCCCTTAATCAATGCAGCAACAGCAGCCGCTTAATTTGCCTTTAAGGAACCGGAGGAATCATGACCAAGTCAGAATTGATAGaa                                                               <  2:1341506/150‑1 (MQ=255)
                                                                                        tgttgAGTTTACTTGACAGATTGCAGGTTTCGTCCTGTAATCAAGCACTAAGGGCGGCTACGGCCGCCCTTAATCAATGCAGCAACAGCAGCCGCTTAATTTGCCTTTAAGGAACCGGAGGAATCATGACCAAGTCAGAATTGATAGaaa                                                              >  2:420727/1‑150 (MQ=255)
                                                                                          ttgAGTTTACTTGACAGATTGCAGGTTTCGTCCTGTAATCAAGCACTAAGGGCGGCTACGGCCGCCCTTAATCAATGCAGCAACAGCAGCCGCTTAATTTGCCTTTAAGGAACCGGAGGAATCATGACCAAGTCAGAATTGATATAAAGa                                                            >  1:22253/1‑150 (MQ=255)
                                                                                          ttgAGTTTACTTGACAGATTGCAGGTTTCGTCCTGTAATCAAGCACTAAGGGCGGCTACGGCCGCCCTTAATCAATGCAGCAACAGCAGCCGCTTAATTTGCCTTTAAGGAACCGGAGGAATCATGACCAAGTCAGAATTGATAGAAAGa                                                            >  1:1251281/1‑150 (MQ=255)
                                                                                          ttgAGTTTACTTGACAGATTGCAGGTTTCGTCCTGTAATCAAGCACTAAGGGCGGCTACGGCCGCCCTTAATCAATGCAGCAACAGCAGCCGCTTAATTTGCCTTTAAGGAACCGGAGGAATCATGACCAAGTCAGAATTGATAGAAAGa                                                            >  2:594724/1‑150 (MQ=255)
                                                                                             aGTTTACTTGACAGATTGCAGGTTTCGTCCTGTAATCAAGCACTAAGGGCGGCTACGGCCGCCCTTAATCAATGCAGCAACAGCAGCCGCTTAATTTGCCTTTAAGGAACCGGAGGAATCATGACCAAGTCAGAATTGATAGAAAGACtt                                                         <  2:441321/150‑1 (MQ=255)
                                                                                             aGTTTACTTGACAGATTGCAGGTTTCGTCCTGTAATCAAGCACTAAGGGCGGCTACGGCCGCCCTTAATCAAGGCAGCACCAGCAGCCGCTTAATTTGCCTTTAAGGAACCGGAGGAATCATGACCAAGTCAGAATTGATAGAAAGACtt                                                         <  2:171253/150‑1 (MQ=255)
                                                                                                   cTTGACAGATTGCAGGTTTCGTCCTGTAATCAAGCACTAAGGGCGGCTACGGCCGCCCTTAATCAATGCAGCAACAGCAGCCGCTTAATTTGCCTTTAAGGAACCGGAGGAATCATGACCAAGTCAGAATTGATAGAAAGACTTGCCAcc                                                   <  2:7019/150‑1 (MQ=255)
                                                                                                   cTTGACAGATTGCAGGTTTCGTCCTGTAATCAAGCACTAAGGGCGGCTACGGCCGCCCTTAATCAATGCAGCAACAGCAGCCGCTTAATTTGCCTTTAAGGAACCGGAGGAATCATGACCAAGTCAGAATTGATAGAAAGACTTGCCAcc                                                   <  1:1650594/150‑1 (MQ=255)
                                                                                                         aGATTGCAGGTTTCGTCCTGTAATCAAGCACTAAGGGCGGCTACGGCCGCCCTTAATCAATGCGGCAACAGCAGCCGCTTAATTTGCCTTTAAGGAACCGGAGGAATCATGACCAAGTCAGAATTGATAGAAAGACTTGCCACCCAGCaa                                             >  1:47200/1‑150 (MQ=255)
                                                                                                           aTTGCAGGTTTCGTCCTGTAATCAAGCACTAAGGGCGGCTACGGCCGCCCTTAATCAATGCAGCAACAGCAGCCGCTTAATTTGCCTTTAAGGAACCGGAGGAATCATGACCAAGTCAGAATTGATAGAAAGACTTGCCACCCAGCAATc                                           <  2:1416424/150‑1 (MQ=255)
                                                                                                            ttGCAGGTTTCGTCCTGTAATCAAGCACTAAGGGCGGCTACGGCCGCCCTTAATCAATGCAGCAACAGCAGCCGCTTAATTTGCCTTTAAGGAACCGGAGGAATCATGACCAAGTCAGAATTGATAGAAAGACTTGCCACCCAGCAATCg                                          <  1:1466912/150‑1 (MQ=255)
                                                                                                               cAGGTTTCGTCCTGTAATCAAGCACTAAGGGCGGCTACGGCCGCCCTTAATCAATGCAGCAACAGCAGCCGCTTAATTTGCCTTTAAGGAACCGGAGGAATCATGACCAAGTCAGAATTGATAGAAAGACTTGCCACCCAGCAATCGcac                                       >  1:651814/1‑150 (MQ=255)
                                                                                                                  gTTTCGTCCTGGAATCAAGCTCTAAGGGCGGCTACGGCCGCCCTTAATCAATGCAGCACCAGCAGCCGCTTAATTTGCCTTTAAGGAACCGGAGGAATCATGACCAAGTCAGAATTGATAGAAAGACTTGCCACCCAGCAATCGCACAtt                                    <  1:929408/150‑1 (MQ=255)
                                                                                                                       gTCCTGTAATCAAGCACTAAGGGCGGCTACGGCCGCCCTTAATCACTGCAGCAACAGCAGCCGCTTAATTTGCCTTTAAGGAACCGGAGGAATCATGACCAAGTCAGAATTGATAGAAAGACTTGCCACCCAGCAATCGCACAGTCCCGc                               <  1:528481/150‑1 (MQ=255)
                                                                                                                       gTCCTGTAATCAAGCACTAAGGGCGGCTACGGCCGCCCTTAATCAATGCAGCAACAGCAGCCGCTTAATTTGCCTTTAAGGAACCGGAGGAATCATGACCAAGTCAGAATTGATAGAAAGACTTGCCACCCAGCAATCGCACATTCCCGc                               >  1:1609112/1‑150 (MQ=255)
                                                                                                                       gTCCTGTAATCAAGCACTAAGGGCGGCTACGGCCGCCCTTAATCAATGCAGCAACAGCAGCCGCTTAATTTGCCTTTAAGGAACCGGAGGAATCATGACCAAGTCAGAATTGATAGAAAGACTTGCCACCCAGCAATCGCACATTCCCGc                               >  1:54598/1‑150 (MQ=255)
                                                                                                                          cTGTAATCAAGCACTAAGGGCGGCTACGGCCGCCCTTAATCAATGCAGCAACAGCAGCCGCTTAATTTGCCTTTAAGGAACCGGAGGAATCATGACCAAGTCAGAATTGATAGAAAGACTTGCCACCCAGCAATCGCACATTCCCGCCaa                            >  2:648546/1‑150 (MQ=255)
                                                                                                                            gTAATCAAGCACTAAGGGCGGCTACGGCCGCCCTTAATCAATGCAGCAACAGCAGCCGCTTAATTTGCCTTTAAGGAACCgga                                                                                             <  1:1110373/83‑1 (MQ=255)
                                                                                                                               aTCAAGCACTAAGGGCGGCTACGGCCGCCCTTAATCAATGCAGCAACAGCAGCCGCTTAATTTGCCTTTAAGGAACCGGAGGAATCATGACCAAGTCAGAATTGATAGAAAGACTTGCCACCCAGCAATCGCACATTCCCGCCAAGACgg                       >  2:472872/1‑150 (MQ=255)
                                                                                                                               aTCAAGCACTAAGGGCGGCTACGGCCGCCCTTAATCAATGCAGCAACAGCAGCCGCTTAATTTGCCTTTAAGGAACCGGAGGAATCATGACCAAGTCAGAATTGATAGAAAGACTTGCCACCCAGCAATCGCACATTCCCGCCAAGACgg                       >  2:472904/1‑150 (MQ=255)
                                                                                                                               aTCAAGCACTAAGGGCGGCTACGGCCGCCCTTAATCAATGCAGCAACAGCAGCCGCTTAATTTGCCTTTAAGGAACCGGAGGAATCATGACCAAGTCAGAATTGATAGAAAGACTTGCCACCCAGCAATCGCACATTCCCGCCAAGACgg                       >  2:229457/1‑150 (MQ=255)
                                                                                                                               aTCAAGCACTAAGGGCGGCTACGGCCGCCCATAATCAATGCAGCAACAGCAGCCGCTTAATTTGCCTTTAAGGAACCGGAGGAATCATGACCAAGTCAGAATTGATAGAAAGACTTGCCACCCAGCAATCGCACATTCCCGCCAAGACgg                       <  1:1292773/150‑1 (MQ=255)
                                                                                                                                 cAATCACTAAGGGCGGCTACGGCCGCCCTTAATCAATGCAGCAACAGCAGCCGCTTAATTTGCCTTTAAGGAACCGGAGGAATCATGACCAAGTCAGAATTGATAGAAAGACTTGCCACCCAGCAATCGCACATTCCCGCCAAGACGGtt                     <  1:1341045/150‑1 (MQ=255)
                                                                                                                                 cAAGCACTAAGGGCGGCTACGGCCGCCCTTAATCAATGCAGCAACAGCAGCCGCTTAATTTGCCTTTAAGGAACCGGAGGAATCATGACCAAGTCAGAATTGATAGAAAGACTTGCCACCCAGCAATCGCACATTCCCGCCAAGACGGtt                     <  1:1571659/150‑1 (MQ=255)
                                                                                                                                 cAAGCACTAAGGGCGGCTACGGCCGCCCTTAATCAATGCAGCAACAGCAGCCGCTTAATTTGCCTTTAAGGAACCGGAGGAATCATGACCAAGTCAGAATTGATAGAAAGACTTGCCACCCAGCAATCGCACATTCCCGCCAAGACGGtt                     <  2:1687023/150‑1 (MQ=255)
                                                                                                                                   ttgattAAGGGCGGCCGTAGCCGCCCTTAGTCAATGCAGCAACAGCAGCCGCTTAATTTGCCTTTAAGGAACCGGAGGAATCATGACCAAGTCAGAATTGATAGAAAGACTTGCCACCCAGCAATCGCACATTCCCGCCAAGACGGTTGa                   >  1:521229/6‑150 (MQ=255)
                                                                                                                                    gCACTAAGGGCGGCTACGGCCGCCCTTAATCAATGCAGCAACAGCAGCCGCTTAATTTGCCTTTAAGGAACCGGAGGAATCATGACCAAGTCAGAATTGATAGAAAGACTTGCCACCCAGCAATCGCACATTCCCGCCAAGACGGTTGaa                  <  2:1648583/150‑1 (MQ=255)
                                                                                                                                        tAAGGGCGGCTACGGCCGCCCTTAATCAATGCAGCAACAGCAGCCGCTTAATTTGCCTTTAAGGAACCGGAGGAATCATGACCAAGTCAGAATTGATAGAAAGACTTGCCACCCAGCAATCGCACATTCCCGCCAAGACGGTTGAAGATg              >  1:509498/1‑150 (MQ=255)
                                                                                                                                          aGGGCGGCTACGGCCGCCCTTAATCAATGCAGCAACAGCAGCCGCTTAATTTGCCTTTAAGGAACCGGAGGAATCATGACCAAGTCAGAATTGATAGAAAGACTTGCCACCCAGCAATCGCACATTCCCGCCAAGACGGTTGAAGATGCa            <  2:1198269/150‑1 (MQ=255)
                                                                                                                                          aGGGCGGCTACGGCCGCCCTTAATCAATGCAGCAACAGCAGCCGCTTAATTTGCCTTTAAGGAACCGGAGGAATCATGACCAAGTCAGAATTGATAGAAAGACTTGCCACCCAGCAATCGCACATTCCCGCCAAGACGGTTGAAGATGCa            <  2:1198284/150‑1 (MQ=255)
                                                                                                                                          aGGGCGGCTACGGCCGCCCTTAATCAATGCAGCAACAGCAGCCGCTTAATTTGCCTTTAAGGAACCGGAGGAATCATGACCAAGTCAGAATTGATAGAAAGACTTGCCACCCAGCAATCGCACATTCCCGCCAAGACGGTTGAAGATGCa            <  2:156848/150‑1 (MQ=255)
                                                                                                                                              cggcTACGGCCGCCCTTAATCAATGCAGCAACAGCAGCCGCTTAATTTGCCTTTAAGGAACCGGAGGAATCATGACCAAGTCAGAATTGATAGAAAGACTTGCCACCCAGCAATCGCACATTCCCGCCAAGACGGTTGAAGATGCAGTaa        >  1:337345/1‑150 (MQ=255)
                                                                                                                                              cggcTACGGCCGCCCTTAATCAATGCAGCAACAGCAGCCGCTTAATTTGCCTTTAAGGAACCGGAGGAATCATGACCAAGTCAGAATTGATAGAAAGACTTGCCACCCAGCAATCGCACATTCCCGCCAAGACGGTTGAAGATGCAGTaa        >  1:337288/1‑150 (MQ=255)
                                                                                                                                              cggcTACGGCCGCCCTTAATCAATGCAGCAACAGCAGCCGCTTAATTTGCCTTTAAGGAACCGGAGGAATCATGACCAAGTCAGAATTGATAGAAAGACTTGCCACCCAGCAATCGCACATTCCCGCCAAGACGGTTGAAGATGCAGTaa        >  1:772571/1‑150 (MQ=255)
                                                                                                                                                    cGGCCGCCCTTAATCAATGCAGCAACAGCAGCCGCTTAATTTGCCTTTAAGGAACCGGAGGAATCATGACCAAGTCAGAATTGATAGAAAGACTTGCCACCCAGCAATCGCACATTCCCGCCAAGACGGTTGAAGATGCAGTAAAAgaga  <  2:1064984/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
AACTTCTCCAACAACGCAATGGCTGAAGCTTTCAAAGCAGCTAAAGGCGAGTAATTCTCTGACTCTTCGGGATTTTTATTCCGAAGTTTGTTGAGTTTACTTGACAGATTGCAGGTTTCGTCCTGTAATCAAGCACTAAGGGCGGCTACGGCCGCCCTTAATCAATGCAGCAACAGCAGCCGCTTAATTTGCCTTTAAGGAACCGGAGGAATCATGACCAAGTCAGAATTGATAGAAAGACTTGCCACCCAGCAATCGCACATTCCCGCCAAGACGGTTGAAGATGCAGTAAAAGAGA  >  NC_000913/963615‑963912

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