Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,210,861 T→G 21.3% 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,8610TG21.3% 114.4 / 3.2 68G27G (GGT→GGG) rpsU30S ribosomal subunit protein S21
Reads supporting (aligned to +/- strand):  ref base T (24/23);  new base G (16/3);  total (42/26)
Fisher's exact test for biased strand distribution p-value = 1.40e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTCACCCTTATAAAAGTCCCTTTCAAAAAAGGCCGCGGTGCTTTACAAAGCAGCAGC‑AATTGCAGTAAAATTCCGCACCATTTTGAAATAAGCTGGCGTTGATGCCAGCGGCAAACCGAATTAATCAAAGGTGAGAGGCACATGCCGGTAATTAAAGTACGTGAAAACGAGCCGTTCGACGTAGCTCTGCGTCGCTTCAAGCGTTCCTGCGAAAAAGCAGGTGTTCT‑GGCGGAAGTTCGTCGTCGTGAG‑TTCTATGAAAAACCGACTACCGAACGTAAGCGCGCTAAAGCTTCTGCAGTGAAACGTCACGCGAAGAAACTGGCTCG‑CGAAAACGCACGCCGCACTCGTCTGTACTAATTCCCCGAGAGCGTTGCTCTCCGATCAGACCGAGTTGTAGTTGTAAGGCCGTGCTTCCGAAAGGAATGCGCGGCTTATTTTCGTTTATGAA  >  NC_000913/3210640‑3211087
                                                                                                                                                                                                                              |                                                                                                                                                                                                                                     
ctcACCCTTATAAAAGTCCCTTTCAAAAAAGGCCGCGGTGCTTTACAAAGCAGCAGC‑AATTGCAGTAAAATTCCGCACCATTTTGAAATAAGCTGGCGTTGATGCCAGCGACAAACCGAATTAATCAAAGGTGAGAGGCACATGCCGGTAATTAAAGTACGTGAAAACAATCCGTTAGACGTAACTCTGAGTCGCTTCAAGCGTTCCTGCGAACAAACAGGGGTTCT‑g                                                                                                                                                                                                                                >  2:61093/1‑228 (MQ=255)
 tcACTCTTATAAAAGTCCCTTTCAAAAAAGGCCGCGGTGCTTTACAAAGCAGCAGC‑AATTGCAGTAAAATTCCGCACCATTTTGAAATAAGCTGGCGTTGATGCCAGCGGCAAACCGAATTAATCAAAGGTAAACGGCCCATGCCGGTAATTAAAACACCTGAAACCAAACCGTTCGACCAAGCCCTGCGGTCCTTCAAGCCTGCCCGCGAAAAAACCGGGGtttgg                                                                                                                                                                                                                                 >  2:694766/1‑224 (MQ=255)
  cACCCTTATAAAAGTCCCTTTCAAAAAAGGCCGCGGTGCTTTACAAAGCAGCAGC‑AATTGCAGTAAAATTCCGCACCATTTTGAAATAAGCTGGCGTTGATGCCAGCGGCAAACCGAATTAATCAAAGGTTAGAGTCACATGCCGGTAATTAAAGTACGTGAAAACAAGCCGTTAGACGTAACTCTGCGTAGCTTCAAGCTTTCCGCGGCAAAAGAAGGTGTTTT‑gcg                                                                                                                                                                                                                              >  1:642754/1‑226 (MQ=255)
   aCCCTTATAAAAGTCCCTTTCAAAAAAGGCCGCGGTGCTTTACAAAGCAGCAGC‑AATTGCAGTAAAATTCCGCACCATTTTGAAATAAGCTGGCGTTGCTGCCAGCGGCAAACCGAATTAATCAAAGGTGAGAGGCAAATGCCGGGAATTAAGGTACGAGAAACCGAGCCGTTCGACGCAGCTCGCCGGCGCTTCAGACGTTCCTTCGTAAAACCGGGTGTTTT‑GCCg                                                                                                                                                                                                                             >  1:391539/1‑228 (MQ=255)
   aCCCTTATAAAAGTCCCTTTCAAAAAAGGCCGCGGTGCTTTACAAAGCAGCAGC‑AATTGCAGTAAAATTCCGCACCATTTTGAAATAAGCTGGCGTTGATGCCAGCGGCAAACCGAATTAATCAAAGGTGAGAGGCACATGCCGGTAATTAAAGTACGTGAAAACGAGCCGTTCGACGTAGCTCTGCGTCGCTTCAAGCGTTCCTGCGAAAAAGCAGGTGTTCT‑GGc                                                                                                                                                                                                                              >  1:422301/1‑227 (MQ=255)
   aCCCTTATAAAAGTCCCTTTCAAAAAAGGCCGCGGTGCTTTACAAAGCAGCAGC‑AATTGCAGTAAAATTCCGCACCATTTTGAAATAAGCTGGCGTTGATGCCAGCGGCAAACCGAATTAATCAAAGGTGAGAGGCACATGCCGGTAATTAAAGTACGTGAAAACAAGCCGTTCGACGTAGCTCTGCGTCGCTTCAAGCGTTCCTGCGAAAAAGCAGGGGTTCT‑gggg                                                                                                                                                                                                                             >  2:489255/1‑226 (MQ=255)
        tataAAAGTCCCTTTCAAAAAAGGCCGCGGTGCTTTACAAAGCAGCAGC‑AATTGCAGTAAAATTCCGCACCATTTTGAAATAAGCTGGCGTTGATGCCAGCGGCAAACCGAATTAATCAAAGGTGAGAGGCACATGCCGGTAATTAAAGTACGTGAAAACGAGCCGTTCGACGTAGCTCTGCGTCGCTTCAAGCGTTCCTGCGAAAAAACAGGGGTTAT‑GGCGGCAgg                                                                                                                                                                                                                        >  2:375476/1‑227 (MQ=255)
        tataAAAGTCCCTTTCAAAAAAGGCCGCGGTGCTTTACAAAGCAGCAGC‑AATTGCAGTAAAATTCCGCACCATTTTGAAATAAGCTGGCGTTGATGCCAGCGGCAAACCGAATTAATCAAAGGTGAGAGGCACATGCCGGTAATTAAAGTACGTGAAAACGAGCCGTTCGACGAAGCTCTGCGTCGCTTCAAGCGTTCCTGCGAAAAAGCAGGTGTGCT‑GGCGGAAGt                                                                                                                                                                                                                        >  1:427007/1‑228 (MQ=255)
        tataAAAGTCCCTTTCAAAAAAGGCCGCGGTGCTTTACAAAGCAGCAGC‑AATTGCAGTAAAATTCCGCACCATTTTGAAATAAGCTGGCGTTGATGCCAGCGGCAAACCGAATTAATCAAAGGTGAGAGGCACATGCCGGTAATTAAAGTACGTGAAAACGAACCGTTCGACTTAGCTCTGCGTCCCTTCAACCGTTCCTGCGAAAAACAAGGGGTTCT‑GGCGGAAGt                                                                                                                                                                                                                        >  2:653326/1‑228 (MQ=255)
        tataAAAGTCCCTTTCAAAAAAGGCCGCGGTGCTTTACAAAGCAGCAGC‑AATTGCAGTAAAATTCCGCACCATTTTGAAATAAGCTGGCGTTGATGCCAGCGGCAAACCGAATTAATCAAAGGTGAGAGGCACATGCCGGTAATTAAAGTACGGGAAAACAAGCCGTTCAACGAAGCCCTGGGGCGCTTCAGGCGTTCCTGCCAAAAAGCAGGGGTTCG‑GGGGGAAGt                                                                                                                                                                                                                        >  2:343395/1‑228 (MQ=255)
         ataAAAGTCCCTTTCAAAAAAGGCCGCGGTGCTTTACAAAGCAGCAGC‑AATTGCAGTAAAATTCCGCACCATTTTGAAATAAGCTGGCGTTGATGCCAGCGGCAAACCGAATTAATCAAAGGTAAGAGGCACATGCCGGTAATTAAAGTACGCGAAAACGAGCCGTTCGACGTAGCTCTGCATCGCTTCAAGCGTTTCTGCTGAAAAGAAGGTGTTCT‑GCCGCAAGtt                                                                                                                                                                                                                       >  1:234339/1‑228 (MQ=255)
          taAAAGTCCCTTTCAAAAAAGGCCGCGGTGCTTTACAAAGCAGCAGC‑AATTGCAGTAAAATTCCGCACCATTTTGAAATAAGCTGGCGTTGATGCCAGCGGCAAACCGAATTAATCAAAGGTGAGAGGCCCATGCCGGTAATTAAAGTACGTGAAAACGAGCCGTTCGCCGTAGCTCTGCGTCGCTTCAAGCGTTCATCCGAAAAAGCAGGGGTTCT‑GCCGGAAGTtc                                                                                                                                                                                                                      >  1:263820/1‑228 (MQ=255)
          taAAAGTCCCTTTCAAAAAAGGCCGCGGTGCTTTACAAAGCAGCAGC‑AATTGCAGTAAAATTCCGCACCATTTTGAAATAAGCTGGCGTTGATGCCAGCGGCAAACCGAATTAATCAAAGGTGAGAGGCACATGCCGGTAATTAAAGTACGTGAAAACGAGCCGTTCGACGTAGCTCTGCGTCGCTTCAAGCGTTCCTGCGAAAAAGCAGGGGTTCT‑GGCGGAAGCtc                                                                                                                                                                                                                      >  1:385472/1‑228 (MQ=255)
           aaaaATCCATTTAAAAAAAGGCCGCGGGGTTTTAAAAAGC‑GGAGCAAATGGCAAGAAAAGTCCGGAACAATTGGAAAAAAGCTGGGGTTGATGCAAGCGGAAAACCGAATAAATCAAAGGTGAGAGGCACATGCCGGTAATTAAAGTACGTGAAAACGAGCCGTTCGACGTAGCTCTGCGTCGCTTCAAGCGTTCCTGCGAAAAAGCAGGTGTTCT‑GGCGGAAGTtcg                                                                                                                                                                                                                     <  2:2848/228‑1 (MQ=255)
             aaGTCCCTTACAAAAACGGCCGCGGTGCTTTACAAAGCAGCAGC‑AATTGACGTAAAATTCCGCCCCATTTTGAAATAAGCTGGCGTTGATGCCAGCGGCAAACCGAATTAATCAAAGGTGAGAGGCACATGCCGGTAATTAAAGTACGTGAAAACGAGCCGTTCGACGTAGCTCTGCGTCGCTTCAAGCGTTCCTGCGAAAAAGCAGGTGTTCT‑GGCGGAAGTtcgtc                                                                                                                                                                                                                   <  1:666349/228‑1 (MQ=255)
               gTCCCTTTCAAAAAAGGCCGCGGTGCTTTACAAAGCAGCAGC‑AATTGCAGTAAAATTCCGCACCATTTTGAATTAACTTGGCGTTGATACCAGCGGCAAACCGAATTAATCAAAGGCGAGAGGCACATGCTGGTAATTAAAGTAAGTGATAACAAGCGGTTCGAAGTAGCTCTTCTTCGCTTCAAACGTTCCTGCGAACAATCAAGAGATCT‑GAGGGAAGttgccg                                                                                                                                                                                                                   >  1:99570/1‑222 (MQ=255)
               gTCCCTTTCAAAAAAGGCCGCGGTGCTTTACAAAGCAGCAGC‑AATTGCAGTAAAATTCCGCACCATTTTGAAATAAGCTGGCGTTGATGCCAGCGGCAAACCGAATTAATCAAAGGTGAGAGGCACATGCCGGTAATTAAAGTACGTGAAAACGAGCCGTTCGACGTAGCTCTGCGTCGCTTCAACCGTTCCTGCGAAAAAGCAGTGGTTCT‑GGCGGAAGTtcctcgt                                                                                                                                                                                                                 >  1:358697/1‑228 (MQ=255)
               gTCCCTTTCAAAAAAGGCCGCGGTGCTTTACAAAGCAGCAGC‑AATTGCAGTAAAATTCCGCACCATTTTGAAATAAGCTGGCGTTGATGCCAGCGGCAAACCGAATTAATCAAAGGTGAGAGGCACATGCCGGCAATTAAAGAACGTGAAAACGAGCCGTTCAACGTAGCTCTGCGTCCCTTTAAGCGTTACCGCGAAAAAGCAGGGTTTCT‑GGCGGAAATtcgtcgc                                                                                                                                                                                                                 >  2:662862/1‑227 (MQ=255)
               gTCCCTTTCAAAAAAGGCCGCGGTGCTTTACAAAGCAGCAGC‑AATTGCAGTAAAATTCCGCACCATTTTGAAATAAGCTGGCGTTGATGCCAGCGGCAAACCGAATTAATCAAAGGTAAGAGGCACATGCCGGCAATTAAAGTACGTGAAAACGAGCCGTTCGACGTAGCTCTGCGTCGCTTCAAGCGTTCCTGCGAAAAAGACGGTGTTCT‑GGCGGAAGGttgtctt                                                                                                                                                                                                                 >  2:496515/1‑226 (MQ=255)
                  tctgtCAAAAAAGGCCGCGGTGCTTTACAAAGCGGCAGC‑AATGGCAGTAAAATACCGAAACATTTTGAAATAAGCTGGCGTTGATTCCAGCGGCAAACCGAATTAGTCAAGGGTGAGGGGCACATGCCGGTAATTAAAGTACGTGAAAACGAGCCGTTCGACGTAGCTCTGCGTCGCTTCAAGCGTTCCTGCGAAAAAGCAGGTGTTCT‑GGCGGAAGTTcgtcgtcgt                                                                                                                                                                                                              <  2:639113/224‑1 (MQ=255)
                  ccTTTCAAAAAAGGCCGCGGTGCTTTACAAAGCAGCAGC‑AATTGCAGTAAAATTCCGCACCATTTTGAAATAAGCTGGCGTTGATGCCAGCGGCAAACCGAATTAATCAAAGGTTAGAGGCCCATGCCGGTAATTAAAGTACGGTAAAACGAGCCGGTCAACGTAACTCTACGTCGCTTCAATCGTTCCTCCGCAAAAGCAGGGGTCTA‑GGCGGAAGttcttctttgt                                                                                                                                                                                                              >  2:412608/1‑223 (MQ=255)
                    tttCAAAAAAGGCCGCGGTGCTTTACAAAGCAGCAGC‑AATTGCAGTAAAATTCCGCACCATTTTGAAATAAGCTGGCGTTGATGCCAGCGGCAAACCGAATTAATCAAAGGTGAGAGGCACATGCCGGTAATTAAAGTACGTGAAAACGACCCGTTCGACGTAGCTCTGCGTCGCTTCAAGCGTTCCTGCGAAAAAGCAGGGTTTCT‑GCCGGAAGTtcgtcgcccggg                                                                                                                                                                                                            >  2:147943/1‑222 (MQ=255)
                     ttaaaaaatGGCCCAGGTGCATTACAAGGAAGCAGC‑AAGAGCGGTAAAATCCCACACTTTTTAGAAATAAGCGAGCGTGGATGCGACCGGCAAACCGAATAAGTCAAAGGTGGAAGGGACATGCCGGAAATTAAAGTACGTGAAAACAAGCCGTACGACGTAGCTCTGGGTCGCTTCAATCGTTCCTGCGAAAAAGAAGGTGTTCT‑GGCGGAAGTTCGTCGTCGTGAg                                                                                                                                                                                                           <  2:99570/225‑1 (MQ=255)
                     ttCAAAAAAGGCCGCGGTGCTTTACAAAGCAGCAGC‑AATTGCAGTAAAATTCCGCACCATTTTGAAATAAGCTGGCGTTGATGCCAGCGGCAAACCGAATTAATCAAAGGAGAGAGGCACATGCCGGTAATTAAAGTACGTGAAAACTAGCCGTTCGACGTAGCTCTGCGTCGCTTCAAGCGTTACTGCGCAAAAGCAAGGGGTCT‑GGCGGAAGTTCGTCGTCGTggg                                                                                                                                                                                                           >  1:672154/1‑226 (MQ=255)
                                ccGGGGTGGTTTACAAAGCAGCAGC‑AAGTGCAGAAAAATTCCGCCCAATTTTGAAATAAGCTGGCGGTTATGCCAGCGGCAAACCAAATTAATCAAAGGTGAGAGGCACATGCCGGTAATTAAAGTACGTGAAACCGAGCCGTTCGACGTAGCTCTGCGTCGCTTCAAGCGTTCCTGCGAAAAAGCAGGTGTTCT‑GGCGGAAGTTCGTCGTCGTGAG‑TTCTATGaaaa                                                                                                                                                                                               <  1:303338/228‑1 (MQ=255)
                                   cgGTGCTTTACAAAGCAGCAGC‑AATTGCAGTAAAATTCCGCACCATTTTGAAATAAGCTGGCGTTGCTCCCAGCGACAAACCGAATTAATCACAGGTGAGAGGCACATGCCGGTAATTAAAGTACGTGAAAACGAACCGTTCGCCGGAGCGCTGCGTCGATTCACGAGTTCCTCCGAAAAATCCGGTGTGCT‑GGCGGAGGATCGTGGT‑GGGAGTTTCGATGAAAAAcc                                                                                                                                                                                            >  2:419631/1‑228 (MQ=255)
                                        ttttAAACAGCAGCAGC‑AATTGCAGTAAAAACCCGCAACATTTAGAAATAAGCGGGGTTTTATGCGAAGGACAAACGCAATAAATCAAAGGGGGGAGGCTACTGCCGGTAATTAAAGTACGTGAAAACGAGCCGTTCGACGTAGCTCTGCGTCGCTTCAAGCGTTCCTGCGAAAAAGCAGGTGTTCT‑GGCGGAAGTTCGTCGTCGTGAG‑TTCTATGAAAAACCGACTa                                                                                                                                                                                       <  2:480795/227‑1 (MQ=255)
                                        cTTTACAAAGCAGCAGC‑AATTGCAGTAAAATTCCGCACCATTTTGAAATAAGCTGGCGTTGATGCCAGCGGCAAACCGAATTAATCAAAGGTGAGAGGCACATGCCGGTAATTAAAGTACGTGAAAACGAGCCGTTCGACGTAGCTCTGCGTCGCTTCAAGCGTTCCTGCGAAAAAGCAGGGGTTCT‑GGCGGAAGTTCGTCGTCGTGAG‑TTCTATGAAAAAccactta                                                                                                                                                                                       >  1:243325/1‑223 (MQ=255)
                                                agcagcagc‑AATTGCAGTAAAATTCCGCACCATTTTGAAATAAGCTGGCGTTGATGCCAGCGGCAAACCGAATTAATCAAAGGTGAGAGGCACATGCCGGTAATTAAAGTACGTGAAAACAAGCCGTTCGACGTAGCTCTGCGTCGCTTCAAGCGTTCCTGCGAAAAAGCAGGTGTTCT‑GGCGGAAGTTCGTCGTCGTGAG‑TTCTATGAAAAACCGACTACCGCACGt                                                                                                                                                                               >  1:508117/1‑228 (MQ=255)
                                                  cagcagc‑aATTGCAGTAAAATTCCGCACCATTTTGAAATAAGCTGGCGTTGATGCCAGCGGCAAACCGAATTAATCAAAGGTGAGAGGCACATGCCGGTAATTAAAGTACGTGAAAACGAGCCGTTCGACGTAGCTCTGCGTCGCTTCAAGCGTTCCTGCGAAAAAGCAGGTGTTCT‑GACGGAAGTTCGTCGCCGTGTG‑TTCTTTGAAAAACGGACCACCGAAccaa                                                                                                                                                                              >  1:23126/1‑224 (MQ=255)
                                                  cagcagc‑aATTGCAGTAAAATTCCGCACCATTTTGAAATAAGCTGGCGTTGATGCCAGCGGCAAACCGAATTAATCAAAGGTGAGAGGCACATGCCGGTAATTAAAGTACGTGAAAACGAGCCGTTCGACGTAGCTCTGCGTCGCTTCAAGCGTTACTGCGAAAAAGCAGGAGTTCT‑AGCGGAACTTCGTCGTCGGGAG‑TTCTATGAAAAACCGCCTACCGAACATaa                                                                                                                                                                             >  1:363051/1‑228 (MQ=255)
                                                          aaTTGCAGTAAAATTCCGCACCATTTTGAAAAAAGCGGGCGTTGATGCCAGCGGCAAACCGAATTAATCAAAGGTGAGAGGCACATGCCGTTAATTAAAGTACGTGAAAACGAGCCGTTCGACGTAGCTCTGCGTCGCTTCAAGCGTTCCTGCGAAAAAGCAGGTGTTCT‑GGCGGAAGTTCGTCGTCGTGAG‑TTCTATGAAAAACCGACTACCGAACGTAAGCGCGCt                                                                                                                                                                      <  1:372063/228‑1 (MQ=255)
                                                           atggCAATAAAATTCCGAACAATTTTGAAATAAGCGGGCGTTGACGCCAACGGCAAACCGAATTAATCAAAGGTGAGAGGCACATGCCGGTAATTAAAGTACGTGAAAACGAGCCGTTCGACGTAGCTCTGCGTCGCTTCAAGCGTTCCTGCGAAAAAGCAGGTGTTCT‑GGCGGAAGTTCGTCGTCGTGAG‑TTCTATGAAAAACCGACTACCGAACGTAAGCGCGCTa                                                                                                                                                                     <  1:652716/225‑1 (MQ=255)
                                                            ttGCAGTAAAATTCCGCACCATTTTGAAATAAGCGTGCGTTGATGCCAGCGGCAAACCGAATTAATCAAAGGTGAGAGGCACATGCCGGTAATTAAAGTACGTGAAAACGAGCCGTTCGACGTAGCTCTGCGTCGCTTCAAGCGTTCCTGCGAAAAAGCAGGTGTTCT‑GGCGGAAGTTCGTCGTCGTGAG‑TTCTATGAAAAACCGACTACCGAACGTAAGCGCGCTaa                                                                                                                                                                    <  1:148238/228‑1 (MQ=255)
                                                                     aaTTCCGCACCATTTTGAAATAAGCTGCCGTTGATGCCAGCGGCAAACCGAATTAATCAAAGGTGAGAGGCACATGCCGGTAATTAAAGTACGTGAAAACGAGCCGTGCGACGTAGCTCTGCGTCGCTTCAAGCGTTCCTGCGAAAAAGCAGGTGTTCT‑GGCGGAAGTTCGTCGTCGTGAG‑TTCTATGAAAAACCGACTACCGAACGTAAGCGCGCTAAAGCTTCTg                                                                                                                                                            <  1:296387/227‑1 (MQ=255)
                                                                            cACCATTTTGAAATAAGCTGGCGTTGATGCCAGCGGCAAACCGAATTAATCAAAGGTGAGAGGCACATGCCGGTAATTAAAGTACGTGAAAACGAGCCGTTCGACGTAGCTCTGCGTCGCTTCAAGCGTTCCTGCGAAAAAGCAGGTGTTCT‑GGCGGAAGTCCGTCGTCGTTAG‑TTCTATGAAAAACCGACTACCGAACGAAAGCGCGCTAAAGCTTCTgacgggaac                                                                                                                                                    >  1:478544/1‑220 (MQ=255)
                                                                               aaTTTGGAAATAAGGTGGCTGGGAAGCCAGCGGAAAACCGAAATAAACAAGGGGGAGGGGCACATGCCGGTAATTAAAGTACGTGAAAACGAGCCGTTCGACGTAGCTCTGCGTCGCTTCAAGCGTTCCTGCGAAAAAGCAGGTGTTCT‑GGCGGAAGTTCGTCGTCGTGAG‑TTCTATGAAAAACCGACTACCGAACGTAAGCGCGCTAAAGCTTCTGCAGTGAAACGt                                                                                                                                                 <  2:90696/227‑1 (MQ=255)
                                                                                 ttttGAAATAAGCTGGCGTTGATGCCAGCGGCAAACCGAATTAATCAAAGGTGAGAGGCACATGCCGGTAATTAAAGGACGTGAAAACGAGCCGTTCGACGTAGCTCTGCGTCGCTTCAAGCGTTCCTGCGAAAAAGCAGGTGTTCT‑GGCGGAAGTTCGTCGTCGTGAG‑TTCTATGAAAAACCGACTACCGAACGTAAGCGCGCGAAAGCTTCCGCAGGCAAACGTCa                                                                                                                                               >  1:379443/1‑228 (MQ=255)
                                                                                           aGCTGGCGTTGATGCCAGCGGCAAACCGAATTAATCAAAGGTGAGAGGCACATGCCGGTAATTAAAGTACGTGAAAACGAGCCGTTCGACGTAGCTCTGCGTCGCTTCAAGCGTTCCTGCGAAAAAGCAGGTGTTCT‑GGCGGAAGTTCGTCGTCGTGAG‑TTCTATGAAAAACCGACTACCGAACGTaa                                                                                                                                                                             >  1:5268/1‑188 (MQ=255)
                                                                                                 ggtgagtGCCAGCGGAAAACACGATAAATCAAGGGGGAGGGGCACTGGCCGGTAATTAAAGTAGGTGAAAACGAGCCGTTCGACGTAGCTCTGCGTCGCTTCAAGCGTTCCTGCGAAAAAGCAGGGGTTCT‑GGCGGAAGTTCGTCGTCGTGAG‑TTCTATGAAAAACCGACTACCGAACGTAAGCGCGCTAAAGCTTCTGCAGTGAAACGTCACGCGAAGAAACTGGCt                                                                                                                               <  2:380262/222‑1 (MQ=255)
                                                                                                                       aaTTAATCAAAGGTGAGAGGCACATGCCGGTAATTAAAGTACGTGAAAACGAGCCGTTCGACGTCGCTCTGCGTCGCTTCAAGCGTTCCTGCGAAAAAGCAGGGGTTCT‑GGCGGACGTTCGTCGTCGTGAG‑TTCTATGAAAAACCGACAACCGAACGAAAGCGCGCTAAAGCTTCTTAAGTGAAACTACACGCGAAGAAACAGGCTCGCCGAAAACACA‑GCCGCAccc                                                                                                        >  2:618530/1‑226 (MQ=255)
                                                                                                                       aaTTAATCAAAGGTGAGAGGCACATGCCGGTAATTAAAGTACGTGAAAACGAGCCGTTCGACGTAGCTCTGCGTCGCTTCAAGCGTTCCTGCGAAAAAGCAGGTGTTCT‑GGCGGAAGTTCGTCGTCGTGAG‑TTCTATGAAAAACCGACTACCGAACGTAAg                                                                                                                                                                            >  1:570531/1‑161 (MQ=255)
                                                                                                                        atgaaTCAAAGGTGAGAGGCACATGCCGGTAATTAAAGTACGTGAAAACGAGCCGTTCGACGTAGCTCTGCGTCGCTTCAAGCGTTCCTGCGAAAAAGCAGGTGTTCT‑GGCGGAAGTTCGTCGTCGTGAG‑TTCTATGAAAAACCGACTACCGAACGTAAgc                                                                                                                                                                           <  2:570531/158‑1 (MQ=255)
                                                                                                                         ttaacaaaaGGTGGAGGGCACATGCCGGTAATTAAAGTATGTGAAAACGAGCCGTTCGACGTAGCTCTGCGTAGCTTCAAGCGTTCCTGCGAAAAAGCAGGTGTTCT‑GGCGGAAGTTCGTCGTCGTGAG‑TTCTATGAAAAACCGACTACCGAACGTAAGCGCGCTAAAGCTTCTGCAGTGAAACGTCACGCGAAGAAACTGGCTCG‑CGAAAACGCACGCCGCACTCGt                                                                                                      <  2:385472/222‑1 (MQ=255)
                                                                                                                           aaTCAAGGGGGAGAGCGACATGCGGGTAATTAAAGTAGTTGAAAAGGAGCGTGTGGACGTAGGTCTGGGGCGCTTCAAGGTTTCCTGCGAAAAAGCAGGTGTTCT‑GGCGGAAGTTCGTCGTCGTGAG‑TTCTATGAAAAACCGACTACCGAACGTAAGCGCGCTAAAGCTTCTGCAGTGAAACGTCACGCGAAGAAACTGGCTCG‑CGAAAACGCACGCCGCACTCGTc                                                                                                     <  2:285829/227‑1 (MQ=255)
                                                                                                                                   ggaagagGCACATGCCGGTAATTAAAGTACGTGAAAACGAGCCGTTCGACGTAGCTCTGCGTCGCTTCAAGCGTTCCTGCGAAAAAGCAGGTGTTCT‑GGCGGAAGTTCGTCGTCGTGAG‑TTCTATGAAAAACCGACTACCGAACGTAAGCGCGCTAAAGCTTCTGCAGTGAAACGTCACGCGAAGAAACTGGCTCG‑CGAAAACGCACGCCGCACTCGTCTGTACTaa                                                                                             <  1:105912/224‑1 (MQ=255)
                                                                                                                                                      aaTTAAAGGTCGGTAAAACGGGCCGTTCGGCGTAGCCCTGCGTTCCTGCAAGCGTTCCTTCGAAAAAGCAGGGGTTCT‑GGCGGGAGTTCGTCGTCGTGAG‑TTCTATGAAAAACCGACTACCGAACGTAAGCGCGCTAAAGCTTCTGCAGTGAAACGTCACGCGAAGAAACTGGCTCG‑CGAAAACGCACGCCGCACTCGTCTGTACTAATTCCCCGAGAGCGTTGctc                                                                          <  2:430552/227‑1 (MQ=255)
                                                                                                                                                            aGTACGTGAAAACGAGCCGTTCGACGTAGCTCTGCGTCGCTTCAAGCGTTCCTGCGAAAAAGCAGGTGTTCT‑GGCGGAAGTTCGTCGTCGTGAG‑TTCTATGAAAAACCGACTACCGAACGTAAGCGCGCTAAAGCTTCTGCAGTGAAACGTCACGCGAAGAAACTGGCTCG‑CGAAAACGCACGCCGCACTCGTCTGTACTAATTCCCCCAGAGCGTTGCTCCCCgcgc                                                                   >  1:514340/1‑225 (MQ=255)
                                                                                                                                                                cgggAAAACGAGCCGTTCGACGTAGCTTTGCGTCCCTTCAAGCGTTCCTGCGAAAAAGCAGGTGTTCT‑GGCGGAAGTTCGTCGTCGTGAG‑TTCTATGAAAAACCGACTACCGAACGTAAGCGCGCTAAAGCTTCTGCAGTGAAACGTCACGCGAAGAAACTGGCTCG‑CGAAAACGCACGCCGCACTCGTCTGTACTAATTCCCCGAGAGCGTTGCTCTCCGATCAGa                                                                <  1:375476/224‑1 (MQ=255)
                                                                                                                                                                cGTGAAAACGAGCCGTTCGACGTAGCTCTGCGTCGCTTCAAGCGTTCCTGCGAAAAAGCAGGTGTTCT‑GGCGGAAGTTCGTCGTCGTGAG‑TTCTATGAAAAACCGACTACCGAACGTAAGCGCGCTAAAGCTTCTGCAGTGAAACGTCACGCGACGAAACTGGCTCG‑CGAAAACGCACGCCCCCCTCGCCTGCCGAAATTACCCCAGAACGTTGCCCTCCGATCaacc                                                               >  2:514247/1‑225 (MQ=255)
                                                                                                                                                                cGTGAAAACGAGCCGTTCGACGTAGCTCTGCGTCGCTTCAAGCGTTCCTGCGAAAAAGCAGGTGTTCT‑GGCGGAAGTTCGTCGTCGTGAG‑TTCTATGAAAAACCGACTACCGAACGTAAGCGCGCTAAAGCTTCTGCAGTGAAACGTCACGCGAAGAACCTGGCTCG‑CAAAAACCCACGCCGCACCCGTCCGTCCTAATTCCCCCACAGCGGTACCCTCCgcacacac                                                               >  2:255865/1‑221 (MQ=255)
                                                                                                                                                                cGTGAAAACGAGCCGTTCGACGTAGCTCTGCGTCGCTTCAAGCGTTCCTGCGAAAAAGCAGGTGTTCT‑GGCGGAAGTTCGTCGTCGTGAG‑TTCTATGAAAAACCGACTACCGAACGTAAGCGCGCTAAAGCTTCTGCAGTGAAACGTCACGCGAAGAAACTGGCTCG‑CGAAAACGCACGCCGCACTCGTCTGTACTAATTCCCCGAGAGCGTTGCTCTCCGATCAGAc                                                               >  1:234487/1‑228 (MQ=255)
                                                                                                                                                                   gAAAACGAGCCGTTCGACGTAGCTCTGCGTCGCTTCAAGCGTTCCTGCGAAAAAGCAGGTGTTCT‑GGCGGAAGTTCGTCGTCGTGAG‑TTCTATGAAAAACCGACTGCCGAACGTAAGCGCGCTAAAGCTTCTGCAGTGAAACGTCACGCGAAGAAACTGGCTCG‑CGAAAACGCAAGCCGCACTCGTCTGCAATAATTCCCCGAGAGGGTTGCTCTCCGATCAGACCgg                                                            >  1:37881/1‑227 (MQ=255)
                                                                                                                                                                   gAAAACGAGCCGTTCGACGTAGCTCTGCGTCGCTTCAAGCGTTCCTGCGAAAAAGCAGGTGTTCT‑GGCGGAAGTTCGTCGTCGTGAG‑TTCTATGAAAAACCGACTACCGAACGTAAGCGCGCTAAAGCt                                                                                                                                                                >  2:135462/1‑129 (MQ=255)
                                                                                                                                                                    aaaaCGAGCCGTTCGACGTAGCTCTGCGTCGCTTCAAGCGTTCCTGCGAAAAAGCAGGTGTTCT‑GGCGGAAGTTCGTCGTCGTGAG‑TTCTATGAAAAACCGACTACCGAACGTAAGCGCGCTAAAGCtt                                                                                                                                                               <  1:135462/129‑1 (MQ=255)
                                                                                                                                                                          atccGTTCGACGTAGCGCTGCTTCGCTGCAAGCGTTCCTGCGAAAAAGCAGGTGTTAT‑GGCGGAAGTTCGTCGTCGTGAT‑TTCTATGAAAAACCGACTACCGAACGTAAGCGCGCTAAAGCTTCTGCAGTGAAACGTCACGCGAAGAAACTGGCTCG‑CGAAAACGCACGCCGCACTCGTCTGTACTAATTCCCCGAGAGCGTTGCTCTCCGATCAGACCGagttgt                                                       <  1:412608/224‑1 (MQ=255)
                                                                                                                                                                              gTTTGACTTTGCTCTGCGCTGCTTCAAGCGTTCCTGCAAAAAAGCAGGTGTTCT‑GGCGTAAGTTCGTCGTCGTGAG‑TTCTATGAAAACCCGACTACCGAACGTAAGCGCGCTAAAGCTTCTGCAGTGAAACGTCACGCGAAGAAACTGGCTCG‑CGAAAACGCACGCCGCACTCGTCTGTACTAATTCCCCGAGAGCGTTGCTCTCCGATCAGACCGagttgtagttg                                                  <  2:150760/227‑1 (MQ=255)
                                                                                                                                                                                        ggttggCGTCCCTTCAAGCGTTCCTGCGAAACAGCAGGTGTTCT‑GGCGGAAGTCCGTCGTCTTGAG‑TTCGAAGAAAACCCGACTACCGAACGTAAGCGCGCTAAAGCTTCTGCAGTGAAACGTCACGCGAAGAAACTGGCTCG‑CGAAAACGCACGCCGCACTCGTCTGTACTAATTCCCCGAGAGCGTTGCTCTCCGATCAGACCGAGTTGTAGTTGTAAGGCCGTg                                        <  2:598827/222‑1 (MQ=255)
                                                                                                                                                                                         ctctGCGTCGCTTCAAGCGTTCCTGCGAAAAAGCAGGTGTTCT‑GGCGGAAGTTCGTCGTCGTGAG‑TTCTATGAAAAACCGACTACCGAACGTAAGCGCGCTAAAGCTTCTGCAGTGAAACGTCACGCGAAGAAACTGGCTCG‑CGAAAACGCACGCCGCACTCGTCTGTACTAATTCCCCGAGAGCGTTGCTCTCCGGTCAGACCGAGTTGTAGTTGTAAGGcctgggt                                      >  1:408866/1‑223 (MQ=255)
                                                                                                                                                                                         ctctGCGTCGCTTCAAGCGTTCCTGCGAAAAAGCAGGTGTTCT‑GGCGGAAGTTCGTCGTCGTGAG‑TTCTATGAAAAACCGACTACCGAACGTAAGCGCGCTAAAGCTTCTGCAGTGAAACGCCACGCGAAGAAACTGGCTCG‑CGAAAACGCACGCCGCAATCGTCTGTACTAATTCCCCGACAGCGTTGCTCTCCGATCAGACCGCGTTGTAGTTGTAAGGcctggct                                      >  1:79138/1‑223 (MQ=255)
                                                                                                                                                                                         ctctGCGTCGCTTCAAGCGTTCCTGCGAAAAAGCAGGGGTTCT‑GGCGGAAGTTCGTCGTCGTGAG‑TTCTATGAAAAACCGACTACCGAACGTAAGCGCGCTAAAGCTTCTGCAGTGAAACGTCACGCGAAGAAACAGGCTCG‑CGAAAACGCACGCCGCCCGCGGCTGTGCTAATTCCCCGAGAGCGGTTCTCCCCGGCCAAACCCAGTTTTAGTTGTAAGGCCGGGcg                                      >  2:419657/1‑227 (MQ=255)
                                                                                                                                                                                             gCGTCGCTTGAAGCGTTCGTGAGAAAAAGCAGGTGTTCT‑GGCGGAAGTTCGTCGTCGTGAG‑TTCTAAGAAAAACCGACTACCGAACGTAAGCGCGCTAAAGCTTCTGCAGTGAAACGTCACGCGAAGAAACTGGCTCG‑CGAAAACGCACGCCGCACTCGTCTGTACTAATTCCCCGAGAGCGTTGCTCTCCGATCAGACCGAGTTGTAGTTGTAAGGCCGTg                                        <  1:264454/222‑1 (MQ=255)
                                                                                                                                                                                              ggTCGTTTCAAGCGTTCCTGCAAAAAAGCAGGTTTTCC‑GGCGGAAGGTCGGCGTCGTGAG‑TTCTAGGAAAAACCGACTACCGAACGTAAGCGCGCTAAAGCTTCTGCAGTGAAACGTCACGCGAAGAAACTGGCTCG‑CGAAAACGCACGCCGCACTCGTCTGTACTAATTCCCCGAGAGCGTTGCTCTCCGATCAGACCGAGTTGTAGTTGTAAGGCCGTGCTTCCGa                                 <  2:408866/227‑1 (MQ=255)
                                                                                                                                                                                               gTCGCTTCAAGCGTTCCTGCGAAAGAGCAGGGATGCT‑GGCGGAAGTTCGTCGTCGTAAG‑TTCTAGGAAAAACCGACTACCGAACGTAAGCGCGCTAAAGCTTCTGCAGTAAAACGTCACGCGAAGAAACTGGCTCG‑CGAAAACGCACGCCGCACTCGTCTGTACTAATTCCCCGAGAGCGTTGCTCTCCGATCAGACCGAGTTGTAGTTGTAAGGCCGTGCTTCCGa                                 <  1:248100/227‑1 (MQ=255)
                                                                                                                                                                                               gTCGCTTCAAGCGTTCCTGCGAAAAAGCAGGTGTTCT‑GGCGGAAGTTCGTCGTCGTGAG‑TTCTATGAAAAACCGACTACCGAACGTAAGCGCGCTAAAGCTTCTGCAGTGAAACGTCACGCGAAGAAACTGGCTCG‑CGAAAACGCACGCCGCACTCGTCTGTACCAATTCCCCCAGAGCGTGGCTCCCCGACCCGCCCCGGTTGTAGTTGTAAGGCCGAGCTTCCGaa                                >  2:41038/1‑228 (MQ=255)
                                                                                                                                                                                                   cTTCAAGCGTTCCTGCGAAAAAGCAGGTGTTCT‑GGCGGAAGTTCGTCGTCGTGAG‑TTCTATGAAAAACCGACTACCGAACGTAAGCGCGCTAAAGCTTCTGCAGTGAAACGTCACGCGAAGAAACTGGCTCG‑CGAAAACGCACGCCGCACTCGTCTGTACTAATTCCCCGAGAGCGTTGCTCTCCGATCAGACCGAGTTGTGGTTGTAAGGCCGTGCTTCCGAAAGGa                            >  1:434442/1‑228 (MQ=255)
                                                                                                                                                                                                        atggTTCCAGCGAAAAAGC‑GGTGTTTTCGGCGGAAGTTCGGCGGAGTGAG‑TTCTATGAAAAACCGACTACCGAACTTAAGCCCGCTAAAGCTTCTGCAGTGAAACGGCACGCGAAGAAACTGGCTCG‑CGAAAACGCACGCCGCACTCGTCTGTACTAATTCCCCGAGAGCGTTGCTCTCCGATCAGACCGAGTTGTAGTTGTAAGGCCGTGCTTCCGAAAGGAATgcg                       <  2:400228/225‑1 (MQ=255)
                                                                                                                                                                                                                 gCGAAAAAGCAGGTGTTCT‑GGCGGAAGTTCGTCGTCGTGAG‑TTCTATGAAAAACCGACTACCGAACGTAAGCGCGCTAAAGCTTCTGCAGTGAAACGTCACGCGAAGAAACTGGCTCG‑CGAAAACGCACGCCGCACTCGTCTGTACTCATTCCCCGAGAGCGTTCCTCTCCGATCGGACCGGGTTGTAGTTGTAAGGCCGTGCTCCCGAAAGGAACTCGCGGCTTAtt              >  1:558931/1‑228 (MQ=255)
                                                                                                                                                                                                                             gtgtTCT‑GGGGGAAGGTCGTCGTGGTGAG‑TTCTATGAAAAACCGACTACCGAACGTAAGCGCGCTAAAGCTTCTGCAGGGAAACGTCACGCGAAGAAACTGGCTCG‑CGAAAACGCACGCCGCACTCGTCTGTACTAATTCCCCGAGAGCGTTGCTCTCCGATCAGACCGAGTTGTAGTTGTAAGGCCGTGCTTCCGAAAGGAATGCGCGGCTTATTTTCGTTTATGaa  <  2:422301/228‑1 (MQ=255)
                                                                                                                                                                                                                              |                                                                                                                                                                                                                                     
CTCACCCTTATAAAAGTCCCTTTCAAAAAAGGCCGCGGTGCTTTACAAAGCAGCAGC‑AATTGCAGTAAAATTCCGCACCATTTTGAAATAAGCTGGCGTTGATGCCAGCGGCAAACCGAATTAATCAAAGGTGAGAGGCACATGCCGGTAATTAAAGTACGTGAAAACGAGCCGTTCGACGTAGCTCTGCGTCGCTTCAAGCGTTCCTGCGAAAAAGCAGGTGTTCT‑GGCGGAAGTTCGTCGTCGTGAG‑TTCTATGAAAAACCGACTACCGAACGTAAGCGCGCTAAAGCTTCTGCAGTGAAACGTCACGCGAAGAAACTGGCTCG‑CGAAAACGCACGCCGCACTCGTCTGTACTAATTCCCCGAGAGCGTTGCTCTCCGATCAGACCGAGTTGTAGTTGTAAGGCCGTGCTTCCGAAAGGAATGCGCGGCTTATTTTCGTTTATGAA  >  NC_000913/3210640‑3211087

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: ‑