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,0010CT96.4% 157.6 / ‑6.0 58S41S (AGC→AGT) rplA50S ribosomal subunit protein L1
Reads supporting (aligned to +/- strand):  ref base C (2/0);  new base T (25/29);  total (28/30)
Fisher's exact test for biased strand distribution p-value = 2.28e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCGCTCAGCTGCAGGAAATCGCGCAGACCAAAGCTGCCGACATGACTGGTGCCGACATTGAAGCGATGACTCGCTCCATCGAAGGTACTGCACGTTCCATGGGCCTGGTAGTGGAGGACTAAGAAA‑TGGCTAAACTGACCAAGCGCATGCGTGTT‑ATCCGCGAGAAAGTTGATGCAACCAAACAGTACGACATCA‑ACGAAGCTATCGCACTGCTGAAAGAGCTGGCGACTGCTAAATTCGTAGAAAGCGTGGACGTAGCTGTTAACCTCGGCATCGACGCTCGTAAATCTGACCAGAACGTACGTGGTGCAACTGTACTGCCGCACGGTACTGGCCGTTCCGTTCGCGTAGCCGTATTTACCCAAGGTGCAAACGCTGAAGCTGCTAAAGCTGCAGGCGCAGAACTGGTAGGTATGGAAGATCTGGCTGACCAGATCAAGAAAGGCGAAATGAACTTTGACGTTGTTATTGCTTCTCCGGATGCAATGCGCGT  >  NC_000913/4178754‑4179246
                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                     
gcgcTCAGCTGCAGGAAAGCGCGCAGACCAAAGCTGCCGACATGACTGGTGCCGACATGGAAGCGATGACGCGCTCCATCGAAGGTACTGCACGTTCCATGGGCCTGGTAGTGGAGGACTAAGAAA‑TGGTTAAACTGACCAAGCGCATGCGTGTT‑ATCCGCGAGAAAGTTGATGCAACCAAACAGTACGACATCA‑ACGAAGCTATCGCACTGCTGAAAGAGCTGGCGACTGCTAAATTCGTAGAAAGtgtg                                                                                                                                                                                                                                                    <  1:395553/251‑1 (MQ=255)
     cAGCTGCAGGAAATCGCGCAGACCAAAGCTGCCGACATGACTGGTGCCGACATTGCAGCGATGACTCGCTCCATCGAAGGTACTGCACGTTCCATGGGCCTGGTAGTGGAGGACTAAGAAA‑TGGTTAAACTGACCAAGCGCATGCGTGTT‑ATCCGCGAGAAAGTTGATGCAAACAAACAGTACGACATCA‑ACGAAGCTATCGCACTGCTGAAAGCGCTGGCGACTGCTAACTTCGAAGAAAGTGTGGACGt                                                                                                                                                                                                                                               >  1:56805/1‑251 (MQ=255)
       gCTGCAGGAAATCGCGCAGACCAAAGCTGCCGACATGACTGGTGCCGACATTGAAGCGATGACTCGCTCCATCGAAGGTACTGCACGTTCCATGGGCCTGGTAGTGGAGGACTAAGAAA‑TGGTTAAACTGACCAAGCGCATGCGTGTT‑ATCCGCGAGAAAGTTGATGCAACCAAACAGTACGACATCA‑ACGAAGCTATCGCACTGCTGAAAGAGCTGGCGACTGCTAAATTCGTAGAAAGTGTGGACGTAg                                                                                                                                                                                                                                             <  1:307080/251‑1 (MQ=255)
                                               ggTGCCGACATTGAAGCGATGACTCGCTCCATCGAAGGTACTGCACGTTCCATGGGCCTGGTAGTGGAGGACTAAGAAA‑TGGTTAAACTGACCAAGCGCATGCGTGTT‑ATCCGCGAGAAAGTTGATGCAACCAAACAGTACGACATCA‑ACGAAGCTATCGCACTGCTGAAAGAGCTGGCGACTGCTAAATTCGTAGAAAGTGTGGACGTAGCTGTTAACCTCGGCATCGACGCTCGTAAATCTGACCAGa                                                                                                                                                                                                      >  1:335803/1‑250 (MQ=255)
                                                gTGCCGACATTGAAGCGATGACTCGCTCCATCGAAGGTACTGCACGTTCCATGGGCCTGGTAGTGGAGGACTAAGAAA‑TGGTTAAACTGACCAAGCGCATGCGTGTT‑ATCCGCGAGAAAGTTGATGCAACCAAACAGTACGACATCA‑ACGAAGCTATCGCACTGCTGAAAGAGCTGGCGACTGCTAAATTCGTAGAAAGTGTGGACGTAGCTGTTAACCTCGGCATCGACGCTCGTAAATCTGACCAGAac                                                                                                                                                                                                    <  1:167661/251‑1 (MQ=255)
                                                      aCATTGCAGCGCTGATTCGCCCCATCACGGGTACCGCATGTTCCATGTGCCTGGTAGTGGAGAACTAAGAGA‑TGGTTAAACTGCCCAAGCGCATGCGTGTT‑ATCCGCGAGAAAGTTTATGCACCCACACAGTCCGACATCA‑TCGCAGCTATCGCACTGCTGAAAGACCTGGCGACTGCTACATCCGTAGCAAGCGTGGACGTAGCTGTAAGCATCGGCATCGACACTCgcacatccaaccagaccgtaagt                                                                                                                                                                                              >  2:46630/1‑229 (MQ=255)
                                                       cATTGAAGCGATGACTCGCTCCATCGAAGGTACTGCACGTTCCATGGGCCTGGTAGTGGAGGACTAAGAAA‑TGGTTAAACTGACCAAGCGCATGCGTGTT‑ATCCGCGAGAAAGTTGATGCAACCAAACAGTACGACATCA‑ACGAAGCTATCGCACTGCTGAAAGAGCTGGCGACTGCTAAATTCGTAGAAAGTGTGGACGTAGCTGTTAACCTCGGCATCGACGCTCGTAAATCTGACCAGAACGTACgtg                                                                                                                                                                                             >  1:88253/1‑251 (MQ=255)
                                                         ttgacgagctGACTTGCTCCGTCGAAGGGACTGCACCAGCCATGATCATGGTAGCGGAGGACTA‑GAAACTGGTTAAAATGACCACGCGCATGCGCGTT‑ATAAGCGAGAAAGATGATGCAACCAAACAGTACGACATCA‑CCGAAGCTATCGCACTGCTGAAAGAGCTGGCGACTGCTAAATTCGCAGAAAGTGTGGACGACGATGTTAACCTCGGCATAGAAGCTAGTACATCTGACCAGAACGTACgtggg                                                                                                                                                                                           >  2:189759/10‑250 (MQ=255)
                                                           gAAGCGATGACTCGCTCCATCGAAGGTACTGCACGTTCCATGGGCCTGGTAGTGGAGGACTAAGAAA‑TGGTTAAACTGACCAAGCGCATGCGTGTT‑ATCCGCGAGAAAGTTGATGCAACCAAACAGTACGACATCA‑ACGAAGCTATCGCACTGCTGAAAGAGCTGGCGACTGCTAAATTCGTAGAAAGTGTGGACGTAGCTGTTAACCTCGGCATCGACGCTCGTAAATCTGACCAGAACGTACgtggtg                                                                                                                                                                                          <  2:41474/250‑1 (MQ=255)
                                                              tctttGACTCGCTCCATCGAAGGTACTGCACGTTTCATGGGCCTGGTAGTGGAGGACTAAGAAA‑TGGTTAAACTGACCAAGCGCATGCGTGTT‑ATCCGCGAGAAAGTTGCTGCACCCAAACAGTACGACCTCA‑ACGAAGCTATCGCACTGCTGAAAGAGCTGGCGACTGCTAAATTCGTAGCAAGTGTGGACGTAGCTGTTAACCTCGGCATCGACGCTCGTCAATCTGTCCCGACAGTACGTGGTGCacc                                                                                                                                                                                      >  2:296845/5‑249 (MQ=255)
                                                                               cGAAGGTACTGCGCGTTCCATGGGCCTGGTAGTGGAGGACTAAGAAA‑TGGTTAAACTGACCAAGCGCATGCGTGTT‑ATCCGCGAGAAAGTTGATGCAACCAAACAGTACGACATCA‑ACGAAGCTATCGCACTGCTGAAAGAGCTGGCGACTGCTAAATTCGTAGAAAGTGTGGACGTAGCTGTTAACCTCGGCATCGACGCTCGTAAATCTGACCAGAACGTACGTGGTGCAACTGTACTGCCGCACGGTa                                                                                                                                                                     <  2:34532/251‑1 (MQ=255)
                                                                                        tGCACGTTCCATGGGCCTGGTAGTGGAGGACTAAGAAA‑TGGTTAAACTGACCAAGCGCATGCGTGTT‑ATCCGCGAGAAAGTTGATGCAACCAAACAGTACGACATCA‑ACGAAGCTATCGCACTGCTGAAAGAGCTGGCGACTGCTAAATTCGTAGAAAGTGTGGACGTAGCTGTTAACCTCGGCATCGACGCTCGTAAATCTGACCAGAACGTACGTGGTGCAACTGTACTGCCGCACGGTACTGGccgtt                                                                                                                                                            >  1:235106/1‑251 (MQ=255)
                                                                                               gccATGGGCCTGGTAGTGGAGGACTAAGAAA‑TGGTTAAACTGACCAAGCGCAGGCGTGTT‑ATCCGCGAGAAAGTTGATGCAACCAAACAGTACGACATCA‑ACGAAGCTATCGCACTGCTGAAAGAGCTGGCGACTGCTAAATTCGTAGAAAGTGTGGACGTAGCTGTTAACCTCGGCATCGACGCTCGTAAATCTGACCAGAACGTACGTGGTGCAACTGTACTGCCGCACGGTACTGTCcggtccgttc                                                                                                                                                      <  2:52115/249‑1 (MQ=255)
                                                                                                  aTGGGCCTGGTAGTGGAGGACTAAGAAA‑TGGTTAAACTGACCAAGCGCATGCGTGTT‑ATCCGCGAGAAAGTTGATGCAACCAAACAGTACGACATCA‑ACGAAGCTATCGCACTGCTGAAAGAGCTGGCGACTGCTAAATTCGTAGAAAGTGTGGACGTAGCTGTTAACCTCGGCATCGACGCTCGTAAATCTGACCAGAACGTACGTGGTGCAACTGTACTGCCGCACGGTACTgg                                                                                                                                                                 <  1:93300/236‑1 (MQ=255)
                                                                                                  aTGGGCCTGGTAGTGGAGGACTAAGAAA‑TGGTTAAACTGACCAAGCGCATGCGTGTT‑ATCCGCGAGAAAGTTGATGCAACCAAACAGTACGACATCA‑ACGAAGCTATCGCACTGCTGAAAGAGCTGGCGACTGCTAAATTCGTAGAAAGTGTGGACGTAGCTGTTAACCTCGGCATCGACGCTCGTAAATCTGACCAGAACGTACGTGGTGCAACTGTACTGCCGCACGGTACTgg                                                                                                                                                                 >  2:93300/1‑236 (MQ=255)
                                                                                                   ggggCCTGGTAGTGGAGGACTAAGAAA‑TGGTTAAACTGACCAAGCGCATGCGTGTT‑ATCCGCGAGAAAGTTGATGCAACCAAACAGTACGACATCA‑ACGAAGCTATCGCACTGCTGAAAGAGCTGGCGACTGCTAAATTCGTAGAAAGTGTGGACGTAGCTGTTAACCTCGGCATCGACGCTCGTAAATCTGACCAGAACGTACGTGGTGCAACTGTACTGCCGCACGGTACTGGCCGTTCCGTTCGCGTa                                                                                                                                                 <  1:317800/250‑1 (MQ=255)
                                                                                                    gggCCTGGTAGTGGAGGACTAAGAAA‑TGGTTAAACTGACCAAGCGCATGCGTGTT‑ATCCGCGAGAAAGTTGATGCAACCAAACAGTACGACATCA‑ACGAAGCTATCGCACTGCTGAAAGAGCTGGCGACTGCCAAATTCGTAGAAAGTGTGGACGTAGCTGTTAACCTCGGCATCGACGCTCGTAAATCTGACCAGAACGTACGTGGTGCAACTGTACTGCCGCACGGTACTGGCCGTTCCGTTCGCGTa                                                                                                                                                 <  2:99531/250‑1 (MQ=255)
                                                                                                             aGTGGAGGACTAAGAAA‑TGGTTAAACTGACCAAGCGCATGCGTGTT‑ATCCGCGAGAAAGTTGATGCAACCAAACAGTACGACATCA‑ACGAAGCTATCGCACTGCTGAAAGAGCTGGCGACTGCTAAATTCGTAGAAAGTGTGGACGTAGCTGTTAACCTCGGCATCGACGCTCGTAAATCTGACCAGAACGTACGTGGTGCAACTGTACTGCCGCACGGTACTGGCCGTTCCGTTCGCGTAGCCGTATTTa                                                                                                                                       <  1:339217/251‑1 (MQ=255)
                                                                                                             aGTGGAGGACTAAGAAA‑GGGTTAAACGGACCAAGCGCATGCGTGTT‑ATCCGCGAGAAAGTTGATGCAACCAAACAGTACGACATCA‑ACGAAGCTATCGCACTGCTGAAAGAGCTGGCGACTGCTAAATCCGTAGAAAGTGTGGACGTAGCTGTTAACCTCGGCATCGACGCTCGTAAATCTGACCAGAACGTACGTGGTGCAACTGTACTGCCGCACGGTACTGGCCGTTCCGTTCGCGTAGCCGTATTTa                                                                                                                                       <  1:217084/251‑1 (MQ=255)
                                                                                                                  aggaCTAAGAAA‑TGGTTAAACTGACCAAGCGAATGCGTGTT‑ATCCGCGAGAAAGTTGATGCAACCAAACAGTACGACATCA‑ACGAAGCTATCGCACTGCTGAAAGAGCTGGCGACTGCTAAATTCGTAGAACGTGTGGACGTAGCTGTTAACCTCGGCATCGACGCTCGTAATTCTGACCAGAACGTACGTGGTGCAACTGTACTGCCTCACGGTACTGGCCGTTCCGTTCGCTTAGCCGTATTTACCCa                                                                                                                                   >  2:329433/1‑250 (MQ=255)
                                                                                                                   ggaCTAAGAAA‑TGGTTAAACTGACCAGGCGCATGCGTTTT‑ATCCGCGAGAAAGTTGATGCAACCAAACAGTACGACATCA‑ACGAAGCTATCGCCCTGCTGAAAGAGCTGGCGACTGCTAAATTCGTAGAAAGTGTGGACGTAGCTGTTAACCGCGGCATCGTCGCTCTTAAATCTGACAAGAGCGTATGTGGTGCAACTGTACGGCCGCACGGGAGTGGCCGTTCCGTTCGCGTTGCCTTATTTACCCAAg                                                                                                                                 <  2:229313/251‑1 (MQ=255)
                                                                                                                      cTAAGAAA‑TGGTTAAACTGACCAAGCGCATGCGTGTT‑ATCCGCGAGAAAGTTGATGCAAACAAACAGTACGACATCACAC‑AAGCTATAGCACTGCTGAAAGAGCTGGCGACTGATAAATTCGTAGAAAGCGTGGACGTAGCTGTTAACCTCGGCATCGACGCTCGTAAATCTGACCAGCACGAACGTGGTGCAACTGAACTGCCGCACGGCACTGGCCGTTCCGTTCGCGTAGCCGTACTTACCCAAGGTg                                                                                                                              >  1:55629/1‑251 (MQ=255)
                                                                                                                        aaGAAA‑TGGTTAAACTGACCAAGCGCATGCGTGTT‑ATCCGCGAGCAAGTTGATGCAACCAAACAGTACGAAATCA‑ACGAAGCTATCGCACTGCTGAAAGAGCTGGCGACTGCTAAATTCGTAGAAAGTGTGGACGTAGCTGTTAACCTCGGCATCGACGCTCGTAACTCTGACCAGAACGTACGTGGTGCAACTGTACTGCCGCACGGCACTGGCCGTTCCGTTCGCGTAGCCGTATTTACCCAAGGTGc                                                                                                                             >  1:47924/1‑250 (MQ=255)
                                                                                                                        aaGAAA‑TGGTTAAACTGACCAAGCGCATGCGTGTT‑ATCCGCGAGAAAGTTGATGCAACCAAACAGTACGACATCA‑ACGAAGCTATCGCACTGCTGAAAGAGCTGGCGACTGCTAAATTCGTAGAAAGTGTGGACGTAGCTGTTAACCTCGGCATCGACGCTCGTAAATCTGACCAGAACGTACGTGGTGCAACTGTACTGCCGCACGGTACTGGCCGTTCCGTTCGCGTAGCCGTATTTACCCACGGTGc                                                                                                                             >  1:162883/1‑250 (MQ=255)
                                                                                                                           aaa‑TGGTTAAACTGACCAAGCGCATGCGTGTT‑ATCCGCGAGAAAGTTGATGCAACCAAACAGTACGACATCA‑ACGAAGCTATCGCACTGCTGAAAGAGCTGGCGACTGCTAAATTCGTAGAAAGTGTGGACGTAGCTGTTAACCTCGGCATCGACGCTCg                                                                                                                                                                                                                    <  1:269615/160‑1 (MQ=255)
                                                                                                                           aaa‑TGGTTAAACTGACCAAGCGCATGCGTGTT‑ATCCGCGAGAAAGTTGATGCAACCAAACAGTACGACATCA‑ACGAAGCTATCGCACTGCTGAAAGAGCTGGCGACTGCTAAATTCGTAGAAAGTGTGGACGTAGCTGTTAACCTCGGCATCGACGCTCg                                                                                                                                                                                                                    >  2:269615/1‑160 (MQ=255)
                                                                                                                           aaa‑TGGTTAAACTGACCAAGCGCATGCGTGTT‑ATCCGCGAGAAAGTTGATGCAACCAAACAGTACGACATCA‑ACGAAGCTATCGCACTGCTGAAAGAGCTGGCGACTGCTAAATTCGTAGAAAGTGTGGACGTAGCTGTTAACCTCGGCATCGACGCTCGTAAATCTGACCAGAACGTACGTGGTGCAACTGTACTGCCGCACGGTACTGGCCGTTCCGTTCGCGTAGCCGTATTTACCCAAGGTGCAAAc                                                                                                                         >  1:193544/1‑251 (MQ=255)
                                                                                                                             a‑TGGTTAAACTGACCAAGCGCATGCGTGTT‑ATCCGCGAGAAAGTTGATGCAACCAAACAGTACGACATCA‑ACGAAGCTATCGCACTGCTGAAAGAGCTGGCGACTGCTAAATTCGTAGAAAGTGTGGACGTAGCTGTTAACCTCGGCATCGACGCTCGTAAATCTGACCAGAACGTACGTGGTGCAACTGTACTGCCGCACGGTACTGGCCGTTCCGTTCGCGTAGCCGTATTTACCCAAGGTGCAAACGc                                                                                                                       <  1:302936/251‑1 (MQ=255)
                                                                                                                                   tAAACTGACCAAGCGCATGCGTGTT‑ATCCGCGAGAAAGTTGATGCAACCAAACAGTACGACATCA‑ACGAAGCTATCGCACTGCTGAAAGAGCTGGCGACTGCTAAATTCGTAGAAAGTGTGGACGTAGCTGTTAACCTCGGCATCGACGCTCGTAAATCTGACCAGAACGTACGTGGTGCAACTGTACTGCCGCACGGTACTGGCCGTTCCGTTCGCGTAGCCGTATTTACCCAAGGTGCAAACGCTGa                                                                                                                    >  1:285693/1‑249 (MQ=255)
                                                                                                                                         gACCAAGCGCATGCGTGTT‑ATCCGCGAGAAAGTTGATGCAACCAAACAGTACGACATCA‑ACGAAGCTATCGCACTGCTGAAAGAGCTGGCGACTGCTAAATTCGTAGAAAGTGTGGACGTAGCTGTTAACCTCGGCATCGACGCTCGTAAATCTGACCAGAACGTACGTGGTGCAACTGTACTGCCGCACGGTACTGGCCGTTCCGTTCGCGTAGCCGTATTTACCCAAGGTGCAAACGCTGAAGCTGCTa                                                                                                            <  2:235106/251‑1 (MQ=255)
                                                                                                                                             aaGCGCATGC‑TGTTTATCCGCGAGAAAGTTGATGCAACCAACCAGTACGACATCA‑ACGAAGCTATCGCACTGCTGAAAGAGCTGGCGACTGCTAAATTCGTAGAAAGTGTGAACGTATCTGTTAACCTCGGCATCGACGCTCGTAAATCTGACCAGAACGTACGTGGTGCAACTGTACTGCCGCACGGTACTGGCCGTTCCGTTCGCGTAGCCGTATTTACCCAAGGTGCAAACGCTGAAGCTGCTaaatc                                                                                                        <  2:212985/251‑3 (MQ=255)
                                                                                                                                                   aTGCGTGTT‑ATCCGCGAGAAAGTTGATGCAACCAAACAGTACGACATCA‑ACGAAGCTATCGCACTGCTGAAAGAGCTGGCGACTGCTACATTCGTAGAAAGTGTGGACGTAGCTGTTAACCTCGGCATCGACGCTCGTAACTCTGACCAGAACGTACGTGGTGCCACTGTACTACCGCACGGTACTGGACGTTCCGTTCGCGTAGCCGTATTTACCCAAGGTGCAAACGCTGAAGCTGCCAAAGCTGCAgg                                                                                                  >  1:392095/1‑251 (MQ=255)
                                                                                                                                                         gtT‑ATCCGCGAGAAAGTTGATGCAACCAAACAGTACGACATCA‑ACGAAGCTATCGCACTGCTGAAAGAGCTGGCGACTGCTAAATTCGTAGAAAGTGTGGACGTAGCTGTTAACCTCGGCATCGACGCTCGTAAATCTGACCAGAACGTACGTGGTGCAACTGTACTGCCGCACGGTACTGGCCGTTCCGTTCGCGTAGCCGTATTTACCCAAGGTGCAACCGCTGAAGCTGCTAAAGCTGCAGgcgc                                                                                               >  2:306320/1‑248 (MQ=255)
                                                                                                                                                           t‑ATCCGCGAGAAAGTTGATGCAACCAAACAGTACGACATCA‑ACGAAGCTACCGCACTGCTGAAAGAGCGGGCGACTGCTAAATTCGTAGAAAGAGTGGACGTGGCTGTTAACCTCGGCAGCGACGCTCGTAAATCTGACCAGAACGGACGTGGTTCAACTGGACTGCCGCACGGTACTGGCCGTTCCGTGCGCGTAGCCGGATTTACCCAAGGTGCAAACGCTGAAGCTTCTAAAGCTGCAGGCTCAGaa                                                                                           <  2:391219/250‑1 (MQ=255)
                                                                                                                                                                   gataaaGTTGATGAAACCAAACAGTACGACATCA‑ACGAAGCTATCGCACTGCTGAAAGAGCTGGCGACTGCTAAATTCGTAGAAAGTGTGGACGTAGCTGTTAACCTCGGCATCGACGCTCGTAAATCTGACCATAACGTACGTGGTGCAACTGTACTGCCGCACGGTACTGGCCGTTCCGTTCGCGTAGCCGTATTTACCCAAGGTGCAAACGCTGAAGCTGCTAAAGCGGCAGGCGCAGAACTggtagg                                                                                   <  2:92777/248‑1 (MQ=255)
                                                                                                                                                                                 aaCCAAACAGTACGACATCA‑ACGAAGCTATCGCACTGCTGAAAGAGCTGGCGACTGCTAAATTCGTAGAAAGTGTGGACGTAGCTGTTAACCTCGGCATCGACGCTCGTAAATCTGACCAGAACGTACGTGGTGCAACTGTACTGCCGCACGGTACTGGCCGTTCCGTTCGCGTAGCCGTATTTCCCCCAGGTGCAAACGCTGAAGCTGCTAACGCTGCAGGCGCAGGACTGGCAGGTATGGAAGAt                                                                         >  2:88790/1‑247 (MQ=255)
                                                                                                                                                                                 aaCCAAACAGTACGACATCA‑ACGAAGCTATCGCACTGCTGAAAGAGCTGGCGACTGCTAAATTCGTAGAAAGTGTGGACGTAGCTGTTAACCTCGGCATCGACGCTCGTAAATCTGACCAGAACGTACGTGGTGCAACTGTACTGCCGCACGGTACTGGCCGTTCCGTTCGCGTAGCCGTATTTACCCAAGGTGCAAACGCTGAAGCTGCTAAAGCTGCAGGCGCAGAACTGGTAGGTATGGAAGAt                                                                         <  1:88790/247‑1 (MQ=255)
                                                                                                                                                                                 aaCCAAACAGTACGACATCA‑ACGAAGCTATCGCACTGCTGAAAGAGCTGGCGACTGCTAAATTCGTAGAAAGTGTGGACGTAGCTGTTAACCTCGGCATCGACGCTCGTAAATCTGACCAGAACGTACGTGGTGCAACTGTACTGCCGCACGGTACTGGCCGTTCCGTTCGCGTAGCCGTATTTACCCAAGGTGCAAACGCTGAAGCTGCTAAAGCTGCAGGCGCAGAACTGGTAGGTATGGAAGATCTgg                                                                     >  1:262113/1‑251 (MQ=255)
                                                                                                                                                                                           tACGACATCA‑ACGAAGCTATCGCACTGCTGAAAGAGCTGGCGACTGCTAAATTCGTAGAAAGTGTGGACGTAGCTGTTAACCTCGGCATCGACGCTCGTAAATCTGACCAGAACGTACGTGGTGCAACTGTACTGCCGCACGGTACTGGCCGTTCCGTTCGCGTAGCCGTATTTACCCAAGGTGCAAACGCTGAAGCTGCTAAAGCTGCCGGCGCAGAACTGGTAGGTATGGAAGATCTGGCTGAcc                                                               >  2:155165/1‑247 (MQ=255)
                                                                                                                                                                                             cGACATCA‑ACGAAGCTATCGCACTGCTGAAAGAGCTGGCGACTGCTAAAGTCGTAGAAAGTGTGGACGTAGCTGTTAACCTCGGCATCGACGCTCGTAAATCTGACCAGAACGTACGTGGTGCAACTGTACTGCCGCACGGTACTGGCCGTTCCGTTCGCGTAGCCGTATTTACCCAAGGTGCAAACGCTGAAGCTGCTAAAGCTGCAGGCGCAGAACTGGTAGGTATGGAAGATCTGGCTGACCAGATCa                                                         <  2:301906/251‑1 (MQ=255)
                                                                                                                                                                                                   cA‑ACGAAGCTATCGCACTGCTGAAAGAGCTGGCGACTGCTAAATTCGTAGAAAGTGTGGACGTAGCTGTTAACCTCGGCATCGACGCTCGTAAATCTGACCAGAGCGTACGTGGTGCAACTGTACTGCCGCACGGTACTGGCCGTTCCGTTCGCGTAGCCGTATTTACCCAAGGTGCAAACGCTGAAGCTGCTAAAGCTGCAGGCGCAGAACTGGTAGGTATGGAAGATCTGGCTGACCAGATCAAGAAAg                                                   <  1:73210/251‑1 (MQ=255)
                                                                                                                                                                                                    a‑aCGAAGCTATCGCACTGCTGAAAGATCTGGCGACTGCTAAATTCGTAGAAAGTGTGGACGTAGCTGTTAACCTCTGCATCGACGCTCGTAAATCTGACCAGAACGGACGGGGTGCAAGTGTACTGCCGCACGGTACTGGCCGTTCCGTTCGCGTAGCCGTATTTACCCAAGGGGCAAACGCTGAAGCTGCTAAAGCTGCAGGCGCAGAACTGGTAGGTATGGAAGATCTGGCTGACCAGATCAAGAAAgg                                                  <  1:391480/251‑1 (MQ=255)
                                                                                                                                                                                                        gAAGCTATCGCACTGCTGAAAGAGCTGGCGACTGCTAAATTCGTAGAAAGTGTGGACGTAGCTGTTAACCTCGGCATCGACGCTCGTAAATCTGACCAGAACGTACGTGGTGCAACTGTACTGCCGCACGGTACTGGCCGTTCCGTTCGCGTAGCCGTATTTACCCAAGGTGCAAACGCTGAAGCTGCTAAAGCTGCAGGCGCAGAACTGGTAGGTATGGAAGATCTGGCTGACCAGATCAAGAAAGGCg                                                >  2:262772/1‑250 (MQ=255)
                                                                                                                                                                                                         aaGCTATCGCACTGCTGAAAGAGCTGGCGACTGCTAAATTCGTAGAAAGTGTGGACGTAGCTGTTAACCTCGGCATCGACGCTCGTAAATCTGACCAGAACGTACGTGGTGCAACTGTACTGCCGCACGGTACTGGCCGTTCCGTTCGCGTAGCCGTATTTACCCAAGGTGCAAACGCTGAAGCTGCTAAAGCTGCAGGCGCAGAACTGGTAGGTATGGAAGATCTGGCTGACCAGATCAAGAAAGGCGaa                                              <  1:296845/251‑1 (MQ=255)
                                                                                                                                                                                                            cTATCGCACTGCTGAAAGAGCTGGCGACTGCTAAATTCGTAGAAAGTGTGGACGTAGCTGTTAACCTCGGCATCGACGCTCGTAAATCTGACCAGAACGTACGTGGTGCAACTGTACTGCCGCACGGTACTGGCCGTTCCGTTCGCGTAGCCGTATTTACCCAAGGTGCAAACGCTGAAGCTGCTAAAGCTGCAGGCGCAGAACTGGTAGGTATGGAAGATCTGGCTGACCAGATCAAGAAAGTCGAAATg                                           >  1:173134/1‑251 (MQ=255)
                                                                                                                                                                                                                  cACTGCTGAAAGAGCTTGCGACTGCTAAATCCGTAGAAAGTGTGGACGTAGCTGTTACCCGCGGCATCGACGCTCGTCAATCTGACCAGAACGTACGTGGTGCAACTGTACTGCCGCTCGGTACTGGCCGTTCCGTGCGCGTAGCCGTATTTACCCAAGGTGCAAACGCTGAAGCTGCTAAAGCTGCAGACGCAGAACTGGTAGGTATGGAAGATCTGGCTGACCAGATCAAGAAATGCGAAACGAACttt                                     <  2:173552/251‑1 (MQ=255)
                                                                                                                                                                                                                   aCTGCTGAAAGAGCTGGCGACTGCTAAATTCGTAGAAAGTGTGGACGTAGCTGTTAACCTCGGCATCGACGCTCGTAAATCTGACCAGAACGTACGTGGTGCAACTGTACTGCCGCACGGTACTGGCCGTTCCGTTCGCGTAGCCGTATTTACCCAAGGTGCAAACGCTGAAGCTGCTAAAGCTGCAGGCGCAGAACTGGTAGGTATGGAAGATCTGGCTGACCAGATCAAGAAAGGCGAAATGAACTTTg                                    >  2:358815/1‑251 (MQ=255)
                                                                                                                                                                                                                            agagCTGGCGACTGCTAAATTCGTAGAAAGTGTGGACGTAGCTGTTAACCTCGGCATCGACGCTCGTAAATCTGACCAGAACGTACGTGGTGCAACTGTACTGCCGCACGGTACTGGCCGTTCCGTTCGCGTAGCCGTATTTACCCAAGGTGCAAACGCTGAAGCTGCTAAAGCTGCAGGCGCAGAACTGGTAGGTATGGAAGATCTGGCTGACCAGATCAAGAAAGGCGAAATGAACTTTGACGTTGTTa                           <  2:86927/251‑1 (MQ=255)
                                                                                                                                                                                                                             gagCTGGCGACTGCTAAATTCGTAGAAAGTGTGGACGTAGCTGTTAACCTCGGCATCGACGCTCGTAAATCTGACCAGAACGTACGTGGTGCAACTGTACTGCCGCACGGTACTGGCCGTTCCGTTCGCGTAGCCGGATTTACCCAAGGTGCAAACGCTGAAGCTGTTAAAGCTGCAGGCGCAGAACTGGTAGGTATGGAAGATCTGGCTGACCAGATCAAGAAAGGCGAAATGAACTTTGACGTTGTTAt                          <  1:252320/251‑1 (MQ=255)
                                                                                                                                                                                                                               gCTGGCGACTGCTAAATTCGTAGAAAGTGTGGACGTAGCTGTTAACCTCGGCATCGACGCGCGTAAATCTGACCAGAACGTACGTGGTGCAACTGTTCTGCCGCACGGTACTGGCCGTTCCGTTCGCGTAGCCGTATTTACCCAAGGTGCAAACGCTGAAGCTGCTAAAGCTGCAGGCGCAGAACAGGTAGGTATGGAAGATCTGGCTGACCAGATCAAGAAAGGCGAAATGAACTTTGACGTTGTTATTg                        <  2:162883/251‑1 (MQ=255)
                                                                                                                                                                                                                                  ggCGACTGCTAAATTCGTAGAAAGTGTGGACGTAGCTGTTAACCTCGGCATCGACGCTCGTAAATCTGACCAGAACGTACGTGGTGCAACTGTACTGCCGCACGGTACTGGCCGTTCCGTTCGCGTAGCCGTATTTACCCAAGGTGCAAACGCTGAAGCTGCTAAAGCTGCAGGCGCAGAACTGGTAGGTATGGAAGATCTGGCTGACCAGATCAAGAAAGGCGAAATGAACTTTGCCGTTGTTATTGCtt                     >  1:328238/1‑251 (MQ=255)
                                                                                                                                                                                                                                  ggCGACTGCTAAATTCGTAGAAAGTGTGGACGTAGCTGTTAACCTCGGCATCGACGCTCGTAAATCTGACCAGAACGTACGTGGTGCAACTGTACTGCCGCACGGTACTGGCCGTTCCGTTCGCGTAGCCGTATTTACCCAAGGTGCAAACGCTGAAGCTGCTAAAGCTGCAGGCGCAGAACTGGTAGGTATGGAAGATCTGGCTGACCAGATCAAGAAAGGCGAAATGAACTTTGACGTTGTTATTGCt                      <  2:338487/250‑1 (MQ=255)
                                                                                                                                                                                                                                     gACTGCTAAATTCGTAGAAAGTGTGGACGTAGCTGTTAACCTCGGCATCGACGCTCGTAAATCTGACCAGAACGTACGTGGTGCAACTGTACTGCCGCACGGTACTGGCCGTTCCGTTCGCGTAGCCGTATTTACCCAAGGTGCAAACGCTGAAGCTGCTAAAGCTGCAGGCGCAGAACTGGTAGGTATGGAAGATCTGGCTGACCAGATCAAGAAAGGCGAAATGAACTTTGACGTTGTTATTGCTtctc                  <  1:262772/251‑1 (MQ=255)
                                                                                                                                                                                                                                        tGCTAAATTCGAAGAAAGTGTGGACGTAGCTGTTAACCTCGGCATCGACGCTCGTAAATCTGACCAGAACGTACGTGGTGCAACTGTACTGCCGCACGGTACTGGCCGTTCCGTTCGCGTAGCCGTATTTACCCAAGGTGCAAACGCTGAAGCTGCTAAAGCTGCAGGCGCAGACCTGGTAGGTATGGAAGATCTCGCTGACCCGATCAAGCAAGGCGAAATGAACTTTGACGTTGTTATTGCTTCTCCgg               >  2:83559/1‑251 (MQ=255)
                                                                                                                                                                                                                                         gCTAAATTCGTAGAAAGTGTGGACGTAGCTGTTAACCTCGGCATCGACGCTCGTAAATCTGACCAGAACGTACGTGGTGCAACTGTACTGCCGCACGGTACTGGCCGTTCCGTTCGCGTAGCCGTATTTACCCAAGGTGCAAACGCTGAAGCTGCTAAAGCTGCAGGCGCAGAACTGGTAGGTATGGAAGATCTGGCTGACCAGATCAAGAAAGGCGAAATGAACTTTGACGTTGTTATTGCTTCTCCGGa              <  2:88253/251‑1 (MQ=255)
                                                                                                                                                                                                                                                  gTAGAAAGAGTGGCCGTAGCTGTTAACCTCGGCATCGACGCTCGTAAATCTGCCCAGAACGTACGTGGTGCAACTGTACTGCCGCACGGTACTGGCCGTTCCGTTCGCGTAGCCGTATTTACCCAAGGTGCAAACGCTGAAGCTGCTAAAGCTGCAGGCGCAGAACTGGTAGGTATGGAAGATCTGGCTGACCAGATCAAGAAAGGCGAAATGAACTTTGACGATGTTATTGCTTCTCCGGATGCAATgcg     >  2:257020/1‑251 (MQ=255)
                                                                                                                                                                                                                                                   tAGAAAGTGTGGACGTAGCTGTTAACCTCGGCATCGACGCTCGTAAATCTGACCAGAACGTACGTGGTGCAACTGTACTGCCGCACGGTACTGGCCGTTCCGTTCGCGTAGCCGTATTTACCCAAGGTGCAAACGCTGAAGCTGCTAAAGCTGCAGGCGCAGAACTGGTAGGTATGGAAGATCTGGCTGACCAGATCAAGAAAGGCGAAATGAACTTTGACGTTGTTATTGCTTCTCCGGATGCAATgcgc    >  1:31622/1‑251 (MQ=255)
                                                                                                                                                                                                                                                    aGAAAGTGTGGACGTAGCTGTTAACCTCGGCAGCGACGCTCGTAAATCTGACCAGAACGTACGTGGTGCAACTGTACTGCCGCACGGTACTGGCCGTTCCGTTCGCGTAGCCGTATTTACCCAAGGTGCAAACGCTGAAGCTGCTAAAGCTGCAGGCGCAGAACTGGTAGGTATGGAAGATCTGGCTGACCAGATCAAGAAAGGCGAAATGAACTTTGACGTTGTTATTGCTTCTCCGGATGCAATGcgcg   <  2:286525/251‑1 (MQ=255)
                                                                                                                                                                                                                                                     gAAAGTGTGGACGTAGCTGTTAACCTCGGCATCGACACTCGTAAATCTGACCAGAACGTACGTGGTGCAACTGTACTGCCGCACGGTACTGGCCGTTCCGTTCGCGTAGCCGTATTTACCCAAGGTGCAAACGCTGACGCTGCTAAAGCTGCAGGCGCAGAACTGGTAGGTATGGAAGATCTGGCTGACCAGATCAAGAAAGGCGAAATGAACTTTGACGTTGTTATTGCTTCTCCGGATGCAATGCGCgt  >  2:292524/1‑251 (MQ=255)
                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                     
GCGCTCAGCTGCAGGAAATCGCGCAGACCAAAGCTGCCGACATGACTGGTGCCGACATTGAAGCGATGACTCGCTCCATCGAAGGTACTGCACGTTCCATGGGCCTGGTAGTGGAGGACTAAGAAA‑TGGCTAAACTGACCAAGCGCATGCGTGTT‑ATCCGCGAGAAAGTTGATGCAACCAAACAGTACGACATCA‑ACGAAGCTATCGCACTGCTGAAAGAGCTGGCGACTGCTAAATTCGTAGAAAGCGTGGACGTAGCTGTTAACCTCGGCATCGACGCTCGTAAATCTGACCAGAACGTACGTGGTGCAACTGTACTGCCGCACGGTACTGGCCGTTCCGTTCGCGTAGCCGTATTTACCCAAGGTGCAAACGCTGAAGCTGCTAAAGCTGCAGGCGCAGAACTGGTAGGTATGGAAGATCTGGCTGACCAGATCAAGAAAGGCGAAATGAACTTTGACGTTGTTATTGCTTCTCCGGATGCAATGCGCGT  >  NC_000913/4178754‑4179246

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