Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,210,859 T→G 19.2% G27G (GGT→GGG)  rpsU → 30S ribosomal subunit protein S21

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,210,8590TG19.2% 128.2 / 2.6 74G27G (GGT→GGG) rpsU30S ribosomal subunit protein S21
Reads supporting (aligned to +/- strand):  ref base T (26/26);  new base G (16/3);  total (45/29)
Fisher's exact test for biased strand distribution p-value = 1.33e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTCACCCTTATAAAAGTCCCTTTCAAAAAAGGCCGCGGTGCTTTACAAAGCAGCAGC‑AATTGCAGTAAAATTCCGCACCATTTTGAAATAAGCTGGCGTTGATGCCAGCGGCAAACCGAATTAATCAAAGGTGAGAGGCACATGCCGGTAATTAAAGTACGTGAAAACGAGCCGTTCGACGTAGCTCTGCGTCGCTTCAAGCGTTCCTGCGAAAAAGCAGGTGTTCT‑GGCGGAAGTTCGTCGTCGTGAG‑TTCTATGAAAAACCGACTACCGAACGTAAGCGCGCTAAAGCTTCTGCAGTGAAACGTCACGCGAAGAAACTGGCTCG‑CGAAAACGCACGCCGC‑ACTCGTCTGTACTAATTCCCCGAGAGCGTTGCTCTCCGATCAGACCGAGTTGTAGTTGTAAGGCCGTGCTTCCGAAAGGAATGCGCGGCTTATTTTCGTTTATGAA  >  NC_000913/3210638‑3211085
                                                                                                                                                                                                                              |                                                                                                                                                                                                                                      
ctcACCCTTATAAAAGTCCCTTTCAAAAAAGGCCGCGGTGCTTTACAAAGCAGCAGC‑AATTGCAGTAAAATTCCGCACCATTTTGAAATAAGCTGGCGTTGATGCCAGCGACAAACCGAATTAATCAAAGGTGAGAGGCACATGCCGGTAATTAAAGTACGTGAAAACAATCCGTTAGACGTAACTCTGAGTCGCTTCAAGCGTTCCTGCGAACAAACAGGGGTTCT‑g                                                                                                                                                                                                                                 >  1:61093/1‑228 (MQ=255)
 tcACTCTTATAAAAGTCCCTTTCAAAAAAGGCCGCGGTGCTTTACAAAGCAGCAGC‑AATTGCAGTAAAATTCCGCACCATTTTGAAATAAGCTGGCGTTGATGCCAGCGGCAAACCGAATTAATCAAAGGTAAACGGCCCATGCCGGTAATTAAAACACCTGAAACCAAACCGTTCGACCAAGCCCTGCGGTCCTTCAAGCCTGCCCGCGAAAAAACCGGGGtttgg                                                                                                                                                                                                                                  >  1:694766/1‑224 (MQ=255)
  cACCCTTATAAAAGTCCCTTTCAAAAAAGGCCGCGGTGCTTTACAAAGCAGCAGC‑AATTGCAGTAAAATTCCGCACCATTTTGAAATAAGCTGGCGTTGATGCCAGCGGCAAACCGAATTAATCAAAGGTTAGAGTCACATGCCGGTAATTAAAGTACGTGAAAACAAGCCGTTAGACGTAACTCTGCGTAGCTTCAAGCTTTCCGCGGCAAAAGAAGGTGTTTT‑gcg                                                                                                                                                                                                                               >  2:642754/1‑226 (MQ=255)
   aCCCTTATAAAAGTCCCTTTCAAAAAAGGCCGCGGTGCTTTACAAAGCAGCAGC‑AATTGCAGTAAAATTCCGCACCATTTTGAAATAAGCTGGCGTTGCTGCCAGCGGCAAACCGAATTAATCAAAGGTGAGAGGCAAATGCCGGGAATTAAGGTACGAGAAACCGAGCCGTTCGACGCAGCTCGCCGGCGCTTCAGACGTTCCTTCGTAAAACCGGGTGTTTT‑GCCg                                                                                                                                                                                                                              >  2:391539/1‑228 (MQ=255)
   aCCCTTATAAAAGTCCCTTTCAAAAAAGGCCGCGGTGCTTTACAAAGCAGCAGC‑AATTGCAGTAAAATTCCGCACCATTTTGAAATAAGCTGGCGTTGATGCCAGCGGCAAACCGAATTAATCAAAGGTGAGAGGCACATGCCGGTAATTAAAGTACGTGAAAACGAGCCGTTCGACGTAGCTCTGCGTCGCTTCAAGCGTTCCTGCGAAAAAGCAGGTGTTCT‑GGc                                                                                                                                                                                                                               >  2:422301/1‑227 (MQ=255)
   aCCCTTATAAAAGTCCCTTTCAAAAAAGGCCGCGGTGCTTTACAAAGCAGCAGC‑AATTGCAGTAAAATTCCGCACCATTTTGAAATAAGCTGGCGTTGATGCCAGCGGCAAACCGAATTAATCAAAGGTGAGAGGCACATGCCGGTAATTAAAGTACGTGAAAACAAGCCGTTCGACGTAGCTCTGCGTCGCTTCAAGCGTTCCTGCGAAAAAGCAGGGGTTCT‑gggg                                                                                                                                                                                                                              >  1:489255/1‑226 (MQ=255)
        tataAAAGTCCCTTTCAAAAAAGGCCGCGGTGCTTTACAAAGCAGCAGC‑AATTGCAGTAAAATTCCGCACCATTTTGAAATAAGCTGGCGTTGATGCCAGCGGCAAACCGAATTAATCAAAGGTGAGAGGCACATGCCGGTAATTAAAGTACGTGAAAACGAGCCGTTCGACGTAGCTCTGCGTCGCTTCAAGCGTTCCTGCGAAAAAACAGGGGTTAT‑GGCGGCAgg                                                                                                                                                                                                                         >  1:375476/1‑227 (MQ=255)
        tataAAAGTCCCTTTCAAAAAAGGCCGCGGTGCTTTACAAAGCAGCAGC‑AATTGCAGTAAAATTCCGCACCATTTTGAAATAAGCTGGCGTTGATGCCAGCGGCAAACCGAATTAATCAAAGGTGAGAGGCACATGCCGGTAATTAAAGTACGTGAAAACGAGCCGTTCGACGAAGCTCTGCGTCGCTTCAAGCGTTCCTGCGAAAAAGCAGGTGTGCT‑GGCGGAAGt                                                                                                                                                                                                                         >  2:427007/1‑228 (MQ=255)
        tataAAAGTCCCTTTCAAAAAAGGCCGCGGTGCTTTACAAAGCAGCAGC‑AATTGCAGTAAAATTCCGCACCATTTTGAAATAAGCTGGCGTTGATGCCAGCGGCAAACCGAATTAATCAAAGGTGAGAGGCACATGCCGGTAATTAAAGTACGTGAAAACGAACCGTTCGACTTAGCTCTGCGTCCCTTCAACCGTTCCTGCGAAAAACAAGGGGTTCT‑GGCGGAAGt                                                                                                                                                                                                                         >  1:653326/1‑228 (MQ=255)
        tataAAAGTCCCTTTCAAAAAAGGCCGCGGTGCTTTACAAAGCAGCAGC‑AATTGCAGTAAAATTCCGCACCATTTTGAAATAAGCTGGCGTTGATGCCAGCGGCAAACCGAATTAATCAAAGGTGAGAGGCACATGCCGGTAATTAAAGTACGGGAAAACAAGCCGTTCAACGAAGCCCTGGGGCGCTTCAGGCGTTCCTGCCAAAAAGCAGGGGTTCG‑GGGGGAAGt                                                                                                                                                                                                                         >  1:343395/1‑228 (MQ=255)
         ataAAAGTCCCTTTCAAAAAAGGCCGCGGTGCTTTACAAAGCAGCAGC‑AATTGCAGTAAAATTCCGCACCATTTTGAAATAAGCTGGCGTTGATGCCAGCGGCAAACCGAATTAATCAAAGGTAAGAGGCACATGCCGGTAATTAAAGTACGCGAAAACGAGCCGTTCGACGTAGCTCTGCATCGCTTCAAGCGTTTCTGCTGAAAAGAAGGTGTTCT‑GCCGCAAGtt                                                                                                                                                                                                                        >  2:234339/1‑228 (MQ=255)
          taAAAGTCCCTTTCAAAAAAGGCCGCGGTGCTTTACAAAGCAGCAGC‑AATTGCAGTAAAATTCCGCACCATTTTGAAATAAGCTGGCGTTGATGCCAGCGGCAAACCGAATTAATCAAAGGTGAGAGGCCCATGCCGGTAATTAAAGTACGTGAAAACGAGCCGTTCGCCGTAGCTCTGCGTCGCTTCAAGCGTTCATCCGAAAAAGCAGGGGTTCT‑GCCGGAAGTtc                                                                                                                                                                                                                       >  2:263820/1‑228 (MQ=255)
          taAAAGTCCCTTTCAAAAAAGGCCGCGGTGCTTTACAAAGCAGCAGC‑AATTGCAGTAAAATTCCGCACCATTTTGAAATAAGCTGGCGTTGATGCCAGCGGCAAACCGAATTAATCAAAGGTGAGAGGCACATGCCGGTAATTAAAGTACGTGAAAACGAGCCGTTCGACGTAGCTCTGCGTCGCTTCAAGCGTTCCTGCGAAAAAGCAGGGGTTCT‑GGCGGAAGCtc                                                                                                                                                                                                                       >  2:385472/1‑228 (MQ=255)
           aaaaATCCATTTAAAAAAAGGCCGCGGGGTTTTAAAAAGC‑GGAGCAAATGGCAAGAAAAGTCCGGAACAATTGGAAAAAAGCTGGGGTTGATGCAAGCGGAAAACCGAATAAATCAAAGGTGAGAGGCACATGCCGGTAATTAAAGTACGTGAAAACGAGCCGTTCGACGTAGCTCTGCGTCGCTTCAAGCGTTCCTGCGAAAAAGCAGGTGTTCT‑GGCGGAAGTtcg                                                                                                                                                                                                                      <  1:2848/228‑1 (MQ=255)
             aaGTCCCTTACAAAAACGGCCGCGGTGCTTTACAAAGCAGCAGC‑AATTGACGTAAAATTCCGCCCCATTTTGAAATAAGCTGGCGTTGATGCCAGCGGCAAACCGAATTAATCAAAGGTGAGAGGCACATGCCGGTAATTAAAGTACGTGAAAACGAGCCGTTCGACGTAGCTCTGCGTCGCTTCAAGCGTTCCTGCGAAAAAGCAGGTGTTCT‑GGCGGAAGTtcgtc                                                                                                                                                                                                                    <  2:666349/228‑1 (MQ=255)
               gTCCCTTTCAAAAAAGGCCGCGGTGCTTTACAAAGCAGCAGC‑AATTGCAGTAAAATTCCGCACCATTTTGAATTAACTTGGCGTTGATACCAGCGGCAAACCGAATTAATCAAAGGCGAGAGGCACATGCTGGTAATTAAAGTAAGTGATAACAAGCGGTTCGAAGTAGCTCTTCTTCGCTTCAAACGTTCCTGCGAACAATCAAGAGATCT‑GAGGGAAGttgccg                                                                                                                                                                                                                    >  2:99570/1‑222 (MQ=255)
               gTCCCTTTCAAAAAAGGCCGCGGTGCTTTACAAAGCAGCAGC‑AATTGCAGTAAAATTCCGCACCATTTTGAAATAAGCTGGCGTTGATGCCAGCGGCAAACCGAATTAATCAAAGGTGAGAGGCACATGCCGGTAATTAAAGTACGTGAAAACGAGCCGTTCGACGTAGCTCTGCGTCGCTTCAACCGTTCCTGCGAAAAAGCAGTGGTTCT‑GGCGGAAGTtcctcgt                                                                                                                                                                                                                  >  2:358697/1‑228 (MQ=255)
               gTCCCTTTCAAAAAAGGCCGCGGTGCTTTACAAAGCAGCAGC‑AATTGCAGTAAAATTCCGCACCATTTTGAAATAAGCTGGCGTTGATGCCAGCGGCAAACCGAATTAATCAAAGGTGAGAGGCACATGCCGGCAATTAAAGAACGTGAAAACGAGCCGTTCAACGTAGCTCTGCGTCCCTTTAAGCGTTACCGCGAAAAAGCAGGGTTTCT‑GGCGGAAATtcgtcgc                                                                                                                                                                                                                  >  1:662862/1‑227 (MQ=255)
               gTCCCTTTCAAAAAAGGCCGCGGTGCTTTACAAAGCAGCAGC‑AATTGCAGTAAAATTCCGCACCATTTTGAAATAAGCTGGCGTTGATGCCAGCGGCAAACCGAATTAATCAAAGGTAAGAGGCACATGCCGGCAATTAAAGTACGTGAAAACGAGCCGTTCGACGTAGCTCTGCGTCGCTTCAAGCGTTCCTGCGAAAAAGACGGTGTTCT‑GGCGGAAGGttgtctt                                                                                                                                                                                                                  >  1:496515/1‑226 (MQ=255)
                  tctgtCAAAAAAGGCCGCGGTGCTTTACAAAGCGGCAGC‑AATGGCAGTAAAATACCGAAACATTTTGAAATAAGCTGGCGTTGATTCCAGCGGCAAACCGAATTAGTCAAGGGTGAGGGGCACATGCCGGTAATTAAAGTACGTGAAAACGAGCCGTTCGACGTAGCTCTGCGTCGCTTCAAGCGTTCCTGCGAAAAAGCAGGTGTTCT‑GGCGGAAGTTcgtcgtcgt                                                                                                                                                                                                               <  1:639113/224‑1 (MQ=255)
                  ccTTTCAAAAAAGGCCGCGGTGCTTTACAAAGCAGCAGC‑AATTGCAGTAAAATTCCGCACCATTTTGAAATAAGCTGGCGTTGATGCCAGCGGCAAACCGAATTAATCAAAGGTTAGAGGCCCATGCCGGTAATTAAAGTACGGTAAAACGAGCCGGTCAACGTAACTCTACGTCGCTTCAATCGTTCCTCCGCAAAAGCAGGGGTCTA‑GGCGGAAGttcttctttgt                                                                                                                                                                                                               >  1:412608/1‑223 (MQ=255)
                    tttCAAAAAAGGCCGCGGTGCTTTACAAAGCAGCAGC‑AATTGCAGTAAAATTCCGCACCATTTTGAAATAAGCTGGCGTTGATGCCAGCGGCAAACCGAATTAATCAAAGGTGAGAGGCACATGCCGGTAATTAAAGTACGTGAAAACGACCCGTTCGACGTAGCTCTGCGTCGCTTCAAGCGTTCCTGCGAAAAAGCAGGGTTTCT‑GCCGGAAGTtcgtcgcccggg                                                                                                                                                                                                             >  1:147943/1‑222 (MQ=255)
                     ttaaaaaatGGCCCAGGTGCATTACAAGGAAGCAGC‑AAGAGCGGTAAAATCCCACACTTTTTAGAAATAAGCGAGCGTGGATGCGACCGGCAAACCGAATAAGTCAAAGGTGGAAGGGACATGCCGGAAATTAAAGTACGTGAAAACAAGCCGTACGACGTAGCTCTGGGTCGCTTCAATCGTTCCTGCGAAAAAGAAGGTGTTCT‑GGCGGAAGTTCGTCGTCGTGAg                                                                                                                                                                                                            <  1:99570/225‑1 (MQ=255)
                     ttCAAAAAAGGCCGCGGTGCTTTACAAAGCAGCAGC‑AATTGCAGTAAAATTCCGCACCATTTTGAAATAAGCTGGCGTTGATGCCAGCGGCAAACCGAATTAATCAAAGGAGAGAGGCACATGCCGGTAATTAAAGTACGTGAAAACTAGCCGTTCGACGTAGCTCTGCGTCGCTTCAAGCGTTACTGCGCAAAAGCAAGGGGTCT‑GGCGGAAGTTCGTCGTCGTggg                                                                                                                                                                                                            >  2:672154/1‑226 (MQ=255)
                                ccGGGGTGGTTTACAAAGCAGCAGC‑AAGTGCAGAAAAATTCCGCCCAATTTTGAAATAAGCTGGCGGTTATGCCAGCGGCAAACCAAATTAATCAAAGGTGAGAGGCACATGCCGGTAATTAAAGTACGTGAAACCGAGCCGTTCGACGTAGCTCTGCGTCGCTTCAAGCGTTCCTGCGAAAAAGCAGGTGTTCT‑GGCGGAAGTTCGTCGTCGTGAG‑TTCTATGaaaa                                                                                                                                                                                                <  2:303338/228‑1 (MQ=255)
                                   cgGTGCTTTACAAAGCAGCAGC‑AATTGCAGTAAAATTCCGCACCATTTTGAAATAAGCTGGCGTTGCTCCCAGCGACAAACCGAATTAATCACAGGTGAGAGGCACATGCCGGTAATTAAAGTACGTGAAAACGAACCGTTCGCCGGAGCGCTGCGTCGATTCACGAGTTCCTCCGAAAAATCCGGTGTGCT‑GGCGGAGGATCGTGGT‑GGGAGTTTCGATGAAAAAcc                                                                                                                                                                                             >  1:419631/1‑228 (MQ=255)
                                        ttttAAACAGCAGCAGC‑AATTGCAGTAAAAACCCGCAACATTTAGAAATAAGCGGGGTTTTATGCGAAGGACAAACGCAATAAATCAAAGGGGGGAGGCTACTGCCGGTAATTAAAGTACGTGAAAACGAGCCGTTCGACGTAGCTCTGCGTCGCTTCAAGCGTTCCTGCGAAAAAGCAGGTGTTCT‑GGCGGAAGTTCGTCGTCGTGAG‑TTCTATGAAAAACCGACTa                                                                                                                                                                                        <  1:480795/227‑1 (MQ=255)
                                        cTTTACAAAGCAGCAGC‑AATTGCAGTAAAATTCCGCACCATTTTGAAATAAGCTGGCGTTGATGCCAGCGGCAAACCGAATTAATCAAAGGTGAGAGGCACATGCCGGTAATTAAAGTACGTGAAAACGAGCCGTTCGACGTAGCTCTGCGTCGCTTCAAGCGTTCCTGCGAAAAAGCAGGGGTTCT‑GGCGGAAGTTCGTCGTCGTGAG‑TTCTATGAAAAAccactta                                                                                                                                                                                        >  2:243325/1‑223 (MQ=255)
                                                agcagcagc‑AATTGCAGTAAAATTCCGCACCATTTTGAAATAAGCTGGCGTTGATGCCAGCGGCAAACCGAATTAATCAAAGGTGAGAGGCACATGCCGGTAATTAAAGTACGTGAAAACAAGCCGTTCGACGTAGCTCTGCGTCGCTTCAAGCGTTCCTGCGAAAAAGCAGGTGTTCT‑GGCGGAAGTTCGTCGTCGTGAG‑TTCTATGAAAAACCGACTACCGCACGt                                                                                                                                                                                >  2:508117/1‑228 (MQ=255)
                                                  cagcagc‑aATTGCAGTAAAATTCCGCACCATTTTGAAATAAGCTGGCGTTGATGCCAGCGGCAAACCGAATTAATCAAAGGTGAGAGGCACATGCCGGTAATTAAAGTACGTGAAAACGAGCCGTTCGACGTAGCTCTGCGTCGCTTCAAGCGTTCCTGCGAAAAAGCAGGTGTTCT‑GACGGAAGTTCGTCGCCGTGTG‑TTCTTTGAAAAACGGACCACCGAAccaa                                                                                                                                                                               >  2:23126/1‑224 (MQ=255)
                                                  cagcagc‑aATTGCAGTAAAATTCCGCACCATTTTGAAATAAGCTGGCGTTGATGCCAGCGGCAAACCGAATTAATCAAAGGTGAGAGGCACATGCCGGTAATTAAAGTACGTGAAAACGAGCCGTTCGACGTAGCTCTGCGTCGCTTCAAGCGTTACTGCGAAAAAGCAGGAGTTCT‑AGCGGAACTTCGTCGTCGGGAG‑TTCTATGAAAAACCGCCTACCGAACATaa                                                                                                                                                                              >  2:363051/1‑228 (MQ=255)
                                                          aaTTGCAGTAAAATTCCGCACCATTTTGAAAAAAGCGGGCGTTGATGCCAGCGGCAAACCGAATTAATCAAAGGTGAGAGGCACATGCCGTTAATTAAAGTACGTGAAAACGAGCCGTTCGACGTAGCTCTGCGTCGCTTCAAGCGTTCCTGCGAAAAAGCAGGTGTTCT‑GGCGGAAGTTCGTCGTCGTGAG‑TTCTATGAAAAACCGACTACCGAACGTAAGCGCGCt                                                                                                                                                                       <  2:372063/228‑1 (MQ=255)
                                                           atggCAATAAAATTCCGAACAATTTTGAAATAAGCGGGCGTTGACGCCAACGGCAAACCGAATTAATCAAAGGTGAGAGGCACATGCCGGTAATTAAAGTACGTGAAAACGAGCCGTTCGACGTAGCTCTGCGTCGCTTCAAGCGTTCCTGCGAAAAAGCAGGTGTTCT‑GGCGGAAGTTCGTCGTCGTGAG‑TTCTATGAAAAACCGACTACCGAACGTAAGCGCGCTa                                                                                                                                                                      <  2:652716/225‑1 (MQ=255)
                                                            ttGCAGTAAAATTCCGCACCATTTTGAAATAAGCGTGCGTTGATGCCAGCGGCAAACCGAATTAATCAAAGGTGAGAGGCACATGCCGGTAATTAAAGTACGTGAAAACGAGCCGTTCGACGTAGCTCTGCGTCGCTTCAAGCGTTCCTGCGAAAAAGCAGGTGTTCT‑GGCGGAAGTTCGTCGTCGTGAG‑TTCTATGAAAAACCGACTACCGAACGTAAGCGCGCTaa                                                                                                                                                                     <  2:148238/228‑1 (MQ=255)
                                                                     aaTTCCGCACCATTTTGAAATAAGCTGCCGTTGATGCCAGCGGCAAACCGAATTAATCAAAGGTGAGAGGCACATGCCGGTAATTAAAGTACGTGAAAACGAGCCGTGCGACGTAGCTCTGCGTCGCTTCAAGCGTTCCTGCGAAAAAGCAGGTGTTCT‑GGCGGAAGTTCGTCGTCGTGAG‑TTCTATGAAAAACCGACTACCGAACGTAAGCGCGCTAAAGCTTCTg                                                                                                                                                             <  2:296387/227‑1 (MQ=255)
                                                                            cACCATTTTGAAATAAGCTGGCGTTGATGCCAGCGGCAAACCGAATTAATCAAAGGTGAGAGGCACATGCCGGTAATTAAAGTACGTGAAAACGAGCCGTTCGACGTAGCTCTGCGTCGCTTCAAGCGTTCCTGCGAAAAAGCAGGTGTTCT‑GGCGGAAGTCCGTCGTCGTTAG‑TTCTATGAAAAACCGACTACCGAACGAAAGCGCGCTAAAGCTTCTgacgggaac                                                                                                                                                     >  2:478544/1‑220 (MQ=255)
                                                                               aaTTTGGAAATAAGGTGGCTGGGAAGCCAGCGGAAAACCGAAATAAACAAGGGGGAGGGGCACATGCCGGTAATTAAAGTACGTGAAAACGAGCCGTTCGACGTAGCTCTGCGTCGCTTCAAGCGTTCCTGCGAAAAAGCAGGTGTTCT‑GGCGGAAGTTCGTCGTCGTGAG‑TTCTATGAAAAACCGACTACCGAACGTAAGCGCGCTAAAGCTTCTGCAGTGAAACGt                                                                                                                                                  <  1:90696/227‑1 (MQ=255)
                                                                                 ttttGAAATAAGCTGGCGTTGATGCCAGCGGCAAACCGAATTAATCAAAGGTGAGAGGCACATGCCGGTAATTAAAGGACGTGAAAACGAGCCGTTCGACGTAGCTCTGCGTCGCTTCAAGCGTTCCTGCGAAAAAGCAGGTGTTCT‑GGCGGAAGTTCGTCGTCGTGAG‑TTCTATGAAAAACCGACTACCGAACGTAAGCGCGCGAAAGCTTCCGCAGGCAAACGTCa                                                                                                                                                >  2:379443/1‑228 (MQ=255)
                                                                                           aGCTGGCGTTGATGCCAGCGGCAAACCGAATTAATCAAAGGTGAGAGGCACATGCCGGTAATTAAAGTACGTGAAAACGAGCCGTTCGACGTAGCTCTGCGTCGCTTCAAGCGTTCCTGCGAAAAAGCAGGTGTTCT‑GGCGGAAGTTCGTCGTCGTGAG‑TTCTATGAAAAACCGACTACCGAACGTaa                                                                                                                                                                              >  2:5268/1‑188 (MQ=255)
                                                                                                 ggtgagtGCCAGCGGAAAACACGATAAATCAAGGGGGAGGGGCACTGGCCGGTAATTAAAGTAGGTGAAAACGAGCCGTTCGACGTAGCTCTGCGTCGCTTCAAGCGTTCCTGCGAAAAAGCAGGGGTTCT‑GGCGGAAGTTCGTCGTCGTGAG‑TTCTATGAAAAACCGACTACCGAACGTAAGCGCGCTAAAGCTTCTGCAGTGAAACGTCACGCGAAGAAACTGGCt                                                                                                                                <  1:380262/222‑1 (MQ=255)
                                                                                                                       aaTTAATCAAAGGTGAGAGGCACATGCCGGTAATTAAAGTACGTGAAAACGAGCCGTTCGACGTCGCTCTGCGTCGCTTCAAGCGTTCCTGCGAAAAAGCAGGGGTTCT‑GGCGGACGTTCGTCGTCGTGAG‑TTCTATGAAAAACCGACAACCGAACGAAAGCGCGCTAAAGCTTCTTAAGTGAAACTACACGCGAAGAAACAGGCTCGCCGAAAACACA‑GCCGC‑Accc                                                                                                        >  1:618530/1‑226 (MQ=255)
                                                                                                                       aaTTAATCAAAGGTGAGAGGCACATGCCGGTAATTAAAGTACGTGAAAACGAGCCGTTCGACGTAGCTCTGCGTCGCTTCAAGCGTTCCTGCGAAAAAGCAGGTGTTCT‑GGCGGAAGTTCGTCGTCGTGAG‑TTCTATGAAAAACCGACTACCGAACGTAAg                                                                                                                                                                             >  2:570531/1‑161 (MQ=255)
                                                                                                                        atgaaTCAAAGGTGAGAGGCACATGCCGGTAATTAAAGTACGTGAAAACGAGCCGTTCGACGTAGCTCTGCGTCGCTTCAAGCGTTCCTGCGAAAAAGCAGGTGTTCT‑GGCGGAAGTTCGTCGTCGTGAG‑TTCTATGAAAAACCGACTACCGAACGTAAgc                                                                                                                                                                            <  1:570531/158‑1 (MQ=255)
                                                                                                                         ttaaccaaggggggggcgCACATGCCGGTAATTAAAGTACGTGAAAGCGAGCCGTTCGACGTAGCTCTGCGTCGCTTCAAGCGTTCCTGCGAAAAAGCAGGTGTTCG‑GGCGGAAGTTCGTCGTCGTGAG‑TTCTATGAAAAACCGACTACCGAACGTAAGCGCGCTAAAGCTTCTGCAGTGAAACGTCACGCGAAGAAACTGGCTCG‑CGAAAACGCACGCCGC‑ACTCGt                                                                                                      <  1:263820/211‑1 (MQ=255)
                                                                                                                         ttaacaaaaGGTGGAGGGCACATGCCGGTAATTAAAGTATGTGAAAACGAGCCGTTCGACGTAGCTCTGCGTAGCTTCAAGCGTTCCTGCGAAAAAGCAGGTGTTCT‑GGCGGAAGTTCGTCGTCGTGAG‑TTCTATGAAAAACCGACTACCGAACGTAAGCGCGCTAAAGCTTCTGCAGTGAAACGTCACGCGAAGAAACTGGCTCG‑CGAAAACGCACGCCGC‑ACTCGt                                                                                                      <  1:385472/222‑1 (MQ=255)
                                                                                                                           aaTCAAGGGGGAGAGCGACATGCGGGTAATTAAAGTAGTTGAAAAGGAGCGTGTGGACGTAGGTCTGGGGCGCTTCAAGGTTTCCTGCGAAAAAGCAGGTGTTCT‑GGCGGAAGTTCGTCGTCGTGAG‑TTCTATGAAAAACCGACTACCGAACGTAAGCGCGCTAAAGCTTCTGCAGTGAAACGTCACGCGAAGAAACTGGCTCG‑CGAAAACGCACGCCGC‑ACTCGTc                                                                                                     <  1:285829/227‑1 (MQ=255)
                                                                                                                                   ggaagagGCACATGCCGGTAATTAAAGTACGTGAAAACGAGCCGTTCGACGTAGCTCTGCGTCGCTTCAAGCGTTCCTGCGAAAAAGCAGGTGTTCT‑GGCGGAAGTTCGTCGTCGTGAG‑TTCTATGAAAAACCGACTACCGAACGTAAGCGCGCTAAAGCTTCTGCAGTGAAACGTCACGCGAAGAAACTGGCTCG‑CGAAAACGCACGCCGC‑ACTCGTCTGTACTaa                                                                                             <  2:105912/224‑1 (MQ=255)
                                                                                                                                                      aaTTAAAGGTCGGTAAAACGGGCCGTTCGGCGTAGCCCTGCGTTCCTGCAAGCGTTCCTTCGAAAAAGCAGGGGTTCT‑GGCGGGAGTTCGTCGTCGTGAG‑TTCTATGAAAAACCGACTACCGAACGTAAGCGCGCTAAAGCTTCTGCAGTGAAACGTCACGCGAAGAAACTGGCTCG‑CGAAAACGCACGCCGC‑ACTCGTCTGTACTAATTCCCCGAGAGCGTTGctc                                                                          <  1:430552/227‑1 (MQ=255)
                                                                                                                                                        ttAAAGTACGTGAAAACGAGCCGTTCGACGTAGCTCTGCGTCGCTTCAAGCGTTCCTGCGAAAAAGCAGGTGTTCT‑GGCGGAAGTTCGTCGTCGTGAG‑TTCTATGAAAAACCGACTACCGAACGTAAGCGCGCTAAAGCTTCTGCAGTGAAACGTCACGTGAATAAACCGGCTCG‑CAGAAACGCACGCCCG‑CCACGACTGTCCTAATCcccccaaggcggtcctaccc                                                                       >  1:601642/1‑212 (MQ=255)
                                                                                                                                                            aGTACGTGAAAACGAGCCGTTCGACGTAGCTCTGCGTCGCTTCAAGCGTTCCTGCGAAAAAGCAGGTGTTCT‑GGCGGAAGTTCGTCGTCGTGAG‑TTCTATGAAAAACCGACTACCGAACGTAAGCGCGCTAAAGCTTCTGCAGTGAAACGTCACGCGAAGAAACTGGCTCG‑CGAAAACGCACGCCGC‑ACTCGTCTGTACTAATTCCCCCAGAGCGTTGCTCCCCgcgc                                                                   >  2:514340/1‑225 (MQ=255)
                                                                                                                                                                cgggAAAACGAGCCGTTCGACGTAGCTTTGCGTCCCTTCAAGCGTTCCTGCGAAAAAGCAGGTGTTCT‑GGCGGAAGTTCGTCGTCGTGAG‑TTCTATGAAAAACCGACTACCGAACGTAAGCGCGCTAAAGCTTCTGCAGTGAAACGTCACGCGAAGAAACTGGCTCG‑CGAAAACGCACGCCGC‑ACTCGTCTGTACTAATTCCCCGAGAGCGTTGCTCTCCGATCAGa                                                                <  2:375476/224‑1 (MQ=255)
                                                                                                                                                                cGTGAAAACGAGCCGTTCGACGTAGCTCTGCGTCGCTTCAAGCGTTCCTGCGAAAAAGCAGGTGTTCT‑GGCGGAAGTTCGTCGTCGTGAG‑TTCTATGAAAAACCGACTACCGAACGTAAGCGCGCTAAAGCTTCTGCAGTGAAACGTCACGCGACGAAACTGGCTCG‑CGAAAACGCACGCCCC‑CCTCGCCTGCCGAAATTACCCCAGAACGTTGCCCTCCGATCaacc                                                               >  1:514247/1‑225 (MQ=255)
                                                                                                                                                                cGTGAAAACGAGCCGTTCGACGTAGCTCTGCGTCGCTTCAAGCGTTCCTGCGAAAAAGCAGGTGTTCT‑GGCGGAAGTTCGTCGTCGTGAG‑TTCTATGAAAAACCGACTACCGAACGTAAGCGCGCTAAAGCTTCTGCAGTGAAACGTCACGCGAAGAACCTGGCTCG‑CAAAAACCCACGCCGC‑ACCCGTCCGTCCTAATTCCCCCACAGCGGTACCCTCCgcacacac                                                               >  1:255865/1‑221 (MQ=255)
                                                                                                                                                                cGTGAAAACGAGCCGTTCGACGTAGCTCTGCGTCGCTTCAAGCGTTCCTGCGAAAAAGCAGGTGTTCT‑GGCGGAAGTTCGTCGTCGTGAG‑TTCTATGAAAAACCGACTACCGAACGTAAGCGCGCTAAAGCTTCTGCAGTGAAACGTCACGCGAAGAAACTGGCTCG‑CGAAAACGCACGCCGC‑ACTCGTCTGTACTAATTCCCCGAGAGCGTTGCTCTCCGATCAGAc                                                               >  2:234487/1‑228 (MQ=255)
                                                                                                                                                                   gAAAACGAGCCGTTCGACGTAGCTCTGCGTCGCTTCAAGCGTTCCTGCGAAAAAGCAGGTGTTCT‑GGCGGAAGTTCGTCGTCGTGAG‑TTCTATGAAAAACCGACTGCCGAACGTAAGCGCGCTAAAGCTTCTGCAGTGAAACGTCACGCGAAGAAACTGGCTCG‑CGAAAACGCAAGCCGC‑ACTCGTCTGCAATAATTCCCCGAGAGGGTTGCTCTCCGATCAGACCgg                                                            >  2:37881/1‑227 (MQ=255)
                                                                                                                                                                   gAAAACGAGCCGTTCGACGTAGCTCTGCGTCGCTTCAAGCGTTCCTGCGAAAAAGCAGGTGTTCT‑GGCGGAAGTTCGTCGTCGTGAG‑TTCTATGAAAAACCGACTACCGAACGTAAGCGCGCTAAAGCt                                                                                                                                                                 >  1:135462/1‑129 (MQ=255)
                                                                                                                                                                    aaaaCGAGCCGTTCGACGTAGCTCTGCGTCGCTTCAAGCGTTCCTGCGAAAAAGCAGGTGTTCT‑GGCGGAAGTTCGTCGTCGTGAG‑TTCTATGAAAAACCGACTACCGAACGTAAGCGCGCTAAAGCtt                                                                                                                                                                <  2:135462/129‑1 (MQ=255)
                                                                                                                                                                          gagctgtgtatgttgcttttggTCGCTTCAAGGGTTCCTGCGAAAAAGCAGGTGTTCT‑GGCGGAAGTTCGTCGTCGTGAG‑TTCTATGAAAAACCGACTACCGAACGTAAGCGCGCTAAAGCTTCTGCAGTGAAACGTCACGCGAAGAAACTGGCTCG‑CGAAAACGCTCGCCGC‑ACTCGTCTGTACTAATTCCCCGAGAGCGTTGCTCTCCGATCAGACCGagttgt                                                       <  2:290737/205‑1 (MQ=255)
                                                                                                                                                                          atccGTTCGACGTAGCGCTGCTTCGCTGCAAGCGTTCCTGCGAAAAAGCAGGTGTTAT‑GGCGGAAGTTCGTCGTCGTGAT‑TTCTATGAAAAACCGACTACCGAACGTAAGCGCGCTAAAGCTTCTGCAGTGAAACGTCACGCGAAGAAACTGGCTCG‑CGAAAACGCACGCCGC‑ACTCGTCTGTACTAATTCCCCGAGAGCGTTGCTCTCCGATCAGACCGagttgt                                                       <  2:412608/224‑1 (MQ=255)
                                                                                                                                                                              gTTTGACTTTGCTCTGCGCTGCTTCAAGCGTTCCTGCAAAAAAGCAGGTGTTCT‑GGCGTAAGTTCGTCGTCGTGAG‑TTCTATGAAAACCCGACTACCGAACGTAAGCGCGCTAAAGCTTCTGCAGTGAAACGTCACGCGAAGAAACTGGCTCG‑CGAAAACGCACGCCGC‑ACTCGTCTGTACTAATTCCCCGAGAGCGTTGCTCTCCGATCAGACCGagttgtagttg                                                  <  1:150760/227‑1 (MQ=255)
                                                                                                                                                                              gTTCGACGTAGCTCTGCGTCGCTTCAAGCGTTCCTGCGAAAAAGCAGGCGTTCT‑GGCGGAAGGTCGCCGTCGTGAG‑TTCTATGAAAAACCGACTACCGAACGTAAGCGCGCTAAAGCTTCTGCAGTAAAACGTCACCCGCAGAAACAGGCCCG‑CGAAAACGCAC‑CCGCAACCCGTCTGGCCCAATTCCCCGGAAGCGTGCCCCTCCGACCAAAcccggtgtgcgtgta                                                 >  1:248100/1‑215 (MQ=255)
                                                                                                                                                                                      tAGCTCTGCGTCGCTTCAAGCGTTCCTGCGAAAAAGCAGGTGTTCT‑GGCGGAAGTTCGTCGTCGTGAG‑TTCTATGAAAAACCGACTACCGAACGTAAGCGCGCTAAAGCTTCTGCAGTGAAACGTCACGCGAAGAAACTGGCTCG‑CGAAAACGCACGCCCC‑ACCCGTCGCGACTAATTCCCCGACGGCGTTGGTCTCCCATACGACCggctgtgctgtgtaaggccgg                                         >  1:146289/1‑208 (MQ=255)
                                                                                                                                                                                        ggttggCGTCCCTTCAAGCGTTCCTGCGAAACAGCAGGTGTTCT‑GGCGGAAGTCCGTCGTCTTGAG‑TTCGAAGAAAACCCGACTACCGAACGTAAGCGCGCTAAAGCTTCTGCAGTGAAACGTCACGCGAAGAAACTGGCTCG‑CGAAAACGCACGCCGC‑ACTCGTCTGTACTAATTCCCCGAGAGCGTTGCTCTCCGATCAGACCGAGTTGTAGTTGTAAGGCCGTg                                        <  1:598827/222‑1 (MQ=255)
                                                                                                                                                                                         ctctGCGTCGCTTCAAGCGTTCCTGCGAAAAAGCAGGTGTTCT‑GGCGGAAGTTCGTCGTCGTGAG‑TTCTATGAAAAACCGACTACCGAACGTAAGCGCGCTAAAGCTTCTGCAGTGAAACGTCACGCGAAGAAACTGGCTCG‑CGAAAACGCACGCCGC‑ACTCGTCTGTACTAATTCCCCGAGAGCGTTGCTCTCCGGTCAGACCGAGTTGTAGTTGTAAGGcctgggt                                      >  2:408866/1‑223 (MQ=255)
                                                                                                                                                                                         ctctGCGTCGCTTCAAGCGTTCCTGCGAAAAAGCAGGTGTTCT‑GGCGGAAGTTCGTCGTCGTGAG‑TTCTATGAAAAACCGACTACCGAACGTAAGCGCGCTAAAGCTTCTGCAGTGAAACGCCACGCGAAGAAACTGGCTCG‑CGAAAACGCACGCCGC‑AATCGTCTGTACTAATTCCCCGACAGCGTTGCTCTCCGATCAGACCGCGTTGTAGTTGTAAGGcctggct                                      >  2:79138/1‑223 (MQ=255)
                                                                                                                                                                                         ctctGCGTCGCTTCAAGCGTTCCTGCGAAAAAGCAGGGGTTCT‑GGCGGAAGTTCGTCGTCGTGAG‑TTCTATGAAAAACCGACTACCGAACGTAAGCGCGCTAAAGCTTCTGCAGTGAAACGTCACGCGAAGAAACAGGCTCG‑CGAAAACGCACGCCGC‑CCGCGGCTGTGCTAATTCCCCGAGAGCGGTTCTCCCCGGCCAAACCCAGTTTTAGTTGTAAGGCCGGGcg                                      >  1:419657/1‑227 (MQ=255)
                                                                                                                                                                                             gCGTCGCTTGAAGCGTTCGTGAGAAAAAGCAGGTGTTCT‑GGCGGAAGTTCGTCGTCGTGAG‑TTCTAAGAAAAACCGACTACCGAACGTAAGCGCGCTAAAGCTTCTGCAGTGAAACGTCACGCGAAGAAACTGGCTCG‑CGAAAACGCACGCCGC‑ACTCGTCTGTACTAATTCCCCGAGAGCGTTGCTCTCCGATCAGACCGAGTTGTAGTTGTAAGGCCGTg                                        <  2:264454/222‑1 (MQ=255)
                                                                                                                                                                                              gtcgcgttaaaggggtccttggAAAAAGAGGGTGGTCT‑GGCGAGAGTCCGCGCTGGGGAG‑GTTTATGAAAAACCGACCACCCAAAGTAAGCCCGCTAAAGCTTCTGCAGTGAAACGTCACGCGAAGAAACTGGCTCG‑CGAAAACGCACGCCGC‑ACTCGTCTGTACTAATTCCCCGAGAGCGTTGCTCTCCGATCAGACCGAGTTGTAGTTGTAAGGCCGTGCTTCCGa                                 <  1:642754/207‑1 (MQ=255)
                                                                                                                                                                                              ggTCGTTTCAAGCGTTCCTGCAAAAAAGCAGGTTTTCC‑GGCGGAAGGTCGGCGTCGTGAG‑TTCTAGGAAAAACCGACTACCGAACGTAAGCGCGCTAAAGCTTCTGCAGTGAAACGTCACGCGAAGAAACTGGCTCG‑CGAAAACGCACGCCGC‑ACTCGTCTGTACTAATTCCCCGAGAGCGTTGCTCTCCGATCAGACCGAGTTGTAGTTGTAAGGCCGTGCTTCCGa                                 <  1:408866/227‑1 (MQ=255)
                                                                                                                                                                                               gTCGCTTCAAGCGTTCCTGCGAAAGAGCAGGGATGCT‑GGCGGAAGTTCGTCGTCGTAAG‑TTCTAGGAAAAACCGACTACCGAACGTAAGCGCGCTAAAGCTTCTGCAGTAAAACGTCACGCGAAGAAACTGGCTCG‑CGAAAACGCACGCCGC‑ACTCGTCTGTACTAATTCCCCGAGAGCGTTGCTCTCCGATCAGACCGAGTTGTAGTTGTAAGGCCGTGCTTCCGa                                 <  2:248100/227‑1 (MQ=255)
                                                                                                                                                                                               gTCGCTTCAAGCGTTCCTGCGAAAAAGCAGGTGTTCT‑GGCGGAAGTTCGTCGTCGTGAG‑TTCTATGAAAAACCGACTACCGAACGTAAGCGCGCTAAAGCTTCTGCAGTGAAACGTCACGCGAAGAAACTGGCTCG‑CGAAAACGCACGCCGC‑ACTCGTCTGTACCAATTCCCCCAGAGCGTGGCTCCCCGACCCGCCCCGGTTGTAGTTGTAAGGCCGAGCTTCCGaa                                >  1:41038/1‑228 (MQ=255)
                                                                                                                                                                                                   cTTCAAGCGTTCCTGCGAAAAAGCAGGTGTTCT‑GGCGGAAGTTCGTCGTCGTGAG‑TTCTATGAAAAACCGACTACCGAACGTAAGCGCGCTAAAGCTTCTGCAGTGAAACGTCACGCGAAGAAACTGGCTCG‑CGAAAACGCACGCCGC‑ACTCGTCTGTACTAATTCCCCGAGAGCGTTGCTCTCCGATCAGACCGAGTTGTGGTTGTAAGGCCGTGCTTCCGAAAGGa                            >  2:434442/1‑228 (MQ=255)
                                                                                                                                                                                                        atggTTCCAGCGAAAAAGC‑GGTGTTTTCGGCGGAAGTTCGGCGGAGTGAG‑TTCTATGAAAAACCGACTACCGAACTTAAGCCCGCTAAAGCTTCTGCAGTGAAACGGCACGCGAAGAAACTGGCTCG‑CGAAAACGCACGCCGC‑ACTCGTCTGTACTAATTCCCCGAGAGCGTTGCTCTCCGATCAGACCGAGTTGTAGTTGTAAGGCCGTGCTTCCGAAAGGAATgcg                       <  1:400228/225‑1 (MQ=255)
                                                                                                                                                                                                                 gCGAAAAAGCAGGTGTTCT‑GGCGGAAGTTCGTCGTCGTGAG‑TTCTATGAAAAACCGACTACCGAACGTAAGCGCGCTAAAGCTTCTGCAGTGAAACGTCACGCGAAGAAACTGGCTCG‑CGAAAACGCACGCCGC‑ACTCGTCTGTACTCATTCCCCGAGAGCGTTCCTCTCCGATCGGACCGGGTTGTAGTTGTAAGGCCGTGCTCCCGAAAGGAACTCGCGGCTTAtt              >  2:558931/1‑228 (MQ=255)
                                                                                                                                                                                                                             gtgtTCT‑GGGGGAAGGTCGTCGTGGTGAG‑TTCTATGAAAAACCGACTACCGAACGTAAGCGCGCTAAAGCTTCTGCAGGGAAACGTCACGCGAAGAAACTGGCTCG‑CGAAAACGCACGCCGC‑ACTCGTCTGTACTAATTCCCCGAGAGCGTTGCTCTCCGATCAGACCGAGTTGTAGTTGTAAGGCCGTGCTTCCGAAAGGAATGCGCGGCTTATTTTCGTTTATGaa  <  1:422301/228‑1 (MQ=255)
                                                                                                                                                                                                                              |                                                                                                                                                                                                                                      
CTCACCCTTATAAAAGTCCCTTTCAAAAAAGGCCGCGGTGCTTTACAAAGCAGCAGC‑AATTGCAGTAAAATTCCGCACCATTTTGAAATAAGCTGGCGTTGATGCCAGCGGCAAACCGAATTAATCAAAGGTGAGAGGCACATGCCGGTAATTAAAGTACGTGAAAACGAGCCGTTCGACGTAGCTCTGCGTCGCTTCAAGCGTTCCTGCGAAAAAGCAGGTGTTCT‑GGCGGAAGTTCGTCGTCGTGAG‑TTCTATGAAAAACCGACTACCGAACGTAAGCGCGCTAAAGCTTCTGCAGTGAAACGTCACGCGAAGAAACTGGCTCG‑CGAAAACGCACGCCGC‑ACTCGTCTGTACTAATTCCCCGAGAGCGTTGCTCTCCGATCAGACCGAGTTGTAGTTGTAAGGCCGTGCTTCCGAAAGGAATGCGCGGCTTATTTTCGTTTATGAA  >  NC_000913/3210638‑3211085

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