Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,179,001 C→T S41S (AGC→AGT rplA → 50S ribosomal subunit protein L1

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,179,0010CT98.6% 228.7 / NA 71S41S (AGC→AGTrplA50S ribosomal subunit protein L1
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base T (32/38);  total (33/38)
Fisher's exact test for biased strand distribution p-value = 4.65e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.80e-01

CGCTCAGCTGCAGGAAATCGCGCAGACCAAAGCTGCCGACATGACTGGTGCCGACATTGAAGCGATGACTCGCTCCATCGAAGGTACTGCACGTTCCATGGGCCTGGTAGTGGAGGACTAAGAAATGGCTAAACTGACCAAGCGCATGCGTGTTATCCGCGAGAAAGTTGATGCAACCAAACAGTACGACATCAACGAAGCTATCGCACTGCTGAAAGAGCTGGCGACTGCTAAATTCGTAGAAAGCGTGGACGTAGCTGTTAACCTCGGCATCGACGCTCGTAAATCTGACCAGAACGTACGTGGTGCAACTGTACTGCCGCACGGTACTGGCCGTTCCGTTCGCGTAGCCGTATTTACCCAAGGTGCAAACGCTGAAGCTGCTAAAGCTGCAGGCGCAGAACTGGTAGGTATGGAAGATCTGGCTGACCAGATCAAGAAAGGCGAAATGAACTTTGACGTTGTTATTGCTTCTCCGGATGCAATGCGCGTT  >  NC_000913/4178755‑4179247
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                      
agcTCAGCTGCAGGAAATCGCGCAGACCAAAGCTGCCGACATGACTGGTGCCGACATTGAAGCGATGACTCGCTCCATCGAAGGTACTGCACGTTCCATGGGCCTGGTAGTGGAGGACTAAGAAATGGTTAAACTGACCAAGCGCATGCGTGTTATCCGCGAGAAAGTTGATGCAACCAAACAGTACGACCTCAACGAAGCTATCGCACTGCTGAAAGCGCTGGCGACTGCTAAATTCGTAGAAAGTGTgg                                                                                                                                                                                                                                                    >  2:18883/2‑251 (MQ=255)
 gcTCAGCTGCAGGAAATCGCGCAGACCAAAGCTGCCGACATGACTGGTGCCGACATTGAAGCGATGACTCGCTCCATCGAAGGTACTGCACGTTCCATGGGCCTGGTAGTGGAGGACTAAGAAATGGTTAAACTGACCAAGCGCATGCGTGTTATCCGCGAGAAAGTTGATGCAACCAAACAGTACGACATCAACGAAGCTATCGCACTGCTGAAAGAGCTGGCGACTGCTAAATTCGTAGAAAGTGTgg                                                                                                                                                                                                                                                    >  1:419610/1‑250 (MQ=255)
 gcTCAGCTGCAGGAAATCGCGCAGACCAAAGCTGCCGACATGACTGGTGCCGACATTGAAGCGATGACTCGCTCCATCGAAGGTACTGCACGTTCCATGGGCCTGGTAGTGGAGGACTAAGAAATGGTTAAACTGACCAAGCGCATGCGTGTTATCCGCGAGAAAGTTGATGCAACCAAACAGTACGACATCAACGAAGCTATCGCACTGCTGAAAGAGCTGGCGACTGCTAAATTCGTAGAAAGTGTgg                                                                                                                                                                                                                                                    >  1:87611/1‑250 (MQ=255)
   tCAGCTGCAGGAAATCGCGCAGACCAAAGCTGCCGACATGACTGGTGCCGACATTGAAGCGATGACTCGCTCCATCGAAGGTACTGCACGTTCCATGGGCCTGGTAGTGGAGGACTAAGAAATGGTTAAACTGACCAAGCGCATGCGTGTTATCCGCGAGAAAGTTGATGCAACCAAACAGTACGACATCAACGAAGCTATCGCACTGCTGAAAGAGCTGGCGACTGCTAAATTCGTAGAAAGTGTGGACg                                                                                                                                                                                                                                                 >  1:303821/1‑251 (MQ=255)
    cAGCTGCAGGAAATCGCGCAGACCAAAGCTGCCGACATGACTGGTGCCGACATTGAAGCGATGACTCGCTCCATCGAAGGTACTGCACGTTCCATGGGCCTGGTAGTGGAGGACTAAGAAATGGTTAAACTGACCAAGCGCATGCGTGTTATCCGCGAGAAAGTTGATGCAACCAAACAGTACGACATCAACGAAGCTATCGCACTGCTGAAAGAGCTGGCGACTGCTAAATTCGTAGAAAGTGTGGAc                                                                                                                                                                                                                                                  <  1:253864/249‑1 (MQ=255)
         gCAGGAAATCGCGCAGACCAAAGCTGCCGACATGACTGGTGCCGACATTGAAGCGATGACTCGCTCCATCGAAGGTACTGCACGTTCCATGGGCCTGGTAGTGGAGGACTAAGAAATGGTTAAACTGACCAAGCGCATGCGTGTTATCCGCGAGAAAGTTGATGCAACCAAACAGTACGACATCAACGAAGCTATCGCACTGCTGAAAGAGCTGGCGACTGCTAAATTCGTAGAAAGTGTGGACGTAGCTg                                                                                                                                                                                                                                           <  2:274445/251‑1 (MQ=255)
           aGGAAATCGCGCAGACCAAAGCTGCCGACATGACTGGTGCCGACATTGAAGCGATGACTCGCTCCATCGAAGGTACTGCACGTTCCATGGGCCTGGTAGTGGAGGACTAAGAAATGGTTAAACTGACCAAGCGCATGCGTGTTATCCGCGAGAAAGTTGATGCAACCAAACAGTACGACATCAACGAAGCTATCGCACTGCTGAAAGAGCTGGCGACTGCTAAATTCGTAGAAAGTGTGGACGTAGCTGtt                                                                                                                                                                                                                                         <  2:205185/251‑1 (MQ=255)
           aGGAAATCGCGCAGACCAAAGCTGCCGACATGACTGGTGCCGACATTGAAGCGATGACTCGCTCCATCGAAGGTACTGCACGTTCCATGGGCCTGGTAGTGGAGGACTAAGAAATGGTTAAACTGACCAAGCGCATGCGTGTTATCCGCGAGAAAGTTGATGCAACCAAACAGTACGACATCAACGAAGCTATCGCACTGCTGAAAGAGCTGGCGACTGCTAAATTCGTAGAAAGTGTGGACGTAGCTGtt                                                                                                                                                                                                                                         <  1:189239/251‑1 (MQ=255)
                   gcgcAGACCAAAGCTGCCGACATGACTGGTGCCGACATTGAAGCGATGACTCGCTCCATCGAAGGTACTGCACGTTCCATGGGCCTGGTAGTGGAGGACTAAGAAATGGTTAAACTGACCAAGCGCATGCGTGTTATCCGCGAGAAAGTTGATGCAACCAAACAGTACGACATCAACGAAGCTATCGCACTGCTGAAAGAGCTGGCGACTGCTAAATTCGTAGAAAGTGTGGACGTAGCTGTTAACCTCgg                                                                                                                                                                                                                                 <  1:58817/251‑1 (MQ=255)
                               gCTGCCGACATGACTGGTGCCGACATTGAAGCGATGACTCGCTCCATCGAAGGTACTGCACGTTCCATGGGCCTGGTAGTGGAGGACTAAGAAATGGTTAAACTGACCAAGCGCATGCGTGTTATCCGCGAGAAAGTTGATGCAACCAAACAGTACGACATCAACGAAGCTATCGCACTGCTGAAAGAGCTGGCGCCTGCTAAATTCGTAGAAAGTGTGGACGTAGCTGTTAACCTCGGCATCGACGCTCg                                                                                                                                                                                                                     >  1:413701/1‑251 (MQ=255)
                                       cATGACTGGTGCCGACATTGAAGCGATGACTCGCTCCATCGAAGGTACTGCACGTTCCATGGGCCTGGTAGTGGAGGACTAAGAAATGGTTAAACTGACCAAGCGCATGCGTGTTATCCGCGAGAAAGTTGATGCAACCAAACAGTACGACATCAACGAAGCTATCGCACTGCTGAAAGAGCTGGCGACTGCTAAATTCGTAGAAAGTGTGGACGTAGCTGTTAACCTCGGCATCGACGCTCGTAAAtatg                                                                                                                                                                                                             <  1:48262/251‑4 (MQ=255)
                                                  ccGACATTGAAGCGATGACTCGCTCCATCGAAGGTACTGCACGTTCCATGGGCCTGGTAGTGGAGGACTAAGAAATGGTTAAACTGACCAAGCGCATGCGGGTTATCCGCGAGAAAGTTGATGCAACCAAACAGTACGACATCAACGAAGCTATCGCACTGCTGAAAGAGCTGGCGACTGCTAAATTCGTAGAAAGTGTGGACGTAGCTGTTAACCGCGGCATCGACGCTCGTAAATCTGACCAGAacgta                                                                                                                                                                                                  <  2:419610/251‑1 (MQ=255)
                                                          gAAGCGATGACTCGCTCCATCGAAGGTACTGCACGTTCCATGGGCCTGGTAGTGGAGGACTAAGAAATGGTTAAACTGACCAAGCGCATGCGTGTTATCCGCGAGAAAGTTGATGCAACCAAACAGTACGACATCAACGAAGCTATCGCACTGCTGAAAGAGCTGGCGACTGCTAAATTCGTAGAAAGTGTGGACGTAGCTGTTAACCTCGGCATCGACGCTCGTAAATCTGACCAGAACGTACGTGGTGc                                                                                                                                                                                          >  2:450465/1‑251 (MQ=255)
                                                              cGATGACTCGCTCCATCGAAGGTACTGCACGTTCCATGGGCCTGGTAGTGGAGGACTAAGAAATGGTTAAACTGACCAAGCGCATGCGTGTTATCCGCGAGAAAGTTGATGCAACCAAACAGTACGACATCAACGAAGCTATCGCACTGCTGAAAGAGCTGGCGACTGCTAAATTCGTAGAAAGTGTGGACGTAGCTGTTAACCTCGGCATCGACGCTCGTAAATCTGACCAGAACGTACGTGGTGCAACt                                                                                                                                                                                      <  1:378709/251‑1 (MQ=255)
                                                              cGATGACTCGCTCCATCGAAGGTACTGCACGTTCCATGGGCCTGGTAGTGGAGGACTAAGAAATGGTTAAACTGACCAAGCGCATGCGTGTTATCCGCGAGAAAGTTGATGCAACCAAACAGTACGACATCAACGAAGCTATCGCACTGCTGAAAGAGCTGGCGACTGCTAAATTCGTAGAAAGTGTGGACGTAGATGTTAACCTCGGCATCGACGCTCGTAAATCTGACCAGAACGTACGTGGTGCAACt                                                                                                                                                                                      >  1:412065/1‑251 (MQ=255)
                                                                aTGACTCGCTCCATCGAAGGTAATGCACGTTCCATGGGCCTGGTAGTGGAGGACTAAGAAATGGTTAAACTGACCAAGCGCATGCGTGTTATCCGCGAGAAAGTTGATGCAACCAAACAGTACGACATCAACGAAGCTATCGCACTGCTGAAAGAGCTGGCGACTGCTAAATTCGTAGAAAGTGTGGACGTAGCTGTTAACCTCGGCATCGACGCTCGTAAATCTGACCAGAACGTACGTGGTGCAACTGt                                                                                                                                                                                    >  2:87559/1‑251 (MQ=255)
                                                                       gCTCCATCGAAGGTACTGCACGTTCCATGGGCCTGGTAGTGGAGGACTAAGAAATGGTTAAACTGACCAAGCGCATGCGTGTTATCCGCGAGAAAGTTGATGCAACCAAACAGTACGACATCAACGAAGCTATCGCACTGCTGAAAGAGCTGGCGACTGCTAAATTCGTAGAAAGTGTGGACGTAGCTGTTAACCTCGGCATCGACGCTCGTAAATCTGACCAGAACGTACGTGGTGCAACTGTACTGCCg                                                                                                                                                                             <  2:164667/251‑1 (MQ=255)
                                                                              cGAAGGTACTGCACGTTCCATGGGCCTGGTAGTGGAGGACTAAGAAATGGTTAAACTGACCAAGCGCATGCGTGTTATCCGCGAGAAAGTTGATGCAACCAAACAGTACGACATCAACGAAGCTATCGCACTGCTGAAAGAGCTGGCGACTGCTAAATTCGTAGAAAGTGTGGACGTAGCTGTTAACCTCGGCATCGACGCTCGTAAATCTGACCAGAACGTACGTGGTGCAACTGTACTGCCGCACGGt                                                                                                                                                                       <  1:235346/250‑1 (MQ=255)
                                                                                     aCTGCACGTTCCATGGGCCTGGTAGTGGAGGACTAAGAAATGGTTAAACTGACCAAGCGCATGCGTGTTATCCGCGAGAAAGTTGATGCAACCAAACAGTACGACATCAACGAAGCTATCGCACTGCTGAAAGAGCTGGCGACTGCTAAATTCGTAGAAAGTGTGGACGTAGCTGTTAACCTCGGCATCGACGCt                                                                                                                                                                                                                       >  2:19981/1‑195 (MQ=255)
                                                                                     aCTGCACGTTCCATGGGCCTGGTAGTGGAGGACTAAGAAATGGTTAAACTGACCAAGCGCATGCGTGTTATCCGCGAGAAAGTTGATGCAACCAAACAGTACGACATCAACGAAGCTATCGCACTGCTGAAAGAGCTGGCGACTGCTAAATTCGTAGAAAGTGTGGACGTAGCTGTTAACCTCGGCATCGACGCt                                                                                                                                                                                                                       <  1:19981/195‑1 (MQ=255)
                                                                                        gCACGTTCCATGGGCCTGGTAGTGGAGGACTAAGAAATGGTTAAACTGACCAAGCGCATGCGTGTTATCCGCGAGAAAGTTGATGCAACCAAACAGTACGACATCAACGAAGCTATCGCACTGCTGAAAGAGCTGGCGACTGCTAAATTCGTAGAAAGTGTGGACGTAGCTGTTAACCTCGGCATCGACGCTCGTAAATCTGACCAGAACGTACGTGGTGCAACTGTACTGCCGCACGGTACTGGccgttc                                                                                                                                                            <  1:229400/251‑1 (MQ=255)
                                                                                                cATGGGCCTGGTAGTGGAGGACTAAGAAATGGTTAAACTGACCAAGCGCATGCGTGTTATCCGCGAGAAAGTTGATGCAACCAAACAGTACGACATCAACGAAGCTATCGCACTGCTGAAAGAGCTGGCGACTGCTAAATTCGTAGAAAGTGTGGACGTAGCTGTTAACCTCGGCATCGACGCTCGTAAATCTGACCAGAACGTACGTGGTGCAACTGTACTGCCGCACGGTACTGGCCGTTCCGTTcgcg                                                                                                                                                    <  1:3402/251‑1 (MQ=255)
                                                                                                     gCCTGGTAGTGGAGGACTAAGAAATGGTGAAACGTACCAAGCGCATGCGGGTTATCCGCGAGAAAGTTGATGCACCCAAACAGTACGACATCAACGAAGCTATCGCACTGCTGAAAGAGCGGGCGACTGCTAAAGTCGTAGAAAGTGTGGACGTAGCTGTTAACCTCGGCATCGACGCTCGTAAATCTGACCAGAACGTACGTGGTGCAACGGTACTGCCGCACGGTACTGGCCGTGCCGTTCGCGTAGcc                                                                                                                                               <  2:413701/251‑1 (MQ=255)
                                                                                                         ggTAGTGGAGGACTAAGAAATGGTTAAACTGACCAAGCGCATGCGTGTTAGCCGCGAGAAAGTTGAGGCAACCAAACAGTACGACATCAACGAAGCTATCGCACTGCTGAAAGAGCTGGCGACGGCTAAATTCGTAGAAAGTGTGGACGTAGCTGTTAACCTCGGCATCGACGCTCGTAAATCGGACCAGAACGGACGTGGTGCAACTGTTCTGCCGCACGGTTCTGGCCGGTCCGTTCGCGTAGCCGt                                                                                                                                             <  2:325649/249‑1 (MQ=255)
                                                                                                          gTAGTGGAGGACTAAGAAATGGTGAAACGGACCAAGCGAATGCGTGTTATCCGCGAGAAAGTTGATGCAACCAAACAGTACGACATCAACGAAGCTATCGCACTGCTGAAAGAGCTGGCGACTGCTAAATTCGTAGAAAGTGTGGACGTAGCTGTTACCCTCGGCATCGACGCTCGTAAATCTGACCAGAACGTACGGGGTGCAACGTTACTGCTGCACGGTACTGGCCGTTCCGTTAGCGTAGCCGTAtt                                                                                                                                          <  2:361736/251‑1 (MQ=255)
                                                                                                           taggggaggaCTAAGAAATGCTTAAAATGACCCAGCGCCTGCGTGTTATCCGCGAGAAAGTTGATGCAACCAAACAGTACGACATCAACGAAGCTATCGCACTGCTGCAAGAGCTGGCGACTGCTAAATTCGTAGAAAGTGTGGACGAAGCTGTTAACCTCGGAATCGACGCTCGTAAATCTGACCAGAACGTACGTGGTGCAACTGTACTGCCGCACGGTTCTGGCCGTTCCGTTCGCGTAGCCGTAttt                                                                                                                                         >  2:168335/5‑251 (MQ=255)
                                                                                                            aGTGGAGGACTAAGAAATGGTTAAACTGACCAAGCGCATGCGTGTTATCCGCGAGAAAGTTGATGCAACCAAACAGTACGACATCAACGAAGCTATCGCACTGCTGAAAGAGCTGGCGACTGCTAAATTCGTAGAAAGTGTGGACGTAGCTGTTAACCTCGGCATCGACGCTCGTAAATCTGACCAGAACGTACGTGGTGCAACTGTACTGCCGCACGGTACTGGCCGTTCCGTTCGCGTAGCCGTAttt                                                                                                                                         >  2:273931/1‑250 (MQ=255)
                                                                                                                  ggaCTAAGAAATGGTTAAACTGACCAAGCGCATGCGTGTTATCCGCGAGAAAGTTGATGCAACCAAACAGTACGACATCAACGAAGCTATCGCACTGCTGAAAGAGCTGGCGACTGCTAAATTCGTAGAAAGTGTGGACGTAGCTGTTAACCTCGGCATCGACGCTCGTAAATCTGACCAGAACGTACGTGGTGCAACTGTACTGCCGCACGGTACTGGCCGTTCCGTTCGCGTAGCCGTATTTACCCAAg                                                                                                                                  <  2:303821/251‑1 (MQ=255)
                                                                                                                      tAAGAAATGGTTAAACTGACCAAGCGCATGCGTGTTATCCGCGAGAAAGTTGATGCAACCAAACAGTACGACATCAACGAAGCTATCGCACTGCTGAAAGAGCTGGCGACTGCTAAATTAGTAGCAAGTGTGGACGTAGCTGTTAACCTCGGCATCGACGCTCGTAAATCTGACCAGAACGTACGTGGTGCAACTGTACTGCCGCACGGTACTGGCCGTTACGATCGCGTAGCCGTAGTTACCCacagggc                                                                                                                              >  2:181076/1‑245 (MQ=255)
                                                                                                                             tggttAAACTGACCACGCGCCTGCGTGTTATCCGAGAGAAAGTTGATGCAACCAAACAGTACGCCATCAACGCAGCTAACGCACCGCTGAAAGAGCTGGCGACTGCTACATTCGCAGCAAGTGTGGACGTAGCTGTTAACCTCGGCATCGACGCTCGTAAATCTGACCAGACCGTACGTGGTGCAACTGTACTGCCGCACGGTACTGGCCGTTCCGTTCGCGTAGCCGTATTTACCCAAGGTGCAAACGCt                                                                                                                       >  1:408556/5‑251 (MQ=255)
                                                                                                                             tggttAAACTGACCAAGCGCATGCGTGTTATCCGCGAGAAAGTTGATGCAACCAAACAGTACGACATCAACGAAGCTATCGCACTGCTGAAAGAGCTGGCGACTGCTAAATTCGTAGAAAGTGTGGACGTAGCTGTTAACCTCGGCATCGACGCTCGTAAATCTGACCAGAACGTACGTGGTGCAACTGTACTGCCGCACGGTACTGGCCGTTCCGTTCGCGTAGCCGTATTTACCCAAGGTGCAAACGCt                                                                                                                       >  1:410428/5‑251 (MQ=255)
                                                                                                                                 tAAACTGACCAAGCGCATGCGTGTTATCCGCGAGAAAGTTGATGCAACCAAACAGTACGACATCAACGAAGCTATCGCACTGCTGAAAGAGCTGGCGACTGCTAAATTCGTAGAAAGTGTGGACGTAGCTGTTAACCTCGGCATCGACGCTCGTAAATCTGACCAGAACGTACGTGGTGCAACTGTACTGCCGCACGGTACTGGCCGTTCCGTTCGCGTAGCCGTATTTACCCAAGGTGCAAACGCTGaa                                                                                                                    <  1:450465/250‑1 (MQ=255)
                                                                                                                                   aaCTGACCAAGCGCATGCGTGTTATCCGCGAGAAAGTTGATGCAACCAAACAGTACGACATCAACGAAGCTATCGCACTGCTGAAAGAGCTGGCGACTGCTAAATTCGTAGAAAGTGTGGACGTAGCTGTTAACCTCGGCATCGACGCTCGTAAATCTGACCAGAACGTACGTGGTGCAACTGTACTGCCGCACGGTACTGGCCGTTCCGTTCGCGTAGCCGTATTTACCCAAGGTGCAAACGCTGAAgct                                                                                                                 <  1:277537/251‑1 (MQ=255)
                                                                                                                                    aCTGACCAAGCGCATGCGTGATATCCGCGAGAAAGTTGATGCAACCAAACAGTACGCCATCAACGAAGCTATCGCACTGCTGAAAGCGCTGGCGACTGCTAAATTCGTAGAAAGTGTGGACGTAGCTGTTAACCTCGGCATCGACGCTCGTAAATCTGACCAGAACGTACGTGGTGCAACTGTACTGCCGCACGGTACTGGCCGTTCCGTTCGCGTAGCCGAATTTACCCACGGTGCAAACGCTGAAgctg                                                                                                                >  2:461963/1‑251 (MQ=255)
                                                                                                                                      tGACCAAGCGCATGCGTGTTATCCGCGAGAAAGTTGATGCAACCAAACAGTACGACATCAACGAAGCTATCGCACTGCTGAAAGAGCTGGCGACTGCTAAATTCGTAGAAAGTGTGGACGTAGCTGTTAACCTCGGCATCGACGCTCGTAAATCTGACCAGAACGTACGTGGTGCAACTGTACTGCCGCACGGTACTGGCCGTTCCGTTCGCGTAGCCGTATTTACCCAAGGGGCAAACGCTGAAgctgct                                                                                                              >  2:273088/1‑251 (MQ=255)
                                                                                                                                      ggtccAAGCGCATGCGTGTTATCCGCGAGAAAGTGGATGCAACCAAACAGGACGACATCAACGAAGCTATCGCACTGCTGAAAGAGCTGGCGACTGCTAAATTCGTAGAAAGTGGTGACGTAGCTGTTAACCTCGGCATCGACGCTCGTAAATCTGTCCAGAACGGACGTGGTGCAACTGTACTGCCGCACGGTACTGGCCGTTCCGTTCGCGTAGCCGTATTTACCCAAGGTGCAAACGCTGAAgctgct                                                                                                              <  1:390114/248‑1 (MQ=255)
                                                                                                                                         ccAAGCGCATGCGTGTTATCCGCGAGAAAGTTGATGCAACCAAACAGTACGACATCAACGAAGCTATCGCACTGCTGAAAGAGCTGGCGACTGCTAAATTCGTAGAAAGTGTGGACGTAGCTGTTAACCTCGGCATCGACGCTCGTAAATCTGACCAGAACGTACGTGGTTCAACTGTACTGCCGCACGGTACTGGCCGTTCCGTTCGCGTAGCCGTATTTACCCAAGGTGCAAACGCTGAAGCTGCTaaa                                                                                                           <  2:386196/251‑1 (MQ=255)
                                                                                                                                          cAAGCGCATGCGTGTTATCCGCGCGAAAGTTGATGCAACAAAACAGTACGACATCAACGAAGCTATCGCACTGCTGAAAGTGCTGGCGACTGCTAAATTCGTAGAAAGTGTGGACGTAGCTGTTAACCTCGGCATCGACGCTCGTAAATCTGACCAGAACGTACGTGGTGCAACTGTACTGCCGCACGGTACTGGCCGTTCCGTTCGCGTAGCCGTATTTACCCAAGGTGCAAACGCTGCAGCTGCTaa                                                                                                            >  2:170183/1‑249 (MQ=255)
                                                                                                                                           aaGCGCATGCGTGTTATCCGCGAGAAAGTCGCTGCAACCAAACAGTACGACATCAACGAAGCTATCGCACTGCTGAAAGAGCTGGCGACTGCTCAATTCGTCGCAAGTGTGGACGAAGCTGTTAACCTCGGCATCGACGCTCGTAAATCTGACCAGCACGTACGTGGTGCAACTGTACTGCCGCACGGTACTGGCCGTTCCGTTCGCGTAGCCGTATTTACCCACGGTGCAAACGCTGAAGCTGCTAAAGc                                                                                                         >  2:461290/1‑251 (MQ=255)
                                                                                                                                            aGCGCAGGCGTGTTATCCGCGAGAAAGTTGATGCAACCAAACAGTACGACATCAACGAAGCTATCGCACTGCTGAAAGAGCTGGCGACTGCTAAATTCGTAGAAAGTGTGGACGTAGCTGTTAACCTCGGCATCGACGCTCGTAAATCTGACCAGAACGTACGTGGTGCAACTGTACTGCCGCACGGTACTGGCCGTTCCGTTCGCGTAGCCGTATTTACCCAAGGTGCAAACGCTGAAGCTGCTAAAGCt                                                                                                        <  1:111467/251‑1 (MQ=255)
                                                                                                                                             gcgcATGCGTGTTATCCGCGAGAAAGTTGATGCAACCAAACAGTACGACATCAACGAAGCTATCGCGCTGCTGAAAGAGCTGGCGACTGCTAAATTCGTAGAAAGTGTGGACGTAGCTGTTAACCTCGGCATCGACGCTCGTAAATCTGACCAGAACGTACGTGGTGCAACTGTACTGCCGCACGGTACTGGCCGTTCCGGGCGCGTAGCCGTATTTACCCAAGGTGCAAACGCGGAAGCTGCTAAAGCTg                                                                                                       <  1:264563/251‑1 (MQ=255)
                                                                                                                                                 aTGCGGGTGATCCGCGAGAAAGTTGATGCAACCAAACAGTACGACATCAACGAAGCTATCGCACTGCTGAAAGAGCTGGCGACTGCTAAATTCGTAGAAAGTGTGGACGTAGCTGTTAACCTCGGCATCGACGCTCGTAAATCTGACCAGAACGTACGTGGTGCAACTGTACTGCCGCACGGTACTGGCCGTTCCGTTCGCGGAGCCGGATTTACCCAAGGTGCAAACGCTGATGCTGTTAAAGCTGCAgg                                                                                                   <  2:183139/251‑1 (MQ=255)
                                                                                                                                                  tGCGTGTTATCCGCGAGAAAGTTGATGCAACCAAACAGTACGACATCAACGAAGCTATCGCACTGCTGAAAGAGCTGGCGCCTGCTAAATTCGTAGAAAGTGTGGACGTAGCTGTTAACCTCGGCATCGACGCTCGTAAATCTGACCAGAACGTACGTGGTGCAACTGTACTGCCGCACGGAACTGGCCGTTCCGTTCGCGTAGCCGTATTTACCCAAGGTGCAAACGCTGAAGCTGCTAACGCTCCAggg                                                                                                  >  2:270413/1‑250 (MQ=255)
                                                                                                                                                   gCGTGTTATCCGCGAGAAAGTTGATGCAACCAAACAGTACGACATCAACGAAGCTATCGCACTGCTGAAAGAGCTGGCGACTGCTAAATTCGTAGAAAGTGTGGACGTAGCTGTTAACCTCGGCATCGACGCTCGTAAATCTGACCAGAACGTACGTGGTGCAACTGTACTGCCGCACGGTACTGGCCGTTCCGTTCGCGTAGCCGTATTTACCCAAGGTGCAAACGCTGAAGCTGCTAAAGCTGCAGgc                                                                                                  <  1:388488/250‑1 (MQ=255)
                                                                                                                                                     gtgtTATCCGCGAGAAAGTTGATGCAACCCACCAGTACGACATCAACGAAGCTATCGCACTGCTGAAAGAGCTGGCGACTGCTAAATTCGTAGAAAGTGTGGACGTAGCTGTTAACCTCGGCATCGACGCTCGTAAATCTGACCAGAACGTACGTGGTGCAACTGTACTGCCGCACGGTACTGGCCGTTCCGTTCGCGTAGCCGTATTTACCCAAGGTGCAAACGCTGAAGCTGCTAAAGCTGCAGgcgc                                                                                                >  2:308896/1‑250 (MQ=255)
                                                                                                                                                     gggtTATCCGCGAGAAAGGTGATGCAACCAAACAGTACGACATCAACGAAGCTATCGCACTGCTGAAAGAGCTGGCGACTGCTAAATTCGTAGAAAGTGTGGACGTAGCTGTTAACCTCGGCATCGACGCTCGTAAATCTGACCAGAACGTACGTGGTGCAACTGTACTGCCGCACGGTACTGGCCGTTCCGTTCGCGTAGCCGTATTTACCCAAGGTGCAAACGCTGAAGCTGCTAAAGCTGCATGCGCa                                                                                               <  2:40069/249‑1 (MQ=255)
                                                                                                                                                         tATCCGCGAGAAAGTTGATGCAACCAAACAGTACGACATCAACGAAGCTATCGCACTGCTGAAAGAGCTGGCGACTGCTAAATTCGTAGAAAGTGTGGACGTAGCTGTTAACCTCGGCATCGACGCTCGTAAATCTGACCAGAACGTACGTGGTGCAACTGTACTGCCGCACGGTACTGGCCGTTCCGTTCGCGTAGCCGTATTTACCCAAGGTGCAAACGCTGAAGCTGCTAAAGCTGCAGGCGCAGaa                                                                                            >  1:420856/1‑250 (MQ=255)
                                                                                                                                                              gcgAGAAAGTTGATGCAACCAAACAGTACGACATCAACGAAGCTATCGCACTGCTGAAAGAGCTGGCGACTGCTAAATTCGTAGAAAGTGTGGACGTAGCTGTTAACCTCGGCATCGACGCTCGTAAATCTGACCAGAACGTACGTGGTGCAACTGTACTGCCGCACGGTACTGGCCGTTCCGTTCGCGTAGCCGTATTTACCCAAGGTGCAAACGCTGAAGCTGCTAAAGCTGCAGGCGCAGAACTggta                                                                                      <  2:87611/251‑1 (MQ=255)
                                                                                                                                                                gagaAAGTTGATGCAACCAAACAGTACGACATCAACGAAGCTATCGCACTGCTGAAAGAGCTGGCGACTGCTAAATTCGTAGAAAGTGTGGACGTAGCTGTTAACCTCGGCATCGACGCTCGTAAATCTGACCAGAACGTACGTGGTGCAACTGTACTGCCGCACGGTACTGGCCGTTCCGTTCGCGTAGCCGTATTTACCCAAGGTGCAAACGCTGAAGCTGCTAAAGCTGCAGGCGCAGAACTggtagg                                                                                    <  1:219478/251‑1 (MQ=255)
                                                                                                                                                                 agaAAATTGCTGAAACCAAACAGTACGAAATCAACGAAGCTATAGCACTGCTGAAAGAGCTGGAGACTGCTAAATTCGTAGAAAGCGTGGACGTAGCTGTTAACCTCGGAATCGAAGCTCGTAAATCTGACCAGAACGTACGTGGTGCAACTGTACTGCCGCACGGTAATGGACGTTCCGTTCGCGTAGCCGTATTTACCCAAGGTGCAAAGGCTGAAGATGCTAAGGATGAAGGAGCAGAACTgggaggt                                                                                   >  2:241814/1‑251 (MQ=255)
                                                                                                                                                                              aaCCAAACAGTACGACATACACGAAGCTATCGCACTGCTGAAAGAGCTGGCGACTGCTAAATTCGTAGAAAGTGTGGACGTCGCTGTTAACCTCGGCATCGACGCTCGTAAATCTGACCAGAACGTACGGGGTGCAACTGTCCTGCCGCACGGTACTGGCCGTTCCGTTCGCGTAGCCGTATTTACCCAAGGAGCAAACGCTGAAGCTGCTAAAGCTACAGGCGCCGAACTGGTAGGTATCGAAGATCTgg                                                                      >  2:166492/1‑251 (MQ=255)
                                                                                                                                                                                    aCAGTACGACATCAACGAAGCTATCGCACTGCTGAAAGAGCTGGCGACTGCTAAATTCGTAGAAAGTGTGGACGTAGCTGTTAACCTCGGCATCGACGCTCGTAAATCTGACCAGAACGTACGTGGTGCAACTGTACTGCCGCACGGTACTGGCCGTTCCGTTCGCGTAGCCGTATTTACCCAAGGTGCAAACGCTGAAGCTGCTAAAGCTGCAGGCGCAGAACTGGTAGGTATGGAAGATCTGGCTGAcc                                                                >  1:336451/1‑251 (MQ=255)
                                                                                                                                                                                      aGTACGACATCAACGAAGCTATCGCACTGCTGAAAGAGCTGGCGACTGCTAACTTCGTAGAAAGTGTGGACGTAGCTGTTAACCTCGGCATCGACGCTCGTAAATCTGACCAGAACGTACGTGGTGCAACTGTACTGCCGCACGGTACTGGCCGTTCCGTTCGCGTAGCCGTATTTACCCAAGGTGCAAACGCTGAAGCTGCTAAAGCTGCAGGCGCAGAACTGGTAGGTATGGAAGATCTGGCTGAcc                                                                >  1:437348/1‑249 (MQ=255)
                                                                                                                                                                                         aCGACATCAACGAAGCTATCGCACTGCTGAAAGAGCTGGCGACTGCTAAATTCGTAGAAAGTGTGGACGTAGCTGTTAACCTCGGCATCGACGCTCGTAAATCTGACCAGAACGTACGTGGTGCAACTGTACTGCCGCACGGTACTGGCCGTTCCGTTCGCGTAGCCGTATTTACCCAAGGTGCAAACGCTGAAGCTGCTAAAGCTGCAGGCGCAGAACTGGTAGGTATGGAAGATCTGGCTGACCAGAt                                                            <  1:95919/250‑1 (MQ=255)
                                                                                                                                                                                         aCGACATCAACGAAGCTATCGCACTGCTGAAAGAGCTGGCGACTGCTAAATTCGTAGAAAGTGTGGACGTAGCTGTTAACCTCGGCATCGACGCTCGTAAATCTGACCAGAACGTACGTGGTGCAACTGTACTGCCGCACGGTACTGGCCGTTCCGTTCGCGTAGCCGTATTTACCCAAGGTGCAAACGCTGAAGCTGCTAAAGCTGCAGGCGCAGAACTGGTAGGTATGGAAGATCTGGCTGACCAGATc                                                           <  2:336451/251‑1 (MQ=255)
                                                                                                                                                                                           gACATCAACGAAGCTATCGCACTGCTGAAAGAGCTGGCGACTGCTAAATTCGTAGAAAGTGTGGACGTAGCTGTTAACCTCGGCATCGACGCTCGTAAATCTGACCAGAACGTACGTGGTGCAACTGTACTGCCGCACGGTACTGGCCGTTCCGTTCGCGTAGCCGTATTTACCCAAGGTGCAAACGCTGAAGCTGCTAAAGCTGCAGGCGCAGAACTGGTAGGTATGGAAGATCTGGCTGACCAGATCa                                                          >  1:2356/1‑250 (MQ=255)
                                                                                                                                                                                                   cGAAGCTAGCGCACGGCTGAAAGAGCGGGCGTCGGCTAAATTCGTAGAAAGTGTGGACGTAGCTGTTAACCTCGGCAGCGACGCTCGTAAATCTGACCAGAACGTACGTGGTGCAACTGTCCTGCCGCACGGTACTGGCCGTTCCGTTCGCGTAGCCGTATTTACCCAAGGTGCAAACGCTGAAGCTGCTAAAGCTGCAGGCGCAGAACTGGTAGGTATGGAAGATCTGGCTGACCAGATCAAGTAAGGCg                                                 <  2:88342/251‑1 (MQ=255)
                                                                                                                                                                                                              cACTGCTGAAAGAGCTGGCGACTGCTAAATTCGTAGAAAGTGTGGACGTAGCTGTTAACCTCGGCATCGACGCTCGTAAATCTGACCAGAACGTACGTGGTGCAACTGTACTGCCGCACGGTACTGGCCGTTCCGTTCGCGTAGCCGTATTTACCCAAGGTGCAAACGCTGAAGCTGCTAAAGCTGCAGGCGCAGAACTGGTAGGTATGGAAGATCTGGCTGACCAGATCAAGAAAGGCGAAATGAACttt                                      >  1:357132/1‑251 (MQ=255)
                                                                                                                                                                                                              cACTGCTGAAAGAGCTGGCGACTGCTAAATTCGTAGAAAGTGTGGACGTAGCTGTTAACCTCGGCATCGACGCTCGTAAATCTGACCAGAACGTACGGGGTGCAACTGTACTGCCGCACGGTACTGGCCGTTCCGTTCGCGTAGCCGTATTTACCCAAGGTGCAAACGCTGAAGCTGCTAAAGCTGCAGGCGCAGAACTGGTAGGTATGGAAGATCTGGCTGACCAGATCAAGAAAGGCGAAATGAAcgtt                                      <  2:413294/251‑4 (MQ=255)
                                                                                                                                                                                                                ctgctgAAAGAGCTGGCGACTGCTAAATTCGTAGAAAGTGTGGACGTAGCTGTTAACCTCGGCATCGACGCTCGTAAATCTGACCAGAACGAACGTGGTGCAACTGTACTGCCGCACGGTACTGGCCGTTCCGTTCGCGTAGCCGTATTTACCCAAGGTGCAAACGCTGAAGCTGCTAAAGCTGCAGGCGCAGAACTGGTAGGTATGGAAGATCTGGCTGACCAGATCCAGAAAGGCGCAATGAACTTTg                                     >  2:136785/1‑250 (MQ=255)
                                                                                                                                                                                                                  gctgAAAGAGCTGGCGACTGCTAAATTCGTAGAAAGTGTGGACGTAGCTGTTAACCTCGGCATCGACGCTCGTAAATCGGACCAGAACGTACGTGGTGCAACTGTACTGCCGCACGGTACTGGCCGTGCCGTTCGCGTAGCCGTATTTACCCAAGGTGCAAACGCTGAAGCTGCTAAAGCTGCAGGCGCAGAACTGGTAGGTATTGAAGATCTGGCTGACCAGATCAAGAAAGGCGAAATGAACTTTGAc                                   <  1:184733/250‑1 (MQ=255)
                                                                                                                                                                                                                   cagAAAGAGCTGGCGACTGCTCAATTCGTAGTAAGTGTGGACGTAGCTGTTAACCTCGGCATCGACGCTCGTAAATCTGACCAGAACGTACGTGGTGCAACCGTACTGACGCACGGTACTGGCCGTTCCGTTCGCGTAGCCGTATTTACCCAAGGTGCAAACGCTGAAGCTGCTAAAGCTGCAGGCGCAGCACTGGAAGGTATGGCAGATCTGGCTGACCAGATCAAGAAAGGCGAAATGAACTTTGACgt                                 >  2:465707/3‑251 (MQ=255)
                                                                                                                                                                                                                          agCTGGCGACTGCTAAATACGTAGAAAGTGTGGACGTAGCTGTTAACCTCGGCATCGACGCTCGTAAATCTGACCAGAACGTACGTGGTGCAACTGTACTGCCGCACGGTACTGGCCGTTCCGTTCGCGTAGCCGTATTTACCCAAGGTGCAAACGCTGAAGCTGCTAAAGCTGCAGGCGCAGAACTGGTAGGTATGGAAGATCTGGCTGACCAGATCAAGAAAGGCGAAATGAACTTTGACGTTGTTAtt                          <  1:409199/251‑1 (MQ=255)
                                                                                                                                                                                                                            cTGGCGACTGCTAAATTCGTAGAAAGTGTGGACGTAGCTGTTAACCTCGGCATCGACGCTCGTAAATCTGACAAGACAGTACGTGGTGCAACTGTACTGCCGCACGGTACTGGCCGTTCCGTTCGCGTAGCCGTATTTACCCAAGGTGCAAACGCTGAAGCTGCTAAAGCTGCAGGCGCAGAACTGGTAGGTATGGAAGATCTGGCTGACCAGATCAAGAAAGGCGAAATGAACTTTGACGTTGTTATTGc                        >  1:432652/1‑251 (MQ=255)
                                                                                                                                                                                                                                 gACTGCTAAATTCGTAGAAAGTGTGGACGTAGCTGTTAACCTCGGCATCGACGCTCGTAAATCTGACCAGAACGTACGTGGTGCAACTGTACTGCCGCACGGTACTGGCCGTTCCGTTCGCGTAGCCGTATTTACCCAAGGTGCAAACGCTGAAGCTGCTAAAGCTGCAGGCGCAGAACTGGTAGGTATGGAAGATCTGGCTGACCAGATCAAGAAAGGCGAAATGAACTTTGACGTTGTTATTGCTtctc                   >  1:36932/1‑251 (MQ=255)
                                                                                                                                                                                                                                    tGCTAAATTCGTAGAAAGTGTGGACGTAGCTGTTAACCTCGGCATCGACGCTCGTAAATCTGACCAGAACGTACGTGGTGCAACTGTACTGCCGCACGGTACTGGCCGTTCCGTTCGCGTAGCCGTATTTACCCAAGGTGCAAACGCTGAAGCTGCTAACGCTGCAGGCGCAGAACTGGTAGGTATGGAAGATCTGGCTGACCAGATCAAGAAAGGCGAAATGAACTTTGACGTTGTTATTGCTTCTCCgg                >  1:465099/1‑251 (MQ=255)
                                                                                                                                                                                                                                     gCTAAATTCGTAGAAAGTGTGGACGTAGCTGTTAACCTCGGCATCGACGCTCGTAAATCTGACCAGAACGTACGTGGTGCAACTGTACTGCCGCACGGTACTGGCCGTTCCGTTCGCGTAGCCGTATTTACCCAAGGTGCAAACGCTGAAGCTGCTAAAGCTGCAGGCGCAGAACTGGTAGGTATGGAAGATCTGGCTGACCAGATCAAGAAAGGCGAAATGAACTTTGACGTTGTTATTGCTTCTCCg                 <  1:87559/249‑1 (MQ=255)
                                                                                                                                                                                                                                          aTTCGTAGAAAGTGTGGACGTAGCTGTTAACCTCGGCATCGACGCTCGTAAATCTGACCAGAACGTACGTGGTGCAACTGTACTGCCGCACGGTACTGGCCGTTCCGTTCGCGTAGCCGTATTTACCCAAGGTGCAAACGCTGAAGCTGCTAAAGCTGCAGGCGCAGAACTGGTAGGTATGGAAGATCTGGCTGACCAGATCAAGAAAGGCGAAATGAACTTTGACGTTGTTATTGCTTCTCCGGATGCaa          <  2:102228/251‑1 (MQ=255)
                                                                                                                                                                                                                                              gTAGAAAGTGTGGACGTAGCTGTTAACCTCGGCATCGACGCTCGTAAATCTGACCAGAACGTACGTGGTGCAACTGTACTGCCGCACGGTACTGGCCGTTCCGTTCGCGTAGCCGTATTTACCCAAGGTGCAAACGCTGAAGCTGCTAAAGCTGCAGGCGCAGAACTGGTAGGTATGGAAGATCTGGCTGACCAGATCAAGAAAGGCGAAATGAACTTTGACGTTGTTATTGCTTCTCCGGATGCAATgcg      <  2:196605/251‑1 (MQ=255)
                                                                                                                                                                                                                                                 gAAAGTGTGGACGTAGCTGTTAACCTCGGCATCGACGCTCGTAAATCTGACCAGAACGTACGTGGTGCAACTGTACTGCCGCACGGTACTGGCCGTTCCGTTCGCGTAGCCGTATTTACCCAAGGTGCAAACGCTGAAGCTGCTAAAGCTGCAGGCGCAGAACTGGTAGGTATGGAAGATCTGGCTGACCAGATCAAGAAAGGCGAAATGAACTTTGACGTTGTTATTGCTTCTCCGGATGCAATGCGCgt   <  1:208945/251‑1 (MQ=255)
                                                                                                                                                                                                                                                  aaaGTGTGGACGTAGCTGTTAACCTCGGCATCGACGCTCGTAAATCTGACCAGAACGTACGTGGTGCAACTGTACTGCCGCACGGTACTGGCCGTTCCGTTCGCGTAGCCGTATTTACCCAAGGTGCAAACGCTGAAGCTGCTAAAGCTGCAGGCGCAGAACTGGTAGGTATGGAAGATCTGGCTGACCAGATCAAGAAAGGCGAAATGAACTTTGACGTTGTTATTGCTTCTCCGGATGCAATGCGCgtt  >  2:370011/1‑251 (MQ=255)
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                      
CGCTCAGCTGCAGGAAATCGCGCAGACCAAAGCTGCCGACATGACTGGTGCCGACATTGAAGCGATGACTCGCTCCATCGAAGGTACTGCACGTTCCATGGGCCTGGTAGTGGAGGACTAAGAAATGGCTAAACTGACCAAGCGCATGCGTGTTATCCGCGAGAAAGTTGATGCAACCAAACAGTACGACATCAACGAAGCTATCGCACTGCTGAAAGAGCTGGCGACTGCTAAATTCGTAGAAAGCGTGGACGTAGCTGTTAACCTCGGCATCGACGCTCGTAAATCTGACCAGAACGTACGTGGTGCAACTGTACTGCCGCACGGTACTGGCCGTTCCGTTCGCGTAGCCGTATTTACCCAAGGTGCAAACGCTGAAGCTGCTAAAGCTGCAGGCGCAGAACTGGTAGGTATGGAAGATCTGGCTGACCAGATCAAGAAAGGCGAAATGAACTTTGACGTTGTTATTGCTTCTCCGGATGCAATGCGCGTT  >  NC_000913/4178755‑4179247

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: