Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,126,518 C→T A121V (GCC→GTC)  yceH → DUF480 domain‑containing protein YceH

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,126,5180CT98.2% 172.9 / ‑4.5 56A121V (GCC→GTC) yceHDUF480 domain‑containing protein YceH
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base T (27/28);  total (28/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.53e-01

CAGAAAACGAACCGTGAACCGGTCATGAATCTGAGTGAATCCGAAGTGCAGGAACAGCTGGATAATCTGGTCAAACGTCATTATCTACGCACAGTGAGCGGTTTTGGTAATCGGGTCACCAAATATGAGCAACGTTTTTGTAATTCAGAATTTGGCGATCTGAAACTGAGCGCAGCGGAAGTGGCGTTAATCACCACGTTGTTATTGCGTGGTGCCCAGACGCCAGGTGAACTGCGCAGCCGCGCCGCGCGAATGTATGAATTCAGCGATATGGCGGAAGTGGAGTCGACGCTGGA‑ACAACTGGCAAATCGCGAAGATGGTCCTTTT‑GTGGTGCGTCTGGCCCGCGAACCGGGTAAACGCGAAAACCGCTACATGCATCTTTTCAGTGGTGAGGTTGAAGATCAGCCGGCGGTGACGGATATGTCGAACGCGGTTGACGGTGATTTACAGGCCCGCGTCGAAGCCCTGGAAATCGAAGTGGCAGAACT  >  NC_000913/1126274‑1126761
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                     
cAGAAAACGAACCGTGAACCGGTCATGAATCTGAGTGAATCCAAAGTGCAGGAACAGCTGGATAATCTGGTCAAACGTCATTATCTACGCACAGTGAGCGGTTTTGGTAATCGGGTCACCAAATATGAGCAACGTTTTTGTAATTCAGAATTTGGCGATCTGAAACTGAGCGCAGCGGAAGTGGCGTTAATCACCACGTTGTTATTGCGTGGTGCCCAGACGCCAGGTGAACTGCGCAGCCGCGTcgcgcg                                                                                                                                                                                                                                                 <  1:37561/251‑1 (MQ=255)
 aGAAAACGAACCGTGAACCGGTCATGAATCTGAGTGAATCCAAAGTGCAGGAACAGCTGGATAATCTGGTCAAACGTCATTATCTACGCACAGTGAGCGGTTTTGGTAATCGGGTCACCAAATATGAGCAACGTTTTTGTAATTCAGAATTTGGCGATCTGAAACTGAGCGCAGCGGAAGTGGCGTTAATCACCACGTTGTTATTGCGTGGTGCCCAGACGCCAGGTGAACTGCGCAGccgcgtcgcg                                                                                                                                                                                                                                                   >  2:334366/1‑248 (MQ=255)
    aaaCGAACCGTGAACCGGTCATGAATCTGAGTGAATCCAAAGTGCAGGAACAGCTGGATAATCTGGTCAAACGTCATTATCTACGCACAGTGAGCGGTTTTGGTAATCGGGTCACCAAATATGAGCAACGTTTTTGTAATTCAGAATTTGGCGATCTGAAACTGAGCGCAGCGGAAGTGGCGTTAATCACCACGTTGTTATTGCGTGGTGCCCAGACGCCAGGTGAACTGCGCAGCCGCGTCGCGCGAATg                                                                                                                                                                                                                                             <  1:652834/251‑1 (MQ=255)
       cGAACCGTGAACCGGTCATGAATCTGAGTGAATCCAAAGTGCAGGAACAGCTGGATAATCTGGTCAAACGTCATTATCTACGCACAGTGAGCGGTTTTGGTAATCGGGTCACCAAATATGAGCAACGTTTTTGTAATTCAGAATTTGGCGATCTGAAACTGAGCGCAGCGGAAGTGGCGTTAATCACCACGTTGTTATTGCGTGGTGCCCAGACGCCAGGTGAACTGCGCAGCCGCGTCGCGCGAATGTAt                                                                                                                                                                                                                                          >  1:291466/1‑251 (MQ=255)
         aaCCGTGAACCGGTCATGAATCTGAGTGAATCCAAAGTGCAGGAACAGCTGGATAATCTGGTCAAACGTCATTATCTACGCACAGTGAGCGGTTTTGGTAATCGGGTCACCAAATATGAGCAACGTTTTTGTAATTCAGAATTTGGCGATCTGAAACTGAGCGCAGCGGAAGTGGCGTTAATCACCACGTTGTTATTGCGTGGTGCCCAGACGCCAGGTGAACTGCGCAGCCGCGTCGCGCGAATGTATGa                                                                                                                                                                                                                                        <  1:302033/251‑1 (MQ=255)
           ccGTGAACCGGTCATGAATCTGAGTGAATCCAAAGTGCAGGAACAGCTGGATAATCTGGTCAAACGTCATTATCTACGCACAGTGAGCGGTTTTGGTAATCGGGTCACCAAATATGAGCAACGTTTTTGTAATTCAGAATTTGGCGATCTGAAACTGAGCGCAGCGGAAGTGGCGTTAATCACCACGTTGTTATTGCGTGGTGCCCAGACGCCAGGTGAACTGCGCAGCCGCGTCGCGCGAATGTATGAAt                                                                                                                                                                                                                                      <  1:329053/251‑1 (MQ=255)
                   cGGTCATGAATCTGAGTGAATCCAAAGTGCAGGAACAGCTGGATAATCTGGTCAAACGTCATTATCTACGCACAGTGAGCGGTTTTGGTAATCGGGTCACCAAATATGAGCAACGTTTTTGTAATTCAGAATTTGGCGATCTGAAACTGAGCGCAGCGGAAGTGGCGTTAATCACCACGTTGTTATTGCGTGGTGCCCAGACGCCAGGTGAACTGCGCAGCCGCGTCGCGCGAATGTATGAATTCAGCGat                                                                                                                                                                                                                              >  1:362192/1‑251 (MQ=255)
                                  gTGAATCCAAAGTGCAGGAACAGCTGGATAATCTGGTCAAACGTCATTATCTACGCACAGTGAGCGGTTTTGGTAATCGGGTCACCAAATATGAGCAACGTTTTTGTAATTCAGAATTTGGCGATCTGAAACTGAGCGCAGCGGAAGTGGCGTTAATCACCACGTTGTTATTGCGTGGTGCCCAGACGCCAGGTGAACTGCGCAGCCGCGTCGCGCGAATGTATGAATTCAGCGATATGGCGGAAGTgg                                                                                                                                                                                                                 >  1:238697/1‑249 (MQ=255)
                                        ccaaagagCAGGAACAGCTGGATAATCTGGTCAAACGTAATTATCTACGCACAGTGAGCGGTTTTGGTAATCGGGTCACCAAATATGTGCAACGTTTTTGTAATTCAGAATTTGGCGATCTGAAACTGAGCGCAGCGGAAGGGGCGTTAATCACCACGTTGTTATTGCGTGGTGCCCAGACGCCAGGAGAACTGCGCAGCCGCGTCGCGCGAATGTATGAATTCAGCGATGTGGCGGagagggggtggagg                                                                                                                                                                                                         >  2:494562/8‑238 (MQ=255)
                                        ccaaaGTGCAGGAACAGCTGGATAATCTGGTCAAACGTCATTATCTACGCACAGTGAGCGGTTTTGGTAATCGGGTCACCAAATATGAGCAACGTTTTTGTAATTCAGAATTTGGCGATCTGAAACTGAGCGCAGCGGAAGTGGCGTTAATCACCACGTTGTTATTGCGTGGTGCCCAGACGCCAGGTGAACTGCGCAGCCGCGTCGCGCGAATGTATGAATTCAGCGATATGGCGGAAGTGGAGTCGACg                                                                                                                                                                                                         >  1:607894/4‑251 (MQ=255)
                                          aaaGTGCAGGAACAGCTGGATAATCTGGTCAAACGTCATTATCTACGCACAGTGAGCGGTTTTGGTAATCGGGTCACCAAATATGAGCAACGTTTTTGTAATTCAGAATTTGGCGATCTGAAACTGAGCGCAGCGGAAGTGGCGTTAATCACCACGTTGTTATTGCGTGGTGCCCAGACGCCAGGTGAACTGCGCAGCCGCGTCGCGCGAATGTATGAATTCAGCGATATGGCGGAAGTGGAGTCGACGCt                                                                                                                                                                                                       >  1:357516/2‑251 (MQ=255)
                                               gCAGGAACAGCTGGATAATCTGGTCAAACGTCATTATCTACGCACAGTGAGCGGTTTTGGTAATCGGGTCACCAAATATGAGCAACGTTTTTGTAATTCAGAATTTGGCGATCTGAAACTGAGCGCAGCGGAAGTGGCGTTAATCACCACGTTGTTATTGCGTGGTGCCCAGACGCCAGGTGAACTGCGCAGCCGCGTCGCGCGAATGTATGAATTCAGCGATATGGCGGCAGGGGAGTCGGCGCCGGa‑a                                                                                                                                                                                                  >  2:312435/1‑250 (MQ=255)
                                                cAGGAACAGCTGGATAATCTGGTCAAACGTCATTATCTACGCACAGTGAGCGGTTTTGGTAATCGGGTCACCAAATATGAGCAACGTTTTTGTAATTCAGAATTTGGCGATCTGAAACTGAGCGCAGCGGAAGTGGCGTTAATCACCACGTTGTTATTGCGTGGTGCCCAGACGCCAGGTGAACTGCGCAGCCGCGTCGCGCGACTGTATGAATTCAGCGATATGGCGGAAggggggacgccactggaac                                                                                                                                                                                                  >  2:660975/1‑232 (MQ=255)
                                                  ggAACAGCTGGATAATCTGGTCAAACGTCATTATCTACGCACAGTGAGCGGTTTTGGTAATCGGGTCACCAAATATGAGCAACGTTTTTGTAATTCAGAATTTGGCGATCTGAAACTGAGCGCAGCGGAAGTGGCGTTAATCACCACGTTGTTATTGCGTGGTGCCCAGACGCCAGGTGAACTGCGCAGCCGCGTCGCGCGAATGTATGAATTCAGCGATAGGGCGGAAGTGGAGTCGACGCTGGa‑acaac                                                                                                                                                                                              <  1:243315/251‑1 (MQ=255)
                                                   gAACAGCTGGATAATCTGGTCAAACGTCATTATCTACGCACAGTGAGCGGTTTTGGTAATCGGGTCACCAAATATGAGCAACGTTTTTGTAATTCAGAATTTGGCGATCTGAAACTGAGCGCAGCGGAAGTGGCGTTAATCACCACGTTGTTATTGCGTGGTGCCCAGACGCCAGGTGAACTGCGCAGCCGCGTCGCGCGAATGTATGAATTCAGCGATATGGCGGAAGTGGAGTCGACGCTGGA‑ACAACt                                                                                                                                                                                             <  1:353657/251‑1 (MQ=255)
                                                      cAGCTGGATAATCTGGTCAAACGTCATTATCTACGCACAGTGAGCGGTTTTGGTAATCGGGTCACCAAATATGAGCAACGTTTTTGTAATTCAGAATTTGGCGATCTGAAACTGAGCGCAGCGGAAGTGGCGTTAATCACCACGTTGTTATTGCGTGGTGCCCAGACGCCAGGGGAACTGCGCAGCCGCGTCGCGCGAGTGGATGAATTCAGCGATATGGCGGAAgggggggcgcacctgggcaaacaggc                                                                                                                                                                                           >  2:147296/1‑226 (MQ=255)
                                                           ggATAATCTGGTCAAACGTCATTATCTACGCACAGTGAGCGGTTTTGGTAATCGGGTCACCAAATATGAGCAACGTTTTTGTAATTCAGAATTTGGCGATCTGAAACTGAGCGCAGCGGAAGTGGCGTTAATCACCACGTTGTTATTGCGTGGTGCCCAGACGCCAGGTGAACTGCGCAGCCGCGTCGCGCGAATGTATGAATTCAGCGATATGGCGGAAGTGGAGTCGACGCTGGA‑ACAACTGGCAAATc                                                                                                                                                                                     <  2:362192/251‑1 (MQ=255)
                                                           ggATAATCTGGTCAAACGTCATTATCTACGCACAGTGAGCGGTTTTGGTAATCGGGTCACCAAATATGAGCAACGTTTTTGTAATTCAGAATTTGGCGATCTGAAACTGAGCGCAGCGGAAGTGGCGTTAATCACCACGTTGTTATTGCGTGGTGCCCAGACGCCAGGTGAACTGCGCAGCCGCGTCGCGCGAATGTATGAATTCAGCGATATGGCGGAAGTGGAGTCGACGCTGGA‑ACAACTGGCAAATc                                                                                                                                                                                     <  1:634179/251‑1 (MQ=255)
                                                           ggATAATCGGGTCAAACGTCATTAGCTACGCACAGTGAGCGGTTTTGGTAATCGGGTCACCAAATATGAGCAACGTTTTTGTAATTCAGAATTTGGCGATCTGAAACTGAGCGCAGCGGAAGTGGCGTTAATCACCACGTTGTTATTGCGTGGTGCCCAGACGCCAGGTGAACTGCGCAGCCGCGTCGCGCGAATGTATGAATTCAGCGATATGGCGGAAGTGGAGTCGACGCTGGA‑ACAACTGGCAAATc                                                                                                                                                                                     <  2:412939/251‑1 (MQ=255)
                                                                aTCTGGTCAAACGTCATTATCTACGCACAGTGAGCGGTTTTGGTAATCGGGTCACCAAATATGAGCAACGTTTTTGTAATTCAGAATTTGGCGATCTGAAACTGAGCGCAGCGGAAGTGGCGTTAATCACCACGTTGTTATTGCGTGGTGCCCAGACGCCAGGTGAACTGCGCAGCCGCGTCGCGCGAATGTATGAATTCAGCGATATGGCGGAAGTGGAGTCGACGCTGGA‑ACAACTGGCAAATCGCGa                                                                                                                                                                                 >  1:241770/1‑250 (MQ=255)
                                                                 tCTGGACAAACGTCATTATCTACGCACAGTGAGCGGTTTTGGTAATCGGGTCACCAAATATGAGCAACGTTTTTGTAATTCAGAATTTGGCGATCTGAAACTGAGCGCAGCGGAAGTGGCGTTAATCACCACGTTGTTATTGCGTGGTGCCCAGCCGCCAGGTGAACTGCGCAGCCGCGTCGCGCGAATGTATGAATTCAGCGATCTGGCGGAAGTGGGGACGACGCTGGA‑ACAACTGGCCAATcgcgg                                                                                                                                                                                 >  2:222993/1‑248 (MQ=255)
                                                                     gTCAAACGTCATTATCTACGCACAGTGAGCGGTTTTGGTAATCGGGTCACCAAATATGAGCAACGTTTTTGTAATTCAGAATTTGGCGATCTGAAACTGAGCGCAGCGGAAGTGGCGTTAATCACCACGTTGTTATTGCGTGGTGCCCAGACGCCAGGTGAACTGCGCAGCCGCGTCGCGCGAATGTATGAATTCAGCGATATGGCGGAAGTGGAGTCGACGCTGGA‑ACAACTGGCAAATCGCGAAGATgg                                                                                                                                                                           <  1:177166/251‑1 (MQ=255)
                                                                         aaCGTCATTATCTACGCACAGGGAGCGGTTTTGGTAATCGGGTCACCAAATATGAGCAACGTTGTTGTAATTCAGAATTTGGCGATCTGAAACTGAGCGCAGCGGAAGTGGCGTTAATCACCACGTTGTTATTGCGTGGTGCCCAGACGCCAGGTGAACTGCGCAGCCGCGTCGCGCGAATGTATGAATGCAGCGATATGGCGGAAGTGGAGTCGACGCTGGA‑ACAACTGGCAAATCGCGAAGATGGTCCt                                                                                                                                                                       <  2:37290/251‑1 (MQ=255)
                                                                              cATTATCTACGCACAGTGAGCGGTTTTGGTAATCGGGTCACCAAATATGAGCAACGTTTTTGTAATTCAGAATTTGGCGATCTGAAACTGAGCGCAGCGGAAGTGGCGTTAATCACCACGTTGTTATTGCGTGGTGCCCAGACGCCAGGTGAACTGCGCAGCCGCGTCGCGCGAATGTATGAATTCAGCGATATGGCGGAAGGGGAGTCGACGCTGGG‑AAAACTGGCAAATCGCGAAGATGGTCCTTTt‑g                                                                                                                                                                  >  2:572443/1‑250 (MQ=255)
                                                                               aTTATCTACGCACAGGGAGCGGTTTTGGTAATCGGGTCACCAAATATGAGCAACGTTTTTGGAATTCAGAATTTGGCGAGCTGAAACTGAGCGCAGCGGAAGTGGCGGTAATCACCACGTTGTTATTGCGTGGTGCCCAGACGCCAGGTGAACTGCGCAGCCGCGTCGCGCGAATGTATGAATTCAGCGATATGGCGGAAGTGGAGTCGACGCTGGA‑ACAACTGGCAAATCGCGAAGATGGTCCTTTt‑gtg                                                                                                                                                                <  2:143653/251‑1 (MQ=255)
                                                                                        gaacaGGGAGAGGTTTTGGTAATCGTGTCACCAAATATGAGCAACGTTTTTGTAATTCAGAATTTGGCGATCTGAAACTGAGCGCAGCGGAAGTGGCGTTAATCACCACGTTGTTATTGCGTGGTGCCCAGACGCCAGGTGAACTGCGCAGCCGCGTCGCGCGAATGTATGAATTCAGCGATATGGCGGAAGTGGAGTCGACGCTGGA‑ACAACTGGCAAATCGCGAAGATGGTCCTTTT‑GTGGTGCGTCTg                                                                                                                                                       <  2:92596/249‑1 (MQ=255)
                                                                                              tGAGCGGTTTTGGTAATCGGGTCACCAAATATGAGCAACGTTTTTGTAATTCAGAATTTGGCGATCTGAAACTGAGCGCAGCGGAAGTGGCGTTAATCACCACGTTGTTATTGCGTGGTGCCCAGACGCCAGGTGAACTGCGCAGCCGCGTCGCGCGAATGTATGAATTCAGCGATATGGCGGAAGGGGAGTCGGCGCTGGG‑ACAACTGGCAAATCGAGGAGATGGGCCTTTT‑TTGGTGCGTCTGGGCcgc                                                                                                                                                 >  2:639242/1‑251 (MQ=255)
                                                                                                  cGGTTTTGGTAATCGGGTCACCAAATATGAGCAACGTTTTTGTAATTCAGAATTTGGCGATCTGAAACTGAGCGCAGCGGAAGTGGCGTTAATCACCACGTTGTTATTGCGTGGTGCCCAGACGCCAGGTGAACTGCGCAGCCGCGTCGCGCGAATGTATGAATTCAGCGATATGGCGGAAGTGGAGTCGACGCTGGA‑ACAACTGGCAAATCGCGAAGATGGTCCTTTT‑GTGGTGCGTCTGGCCCGCGAAc                                                                                                                                             <  1:429706/251‑1 (MQ=255)
                                                                                                      tttGGTACACGGGTCACCAAATATGAGCAACGTTTTTGTAATTCAGAATTTGGCGATCCGAAACTGAGCGCAGCGGAAGTGGCGGTAATCACCACGTTGTTATTGCGTGGTGCCCAGACGCCAGGTGAACTGCGCAGCCGCGTCGCGCGAATGTATGAATTCAGCGATATGGCGGAAGTGGAGTCGGCGCTGGA‑ACAACTGGCAAATCGCGAAGATGGTCCTTTT‑GTGGTGCGTCTGGCCCGCGAACCgg                                                                                                                                          <  2:328012/250‑1 (MQ=255)
                                                                                                               cGGGTCACCAAATATGAGCAACGTTTTTGTAATTCAGAATTTGGCGATCTGAAACTGAGCGCAGCGGAAGTGGCGTTAATCACCACGTTGTTATTGCGTGGTGCCCAGACGCCAGGTGAACTGCGCAGCCGCGTCGCGCGAATGTATGAATTCAGCGATATGGCGGAAGTGGAGTCGACGCTGGA‑ACAACTGGCAAATCGCGAAGATGGTCCTTTT‑GTGGTGCGTCTGGCCCGCGAACCGGGTAAAcgcg                                                                                                                                 >  1:115152/1‑250 (MQ=255)
                                                                                                                        aaaTATGAGCAACGTTTTTGTAATTCAGAATTTGGCGATCTGAAACTAAGCGCAGCGGAAGTGGCGTTAATCACCACGTTGTTATTGCGTGGTGCCCAGACGCCAGGTGAACTGCGCAGCCGCGTCGCGCGAATGTATGAATTCAGCGATATGGCGGAAGAGGAGTCGCCGCTGGG‑ACAAACGGCAAATCGCGAAGATGGGCCTTTT‑TTGGTGCGTCTGGGCCGCGAACCGGGTAAAAGCGGAAAACGCt                                                                                                                        >  2:639531/1‑250 (MQ=255)
                                                                                                                              gAGCAACGTTTTTGTAATTCAGAATTTGGCGATCTGAAACTGAGCGCAGCGGAAGTGGCGTTAATCACCACGTTGTTATTGCGTGGTGCCCAGACGCCAGGTGAACTGCGCAGCCGCGTCGCGCGAATGTATGAATTCAGCGATATGGCGGAAGTGGAGTCGACGCTGGA‑ACAACTGGCAAATCGCGAAGATGGTCCTTTT‑GTGGTGCGTCTGGCCCGCGAACCGGGTAAACGCGAAAACCGCTACATGCa                                                                                                                 <  1:222993/251‑1 (MQ=255)
                                                                                                                                        tttGTAATTCAGAATTTGGCGATCTGAAACTGAGCGCAGCGGAAGTGGCGTTAATCACCACGTTGTTATTGCGTGGTGCCCAGACGCCAGGTGAACTGCGCAGCCGCGTCGCGCGAATGTATGAATTCAGCGATATGGCGGAAGTGGAGGCGGCGCTGGG‑ACAAATGGCAAATCGCGAAGATGGTCCTTTT‑GTGGTGCGTCTGGCCCGCGAACCGGGTAAACGCGAAAACCGCTACATGCATCTTTTCAgg                                                                                                       >  1:459047/1‑250 (MQ=255)
                                                                                                                                          tGTAATTCAGAATTTGGCGATCTGAAACTGAGCGCAGCGGAAGTGGCGTTAATCACCACGTTGTTATTGCGTGGTGCCCAGACGCCAGGTGAACTGCGCAGCCGCGTCGCGCGAATGTATGAATTCAGCGATATGGCGGAAGTGGAGTCGACGCTGGA‑ACAACTGGCAAATCGCGAAGATGGTCCTTTT‑GTGGTGCGTCTGGCCCGCGAACCGGGTAAACGCGAAAACCGCTACATGCATCTTTTCAgtgg                                                                                                     >  1:621538/1‑251 (MQ=255)
                                                                                                                                          tGTAATTCAGAATTTGGCGATCTGAAACTGAGCGCAGCGGAAGGGGCGTTAATCACCACGTTGTGATGGCGTGGTGCCCAGACGCCAGGTGAACGGCGCAGCCGCGTCGCGCGAATGTATGAATTCAGCGATATGGCGGAAGTGGAGTCGACGCTGGA‑ACAACTGGCAAATCGCGAAGATGGTCCTTTT‑GTGGTGCGTCTGGCCCGCGAACCGGGTAAACGCGAAAACCGCTACATGCACCTTTTCAgtg                                                                                                      <  2:621538/250‑1 (MQ=255)
                                                                                                                                            tAATTCAGAATTTGGCGATCTGAAACTGAGCGCAGCGGAAGTGGCGTTAATCACCACGTTGTTATTGCGTGGTGCCCAGACGCCAGGTGAACTGCGCAGCCGCGTCGCGCGAATGTATGAATTCAGCGATATGGCGGAAGTGGAGTCGACGCTGGA‑ACAACTGGCAAATCGCGAAGATGGTCCTTTT‑GTGGTGCGTCTGGCCCGCGAACCGGGTAAACGCGAAAACCGCTACATGCATCTTTTCAgtggtg                                                                                                   >  1:182273/1‑251 (MQ=255)
                                                                                                                                              aTTCAGAATTTGGCGATCTGAAACTGAGCGCAGCGGAAGTGGCGTTAGTCAAAACGCTGTTATTGCGTGGAGCCCAGACGCCAGGAGAACTGCGCAGCCGCGTCGCGCGACTGTATGAATTCAGCGATATGGCGGAAGTGGAGTCGGCACTGGA‑ACCACTGGAAAATCGCGAAGATGGTTCTTTT‑GTGTAGCGTCTGGCCCGCGAACCGGGAAAACGCGAAAACCGCTACATGCATCTTTTCAggtgtg                                                                                                   >  1:412953/1‑244 (MQ=255)
                                                                                                                                                  aGAATTTGGCGATCTGAAACTGAGCGCAGCGGAAGTGGCGTTAATCACCACGTTGTTATTGCGTGGTGCCCAGACGCCAGGTGAACTGCGCAGCCGCGTCGCGCGAATGTATGAATTCAGCGATATGGCGGAAGTGGAGTCGACGCTGGA‑ACAACTGGCAAATCGCGAAGATGGTCCTTTT‑GTGGTGCGTCTGGCCCGCGAACCGGGTAAACGCGAAAACCGCTACATGCCTCTTTTCAGTGGTGAGGTTg                                                                                             >  1:44830/1‑251 (MQ=255)
                                                                                                                                                     aTTTGGCGATCTGAAACTGAGCGCAGCGGAAGTGGCGTTAATCACCACGTTGTTATTGCGTGGTGCCCAGACGCCAGGTGAACTGCGCAGCCGCGCCGCGCGAATGTATGAATTCAGCGATATGGCGGAAGTGGAGTCGACGCTGGA‑ACAACTGGCAAATCGCGAAGATGGTCCTTTT‑GTGGTGCGTCTGGCCCGCGAACCGGGTAAACGCGAAAACCGCTACATGCATCTTTTCAGTGGTGAGGTTGa                                                                                            >  1:67016/1‑249 (MQ=255)
                                                                                                                                                           cGATCTGAAACTGAGCGCAGCGGAAGTGGCGTTAATCACCACGGTGATATTGCGTGGTGCCCAGACGACAGGTGAACTGCGCAGACGCGTCGCGCGAATGTATGAATTCAGCGATATGGCGGAAGTGGGGTCGGCGGTCGAGACA‑CTGGCCAAACGCGAATAGGGTCCTTTT‑GTGGGGGGTCTGGCCCGGGAAACGGGGTAACGCGCCAACCACGAAATGCGTCTTTTCAATGGAGAgggtgctgagcagc                                                                                    >  2:73738/1‑239 (MQ=255)
                                                                                                                                                             aTCTGAAACTGAGCGCAGCGGAAGTGGCGTTAATCACCACGTTGTTATTGCGTGGTGCCCAGACGCCAGGTGAACTGCGCAGCCGCGTCGCGCGAATGTATGAATTCAGCGATATGGCGGAAGTGGAGTCGACGCTGGA‑ACAACTGGCAAATCGCGAAGATGGTCCTTTT‑GTGGTGCGTCTGGCCCGCGAACCGGGTAAACGCGAAAACCGCTACATGCATCTTTTCAGTGGTGGGGTTGAAGATCAGCcg                                                                                  >  1:561599/1‑251 (MQ=255)
                                                                                                                                                                       gAGCGCAGCGGAAGTGGCGTTAATCACCACGTTGTTATTGCGTGGTGCCCAGACGCCAGGTGAACTGCGCAGCCGCGTCGCGCGAATGTATGAATTCAGCGATATGGCGGAAGTGGAGTCGACGCTGGA‑ACAACTGGCAAATCGCGAAGATGGTCCTTTT‑GTGGTGCGTCTGGCCCGCGAACCGGGTAAACGCGAAAACCGCTACATGCATCTTTTCAGTGGTGAGGTTGAAGATCAGCCGGCGGTGAc                                                                          <  1:639242/249‑1 (MQ=255)
                                                                                                                                                                        aGCGCAGCGGAAGTGGCGTTAATCACCACGTTGTTATTGCGTGGTGCCCAGACGCCAGGTGAACTGCGCAGCCGCGTCGCGCGAATGTATGAATTCAGCGATATGGCGGAAGTGGAGTCGACGCTGGA‑ACAACTGGCAAATCGCGAAGATGGTCCTTTT‑GTGGTGCGTCTGGCCCGCGAACCGGGTAAACGCGAAAACCGCTACATGCATCTTTTCAGTGGTGAGGTTGAAGATCAGCCGGCGGTGACGGa                                                                       <  1:639531/251‑1 (MQ=255)
                                                                                                                                                                            cAGCGGAAGTGGCGTTAATCACCACGTTGTTATTGCGTGGCGCAAAGACGCAAGGTGAACTGCGCAGCAGCGTCGCGCGAATGGATGAATTCAGCGATATGGCGGAAGGGGAGACGAAGCAGGA‑ACAACTGGAAAAACGCGACGATGGTCCTTTTTATGG‑GCGTCTGGACCACGAACCGCGTACACGAGCAAACCTCGACATGCCTTCTTTCCGTGGTGAGGGTGAAGAtacatcggccgggtgggaaagt                                                                   >  2:237381/1‑230 (MQ=255)
                                                                                                                                                                               cGGAAGTGGCGTTAATCACCACGTTGTTATTGCGTGGTGCCCAGACGCCAGGTGAACTGCGCAGCCGCGTCGCGCGAATGTATGAATTCAGCGATATGGCGGAAGTGGAGTCGACGCTGGA‑ACAACTGGCAAATCGCGAAGATGGTCCTTTT‑GTGGTGCGTCTGGCCCGCGAACCGGGTAAACGCGAAAACCGCTACATGCATCTTTTCAGTGGTGAGGTTGAAGATCAGCCGGCGGTGACGGATATGTCg                                                                <  1:312435/251‑1 (MQ=255)
                                                                                                                                                                                ggAAGTGGCGTTAATCACCACGTTGTTATTGCGTGGTGCCCAGACGCCAGGTGAACTGCGCAGCCGCGTCGCGCGAATGTATGAATTCAGCGATATGGCGGAAGTGGAGTCGACGCTGGA‑ACAACTGGCAAATCGCGAAGATGGTCCTTTT‑GGGGTGCGTCTGGCCCGCGAACCGGGTAAACGCGAAAACCGCTACATGCATCTTTTCAGTGGTGAGGTTGAAGATCAGCCGGCGGTGCCGGATATGTCg                                                                >  1:530442/1‑250 (MQ=255)
                                                                                                                                                                                  aaGTGGCGTTAATCACCACGTTGTTATTGCGTGGTGCCCAGACGCCAGGTGAACTGCGCAGCCGCGTCGCGCGAATGTATGAATTCAGCGATATGGCGGAAGTGGAGTCGACGCTGGA‑ACAACTGGCAAATCGCGAAGATGGTCCTTTT‑GTGGTGCGTCTGGCCCGCGAACCGGGTAAACGCGAAAACCGCTACATGCATCTTTTCAGTGGTGAGGTTGAAGATCAGCCGGCGGTGACGGATATGTCGAAc                                                             <  2:241770/251‑1 (MQ=255)
                                                                                                                                                                                                     gttggtagtGCGTGGTGCCCAGACGCCAGGTGAACTGCGCAGCCGCGTCGCGCGAATGTATGAATTCAGCGATATGGCGGAAGTGGAGTCGACGCTGGA‑ACAACTGGCAAATCGCGAAGATGGTCCTTTT‑GTGGTGCGTCTGGCCCGCGAACCGGGTAAACGCGAAAACCGCTACATGCATCTTTTCAGTGGTGAGGTTGAAGATCAGCCGGCGGTGACGGATATGTCGAACGCGGTTGACGGTGATTTAc                                          <  2:238697/243‑1 (MQ=255)
                                                                                                                                                                                                        gttATTGCGTGGTGCCCAGACGCCAGGTGAACTGCGCAGCCGCGTCGCGCGAATGTATGAATTCAGCGATATGGCGGAAGTGGAGTCGACGCTGGA‑ACAACTGGCAAATCGCGAAGATGGTCCTTTT‑GTGGTGCGTCTGGCCCGCGAACCGGGTAAACGCGAAAACCGCTACATGCATCTTTTCAGTGGTGAGGTTGAAGATCAGCCGGCGGGGACGGGTATGTCGAACGCGGTTGACGGGGATTTACAgg                                       >  1:310732/1‑251 (MQ=255)
                                                                                                                                                                                                          tATTGCGTGGTGCCCAGACGCCAGGTGAACTGCGCAGCCGCGTCGCGCGAATGTATGAATTCAGCGATATGGCGGAAGTGGAGTCGACGCTGGA‑ACAACTGGCAAATCGCGAAGATGGTCCTTTT‑GTGGTGCGTCTg                                                                                                                                                       >  2:106894/1‑137 (MQ=255)
                                                                                                                                                                                                          tATTGCGTGGTGCCCAGACGCCAGGTGAACTGCGCAGCCGCGTCGCGCGAATGTATGAATTCAGCGATATGGCGGAAGTGGAGTCGACGCTGGA‑ACAACTGGCAAATCGCGAAGATGGTCCTTTT‑GTGGTGCGTCTg                                                                                                                                                       <  1:106894/137‑1 (MQ=255)
                                                                                                                                                                                                              gcggggtgCCCAGGCGCCAGGTTAACTGCGCAGCCGCGTCGCGCGAATGTATGAATTCAGCGATATGGCGGAAGTGGGGGCGACGCTTGA‑GCAACTGGCAAATCGCGAAGATTGTCCTTTT‑GTGGTGCGTCTGGCCCGCGAACCGGGGAAACGCGAAAACCGCTACATGCATCTTTTCAGTGGTGAGGTTGAAGATCAGCCGGCGGTGACGGATATGTCGAACGCGGTTGACGGTGATTTACAGGCCcgcg                                 <  2:561599/247‑1 (MQ=255)
                                                                                                                                                                                                                 gggtgCCCAGCCGCCAGGTGAACGGCGCAGCCGCGTCGCGCGAATGTATGAATTCAGCGATATGGCGGAAGTGGAGTCGACGCTGGA‑ACAACTGGCAAATCGAGAAGATGGTCCTTTT‑GAGGTGCGTCTGGCCCGCGAACCGGGTAAACGCGAAAACCGCTACATGCATCTTTTCAGTGGTGAGGTTGAAGATCAGCCGGCGGTGACGGATATGTCGAACGCGGTTGACGGTGATTTACAGGCCCGCGTCg                              <  2:235451/250‑1 (MQ=255)
                                                                                                                                                                                                                           aCGCCAGGTGAACTGCGCAGCCGCGTCGCGCGAATGTATGAATTCAGCGATATGGCGGAAGTGGAGGCGACGCTGGA‑ACAACTGGCAAATCGCGAAGATGGTCCTTTT‑GTGGTGCGTCTGGCCCGCGAACCGGGTAAACGCGAAAACCGCTACATGCATCTTTTCAGTGGTGAGGTTGAAGATCAGCCGGCGGTGACGGATATGTCGAACGCGGTTGACGGTGATTTACAGGCCCGCGTCGAAGCCCTGGa                    <  1:494562/251‑1 (MQ=255)
                                                                                                                                                                                                                               cAGGTGAACTGCGCAGCCGCGTCGCGCGAATGTATGAATTCAGCGATATGGCGGAAGTGGAGTCGACGCTGGA‑ACAACTGGCAAATCGCGAAGATGGTCCTTTT‑GTGGTGCGTCTGGCCCGCGAACCGGGTAAACGCGAAAACCGCTACATGCATCTTTTCAGTGGTGAGGTTGAAGATCAGCCGGCGGTGACGGATATGTCGAACGCGGTTGACGGTGATTTACAGGCCCGCGTCGAAGCCCTGGAAATc                >  1:37916/1‑251 (MQ=255)
                                                                                                                                                                                                                                           gcAGCCGCGTCGCGCGAATGTATGAATTCAGCGATATGGCGGAAGTGGAGTCGACGCTGGA‑ACAACTGGCAAATCGCGAAGATGGTCCTTTT‑GTGGTGCGTCTGGCCCGCGAACCGGGTAAACGCGAAAACCGCTACATGCATCTTTTCAGTGGTGAGGTTGAAGATCAGCCGGCGGTGACGGATATGTCGAACGCGGTTGACGGTGATTTACAGGCCCGCGTCGAAGCCCTGGAAATCGAAGTGGCAGaa    <  1:572443/251‑1 (MQ=255)
                                                                                                                                                                                                                                             aGCCGCGTCGCGCGAATGTATGAATTCAGCGATATGGCGGAAGTGGAGTCGACGCTGGA‑ACAACTGGCAAATCGCGAAGATGGTCCTTTT‑GTGGTGCGTCTGGCCCGCGAACCGGGTAAACGCGAAAACCGCTACATGCATCTTTTCAGTGGTGAGGTTGAAGATCAGCCGGCGGTGACGGATATGTCGAACGCGGTTGACGGTGATTTACAGGCCCGCGTCGAAGCCCTGGAAATCGAAGTGGCAGAACt  <  1:73738/251‑1 (MQ=255)
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                     
CAGAAAACGAACCGTGAACCGGTCATGAATCTGAGTGAATCCGAAGTGCAGGAACAGCTGGATAATCTGGTCAAACGTCATTATCTACGCACAGTGAGCGGTTTTGGTAATCGGGTCACCAAATATGAGCAACGTTTTTGTAATTCAGAATTTGGCGATCTGAAACTGAGCGCAGCGGAAGTGGCGTTAATCACCACGTTGTTATTGCGTGGTGCCCAGACGCCAGGTGAACTGCGCAGCCGCGCCGCGCGAATGTATGAATTCAGCGATATGGCGGAAGTGGAGTCGACGCTGGA‑ACAACTGGCAAATCGCGAAGATGGTCCTTTT‑GTGGTGCGTCTGGCCCGCGAACCGGGTAAACGCGAAAACCGCTACATGCATCTTTTCAGTGGTGAGGTTGAAGATCAGCCGGCGGTGACGGATATGTCGAACGCGGTTGACGGTGATTTACAGGCCCGCGTCGAAGCCCTGGAAATCGAAGTGGCAGAACT  >  NC_000913/1126274‑1126761

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