Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,344,935 G→A T208T (ACC→ACT yjdF ← DUF2238 domain‑containing inner membrane protein YjdF

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,344,9350GA98.0% 154.6 / ‑4.4 50T208T (ACC→ACTyjdFDUF2238 domain‑containing inner membrane protein YjdF
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base A (33/15);  total (35/15)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.32e-01

GATCCAGATCCATCTGCTGGATAAATATCGCAAAGTACCGCCCGAGCCTGTTCATGCTGATATTCCGGTCGGTGCAGTGAGTGATGTGAAAGCCTGACGTGGACAATAGCTAACGAAAATGGCCTGATGCGATATGTTTATCAGGCCACGTCTTAATGTAATGCTTTGGATTTGCAAATTTTGTAGACCGGATAAGGAATTCACGCCGCATCCGGCATCAACAAAGCGCAAGTTGTTATCCGGTTATCAAGCCAAAGCGCCGTAGCTGGCGGCAATGAAAGCGAGCGAGGAATATCACCGTCGTTAATGCGCCAAGCAGCGCGCAGAACATATCGGATTGCGTGTCCCACTGGTCGCCCTGGGTGCCGAGAAAATCATCCGCTCCCTGACCCATTGCCAGCGCTGCCCACCACTCGATTAATTCATACATCGCGCTTATCGCCAGCGCCACGCAGCAGACCAGGAACGCCACCATTTTACGTCCGCGC  >  NC_000913/4344694‑4345181
                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                      
gatccagatccaTCTGCTGGATAAATATCGCAAAGTACCGCCCGAGCCTGTTCATGCTGATATTCCGGTCGGTGCAGTGAGTGATGTGAAAGCCTGACGTGGACAATAGCTAACGAAAATGGCCTGATGCGATATGTTTATCAGGCCACGTCTTAATGTAATGCTTTGGATTTGCAAATTTTGTAGACCGGATAAGGAATTCACGCCGCATCCGGCATCAACAAAGCGCAAGTTGTTATCCAGTTATCa                                                                                                                                                                                                                                                 >  1:144399/1‑249 (MQ=255)
            tCTGCTGGATAAATATCGCAAAGTACCGCCCGAGCCTGTTCATGCTGATATTCCGGTCGGTGCAGTGAGTGATGTGAAAGCCTGACGTGGACAATAGCTAACGAAAATGGCCTGATGCGATATGTTTATCAGGCCACGTCTTAATGTAATGCTTTGGATTTGCAAATTTTGTAGACCGGATAAGGAATTCACGCCGCATCCGGCATCAACAAAGCGCAAGTTGTTATCCAGTTATCAAGCCAAAGCGCCGt                                                                                                                                                                                                                                   >  1:60875/1‑251 (MQ=255)
             ctgctgGATAAATATCGCAAAGTACCGCCCGAGCCTGTTCATGCTGATATTCCGGTCGGTGCAGTGAGTGATGTGAAAGCCTGACGTGGACAATAGCTAACGAAAATGGCCTGATGCGATATGTTTATCAGGCCACGTCTTAATGTAATGCTTTGGATTTGCAACTTTTGTAGACCGGATAAGGAATTCACGCCGCATCCGGCATCAACAAAGCGCAAGTTGTTATCCAGTTATCAACCCAAAGCGCCGt                                                                                                                                                                                                                                   >  1:136247/1‑250 (MQ=255)
               gctgGATAAATATCGCACAGTACCGCCCGAGCCTGTTCATGCTGATGTTCCGGTCGGTGCAGTGAGTGATGTGAAAGCCTGACGTGGACAATAGCTAACGACAATGGCCTGATGCGATATGTTTATCAGGCAACGTCTTACTGTAATGCTTTGGATTTGCAAATTTTGTAGACCGGATAAGGAATTCACGCCGCATCAGGCATCAACAAAGCGCAAGTTGTTATCCCGTTATCAAGCCAAAGCGCCGTAGc                                                                                                                                                                                                                                >  2:91424/1‑251 (MQ=255)
               gctgGATAAATATCGCAAAGTACCGCCCGAGCCTGTTCATGCTGATATTCCGGTCGGTGCAGTGAGTGATGTGAAAGCCTGACGTGGACAATAGCTAACGAAAATGGCCTGATGCGATATGTTTATCAGGCCACGTCTTAATGTAATGCTTTGGATTTGCAAATTTTGTAGACCGGATAAGGAATTCACGCCGCATCCGGCATCAACAAAGCGCAAGTTGTTATCCAGTTATCAAGCCAAAGCGCCGTAGc                                                                                                                                                                                                                                >  1:118028/1‑251 (MQ=255)
                      aaaTATCGCAAAGTACCGCCCGAGCCTGTTCATGCTGATATTCCGGTCGGTGCAGTGAGTGATGTGAAAGCCTGACGTGGACAATAGCTAACGAAAATGGCCTGATGCGATATGTTTATCAGGCCACGTCTTAATGTAATGCTTTGGATTTGCAAATTTTGTAGACCGGATAAGGAATTCACGCCGCATCCGGCATCAACAAAGCGCAAGTTGTTATCCAGTTATCAAGCCAAAGCGCCGTAGCTggcggc                                                                                                                                                                                                                         >  1:368336/1‑251 (MQ=255)
                             gCAAAGTACCGCCCGAGCCTGTTCATGCTGATATTCCGGGCGGTGCAGTGAGTGATGTGAAAGCCTGACGTGGACAATAGCTAACGAAAAGGGCCTTATGCGATATGTTTATCAGGCCACGTCTTAATGTAATGCTTTGGATTTGCAAATTTTGTAGACCGGATAAGGAATTCACGCCGCATCCGGCAGCAACAAAGCGCAAGTTGTTATCCAGTTATCAAGCCAAAGCGCCGTAGCTGGCGGCAATGaaa                                                                                                                                                                                                                  <  1:51780/251‑1 (MQ=255)
                                   tACCGCCCGAGCCTGTTCATGCTGATATTCCGGTCGGTGCAGTGAGTGATGTGAAAGCCTGACGTGGACAATAGCTAACGAAAATGGCCTGATGCGATATGTTTATCAGGCCACGTCTTAATGTAATGCTTTGGATTTGCAAATTTTGTAGACCGGATAAGGAATTCACGCCGCATCCGGCATCAACAAAGCGCAAGTTGTTATCCAGTTATCAAGCCAAAGCGCCGTAGCTGGCGGCAATGAAagcgagc                                                                                                                                                                                                            >  1:297666/1‑251 (MQ=255)
                                    aCCGCCCGAGCCTGTTCATGCTGATATTCCGGTCGGTGCAGTGAGTGATGTGAAAGCCTGACGTGGACAATAGCTAACGAAAATGGCCTGATGCGATATGTTTATCAGGCCACGTCTTAATGTAATGCTTTGGATTTGCAAATTTTGTAGACCGGATAAGGAATTCACGCCGCATCCGGCATCAACAAAGCGCAAGTTGTTATCCAGTTATCAAGCCAAAGTGCCGTAGCTGGCGGCAATGAAagcgagcg                                                                                                                                                                                                           >  1:337390/1‑251 (MQ=255)
                                      cGCCCGAGCCTGTTCATGCTGATATTCCGGTCGGTGCAGTGAGTGATGTGAAAGCCTGACGTGGACAATAGCTAACGAAAATGGCCTGATGCGATATGTTTATCAGGCCACGTCTTAATGTAATGCTTTGGATTTGCAAATTTTGTAGACCGGATAAGGAATTCACGCCGCATCCGGCATCAACAAAGCGCAAGTTGTTATCCAGTTATCAAGCCAAAGCGCCGTAGCTGGCGGCAATGAAagcgagcg                                                                                                                                                                                                           >  1:286636/1‑249 (MQ=255)
                                                                   gtcggtGCAGTGAGTGATGTGAAAGCCTGACGTGGACAATAGCTAACGAAAATGGCCTGATGCGATATGTTTATCAGGCCACGTCTTAATGTAATGCTTTGGATTTGCAAATTTTGTAGACCGGATAAGGAATTCACGCCGCATCCGGCATCAACAAAGCGCAAGTTGTTATCCAGTTATCAAGCCAAAGCGCCGTAGCTGGCGGCAATGAAAGCGAGCGAGGAATATCACCGTCGTTAATGCGCCAagc                                                                                                                                                                             >  1:120656/1‑250 (MQ=255)
                                                                          cAGTGAGTGATGTGAAAGCCTGACGTGGACAATAGCTAACGAAAATGGCCTGATGCGATATGTTTATCAGGCCACGTCTTAATGTAATGCTTTGGATTTGCAAATTTTGTAGACCGGATAAGGAATTCACGCCGCATCCGGCATCAACAAAGCGCAAGTTGTTATCCAGTTATCAAGCCAAAGCGCCGTAGCTGGCGGCAATGAAAGCGAGCGAGGAATATCACCGTCGTTAATGCGCCAAGCAgcgcgc                                                                                                                                                                      >  1:105182/1‑250 (MQ=255)
                                                                                            ccTGACGTGGACAATAGCTAACGAAAATGGCCTGATGCGATATGTTTATCAGGCCACGTCTTAATGTAATGCTTTGGATTTGCAAATTTTGTAGACCGGATAAGGAATTCACGCCGCATCCGGCATCAACAAAGCGCAAGTTGTTATCCAGTTATCAAGCCAAAGCGCCGTAGCTGGCGGCAATGAAAGCGAGCGAGGAATATCACCGTCGTTAATGCGCCAAGCAGCGCGCAGAACATATCGGATTGCgt                                                                                                                                                   >  2:379573/1‑251 (MQ=255)
                                                                                                    gggtataaGCTACCGAAAATGGCCGGATTGGATAGTGGTATCAGGCCACTTCGTAATGTAATTCTTTGCATTTGCAAATTTTTTAGACCGGAAAAGGAATTCACTCCGCAGCCGGCATCAACAAAGCGCAAGGTGTGATCCATTTATCAAGCCAAAGCGCCGTAGCTGGCGGCAAGGAAAGCGAGCGAGGAATATCTCCGTCGTTAATGCGCCAAGCAGAGCGCAGAACATATCGGATTGCGTGTCCCAc                                                                                                                                            <  2:477975/243‑1 (MQ=255)
                                                                                                    ggACAATAGCTAACGAAAATGGCCTGATGCGATATGTTTATCAGGCAACGTCTTAATGTAATGCTTTGGATTTGCAAATTTTGTAGACCGGATAAGGAATTCACGCCGCATCCGGCATCAACAAAGCGCAAGTTGTTATCCGGTTATCAAGCCAAAGCGCCGTAGCTGGCGGCAATGAAAGCGAGCGAGGAATATCACCGTCGTTCATGCGCCAAGCAGCGCGCAGAACATATCGGATTGCGTGTCCCACt                                                                                                                                           >  2:144438/1‑251 (MQ=255)
                                                                                                     gACAATAGCTAACGAAAATGGCCTGATGCGATATGTTTATCAGGCCACGTCTTAATGTAATGCTTTGGATTTGCAAATTTTGTAGACCGGATAAGGAATTCACGCCGCATCCGGCATCAACAAAGCGCAAGTTGTTATCCAGTTATCAAGCCAAAGCGCCGTAGCTGGCGGCAATGAAAGCGAGCGAGGAATATCACCGTCGTTAATGCGCCAAGCAGCGCGCAGAACATATCGGATTGCGTGTCCCACTg                                                                                                                                          <  2:120656/251‑1 (MQ=255)
                                                                                                      aCAATAGCTAACGAAAATGGCCTGATGCGATATGTTTATCAGGCCACGTCTTAATGTAATGCTTTGGATTTGCAAATTTTGTAGACCGGATAAGGAATTCACGCCGCATCCGGCATCAACAAAGCGCAAGTTGTTATCCAGTTATCAAGCCAAAGCGCCGTAGCTGGCGGCAATGAAAGCGAGCGAGGAATATCACCGTCGTTAATGCGCCAAGCAgcgc                                                                                                                                                                        <  1:32733/220‑1 (MQ=255)
                                                                                                      aCAATAGCTAACGAAAATGGCCTGATGCGATATGTTTATCAGGCCACGTCTTAATGTAATGCTTTGGATTTGCAAATTTTGTAGACCGGATAAGGAATTCACGCCGCATCCGGCATCAACAAAGCGCAAGTTGTTATCCAGTTATCAAGCCAAAGCGCCGTAGCTGGCGGCAATGAAAGCGAGCGAGGAATATCACCGTCGTTAATGCGCCAAGCAgcgc                                                                                                                                                                        >  2:32733/1‑220 (MQ=255)
                                                                                                       cactAGCTAATGAAAATGGCCTGAGGCGATATGTTTCTCAGGCAAAGTCTTAATGGAGTGCTTTGGCTTTGCAAATTTTGTAGACCGGATAAGGAATTCACGCCGCATCAGGCATCAACAAAGCGCAAGATGTAATCCAGTTACCAAGCCACAGCGCCGAAGCTGGCGGCAATGAAGGCGAGCGAGGAATATCACCGGCGTTACTGCGCCCAGCACAGCGCAGAAAATATCGGATTGAGGGTccccggggg                                                                                                                                        >  2:27631/4‑245 (MQ=255)
                                                                                                        aaTAGCTAACGAAAATGGCCTGATGCGATATGTTTATCAGGCCACGTCTTAATGTAATGCTTTGGATTTGCAAATTTTGTAGACCGGATAAGGAATTCACGCCGCATCCGGCATCAACAAAGCGCAAGTTGTTATCCAGTTATCAAGCCAAAGCGCCGTAGCTGGCGGCAATGCAAGCGAGCGAGGAATATCACCGTCGTTAATGCGCCAAGCAGCGCGCAGAAAATATCGGATTGCGTGTCCCACTGGTc                                                                                                                                       >  2:296215/1‑251 (MQ=255)
                                                                                                                   aaaaTGGCCTGATGCGATATGTTTATCAGGCCCCGTCTTAATGTAATGCTTTGGATTTGCACATTTTGTAGACCGGCTAAGGACTTCACGCCGCATCCGGCATCAACAAAGCGCAAGTTGTTATCCAGTTATCAAGCCAAAGCGCCGTAGCTGGCGGCAATGAAAGCGAGCGAGGAATATCACCGTCGTTAATGCGCCAAGCAGCGCGCAGAACATATCGGCTTGCGTGTCCCACTCGTCGCCCTGGGTGc                                                                                                                            >  2:254551/1‑251 (MQ=255)
                                                                                                                                     atGTTTATCAGGCCACGTCTTAATGTAATGCTTTGGATTTGCAAATTTTGTAGCCCGGATAAGGACTTCACGCCGCATCCGGCATCAACAAAGCGCAAGTTGTTATCCAGTTATCAAGCCAACGCGCCGTAGCTGGCGGCAATGCCAGCGAGAGAGGAATATCACCGTCGTTAATGCGCCCAGCACCGCGCCGAACATATCGGATTGCGTGTCCCACTTGACGCCATGGGTGCCGAGAAAATCATCAGCTc                                                                                                          >  2:478197/1‑251 (MQ=255)
                                                                                                                                        gTTTATCAGGCCACGTCTTAATGTAATGCTTTGGATTTGCAAATTTTGTAGACCGGATAAGGAATTCACGCCGCATCCGGCATCAACAAAGCGCAAGTTGTTATCCAGTTATCAAGCCAAAGCGCCGTAGCTGGCGGCAATGAAAGCGAGCGAGGAATATCACCGTCGTTAATGCGCCAAGCAGCGCGCAGAACATATCGGATTGCGTGTCCCACTGGTCGCCCTGGGTGCCGAGAAAATCATCCGCTccc                                                                                                        <  1:356953/251‑1 (MQ=255)
                                                                                                                                          ttATCAGGCCACGTCTTAATGTAATGCTTTGGATTTGCAAATTTTGTAGACCGGATAAGGAATTCACGCCGCATCCGGCATCAACAAAGCGCAAGTTGTTATCCAGTTATCAAGCCAAAGCGCCGTAGCTGGCGGCAATGAAAGCGAGCGAGGAATATCACCGTCGTTAATGCGCCAAGCAGCGCGCAGAACATATCGGATTGCGTGTCCCACTGGTCGCCCTGGGTGCCGAGAAAATCATCCGCTCCCt                                                                                                       <  1:296215/250‑1 (MQ=255)
                                                                                                                                             tCAGGCCACGTCTTAATGTAATGCTTTGGATTTGCAAATTTTGTAGACCGGATAAGGAATTCACGCCGCATCCGGCATCAACAAAGCGCAAGTTGTTATCCAGTTATCAAGCCAAAGCGCCGTAGCTGGCGGCAATGAAAGCGAGCGAGGAATATCACCGTCGTTAATGCGCCAAGCAGCGCGCAGAACATATCGGATTGCGTGTCCCACTGGTCGCCCTGGGTGCCGAGAAAATCATCCGCTCCCTGAcc                                                                                                   >  1:311096/1‑251 (MQ=255)
                                                                                                                                                 gCCACGTCTTAATGTAATGCTTTAGATTTGCAAATTTTGTAGACCGGCTAAGGAATTCACGCCGCATCCGGCATCAACAAAGCGCAAGTTGTTATCCAGTTATCAAGCCAAAGCGCCGTAGCTGGCGGCAATGAAAGCGAGCGAGGAATATCACCGTCGTTAATGCGCCAAGCAGCGCGCGGAACATATCGGATTGCGTGTCCCACTGGTCGCCCCGGGTACCAGGAAAATCATCCGCTCCCTGAcccttt                                                                                               >  2:193641/1‑248 (MQ=255)
                                                                                                                                                 cccACGTCTTAATGTAATGCTTTGGATTTGCAAATTTTGTAGACCGTATAAGGAATTCACGCCGCATCCGGCATCAACAAAGCGCAAGTTGTTATCCAGTTATCAAGCCAAAGCGCCGTAGCTGGCGGCAATGAAAGCGAGCGAGGAATATCACCGTCGTTAATGCGCCAAGCAGCGCGCAGAACATATCGGATTGCGTGTCCCACTGGTCGCCCTGGGTGCCGAGAAAATCATCCGCTCCCTGACCCAtt                                                                                               <  2:311096/250‑1 (MQ=255)
                                                                                                                                                      gTCTGAATGTAATGCTTTGGATTTGCAAATTTTGTAGACCGGATAAGGAATTCACGCCGCATCCGGCATCAACAAAGCGCAAGTTGTTAGCCAGTTATCAAGCCAAAGCGCCGTAGCTGGCGGCAATGAAAGCGAGCGAGGAATATCACCGTCGTTAATGCGCCAAGCAGCGCGCAGAACATATCGGATTGCGTGTCCCACTGGTCGCCCTGGGGGCCGAGAAAATCATCCGCTCCCTGACCCATTGCCAg                                                                                          <  2:136247/251‑1 (MQ=255)
                                                                                                                                                           aatgtaatgCTTTGGATTTGCAAATTTTGTAGACCGGATAAGGAATTCACGCCGCATCCGGCCTCAACAAAGCGCAAGTTGTTATCCAGTTATCAAGCCAAAGCGCCGTAGCTGGCGGCAATGAAAGCGAGCGAGGAATATCACCGTCGTTAATGCGCCAAGCAGCGCGCAGAACATATCGGATTGCGTGTCCCACTGGTCGCCCTGGGTGCCGAGAAAATCATCCGCTCCCTGACCCATTGCCAGCGCTg                                                                                     <  2:489641/251‑1 (MQ=255)
                                                                                                                                                           aatgtaatgCTTTGGATTTGCAAATTTTGTAGACCGGATAAGGAATTCACGCCGCATCCGGCATCAACAAAGCGCAAGTTGTTATCCAGTTATCAAGCCAAAGCGCCGTAGCTGGCGGCAATGAAAGCGAGCGAGGAATATCACCGTCGTTAATGCGCCAAGCAGCGCGCAGAACATATCGGATTGCGTGTCCCACTGGTCGCCCTGGGTGCCGAGAAAATCATCCGCTCCCTGACCCATTGCCAGCGCTg                                                                                     >  1:263974/1‑251 (MQ=255)
                                                                                                                                                              gtaatgATTTGGATTTGCAAATTTTGTAGACCGGATAAGGAATACACGCCGCATCCGGCATCAACAAAGCGCAAGTTGTTATCCAGTTATCAAGCCAAAGCGCCGTAGCTGGCGGCAATGAAAGCGAGCGAGGAATATCACCGTCGTTAATGCGCCAAGCAGCGCGCAGAACATATCGGATTGCGTGTCCCACTGGTCGCCCTGGGTGCCGAGAAAATAATCCGCTCCCTGACCCATTGCCAGAGCTGccc                                                                                  >  1:190693/1‑251 (MQ=255)
                                                                                                                                                                          aTTTGCAAATTTTGTAGACCGGATAAGGAATTCACGCCGCATCCGGCATCAACAAAGCGCAAGTTGTTATACAGTTATCAAGCCAAAGCGCCGTAGCTGGCGGCAATGAAAGCGAGCGAGGAATATCACCGTAGTTAATGCGCCTAGCAGCACGCAGAACATATCGGATTGCGTGTCCCACTGGTCGCCCTGGGTGCAGAGAAAATCATCCGCTCCCTGACCCTTTGCACGCGCTGccccccccgcatta                                                                       >  2:54978/1‑239 (MQ=255)
                                                                                                                                                                            ttGCAAATTTTGTAGACCGGATAAGGTATTCACGCCGCATCCGGCATCAACAAAGCGCAAGTTGTTATCCAGTTATCAAGCCAAAGCGCCGTAGCTGGCGGCAATGAAAGCGAGCGAGGAATATCACCGTCGTTAATGCGCCAAGCAGCGCGCAGAACATATCGGATTGCGTGTCCCACTGGTCGCCCTGGGTGCCGAGAAAATCATCCGCTCCCTGACCCATTGCCAGCGCTGCCCACCACTCGATTAAt                                                                    >  1:194470/1‑251 (MQ=255)
                                                                                                                                                                             tGCAAATTTTGTAGACCGGATAAGGAATCCACGCCGCATCCGGCATCAACAAAGCGCAAGTTGTTATCCAGTTATCATGCCAAAGCGCCGTAGCTGGCGGCAATGAAAGCGAGCGATGAATATCACCGTCGTTAATGCGCCAAGCAGCGCGCAGAACATATCGGATTGCGTGTCCCACTGGGCGCCCTGGGTGCCGAGAAAATCATCCGCTCCCTGACCCATTGCCAGCGCTGCCCACCACTCGATTAAt                                                                    <  1:314576/250‑1 (MQ=255)
                                                                                                                                                                                     ttGTAGACCGGATAAGGAATTCACGCCGCATCCGGCATCAACAAAGCGCAAGTTGTTATCCAGTTATCAAGCCAAAGCGCCGTAGCTGGCGGCAATGAAAGCGAGCGAGGAATATCACCGTCGTTAATGCGCCAAGCAGCGCGCAGAACATATCGGATTGCGTGTCCCACTGGTCGCCCTGGGTGCCGAGAAAATCATCCGCTCCCTGACCCATTGCCAGCGCTGCCCACCACTCGATTAATTCATAAATc                                                           >  2:381547/1‑251 (MQ=255)
                                                                                                                                                                                     ttGTAGACCGGATAAGGAATTCACGCCGCATCCGGCATCAACAAAACGCAAGTTGTTATCCAGTTATCAAGCCAAAGCGCCGTAGCTGGCGGCAATGAAAGCGCGCGAGGAATATCACCGTCGTTAATACGCCAAGCAGCGCGCAGAACATATCGGATTGCGTGTCCCACTGGTCGCCCTGGGGGCCGAGAAAATCATCCGCTCCCTGACCCATTGCCAGCGCTGCCCACCACTCGCTTAATTCATAcctc                                                           >  2:309910/1‑248 (MQ=255)
                                                                                                                                                                                                  aagtaaTTCACGCCGCAGCCGGCATCAACAAAGCGCAAGTTGTTATCCAGTTATCAAGCCAAAGCGCCGTAGCTGGCGGCAATGAAAGCGAGCGAGGAATATCACCGTCGTTAATGCGCCAAGCAGCGCGCAGAACATATCGGATTGCGTGTCCCACTGGTCGCTCTGGGTGCCGAGAAAATCATCCGCTCCCTGACCCATTGCCAGCGCTGCCCACCACTCGATTAATTCATACATCGCGCTTATCGcc                                               <  2:268381/246‑1 (MQ=255)
                                                                                                                                                                                                  aaGGAATTCACGCCGCATCCGGCATCAACAAAGCGCAAGTTGTTATCCAGTTATCAAGCCAAAGCGCCGTAGCTGGCGGCAATGAAAGCGAGCGAGGAATATCACCGTCGTTAATGCGCCAAGCAGCGCGCAGAACATATCGGATTGCGTGTCCCACTGGTCGCCCTGGGTGCCGAGAAAATCATCCGCTCCCTGACCCATTGCCAGCGCTGCCCACCACTCGATTAATTCATACATCGCGCTTATCGcc                                               <  1:175174/250‑1 (MQ=255)
                                                                                                                                                                                                   aGGAATTCACGCCGCATCCGGCATCAACAAAGCGCAAGTTGTTATCCAGTTATCAAGCCAAAGCGCCGTAGCTGGCGGCAATGAAAGCGAGCGAGGAATATCACCGTCGTTAATGCGCCAAGCAGCGCGCAGAACATATCGGATTGCGTGTCCCACTGGTCGCCCTGGGTGCCGAGAAAATCATCCGCTCCCTGACCCATTGCCAGCGCTGCCCACCACTCGCTTAATTAATACATCGCGCTTATCGca                                               >  1:345098/1‑248 (MQ=255)
                                                                                                                                                                                                       aTTCACGCCGCATCCGGCATCAACAAAGCGCAAGTTGTTATCCAGTTATCAAGCCAAAGCGCCGTAGCTGGCGGCAATGAAAGCGAGCGAGGAATATCACCGTCGTTAATGCGCCAAGCAGCGCGCAGAACATATCGGATTGCGTGTCCCACTGGTCGCCCTGGGTGCCGAGAAAATCATCCGCTCCCTGACCCATTGCCAGCGCTGCCCACCACTCGATTAATTCATACATCGCGCTTATCGCCAGCGcc                                         >  1:170937/1‑251 (MQ=255)
                                                                                                                                                                                                                 cATCCGGCATCAACAAAGCGCAAGTTGTTATCCAGTTATCAAGCCAAAGCGCCGTAGCTGGCGGCAATGAAAGCGAGCGAGGAATATCACCGTCGTTAATGCGCCAAGCAGCGCGCAGAACATATCGGATTGCGTGTCCCACTGGTCGCCCTGGGTGCCGAGAAAATCATCCGCTCCCTGACCCATTGCCAGCGCTGCCCACCACTCGATTAATTCATACATCGCGCTTATCg                                                 >  1:412149/1‑233 (MQ=255)
                                                                                                                                                                                                                 cATCCGGCATCAACAAAGCGCAAGTTGTTATCCAGTTATCAAGCCAAAGCGCCGTAGCTGGCGGCAATGAAAGCGAGCGAGGAATATCACCGTCGTTAATGCGCCAAGCAGCGCGCAGAACATATCGGATTGCGTGTCCCACTGGTCGCCCTGGGTGCCGAGAAAATCATCCGCTCCCTGACCCATTGCCAGCGCTGCCCACCACTCGATTAATTCATACATCGCGCTTATCg                                                 <  2:412149/233‑1 (MQ=255)
                                                                                                                                                                                                                     cGGCATCAACAAAGCGCAAGTTGTTATCCAGTTATCAAGCCAAAGCGCCGTAGCTGGCGGCAATGAAAGCGAGCGAGGAATATCACCGTCGTTAATGCGCCAAGCAGCGCGCAGAACATATCGGATTGCGTGTCCCACTGGTCGCCCTGGGTGCCGAGAAAATCATCCGCTCCCTGACCCATTGCCAGCGCTGCCCACCACTCGATTAATTCATACATCGCGCTTATCGCCAGCGCCCCGCAGCAGAcc                             >  2:501190/1‑249 (MQ=255)
                                                                                                                                                                                                                              caaAGCGCAAGTTGTTATCCAGTTATCAAGCCAAAGCGCCGTAGCTGGCGGCAATGAAAGCGAGCGAGGAATATCACCGTCGTTAATGCGCCAAGCAGCGCGCAGAACATATCGGATTGCGTGTCCCACTGGTCGCCCTGGGTGCCGAGAAAATCATCCGCTCCCTGACCCATTGCCAGCGCTGCCCACCACTCGATTAATTCATACATCGCGCTTATCGCCAGCGCCACGCAGCAGACCAGGAACGccac                  >  2:122753/1‑251 (MQ=255)
                                                                                                                                                                                                                                           gttATCCAGTTATCAAGCCAAAGCGCCGTAGCTGGCGGCAATGAAAGCGAGCGAGGAATATCACCGTCGTTAATGCGCCAAGCAGCGCGCAGAACATATCGGATTGCGTGTCCCACTGGTCGCCCTGGGTGCCGAGAAAATCATCCGCTCCCTGACCCATTGCCAGCGCTGCCCACCACTCGCTTAATTCATACATCGCGCTTATCGCCAGCGCCACGCAGCAGACCAGGAACGccacagatccgcacacc     >  2:123566/1‑238 (MQ=255)
                                                                                                                                                                                                                                           gttATCCAGTTATCAAGCCAAAGCGCCGTAGCTGGCGGCAATGAAAGCGAGCGAGGAATATCACCGTCGTTAATGCGCCAAGCAGCGCGCAGAACATATCGGATTGCGTGTCCCACTGGTCGCCCTGGGTGCCGAGAAAATCATCCGCTCCCTGACCCATTGCCAGCGCTGCCCACCACTCGATTAATTCATACATCGCGCTTATCGCCAGCGCCACGCAGCAGACCAGGAACGccac                  <  1:123566/238‑1 (MQ=255)
                                                                                                                                                                                                                                             tATCCAGTTATCAAGCCAAAGCGCCGTAGCTGGCGGCAATGAAAGCGAGCGAGGAATATCACCGTCGTTAATGCGCCAAGCAGCGCGCAGAACATATCTGATTGCGTGTCCCACTGGTCGCCCTGGGTGCCGAGAAAATCATCCGCTCCCTGACCCATTGCCAGCGCTGCCCACCACTCGATTAATTCATACATCGCGCTTATCGCCAGCGCCACGCAGCAGACCAGGAACGccacca                <  1:32097/238‑1 (MQ=255)
                                                                                                                                                                                                                                             tATCCAGTTATAAAGCCAAAGCGCCGTAGCTGGCGGCCATGAAAGCGAGCGAGGAATATCACCGTCGTTAATGCGCCAAGCAGAGCGCAGAACATATCGGATTGCGTGTCCCACTGGTCGCCCTGGGTGCCGAGAAACTCATCCGCTCCCTGACCCATTGCCAGCGCTGCCCACCACTCGATTAATTCATACATCGCGCTTATCGCCAGCGCCACGCAGCAGACCAGGAACGccacca                >  2:32097/1‑238 (MQ=255)
                                                                                                                                                                                                                                              aTCCAGTTATCAAGCCAAAGCGCCGTAGCTGGCGGCAATGAAAGCGAGCGAGGAATATCACCGTCGTTAATGCGCCAAGCAGCGCGCAGAACATATCGGATTGCGTGTCCCACTGGTCGCCCTGGGTGCCGAGAAAATCATCCGCTCCCTGACCCATTGCCAGCGCTGCCCACCACTCGATTAATTCATACATCGCGCTTATCGCCATCGCCACGCAGCAGACCAGGAACGCCACCATTTTACGTCCgcgc  >  2:479665/1‑251 (MQ=255)
                                                                                                                                                                                                                                              aTCCAGTTATCAAGCCAAAGCGCCGTAGCTGGCGGCAATGAAAGCGAGCGAGGAATATCACCGTCGTTAATGCGCCAAGCAGCGCGCAGAACATATCGGATTGCGTGTCCCACTGGTCGCCCTGGGTGCCGAGAAAATCATCCGCTCCCTGACCCATTGCCAGCGCTGCCCACCACTCGATTAATTCATACATCGCGCTTATCGCCAGCGCCACGCAGCAGACCAGGAACGCCACCATTTTACGTCCgcgc  >  1:160628/1‑251 (MQ=255)
                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                      
GATCCAGATCCATCTGCTGGATAAATATCGCAAAGTACCGCCCGAGCCTGTTCATGCTGATATTCCGGTCGGTGCAGTGAGTGATGTGAAAGCCTGACGTGGACAATAGCTAACGAAAATGGCCTGATGCGATATGTTTATCAGGCCACGTCTTAATGTAATGCTTTGGATTTGCAAATTTTGTAGACCGGATAAGGAATTCACGCCGCATCCGGCATCAACAAAGCGCAAGTTGTTATCCGGTTATCAAGCCAAAGCGCCGTAGCTGGCGGCAATGAAAGCGAGCGAGGAATATCACCGTCGTTAATGCGCCAAGCAGCGCGCAGAACATATCGGATTGCGTGTCCCACTGGTCGCCCTGGGTGCCGAGAAAATCATCCGCTCCCTGACCCATTGCCAGCGCTGCCCACCACTCGATTAATTCATACATCGCGCTTATCGCCAGCGCCACGCAGCAGACCAGGAACGCCACCATTTTACGTCCGCGC  >  NC_000913/4344694‑4345181

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