Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,146,384 A→G V71A (GTA→GCA)  yceF ← m(7)GTP pyrophosphatase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,146,3840AG93.1% 163.9 / 2.1 58V71A (GTA→GCA) yceFm(7)GTP pyrophosphatase
Reads supporting (aligned to +/- strand):  ref base A (3/1);  new base G (26/28);  total (29/29)
Fisher's exact test for biased strand distribution p-value = 6.11e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.19e-02

CTCTTAAAGCTACCCGCGCAGTGCAGGGGATGCTCTTTACGCACGTAATTATCAATCTCCGCCTCGCTCAGGTGACGGAAATGGACGTCAAAAGGCTCCACTTCTGTTTGCAGATGCCCATTCGCCGAATTAAACAACGCCAGTCCAGTATAGAAGGTGACGATATTGCCGCTGGCTTTGCGTAATTGCAGACGGGCATTTTCTTCCGTTAACGGTTTACCGGTGATTTCGCCATCAAGAACACATACCTGGTCTGAACCAATAATTAAATGATCCGGATAACGTGACGCCAGAGATTGCGCTTTTTCTTGTGCCAGTCGAAGCACCAACTGTCGCGGTGATTCGTCGCTGCGTGGGGTTTCGTCGACCTCTGGTGCTGCACATTCAAAAGAGATTTGCAGTTTTTCCAGAAGAGCGCGACGCCATGGCGAGGTGGAGGCTAAAATAAGTTTAGGCATATTTTTTTCCATCAGATATAGCGTAT  >  NC_000913/1146138‑1146621
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                             
ctctTAAAGCTACCCGCGCAGTGCAGGGGATGCTCTTTACGCACGTAATTATCAATCTCCGCCTCGCTCAGGTGACGGAAATGGACGTCAAAAGGCTCCACTTCTGTTTGCAGATGCCCATTCGCCGAATTAAACAACGCCAGTCCAGTATAGAAGGTGACGATATTGCCGCTGGCTTTGCGTAATTGCAGACGGGCATTTTCTTCCGTTAACGGTTTACCGGTGATTTCGCCATCAAGAACACATGCCTg                                                                                                                                                                                                                                           >  1:261988/1‑251 (MQ=255)
             ccGCGCAGTGCAGGGGATGCTCTTTACGCACGTAATTATCAATCTCCGCCTCGCTCAGGTGACGGAAATGGACGTCAAAAGGCTCCACTTCTGTTTGCAGATGCCCATTCGCCGAATTAAACAACGCCAGTCCAGTATAGAAGGTGACGATATTGCCGCTGGCTTTGCGTAATTGCAGACGGGCATTTTCTTCCGTTAACGGTTTACCGGTGATTTCGCCATCAAGAACACATGCCTGGACTGAACCaat                                                                                                                                                                                                                               >  1:197366/1‑250 (MQ=255)
                     ggCAGGGGATGCTCTTTACGCACGTAATTTTCAATCTCCGCCTCGCTCAGGTGACGGAAATGGACGTCAAAAGGCTCCACTTCTGTTTGCAGATGCCCATTCGCCGAATTAAACAACGCCAGTCCAGTATAGAAGGTGACGATATTGCCGCTGGCTTTGCGTAATTGCAGACGGGCATTTTCTTCCGTTAACGGTTTACCGGTGATTTCGCCATCAAGAACACATGCCTGGTCTGAACCAATAATTAAATg                                                                                                                                                                                                                      <  2:14240/250‑1 (MQ=255)
                        aGGGGATGCTCTTTACGCACGTAATTATCAATCTCCGCCTCGCTCAGGTGACGGAAATGGACGTCAAAAGGCTCCACTTCTGTTTGCAGATGCCCATTCGCCGAATTAAACAACGCCAGTCCAGTATAGAAGGTGACGATATTGCCGCTGGCTTTGCGTAATTGCAGACGGGCATTTTCTTCCGTTAACGGTTTACCGGTGATTTCGCCATCAAGAACACATGACTGGTCTGCACCACTAATTAAATGATc                                                                                                                                                                                                                   >  2:336446/1‑251 (MQ=255)
                           ggATGCTCTTTACGCACGTAATTATCAATCTCCGCCTCGCTCAGGTGACGGAAATGGACGTCAAAAGGCTCCACTTCTGTTTGCAGATGCCCATTCGCCGAATTAAACAACGCCAGTCCAGTATAGAAGGTGACGATATTGCCGCTGGCTTTGCGTAATTGCAGACGGGCATTTTCTTCCGTTAACGGTTTACCGGTGATTTCGCCATCAAGAACACATGCCTGGTCTGAACCAATAATTAAATGATCCgg                                                                                                                                                                                                                >  1:377792/1‑251 (MQ=255)
                                    ttACGCACGTAATTATCAATCTCCGCCTCGCTCAGGTGACGGAAATGGACGTCAAAAGGCTCCACTTCTGTTTGCAGATGCCCATGCGCCGAATTAAACAACGCCAGTCCAGTAGAGAAGGTGACGATATTGCCGCTGGCTTTGCGTAATTGCAGACGGGCATTTTCTTCCGTTAACGGTTTACCGGTGATTTCGCCATCAAGAACACATGCCTGTTCTGAACCAATAATTAAATGATCCGGATAACg                                                                                                                                                                                                          <  1:423825/248‑1 (MQ=255)
                                        gCACGTAATTATCAATCTCCGCCTCGCTCAGGTGACGGAAATGGACGTCAAAAGGCTCCACTTCTGTTTGCAGATGCCCATTCGCCGAATTAAACAACGCCAGTCCAGTATAGAAGGTGACGATATTGCCACTGGCTTTGCGTAATTGCAGACGGGCATTTTCTTCCGTTAACGGTTTACCGGTGATTTCGCCATCAAGAACACATACCTGGTCTGAACCAATAATTAAATGATCCGGATAACGTGACGcc                                                                                                                                                                                                   >  1:162342/1‑251 (MQ=255)
                                              aaTTATCAATCTCCGCCTCGCTCAGGTGACGGAAATGGACGTCAAAAGGCTCCACTTCTGTTTGCAGATGCCCATTCGCCGAATTAAACAACGCCAGTCCAGTATAGAAGGTGACGATATTGCCGCTGGCTTTGCGTAATTGCAGACGGGCATTTTCTTCCGTTAACGGTTTACCGGTGATTTCGCCATCAAGAACACATGCCTGGTCTGAACCAATAATTAAATGATCCGGATAACGTGACGCCag                                                                                                                                                                                                 >  1:218687/1‑247 (MQ=255)
                                                ttATCAATCTCCGCCTCGCTCAGGTGACGGAAATGGACGTCAAAAGGCTCCACTTCTGTTTGCAGATGCCCATTCGCCGAATTAAACAACGCCAGTCCAGTATAGAAGGTGACGATATTGCCGCTGGCTTTGCGTAATTGCAGACGGGCATTTTCTTCCGTTAACGGTTTACCGGTGATTTCGCCATCAAGAACACATGCCTGGTCTGAACCAATAATTAAATGATCCGGATAACGTGACGCCAGAGATTg                                                                                                                                                                                           >  1:117772/1‑251 (MQ=255)
                                                 tATCAATCTCCGCCGCGCTCAGGTGACGGAAATGGACGTCAAAAGGCTCCACTTCTGTTTGCAGATGCCCATTCGCCGAATTAAACAACGCCAGTCCAGTATAGAAGGTGACGATATTGCCGCTGGCTTTGCGTAATTGCAGACGGGCATTTTCTTCCGTTAACGGTTTACCGGTGATTTCGCCATCAAGAACACATGCCTGGTCTGAACCAATAATTAAATGATCCGGATAACGTGACGCCAGAGATTg                                                                                                                                                                                           <  2:355716/250‑1 (MQ=255)
                                                  aTCAATCTCCTCCTCTCTCGGGTGACGGAAATGTACGTCAAAAGGCCCCACTTCTGTTTGCAGATGCCCAATCGCCGAATTAAACAACGACAGTCCAGTATAGAAGGTGTCGATATTGCCGCTGGCTTTGCGTAATTGCAGACGGGCATTTTCTTCCGTTAACGGTTTACCGGTGATTTCGCCATCAATAACACATGCCTGGTCTGAACCAATAAGTAAATGATCCGGATATCGTGACGCCagcgactgc                                                                                                                                                                                          <  2:204535/250‑8 (MQ=255)
                                                          cccccccgttCAGGTGCCGGAAAAGGCCGTCAAAAGTCTCCACTTCTGTTTGCAGATGCCCCTTCGCCGAATTAACCAACGCCAGTCCAGTATAGAAGGTGACGATATTGCCGCTGGCTTTGCGTAATTGCCGTCGGGCATTTTCTTCCGTTAACGGTTTACCGGTGATTTCGCCATCAAGAACACATGCCTTGTCTGAACCAATAATTAAATGATCCTGATAACGTGACGCCAGAGATTGCGCttttt                                                                                                                                                                                   <  2:271341/240‑1 (MQ=255)
                                                            gCCTCGCTCAGGTGACGGAAATGGACGTCAAAAGGCTCCACTTCTGTTTGCAGATGCCCATTCGCCGAATTAAACAACGCCAGTCCAGTATAGAAGGTGACGATATTGCCGCTGGCTTTGCGTAATTGCAGACGGGCATTTTCTTCCGTTAACGGTTTACCGGTGATTTCGCCATCAAGAACACATGCCTGGTCTGAACCAATAATTAAATGATCCGGATAACGTGACGCCAGAGATTGCGCTTTTTCTtg                                                                                                                                                                               <  1:171040/251‑1 (MQ=255)
                                                              cTCGCTCAGGTGACGGAAATGGACGTCAAAAGGCTCCACTTCTGTTTGCAGATGCCCATTCGCCGAATTAAACAACGCCAGTCCAGTATAGAAGGTGACGATATTGCCGATGGCTTTGCGTAATTGCAGACGGGCATTTTCTTCCGTTAACGGTTTACCGGTGATTTCGCCATCAAGAACACATGCCTGGTCTGAACCAATAATTAAATGATCCGGATAACGTGACGCCAGAGATTGCGCTTTTTCTtgtg                                                                                                                                                                             >  2:438209/1‑251 (MQ=255)
                                                                 gCTCAGGTGACGGAAAGGGACGTAAAAAAGCTCCACTTCTGTTTGCAGATGCCCATTCGCCGAATTAAACAACGCCAGTCTAGTATAGACGGTGACGATATTGCCGCTGGCTTTGCGTAATTGCAGACGGGCATTTTCTTCCTTTAACGGTTTACCGGTGATTTCGCCATCAAGAACAAATGCCTGGGCTGAACCAATAATTAAATGATCCTGTTACCGTGACTCCAGAGCTTGCGCTTTTTCTTGTGcc                                                                                                                                                                           >  1:260080/1‑250 (MQ=255)
                                                                    cAGGTGACGGAAATGGACGTCAAAAGGCTCCACTTCTGTTTGCAGATGCCCATTCGCCGAATTAAACAACGCCAGTCCAGTATAGAAGGTGACGATATTGCCGCTGGCTTTGCGTAATTGCAGACGGGCATTTTCTTCCGTTAACGGTTTACCGGTGATTTCGCCATCAAGAACACATGCCTGGTCTGAACCAATAATTAAATGATCCGGATAACGTGACGCCAGAGATTGCGCTTTTTCTTGTGCCAg                                                                                                                                                                         <  1:260548/249‑1 (MQ=255)
                                                                         gACGGAAATGGACGTCAAAAGGCTCCACTTCTGTTTGCAGATGCCCATTCGCCGAATTAAACAACGCCAGTCCAGTATAGAAGGTGACGATATTGCCGCTGGCTTTGCGTAATTGCAGACGGGCATTTTCTTCCGTTAACGGTTTACCGGTGATTTCGCCATCAAGAACACATGCCTGGTCTGAACCAATAATTAAATGATCCGGATAACGTGACGCCAGAGATTGCGCTTTTTCTTGTGCCAGTCGAAg                                                                                                                                                                   <  1:336446/250‑1 (MQ=255)
                                                                         gACGGAAATGGACGTCAAAAGGCTCCACTTCTGTTTGCAGATGCCCATTCGCCGAATTAAACAACGCCAGTCCAGTATAGAAGGTGACGATATTGCCGCTGGCTTTGCGTAATTGCAGACGGGCATTTTCTTCCGTTAACGGTTTACCGGTGATTTCGCCATCAAGAACACATGCCTGGTCTGAACCAATAATTAAATGATCCGGATAACGTGACGCCAGAGATTGCGCTTTTTCTTGTGCCAGTCGAAg                                                                                                                                                                   <  1:422445/250‑1 (MQ=255)
                                                                            ggAAATGGACGTCAAAAGGCTCCACTTCTGTTTGCAGATGCCCATTCGCCGAATTAAACAACGCCAGTCCAGTATAGAAGGTGACGATATTGCCGCTGGCTTTGCGTAATTGCAGACGGGCATTTTCTTCCGTTAACGGTTTACCGGTGATTTCGCCATCAAGAACACATGCCTGGTCTGAACCAATAATTAAATGATCCGGATAACGTGACGCCaga                                                                                                                                                                                                <  2:314977/218‑1 (MQ=255)
                                                                                aTGGACGTCAAAAGGCTCCACTTCTGTTTGCAGATGCCCATTCGCCGAATTAAACAACGCCAGTCCAGTATAGAAGGTGACGATATTGCCGCTGGCTTTGCGTAATTGCAGACGGGCATTTTCTTCCGTTAACGGTTTACCGGTGATTTCGCCATCAAGAACACATGCCTGGTCTGAACCAATAATTAAATGATCCGGATAACGTGACGCCGGAGATTGCGCTTTTTCTTGTGCCAGTCGAAGCACCAACt                                                                                                                                                           >  2:157013/1‑251 (MQ=255)
                                                                                        cAAAAGGCTCCACTTCTGTTTGCAGATGCCCATTCGCCGAATTAAACAACGCCAGTCCAGTATAGAAGGTGACGATATTGCCGCTGGCTTTGCGTAATTGCAGACGGGCATTTTCTTCCGTTAACGGTTTACCGGTGATTTCGCCATCAAGAACACATGCCTGGTCTGAACCAATAATTAAATGATCCGGATAACGTGACGCCAGAGATTGCGCTTTTTCTTGTGCCAGTCGAAGCACCAACTGTCGCGGt                                                                                                                                                   <  1:419365/251‑1 (MQ=255)
                                                                                           aaGGCTCCCATTCTGTTTGCAGATACCCATTCCCCGAATTCAACACCGCCAGTCCCGAATAGAAAGTGACGATCTTGCCGCTGGCTTTGCGTAATTGCAGACGAGCATATTCTTCCGTTAACGGGTTACCGGAGATTTCGCCATCCAGAAAACATGCCTGGTCTGAACCAATAATTAAATGATCCGGAAAACGTGACACCAGAGCTTGCGCTTTTTCTTGTGCCAGTCGAACCACAAAcgatcgcagggat                                                                                                                                                >  2:259273/1‑239 (MQ=255)
                                                                                           aaGGCTCCACTTCTGTTTGCAGATGCCCATTCGCCGAATTAAACAACGCCAGTCCAGTATAGAAGGTGACGATCTTGCCGCTGGCTTTGCGTAATTGCAGACGGGCATTTTCTTCCGTTAACGGTTTACCGGTGATTTCGCCATCAAGAACACATGCCTGGTCTGAACCAATAATTAAATGATCCGGCTAACGTGACGCCAGAGATTGCGCTTTTTCTTGTGCCAGTCGAAGCACCAACTGTCGCGTTGAt                                                                                                                                                >  2:382971/1‑251 (MQ=255)
                                                                                             ggCTCCACTTCTGTTTGCAGATGCCCATTCGCCGAATTAAACAACGCCAGTCCAGTATAGAAGGTGACGATATTGCCGCTGGCTTTGCGTAATTGCAGACGGGCATTTTCTTCCGTTAACGGTTTACCGGTGATTTCGCCATCAAGAACACATGCCTGGTCTGAACCAATAATTAAATGATCCGGATAACGTGACGCCAGAGATTGCGCTTTTTCTTGTGCCAGTCGAAGCACCAACTGTCGCGGTGATtc                                                                                                                                              <  1:157013/251‑1 (MQ=255)
                                                                                                tCCACTTCTGTTTGCAGATGCCCATTCGCCGAATTAAACAACGCCAGTCCAGTATAGAAGGTGACGATATTGCCGCTGGCTTTGCGTAATTGCAGACGGGCATTTTCTTCCGTTAACGGTTTACCGGTGATTTCGCCATCAAGAACACATACCTGGTCTGAACCAATAATTAAATGATCCGGATAACGTGACGCCAGAGATTGCGCTTTTTCTTGTGCCAGTCGAAGCACCAACTGTCGCGGTGCTtcgtc                                                                                                                                           >  2:128635/1‑251 (MQ=255)
                                                                                                    cTTCTGTTTGCAGATGCCCATTCGCCGAATTAAACAACGCCAGTCCAGTATAGAAGGTGACGATATTGCCGCTGGCTTTGCGTAATTGCAGACGGGCATTTTCTTCCGTTAACGGTTTACCGGTGATTTCGCCATCAAGAACACATGCCTGGTCTGAACCAATAATTAAATGATCCGGATAACGTGACGCCAGAGATTGCGCTTTTTCTTGTGCCAGTCGAAGCACCAACTGTCGCGGTGATTCGTCGCTg                                                                                                                                       >  1:70440/1‑251 (MQ=255)
                                                                                                      tCTGTTTGCAGATGCCCATTCGCCGAATTAAACAACGCCAGTCCAGTATAGAAGGTGACGATATTGCCGCTGGCTTTGCGTAATTGCAGACGGGCATTTTCTTCCGTTAACGGTTTACCGGTGATTTCGCCATCAAGAACACATGCCTGGTCTGAACCAATAATTAAATGATCCGGATAACGTGACGCCAGAGATTGCGCTTTTTCTTGTGCCAGTCGAAGCACCAACTGTCGCGGTGATTCGTCGCTGCg                                                                                                                                     >  1:477341/1‑251 (MQ=255)
                                                                                                        tGTTTGCAGATGCCCATTCGCCGAATTAAACAACGCCAGTCCAGTATAGAAGGTGACGATATTGCCGCTGGCTTTGCGTAATTGCAGACGGGCATTTTCTTCCGTTAACGGTTTACCGGTGATTTCGCCATCAAGAACACATACCTGGTCTGAACCAATAATTAAATGATCCGGATAACGTGACGCCAGAGATTGCGCTTTTTCTTGTGCCAGTCGAAGCACCAACTGTCGCGGTGATTCGTCGCTGCGTg                                                                                                                                   >  1:118761/1‑251 (MQ=255)
                                                                                                         gTTTGCAGATGCCCATTCGCCGAATTAAACAACGCCAGTCCAGTATAGAAGGTGACGATATTGCCGCTGGCTTTGCGTAATTGCAGACGGGCATTTTCTTCCGTTAACGGTTTACCGGTGATTTCGCCATCAAGAACACATGCCTGGTCTGAACCAATAATTAAATGATCCGGATAACGTGACGCCAGAGATTGCGCTTTTTCTTGTGCCAGTCGAAGCACCAACTGTCGCGGTGATTCGTCGCTGCGTgg                                                                                                                                  >  2:436507/1‑251 (MQ=255)
                                                                                                            tGCAGATGCCCGTTCGCCGAATTAAACAACGCCAGTCCAGAATCGAAGGTGACGATCTTGACGCTGGCTTTGAGTAATTGCAGACAGGCATTTTCTTCCTTTAACGGTTTACCGATGATTTCGCCATCCAGAACACATGCCAGGTCTGAACCAATAATTAAATGATCCGGATAACGCGACGCCCGAGATTGCGCTTTTTCTTGTGCCAGTCGAAGCACACACTGTAGCGGTGATTTGTATCTGCGTggggg                                                                                                                               >  1:439320/1‑250 (MQ=255)
                                                                                                            tGCAGATGCCCATTCGCCGAATTAAACAACGCCAGTCCAGTATAGAAGGTGACGATATTGCCGCTGGCTTTGCGTAATTGCAGACGGGCATTTTCTTCCGTTAACGGTTTACCGGTGATTTCGCCATCAAGAACACATGCCTGGTCTGAACCAATAATTAAATGATCCGGATAACGTGACGCCAGAGATTGCGCTTTTTCTTGTGCCAGTCGAAGCACCAACTGTCGCGGTGAGTCGTCGCTGCGTgggg                                                                                                                                <  1:413360/250‑1 (MQ=255)
                                                                                                                 aTGCCCATTCGCCGAATTAAACAACGCCAGTCCAGTATAGAAGGTGACGATATTGCCGCTGGCTTTGCGTAATTGCAGACGGGCATTTTCTTCCGTTAACGGTTTACCGGTGATTTCGCCATCAAGAACACATGCCTGGTCTGAACCAATAATTAAATGATCCGGATAACGTGACGCCAGAGATTGCGCTTTTTCTTGTGCCAGTCGAAGCACCAACTGTCGCGGTGATTCGTCGCTGCGTGGGGTTtcgt                                                                                                                          <  2:1155/251‑1 (MQ=255)
                                                                                                                     ccATTCGCCGAATTAAACAACGCCAGTCCAGTATAGAAGGTGACGATATTGCCGCTGGCTTTGCGTAATTGCAGACGGGCATTTTCTTCCGTTAACGGTTTACCGGTGATTTCGCCATCAAGAACACATGCCTGGTCTGAACCAATAATTAAATGATCCGGATAACGTGACGCCAGAGATTGCGCTTTTTCTTGTGCCAGTCGAAGCACCAACTGTCGCGGTGATTCGTCGCTGCGTGGGTTTTCGTCGAc                                                                                                                      <  1:382971/251‑1 (MQ=255)
                                                                                                                            ccGAATTAAACAACGCCAGTCCAGTATAGAAGGTGACGATATTGCCGCTGGCTTTGCGTAATTGCAGACGGGCATTTTCTTCCGTTAACGGTTTACCGGTGATTTCGCCATCAAGAACACATACCTGGTCTGAACCAATAATTAAATGATCCGGATAACGTGACGCCAGAGATTGCGCTTTTTCTTGTGCCAGTCGAAGCACCAACTGTCGCGGTGATTCGTCGCTGCGTGGGGTTTCGTCGACCTCTGGt                                                                                                               <  1:128635/251‑1 (MQ=255)
                                                                                                                             cGAATTACACAACGCCAGTCCAGTATAGAAGGTGACGATATTGCCGCTGGCTTTGCGTAATTGCAGACGGGCATTTTCTTCCGTTAACGGTTTACCGGTGATTTCGCCATCAAGAACACATGCCTGGTCTGAACCAATAATTAAATGATCCGGATAACGTGACGCCAGAGATTGCGCTTTTTCTTGTGCCAGTCGAAGCACCAACTGTCGCGGTGATTCGTCGCTGCGGGGGGTTTCGTCGACCTTTGGtg                                                                                                              >  2:503594/1‑251 (MQ=255)
                                                                                                                               aaTTAAACAACGCCAGTACAGTATAGAAGGTGACGATATTGCCGCTGGCTTTGCGTAATTGCAGACGGGCATTTTCTTCCGTTAACGGTTTACCGGTGATTTCGCCATCAAGAACACATGCCTGGTCTGAACCAATAATTAAATGATCCGGATAACGTGACGCCAGAGATTGCGCTTTTTCTTGTGCCAGTCGAAGCACCAACTGTCGCGGTGATTCGTCGCTGCGTGGGGTTTCGTCGACCTCTGGtgcg                                                                                                            >  2:297358/1‑250 (MQ=255)
                                                                                                                                  tAAACAACGCCAGTCCAGTATAGAAGGTGACGATATTGCCGCTGGCTTTGCGTAATTGCAGACGGGCATTTTCTTCCGTTAACGGTTTACCGGTGATTTCGCCATCAAGAACACATGCCTGGTCTGAACCAATAATTAAATGATCCGGATAACGTGACGCCAGAGATTGCGCTTTTTCTTGTGCCAGTCGAAGCACCAACTGTCGCGGTGATTCGTCGCTGCGTGGGGTTTCGTCGACCTCTGGtgctgc                                                                                                          >  1:40731/1‑250 (MQ=255)
                                                                                                                                    aacaacGCCAGTCCAGTATAGAAGGGGACGATATTGCCGCTGGCTTTGCGTAATTGCAGACGGGCATTTTCTTCCGTTAACGGTTTACCGGTGATTTCGCCATCAAGAACACATGCCTGGTCTGAACCAATAATTAAATGATCCGGATAAAGTGACGCCAGAGATTGCGCTTTTTCTTGTGACAGTCGAAGCACCAACTGTCGCGGTGTTTAGTCGCTGCGTGGGGTTTCGTCGACCTCTGGGGCTGcac                                                                                                        >  2:434834/1‑250 (MQ=255)
                                                                                                                                     acaacGCCAGTCCAGTATAGAAGGTGACTATATTGCCGCTGGCTTTGCGTAATTGCAGACGGGCATTTTCTTCCGTTAACGGTTTACCGGTGATTTCGCCATCAAGAACACATGCCTTGTCTGAACCAATAATTAAATGATCCGGATAACGTGACGCCAGAGATTTCGCTTTTTCTTGTGCCATTCGAAGCACCAACTGTCGCGGTGATTCGTCGCTGCGTGGGGTTTCGTCGACCTCTGGTGCTGCACAt                                                                                                      <  1:276819/251‑1 (MQ=255)
                                                                                                                                            cagtccagtATAGAAGGTGACGATATTGCCGCTGGCTTTGCGTAATTGCAGACGGGCATTTTCTTCCGTTAACGGTTTACCGGTGATTTCGCCATCAAGAACACATGCCTGGTCTGAACCAATAATTAAATGATCCGGATAACGTGACGCCAGAGATTGCGCTTTTTCTTGTGCCAGTCGAAGCACCAACTGTCGCGGTGATTCGTCGCTGCGTGGGGTTTCGTCGACCTCTGGTGCTGCACATTCAAAag                                                                                               <  2:58569/251‑1 (MQ=255)
                                                                                                                                             agtccagtATAGAAGGTGACGATATTGCCGCTGGCTTTGCGTAATTGCAGACGGGCATTTTCTTCCGTTAACGGTTTACCGGTGATTTCGCCATCAAGAACACATGCCTGGTCTGAACCAATAATTAAATGATCCGGATAACGTGACGCCAGAGATTGCGCTTTTTCTTGTGCCAGTCGAAGCACCAACTGTCGCGGTGATTCGTCGCTGCGTGGGGTTTCGTCGACCTCTGGTGCTGCACATTCAAAaga                                                                                              <  2:371993/251‑1 (MQ=255)
                                                                                                                                                         aaGGTGACGATATTGCCGCTGGCTTTGCGTAATTGCAGACGGGCATTTTCTTCCGTTAACGGTTTACCGGTGATTTCGCCATCAAGAACACATGCCTGGTCTGAACCAATAATTAAATGATCCGGATAACGTGACGCCAGAGATTGCGCTTTTTCTTGTGCCAGTCGAAGCACCAACTGTCGCGGTGATTCGTCGCTGCGTGGGGTTTCGTCGACCTCTGGTGCTGCACATTCAAAAGAGATTTGCAGttt                                                                                  >  2:358405/1‑251 (MQ=255)
                                                                                                                                                            gTGACGATATTGCCGCTGGCTTTGCGTAATTGCAGACGGGCATTTTCTTCCGTTAACGGTTTACCGGTGATTTAGCCATCAAGAACACATGCCAGGTCTGAACCAATAATTAAATGATCCGGATAACGTGACGCCAGAGATTGCGCTTTTTCTTGTGCCAGTCGAAGCACCAACTGTCGCGGTGATTCGTCGCTGCGTGGGGTTTCGTCGACCTCTGGTGCTGCACATTCAAACGAGATTTGCAGTTTTTc                                                                               >  1:472408/1‑251 (MQ=255)
                                                                                                                                                              gACGATAGTGCCGCTGGCTTTGCGTAATTGCAGACGGGCATTTTCTTCCGTTAACGGTTTACCGGTGATTTCGCCATCAAGAACACATGCCTGGTCTGAACCAATAATTAAATGATCCGGATAACGTGACGCCAGAGATTGCGCTTTTTCTTGTGCCAGTCGAAGCACCAACTGTCGCGGTGATTCGTCGCTGCGTGGGGTTTCGTCGACCTCTGGTGCTGCACATTCAAAAGAGATTTGCAGTTTTTcc                                                                              <  2:83147/250‑1 (MQ=255)
                                                                                                                                                                    ccgacgcttcttccgatctAATTGCAGACGGGCATTTTCTTCCGTTAACGGTTTACCGGTGATTTCGCCATCAAGAACACATGCCTGGTCTGAACCAATAATTAAATGATCCGGATAACGTGACGCCAGAGATTGCGCTTTTTCTTGTGCCAGTCGAAGCACCAACTGTCGCGGTGATTCGTCGCTGCGTGGGGTTTCGTCGACCTCTGGTGCTGCACATTCAAAAGAGATTTGCAGTTTTTCCAGAagag                                                                       <  2:443856/233‑1 (MQ=255)
                                                                                                                                                                           atGGCTTTGCGTAATTGCAGACGGGCATTTTCTTCCGTTAACGGTTTACCGGTGATTTCGCCATCAAGAACACATGCCTGGTCTGAACCAATAATTAAATGATCCGGATAACGTGACGCCAGAGATTGCGCTTTTTCTTGTGCCAGTCGAAGCACCAACTGTCGCGGTGATTCGTCGCTGCGTGGGGTTTCGTCGACCTCTGGTGCTGCACATTCAAAAGAGATTTGCAGTTTTTCCAGAAGAGCGCGACg                                                                <  1:438209/250‑1 (MQ=255)
                                                                                                                                                                                     gTAATTGCAGACGGGCATTTTCTTCCGTTAACGGTTTACCGGTGATTTCGCCATCAAGAACACATGCCTGGTCTGAACCAATAATTAAATGATCCGGATAACGTGACGCCAGAGATTGCGCTTTTTCTTGTGCCAGTCGAAGCACCAACTGTCGCGGTGATTCGTCGCTGCGTGGGGTTTCGTCGACCTCTGGTGCTGCACATTCAAAAGAGATTTGCAGTTTTTCCAGAAGAGCGCGACGCCATGGCg                                                        >  2:370218/1‑249 (MQ=255)
                                                                                                                                                                                       aaTTGCAGACGGGCATTTTCTTCCGTTAACGGTTTACCGGTGATTTCGCCATCAAGAACACATGCCTGGTCTGAAACAATAATTAAATGATCCGGATAACGTGACGCCAGAGATTGCGCTTTTTCTTGTGCCAGTCGAAGCACCAACTGTCGCGGTGATTCGTCGCTGCGTGGGGTTTCGTCGACCTCTGGTGCTGCACATTCAAAAGAGATTTGCAGTTTTTCCAGAagag                                                                       >  1:443856/1‑232 (MQ=255)
                                                                                                                                                                                          tGCAGACGGGCATTTTCTTCCGTTAACGGTTTACCGGTGATTTCGCCATCAAGAACACATGCCTGGTCTGAACCAATAATTAAATGATCCGGATAACGTGACGCCAGAGATTGCGCTTTTTCTTGTGCCAGTCGAAGCACCAACTGTCGCGGTGATTCGTCGCTGCGTGGTGTTTCGTCGACCTCTGGTGCTGCACATTCAAAAGAGATTTGCAGTTTTTCCAGAAGAGCGCGACGCCATGGCGAGGTgg                                                  <  1:22338/250‑1 (MQ=255)
                                                                                                                                                                                                   gCATTTTCTTCCGTTAACGGTTTACCGGTGATTTCGCCATCAAGAACACATGCCTGGTCTGAACCAATAATTAAATGATCCGGATAACGTGACGCCAGAGATTGCGCTTTTTCTTGTGCCAGTCGAAGCACCAACTGTCGCGGTGATTCGTCGCTGCGTGGGGTTTCGTCGACCTCTGGTGCTGCACATTCAAAAGAGATTTGCAGTTTTTCCAGAAGAGCGCGACGCCATGGCGAGGTGGAGGCTAAAAt                                        <  1:358405/251‑1 (MQ=255)
                                                                                                                                                                                                    cATTTTCTTCCGTTGACGATTTACCGTTGATTTCGCCATCAAGAACATATGCCTGGTCTGAACCAATAATTAGATGTTCCGGATACATTGACGCCAGAGCTTGCGCTTTTTCTTGTGCCCGTCGAAGCACCAACTGTCGCGGTGATTCGTCGCTGCGTGGTGTTTCGTCGACCTCTGGTGCTGCACATTCAAAAGAGATTTGGAGTTTTTCCCGAAGAGCGGGACGCCAAGGCGAGGTGAAGGCtcaatt                                        >  2:318513/1‑245 (MQ=255)
                                                                                                                                                                                                       tttCTTCCGTTAACGGTTTACCGGTGATTTCGCCATCAAGAACACATGCCTGGTCTGAACCAATAATTAAATGATCCGGATAACGTGACGCCAGAGATTGCGCTTTTTCTTGTGCCAGTCGAAGCACCAACTGTCGCGGTGATTCGTCGCTGCGTGGGGTTTCGTCGACCTCTGGTGCTGCACATTCAAAAGAGATTTGCAGTTTTTCCAGAAGAGCGCGACGCCATGGCGAGGTGGAGGCTAAAATAAg                                     <  1:503594/250‑1 (MQ=255)
                                                                                                                                                                                                            tcCGGTAACGGTTCACCGGTGAATTCTCCATCAAGAACTCAGGCCTGGGCTGAACCAATAATTAAATGATCCGGATAACGTGACGCCAGAGATTGCGCTTTTTCTTGTGCCAGTCGAAGCACCAACTGTCGCGGTGATTCGTCGCTGCGTGGGGTTTCGTCGACCTCTGGTGCTGCTCATTCAAAAGAGATTTGCAGTTTTTCCAGAAGAGCGCGACGCCATGGCGAGGTGGAGGCTAAAATAAGTTTAg                                <  2:378885/250‑1 (MQ=255)
                                                                                                                                                                                                                          accgtggATTGCGCCATCAAGAACACAGGCCTGGTCTGAACCAATAATTAAATGATCCGGATAACGTGACGCCAGAGATTGCGCTTTTTCTTGTGCCAGTCGAAGCACCAACTGTCGCGGTGATTCGTCGCTGCGTGGGGTTTCGTCGGCCTCTGGTGCTGCACATTCAAAAGAGATTTGCAGTTTTTCCAGAAGAGCGCTACGCCATGGCGAGGTGGAGGCTAAAATAAGTTTAGGCATATTTTTTTCCa                 <  2:472408/245‑1 (MQ=255)
                                                                                                                                                                                                                                   ttCGCCATCAAGAACACATGCCTGGTCTGAACCAATAATTAAATGATCCGGATAACGTGACGCCAGAGATTGCGCTTTTTCTTGTGCCAGTCGAAGCACCAACTGTCGCGGTGATTCGTCGCTGCGTGGGGTTTCGTCGACCTCTGGTGCTGCACATTCAAAAGAGATTTGCAGTTTTTCCAGAAGAGCGCGACGCCATGGCGAGGTGGAGGCTAAAATAAGTTTAGGCATATTTTTTTCCATCAGatatt        >  2:157433/1‑250 (MQ=255)
                                                                                                                                                                                                                                    tCGCCATCAAGAACACATGCCTGGTCTGAACCAATAATTAAATGATCCGGATAACGTGACGCCAGAGATTGCGCTTTTTCTTGTGCCAGTCGAAGCACCAACTGTCGCGGTGATTCGTCGCTGCGTGGGGTTTCGTCGACCTCTGGTGCTGCACATTCAAAAGAGATTTGCAGTTTTTCCAGAAGAGCGCGACGCCATGGCGAGGTGGAGGCTAAAATAAGTTTAGGCATATTTTTTTCCATCAGAtata        <  1:436507/250‑1 (MQ=255)
                                                                                                                                                                                                                                       ccATCCAGAACACATGCCTGGTCTGAACCAATAATTAAATGATCCGGCTAACGTGACGCCAGAGATTGCGCTTTTTCTTGTGCCAGTCGAAGCACCAACTGTCGCGGTGATTCGTCGCTGCGTGGGGTTTCGTCGACCTCTGGTGCTGCACATTCAAAAGAGATTTGCAGTTTTTCCAGAAGAGCGCGACGCCATGGCGAGGTGGAGGATCATAGAAGTTTAGGCATCTTTTTTTCCATAAGTTATAGCGt    >  2:11657/1‑251 (MQ=255)
                                                                                                                                                                                                                                         aTCAAGAACACATGCCTGGTCTGAACCAATAATTAAATGATCCGGATAACGTGACGCCAGAGATTGCGCTTTTTCTTGTGCCAGTCGAAGCACCAACTGTCGCGGTGATTCGTCGCTGCGTGGGGTTTCGTCGACCTCTGGTGCTGCACATTCAAAAGAGATTTGCAGTTTTTCCAGAAGAGCGCGACGCCATGGCGAGGTGGAGGCTAAAATAAGTTTAGGCATATTTTTTTCCATCAGATATAGCGTAt  <  2:477341/251‑1 (MQ=255)
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                             
CTCTTAAAGCTACCCGCGCAGTGCAGGGGATGCTCTTTACGCACGTAATTATCAATCTCCGCCTCGCTCAGGTGACGGAAATGGACGTCAAAAGGCTCCACTTCTGTTTGCAGATGCCCATTCGCCGAATTAAACAACGCCAGTCCAGTATAGAAGGTGACGATATTGCCGCTGGCTTTGCGTAATTGCAGACGGGCATTTTCTTCCGTTAACGGTTTACCGGTGATTTCGCCATCAAGAACACATACCTGGTCTGAACCAATAATTAAATGATCCGGATAACGTGACGCCAGAGATTGCGCTTTTTCTTGTGCCAGTCGAAGCACCAACTGTCGCGGTGATTCGTCGCTGCGTGGGGTTTCGTCGACCTCTGGTGCTGCACATTCAAAAGAGATTTGCAGTTTTTCCAGAAGAGCGCGACGCCATGGCGAGGTGGAGGCTAAAATAAGTTTAGGCATATTTTTTTCCATCAGATATAGCGTAT  >  NC_000913/1146138‑1146621

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑