Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,127,334 G→A G177S (GGC→AGC)  yceM → putative oxidoreductase YceM

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,127,3340GA91.1% 118.5 / ‑1.1 48G177S (GGC→AGC) yceMputative oxidoreductase YceM
Reads supporting (aligned to +/- strand):  ref base G (2/2);  new base A (19/22);  total (23/25)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.84e-01

GCGCGATGCTGAACGGCTGGTGGAACTGGCGGCGCGGAAAAAACTGACGTTGATGGTCGGTTTTAACCGTCGTTTCGCACCACTCTACGGTGAGTTAAAAACGCAACTCGCCACCGCAGCCTCGCTAAGAATGGATAAACATCGTAGCAATAGCGTCGGGCCACACGATCTTTATTTCACGTTGCTGGATGATTATCTGCATGTGGTGGATACCGCGCTGTGGTTGTCGGGCGGCAAAGCCTCTCTGGATGGCGGTACGCTACTGACTAACGACGCTGGCGAAATGCTGTTTGCCGAGCACCATTTTTCGGCTGGTCCTTTGCAGATCACCACCTGTATGCATCGCCGTGCCGGAAGTCAGCGTGAAACCGTGCAGGCCGTGACTGACGGTGCGCTCATCGACATTACGGATATGCGCGAATGGCGTGAGGAGCGCGGGCAGGGCGTAGTGCATAAACCGATTCCTGGTTG  >  NC_000913/1127105‑1127575
                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                 
gcgcGATGCTGAACGGCTGGTGGAACCGGCGGCGCGGAAAAAACTGACGGTGAGGGGCGGTTTTAACCGTCGTTTCGCACCACTCTACGGTGAGTTAAAAACGCAACTCGCCACCGCAGCCTCGCTAAGAATGGATAAACATCGTAGCAATAGCGTCGGGCCACACGATCTTTATTTCACGTTGCTGGATGATTATCTGCATGTGGTGGAGACCGCGCTGTGGTTGTCGAGCGACAAAGCCTCTCTGGATg                                                                                                                                                                                                                              <  2:434813/251‑1 (MQ=255)
      tGCTGAACGGCTGGTGGAACTGGCGGCGCGGAAAAAACTGACGTTGATGGTCGGTTTTAACCGTCGTTTCGCACCACTCTACGGTGAGTTAAAAACGCAACTCGCCACCGCAGCCTCGCTAAGAATGGATAAACATCGTAGCAATAGCGTCGGGCCACACGATCTTTATTTCACGTTGCTGGATGATTATCTGCATGTGGTGGATACCGCGCTGTGGTTGTCGAGCGACAAAGCCTCTCTGGCTGGCGGt                                                                                                                                                                                                                         >  1:296467/1‑250 (MQ=255)
                    ggtaaCTGGCGGCGCGGAAAAAACTGACGTTGATGGTCGGTTTTAACCGTCGTTTCGCACCACTCTACGGTGAGTTAAAAACGCAACTCGCCACCGCAGCCTCGCTAAGAATGGATAAACAGCGTAGCAATAGCGTCGGGCCACACGATCTTTATTTCACGTTGCTGGATGATTATCTGCATGTGGTGGATACCGCGCTGTGGTTGTCGAGCGACAAAGCCTCTCTGGATGGCGGTACGCTACTGACTAac                                                                                                                                                                                                          <  2:121886/248‑1 (MQ=255)
                       aaCTGGCGGCGCGGAAAAAACTGACGTTGATGGTCGGTTTTAACCGTCGTTTCGCACCACTCTACGGTGAGTTAAAAACGCAACTCGCCACCGCAGCCTCGCTAAGAATGGATAAACATCGTAGCAATAGCGTCGGGCCACACGATCTTTATTTCACGTTGCTGGATGATTATCTGCATGTGGTGGATACCGCGCTGTGGTTGTCGAGCGACAAAGCCTCTCTGGATGGCGGTACGCTACTGACTAacgac                                                                                                                                                                                                       <  1:44435/251‑1 (MQ=255)
                                  cgGAAAAAACTGACGTTGATGGTCGGTTTTAACCGTCGTTTCGCACCACTCTACGGTGAGTTAAAAACGCAACTCGCCACCGCAGCCTCGCTAAGAATGGATAAACATCGTAGCAATAGCGTCGGGCCACACGATCTTTATTTCACGTTGCTGGATGATTATCTGCATGTGGTGGATACCGCGCTGTGGTTGTCGAGCGACAAAGCCTCTCTGGATGGCGGTACGCTACTGACTAACGACGCTGGCGAAAt                                                                                                                                                                                            <  1:6489/251‑1 (MQ=255)
                                     aaaaaaCTGACGTTGATGGTCGTTTTTAACCGTCGTTTCGCACCACTCTACGGTGAGTTAAAAACTCAACTCGCCACCGCAGCCTCGCTAAGAATGGATAAACATCGTAGCAATAGCGTCGGGCCACACGATCTTTATTGCACGTTGCTGGATGATTATCTGCATGTGGTGGATACCGCGCTGTGGTTGTCGAGCGACAAATCCTCTCTGGATGGCGGTACGCTACTGACTATCGACGCTGGCGAAATGCt                                                                                                                                                                                         <  1:294675/251‑1 (MQ=255)
                                                         cGGTTTTAACCGTCGTTTCGCACCTCTCTACGGTGAGTTAAAAACGCAACTCGCCACCGCAGCCTCGCTAAGAATGGATAAACATCGTAGCAATAGCGTCGTGCCACACGATCTTTATTTCACGTTGCTGGATGATTATCTGCATGTGGTGGATACCGCGCTGTGGTTGTCGAGCGACAAAGCCTCTCTGGATGGCGGTACGCTACTGACTAACGACGCTGGCGAAATGCTGTTTGCCGAGCACCAtttt                                                                                                                                                                      <  1:355041/250‑1 (MQ=255)
                                                                              accacTCTACGGTGAGTTAAAAACGCAACTCGCCACCGCAGCCTCGCTAAGAATGGATAAACATCGTAGCAATAGCGTCGGGCCACACGATCTTTATTTCACGTTGCTGGATGATTATCTGCATGTGGTGGATACCGCGCTGTGGTTGTCGAGCGACAAAGCCTCTCTGGATGGCGGTACGCTACTGACTAACGACGCTGGCGAAATGCTGTTTGCCGAGCACCATTTTTCGGCTGGTCCTTTGCAGATcc                                                                                                                                                >  2:8776/1‑250 (MQ=255)
                                                                                      aCGGTGAGTTAACAAAGCAACTCACAACCGCAGCCTCGCTAAGAATGGATACAAATCGTCGCAATAGCGTCGGGCCACACGATCTTTATTTAACGTTGCTGGATGATTCTCTGCATGTGGTGGATACCGCGCTGCGGTTGTCGGGCGACAAAGCCTATCTGGATCGCGGTACGCTACTGCCTCACGACGCGGGCGAAATGCTGTTTGCCGAGCCCCATTTTTCGGCTGGCCCTTTGCATCTCACCACCTGt                                                                                                                                        >  2:24803/1‑251 (MQ=255)
                                                                                      aCGGTGAGTTAAAAACGCAACTCTCCACCTCAGCCGCGCTAAGAATGGATAAACATCGTAGCAATAGCGTCGGGCCACACGATCTTTATTTCACGTTGCTGGATGATTATCTGCATGTGGTGGATACCGCGCTGTGGTTGTCGAGCGACAAAGCCTCTCTGGATGGCGGTACGCTACTGACTAACGACGCTGGCGAAATGCTGTTTGCCGAGCACCATTTTTCGGCTGGTCCTTTGCAGATCACCGCCTGt                                                                                                                                        <  2:153944/251‑1 (MQ=255)
                                                                                      aCGGTGAGTTAAAAACGCAACTCGCCACCGCAGCCTCGCTAAGAATGGATAAACATCGTAGCAATAGCGTCGGGCCACACGATCTTTATTTCACGTTGCTGGATGATTATCTGCATGTGGTGGATACCGCGCTGTGGTTGTCGAGCGACAAAGCCTCTCTGGATGGCGGTACGCTACTGACTAACGACGCTGGCGAAATGCTGTTTGCCGAGCACCATTTTTCGGCTGGTCCTTTGCAGATCaccacc                                                                                                                                           <  1:305845/248‑1 (MQ=255)
                                                                                         gTGATTTAAAAACGCCACCCGCCTCCGCAGCCTCGCCAAGAATGGATAAACATCGTAGCAATAGCGTCGGGCCACACGATCTTTATGTCACGTTGCTGGATGATTATCTGCATGTGGTGGATACCGCGCTGTGGTTGTCGAGCGACAAAGCCTCTCTGGATGGCGGTACGCTACTGACTAACGACGCTGGCGAAATGCTGTTTGCCGAGCACCATTTTTCGGCTGGTCCTTTGCAGATCACCACCTGTATg                                                                                                                                     <  1:267442/251‑1 (MQ=255)
                                                                                                            cGCCACCGCAGCCTCTCTAAGAATGGATAAACATCGTAGCAATAGCGTCGGGCCACACGATCTTTATTTCACGTTGCTGGATGATTATCTGCATGTGGTGGATACCGCGCTGTGGTTGTCGAGCGACAAAGCCTCTCTGGATGGCGGTACGCTACTGACTAACGACGCTGGCGAAATGCTGTTTGCCGAGCACCATTTTTCGGCTGGTCCTGTGCAGATCACCACCTGTATGCATCGCCGTGCCGGAAg                                                                                                                    <  2:182686/249‑1 (MQ=255)
                                                                                                             gCCACCGCAGCCTCGCTAAGAATGGATAAACATCGTAGCAATAGCGTCGGGCCCCACGCTCTTTATTTCACGTTGCTGGATGATTATCTGCATGTGGTGGATACCGCGCTGTGGTTGTCGCGCGACAAAGCCTCTCTGGATGGCGGTACGCTACTGACTAACGACGCTGGCGAAATGCTGTTTGCCGAGCACCATTTTTCGGCTGGTCCTTTTCAGATCACCACCTGTAGGCCTCGCCGTGCCGGGAGTc                                                                                                                  >  2:315932/1‑250 (MQ=255)
                                                                                                             gCCACCGCAGCCTCGCTAAGAATGGATAAACATCGTAGCAATAGCGTCGGGCCACACGATCTTTATTTCACGTTGCTGGATGATTATCTGCATGTGGTGGATACCGCGCTGTGGTTGTCGAGCGACAAAGCCTCTCTGGATGGCGGTACGCTACTGACTAACGACGCTGGCGAAATGCTGTTTGCCGAGCACCATTTTTCGGCTGGTCCTTTGCAGATCACCACCTGTATGCATCGCCGTGCCGGACGTc                                                                                                                  >  1:49498/1‑250 (MQ=255)
                                                                                                             cccACCGCAGCCTCGCGAAGAATGGATAAACATCGGAGCAATAGCGTCGGGCCACACGATCTTTATTGCACGTTGCTGGATTATTATCTGCATGTGGTGGATCCCGCGCTGTGGTTGTCGAGCGACAAAGCCTCTCTGGATGGTGGTACGCTACTGACTAACGACGCTGGCGAAATGCTGTTTTCCGTGCACCATTTTTCGGCTGGGCCTTTGCAGATCACCACCTGTATGCATCGCCGTGCCGGAAGTCa                                                                                                                 <  1:239731/250‑1 (MQ=255)
                                                                                                                cccGCAGCCTCGCTAAGAATGGATAAACATCGTAGCAATAGCGTCGGGCCACACGATCTTTATTTCACGTTGCTGGATGATTATCTGCATGTGGTGGATACCGCGCTGTGGTTGTCGAGCGACAAAGCCTCTCTGGATGGCGGTACGCTACTGACTAACGACGCTGGCGAAATGCTGTTTGCCGAGCACCATTTTTCGGCTGGTCCTTTGCAGATCACCACCTGTATGCATCGCCGTGCCGGAAGTCAGCg                                                                                                              <  2:489513/250‑1 (MQ=255)
                                                                                                                 ccGCAGCCTCGCTAAGAATGGATAAACATCGTAGCAATAGCGTCGGGCCACACGATCTTTATTTCACGTTGCTGGATGATTATCTGCATGTGGTGGATACCGCGCTGTGGTTGTCGAGCGACAAAGCCTCTCTGGATGGCGGTACGCTACTGACTAACGACGCTGGCGAAATGCTGTTTGCCGAGCACCATTTTTCGGCTGGTCCTTTGCAGATCACCACCTGTATGCATCGCCGTGCCGGAAGTCAGCGt                                                                                                             >  1:388186/1‑251 (MQ=255)
                                                                                                                          cGCTAAGAATGGATAAACATCGTAGCAATAGCGTCGGGCCACACGATCTTTATTTCACGTTGCTGGATGATTATCTGCATGTGGTGGATACCGCGCTGTGGTTGTCGAGCGACAAAGCCTCTCTGGATGGCGGTACGCTACTGACTAACGACGCTGGCGAAATGCTGTTTGCCGAGCACCATTTTTCGGCTGGTCCTTTGCAGATCACCACCTGTATGCATCGCCGTGCCGGAAGTCAGCGTGAAACCGTg                                                                                                    <  1:163404/251‑1 (MQ=255)
                                                                                                                              aaGAATGGATAAACATCGTAGCAATAGCGTCGGGCCACACGATCTTTATTTCACGTTGCTGGATGATTATCTGCATGTGGTGGATACCGCGCTGTGGTTGTCGAGCGACAAAGCCTCTCTGGATGGCGGTACGCTACTGACTAACGACGCTGGCGAAATGCTGTTTGCCGAGCACCATTTTTCGGCTGGTCCTTTGCAGATCACCACCTGTATGCATCGCCGTGCCGGAAGTCCGCGTGAAAACGTGCAgg                                                                                                >  1:347112/1‑251 (MQ=255)
                                                                                                                                   tGGATAAACATCGTAGCAATAGCGTCGGGCCACACGATCTTTATTTCACGTTGCTGGATGATTATCTGCATGTGGTGGATACCGCGCTGTGGTTGTCGAGCGACAAAGCCTCTCTGGATGGCGGTACGCTACTGACTAACGACGCTGGCGAAATGCTGTTTGCCGAGCACCATTTTTCGGCTGGTCCTTTGCAGATCACCACCTGTATGCATCGCCGTGCCGGAAGTCAGCGTGAAACCGTGCAGGCCGtg                                                                                           >  2:441461/1‑251 (MQ=255)
                                                                                                                                       tAAACATCGTAGCAATAGCGTCGGGCCACACGCTCTTTATTTCACGTTGCTGGATGATTATCTGCATGTGGTGGATACCGCGCTGTGGTTGTCGCGCGACAAAGCCTCACTGGATGGCGGTACGCTACTGACTAACGACGCTGGCGAAATGCTGTTTGCCGAGCACCATTTTTCGGCTGGACCTTTGCAGATACCCACCTGGATGCATCGCCGTGCCGGCAGTCAGCGTGAAACAGAGCAGGCCGcgactg                                                                                       >  2:91699/1‑251 (MQ=255)
                                                                                                                                              cGTAGCAATAGCGTCGGGCCACACGATCTTTATTTCACGTTGCTGGATGATTATCTGCATGTGGTGGATACCGCGCTGTGGTTGTCGAGCGACAAAGCCTCTCTGGATGGCGGTACGCTACTGACTAACGACGCTGGCGAAATGCTGTTTGCCGAGCACCATTTTTCGGCTGGTCCTTTGCAGAACACCACCTGTATGCATCGCCGTGCCGGAAGTAAGCGCGAAACCGTGCCGGCCGTGACTGACGGTgc                                                                                >  2:353519/1‑251 (MQ=255)
                                                                                                                                                tAGCAATAGCGTCGGGCCACACGATCTTTATTTCACGTTGCTGGATGATTATCTGCATGTGGTGGATACCGCGCTGTGGTTGTCGAGCGACAAAGCCTCTCTGGATGGCGGTACGCTACTGACTAACGACGCTGGCGAAATGCTGTTTGCCGAGCACCATTTTTCGGCTGGTCCTTTGCAGATCACCACCTGTATGCATCGCCGTGCCGGAAGTCAGCGTGAAACCGTGCAGGCCGTGACTGACGGTgcgc                                                                              >  2:29308/1‑251 (MQ=255)
                                                                                                                                                   cAATAGCGTCGGGCCACACGATCTTTATTTCACGTTGCTGGATGATTATCTGCATGTGGTGGATACCGCGCTGTGGTTGTCGAGCGACAAAGCCTCTCTGGATGGCGGTACGCTACTGACTAACGACGCTGGCGAAATGCTGTTTGCCGAGCACCATTTTTCGGCTGGTCCTTTGCAGATCACCACCTGTATGCATCGCCGTGCCGGAAGTCAGCGTGAAACCGTGCAGGCCGTGACTGACGGTGCGCTc                                                                            >  1:505338/1‑250 (MQ=255)
                                                                                                                                                              ggCCACACGATCTTTATTTCACGTTGCTGGATGATTATCTGCATGGGGTGGATACCGCGCTGTGGTTGTCGAGCGACAAAGCCTCTCTGGATGGCGGTACGCTACTGACTAACGACGCTGGCGAAATGCTGTTTGCCGAGCACCATTTTTCGGCTGGTCCTTTGCAGATCACCACCTGTATGAATCGCCGGGCCGGAAGTCAGCGTGAAACAGTGCCGGCCGTGACTGACGTGGCGCTCATAGGCATTACg                                                                >  2:413339/1‑251 (MQ=255)
                                                                                                                                                                      gctCTTTATTTCACGTTGCTGGATGATTATCTGCATGTGGTGGATACCGCGCTGTGGTTGTCGAGCGACAAAGCCTCTCTGGATGGCGGTACGCTACTGACTAACGACGCTGGCGAAATGCTGTTTGCCGAGCACCATTTTTCGGCTGGTCCTTTGCAGATCACCACCTGTATGCATCGCCGTGCCGGAAGTCAGCGTGAAACCGTGCAGGCCGTGACTGACGGTGCGCCCATCGACATTACggctaggcg                                                        >  2:282239/3‑244 (MQ=255)
                                                                                                                                                                      gATCTTTATTTCACGTTGCTGGATGATTATCTGCATGTGGTGGATACCGCGCTGTGGTTGTCGAGCGACAAAGCCTCTCTGGATGGCGGTCCGCTACTGACTAACGACGCTGGCGAAATGCTGTTTGCCGAGCACCATTTTTCGGCTGGTCCTTTGCAGATCACCCCCTGTATGCATCGCCGTGCCGGATGTCAGCGTTAAACCGTTCAGGCCGTGACTGACGGTGCGCTCATCTACATTACGGATTTgcg                                                        >  2:188750/1‑251 (MQ=255)
                                                                                                                                                                          tttATTTCACGTTGCTGGATGATTATCTGCATGTGGTGGATACCGCGCTGTGGTTGTCGAGCGACAAAGCCTCTCTGGATGGCGGTACGCTACTGACTAACGACGCTGGCGAAATGCTGTTTGCCGAGCACCATTTTTCGGCTGGTCCTTTGCAGATCACCACCTGTATGCATCGCCGTGCCGGAAGTCAGCGTGAAACCGTGCAGGCCGTTACTGACGGTGCGCTCATCGACATTACGGATATGCGCGa                                                     >  1:370815/1‑250 (MQ=255)
                                                                                                                                                                           ttattaCACGTTTCTGGATGGTTATCTTCATGTGGTGGATACCGCGCTGGGGTTTGCGAGCGACAAATCCTCTCTGGATGGCGGTACGCTACGGACTAACGACGCTGGCGAAATGCTGTTTGCTGAGCACCATTTTTCGGCTGGTCCTTTGCAGATCACCACCTGTATGCTTCGCCGTGCCGGAAGTCAGCGTGAAACCGTGCAGGCCGTGACTGACGGTGCGCTCATCGACATTACGGATATgcgc                                                       <  2:309352/247‑1 (MQ=255)
                                                                                                                                                                           ttATTTCACGTTGCTGGATGATTATCTGCATGTGGTGGATACCGCGCTGTGGTTGTCGGGCGGCAAAGCCTCTCTGGATGGCGGTACGCTACTGACTAACGACGCTGGCGAAATGCTGTTTGCCGAGCACCa                                                                                                                                                                          <  2:237785/132‑1 (MQ=255)
                                                                                                                                                                           ttATTTCACGTTGCTGGATGATTATCTGCATGTGGTGGATACCGCGCTGTGGTTGTCGGGCGGCAAAGCCTCTCTGGATGGCGGTACGCTACTGACTAACGACGCTGGCGAAATGCTGTTTGCCGAGCACCa                                                                                                                                                                          >  1:237785/1‑132 (MQ=255)
                                                                                                                                                                           ttATTTCACGTTGCTGGATGATTATCTGCATGTGGTGGATACCGCGCTGTGGTTGTCGAGCGACAAAGCCTCTCTGGATGGCGGTACGCTACTGACTAACGACGCTGGCGAAATGCTGTTTGCCGAGCACCATTTTTCGGCTGGTCCTTTGCAGATCACCACCTGTATGCATCGCCGTGCCGGAAGTCAGCGTGAAACCGTGCAGGCCGTGACTGACGGTGCGCTCATCGACATTACGGATATgc                                                         <  1:8776/245‑1 (MQ=255)
                                                                                                                                                                             aTTTCACGTTGCTGGATGATTATCTGCATGTGGTGGATCCCGCGCTGTGGTTGTCGAGCGACAAAGCCTCTCTGGATGGCGGTACGCTACTGACTAACGACGCTGGCGAAATGCTGTTTGCCGAGCACCATTTTTCGGCTGGTCCTTTGCAGATCACCACCTGTATGCATCGCCGTGCCGGAAGTCAGCGTGAAACCGTGCAGGCCGTGACTGACGGTGCGCCCATCGACATTACGGATATGCGCGAATg                                                  <  2:319717/250‑1 (MQ=255)
                                                                                                                                                                                tCACGTTGCTGGATGATTATCTGCATGTGGTGGATTCCGCGCTGTGGTTGTCGAGCGACAAAGCCTCTCTGGATGGCGGTACGCTACTGACTAACGACGCTGGCGAAATGCTGTTTGCCGAGCACCATTTTTCGGCTGGTCCTTTGCAGATCACCCCCTGTATGCATCGCCGTGCCGGAAGTCAGCGTGAAACCGTGCAGGCCGTGACTGACGGTGCGCTCATCGACATTACGGATATGCGCGAATGGCg                                               <  1:40611/250‑1 (MQ=255)
                                                                                                                                                                                        cTGGATGATTATCTGCATGTGGTGGATACCGCGCTGTGGTTGTCGAGCGACAAAGCCTCTCTGGATGGCGGTACGCTACTGACTAACGACGCTGGCGAAATGCTGTTTGCCGAGCACCATTTTTCGGCTGGTCCTTTGCAGATCACCCCCTGTATGCATCGCCGTGCCGGCAGTCAGCGTGAAACCGTGCAGGCCGTGCATGACAGTGCGCTCATCAGCATTACGGATATGCGCGAATGGCGTGAGGGgcg                                      >  2:34707/1‑251 (MQ=255)
                                                                                                                                                                                        cTGGATGATTATCTGCATGTGGTGGATACCGCGCTGTGGTTGTCGAGCGACAAAGCCTCTCTGGATGGCGGTACGCTACTGACTAACGACGCTGGCGAAATGCTGTTTGCCGAGCACCATTTTTCGGCTGGTCCTTTGCAGATCACCACCTGTATGCATCGCCGTGCCGGAAGTCAGCGTGAAACCGTGCAGGCCGTGACTGACGGTGCGCTCATCGACATTACGGATATGCGCGAATGGCGAGAGGAgcg                                      >  1:148385/1‑251 (MQ=255)
                                                                                                                                                                                          ggATGATTATCTGCATGTGGTGGATACCGCGCTGTGGTTGTCGAGCGACAAAGCCTCTCTGGATGGCGGTACGCTACTGACTAACGACGCTGGCGAAATGCTGTTTGCCGAGCACCATTTTTCGGCTGGTCCTTTGCAGATCACCACCTGTATGCATCGCCGTGCCGGAAGTCAGCGTGAAACCGTGCAGGCCGTGACTGACGGTGCGCTCATCGACATTACGGATATGCGCGAATGGCGTGAGGAgcg                                      <  1:32197/249‑1 (MQ=255)
                                                                                                                                                                                           gatgatTATCTGCATGTGGTGGATACCGCGCTGTGGTTGTCGAGCGACAAAGCCTCTCTGGATGGCGGTACGCTACTGACTAACGACGCTGGCGAAATGCTGTTTGCCGCGCACCATTTTTCGGCTGGTCCTTTGCAGATCACCACCTGTATGCATCGCCGTGCCGGAAGTCAGCGTGAAACCGTGCAGGCCGTGACTGACGGTGCGCTCATCGCCATTACGGATATGCGCGAATGGCGTGAGGAGCGCgg                                   >  2:198012/1‑251 (MQ=255)
                                                                                                                                                                                                      gCATGTGGTGGATACCGCGCTGTGGTTGTCGGGCGACAATGCCTCTCGGTATGGCGGGACGCTACGGACTAACGACTCTGGCGAAATGCTGTTCGCCGAGCACCATTGTGCGGCTGGTCCTTTGCAGAGCACCACCTGTATGCATCGTCGTGCCGGAAGTCAGCGGGAAACCGTGCAGGCCGTGACTGATGGGGCGCTCATCGACATTTCGGATATGCGCGAATGGCTTGAGGAGCGCTGGCAGGGCg                           <  2:494191/248‑1 (MQ=255)
                                                                                                                                                                                                      gCATGGGGGGGAAACCGCCCTGTGGTTGTCGAGCGACAAAGCCTCTCTGGATGGCGGTACGCTACTGACTAACGACGCTGGCGAAATGCTGTTTGCCGAGCACCATTTTTCGGCTGGTCCTTTGCAGATCACCACCTGTGTGCATCGCCGTGCCGGAAGTCAGCGTGAAACCGTGCAGGCCGTGACTGACGGTGCGCTCAGCGACATTACGGATATGCGCGAATGGCGTGAGGAGCGCGGGCAGGGCg                           <  2:494196/248‑1 (MQ=255)
                                                                                                                                                                                                       cattgggtggATACCGCGCTGTGGTTGGCGAGCGACAAAGCCTCTCTGGATGGCGGTTCGCTACTGACGAACGACGCTGGCGAAATGCTGTTTGCCGAGCACCATTTTTCGGCTGGTCCTTTGCAGATCGCCGCCTGTATGCATCGCCGTGCCGGAAGTCAGCGGGAAACCGTGCAGGCCGTGACTGACGGTGCGCTCATCGACATTACGGATATGCGCGAATGGCGTGAGGAGCGCGGGCAGGGCGTAGt                       <  2:106865/246‑1 (MQ=255)
                                                                                                                                                                                                           tggtggATACCGCGCTGTGGTTGTCGAGCGACAAAGCCTCTCTGGATGGCGGTACGCTACTGACTAACGACGCTGGCGAAATGCTGTTTGCCGAGCACCATTTTTCGGCTGGTCCTTTGCAGATCACCACCTGTATGCATCGCCGTGCCGGAAGTCAGCGTGAAACCGTGCAGGCCGTGACTGACGGGGCGCTCATCGACATTACGGATATGCGCGAAGGGCGTGGGGAGCGCGGGCAGGGCGTAGTGCAt                   >  2:385524/1‑251 (MQ=255)
                                                                                                                                                                                                               ggATACCGCGCTGTGGTTGTCGAGCGACAAAGCCTCTCTGGATGGCGGTACGCTACTGACTAACGACGCTGGCGAAATGCTGTTTGCCGAGCACCATTTTTCGGCTGGTCCTTTGCAGATCACCACCTGTATGCATCGCCGTGCCGGAAGTCAGCGTGAAACCGTGCAGGCCGTGACTGACGGTGCGCTCATCGACATTACGGATATGCGCGAATGGCGTGAGGAGCGCGGGCAGGGCGTAGTGCATaaa                <  1:353519/250‑1 (MQ=255)
                                                                                                                                                                                                                    ccGCGCTGTGGTTGTCGAGCGACAAAGCCTCTCTGGATGGCGGTACGCTACTGACTAACGACGCTGGCGAAATGCTGTTTGCCGAGCACCATTTTTCGGCTGGTCCTTTGCAGATCACCACCTGTATGCATCGCCGTGCCGGAAGTCAGCGTGAAACCGTGCAGGCCGTGACTGACGGTGCGCTCATCGACATTACGGATATGCGCGAATGGCGTGAGGAGCGCGGGCAGGGCGTAGTGCATAAACCGAtt          >  1:92471/1‑251 (MQ=255)
                                                                                                                                                                                                                     cgcgCTGGGGTTGTCGAGCGACAAAGCCTCTCTGGATGGCGGTACGATACTGACTAACGACGCTGGCGAAATGCTGTTTGCCGAGCACCATTTTTCGGCTGGTCCTTTGCAGATCCCCACCTGTATGCATCGCCGTGCCGGAAGTCAGCGTGAAACCGTGCAGGCCGTGACTGACGGTGCGCTCATCGGCATTACGGATATGCGCGAATGGCGGGCGGAGGGCGGGCAGGGGGGGGTGAATAAACCGAttt         >  2:40591/1‑250 (MQ=255)
                                                                                                                                                                                                                           gggggtGTCGAGCGACAAAGCCTCTCTGGATGGCGGTACGCTACTGACTAACGCCGCTGGCGAAATGCTGTTTGCCGAGTTCCATTTTTCGGCTGGTCCTTTGCAGATCACCACCTGTATGCATCGCCGGGCCGGAAGTCAGCGTGAAACCGTGCAGGCCGTGACTGACGGTGCGCTCATCGACATTACGGATATGCGCGAATGGCGTGAGGAGCGCGGGCAGGGCGTAGTGCATAAACCGATTCCTgg     <  2:296467/244‑1 (MQ=255)
                                                                                                                                                                                                                             gtttGTCGTGCGACTAAGCCTCTCTGGATGGCGGTACGCTACGGACTATCGACGCTGACGCACTGCTGTTTGCCTAGCACCATTTTTCGGCTTGTCCTTTGCATATCCCCACCTGTATGCATCGCTGTGCCGGAATTCAGCGTGAAACCGTGCAGGCCGTGACTGACGGTGCGTTCATCGACAGTACGGATATGCGCGAATGTCGTGAGGAGCGCGGGCAGGGCGTAGTGCATAAACCGAGTCCTGGTTg  <  2:458552/248‑1 (MQ=255)
                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                 
GCGCGATGCTGAACGGCTGGTGGAACTGGCGGCGCGGAAAAAACTGACGTTGATGGTCGGTTTTAACCGTCGTTTCGCACCACTCTACGGTGAGTTAAAAACGCAACTCGCCACCGCAGCCTCGCTAAGAATGGATAAACATCGTAGCAATAGCGTCGGGCCACACGATCTTTATTTCACGTTGCTGGATGATTATCTGCATGTGGTGGATACCGCGCTGTGGTTGTCGGGCGGCAAAGCCTCTCTGGATGGCGGTACGCTACTGACTAACGACGCTGGCGAAATGCTGTTTGCCGAGCACCATTTTTCGGCTGGTCCTTTGCAGATCACCACCTGTATGCATCGCCGTGCCGGAAGTCAGCGTGAAACCGTGCAGGCCGTGACTGACGGTGCGCTCATCGACATTACGGATATGCGCGAATGGCGTGAGGAGCGCGGGCAGGGCGTAGTGCATAAACCGATTCCTGGTTG  >  NC_000913/1127105‑1127575

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