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→GGGrpsU30S 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

CTCACCCTTATAAAAGTCCCTTTCAAAAAAGGCCGCGGTGCTTTACAAAGCAGCAGCAATTGCAGTAAAATTCCGCACCATTTTGAAATAAGCTGGCGTTGATGCCAGCGGCAAACCGAATTAATCAAAGGTGAGAGGCACATGCCGGTAATTAAAGTACGTGAAAACGAGCCGTTCGACGTAGCTCTGCGTCGCTTCAAGCGTTCCTGCGAAAAAGCAGGTGTTCTGGCGGAAGTTCGTCGTCGTGAGTTCTATGAAAAACCGACTACCGAACGTAAGCGCGCTAAAGCTTCTGCAGTGAAACGTCACGCGAAGAAACTGGCTCGCGAAAACGCACGCCGCACTCGTCTGTACTAATTCCCCGAGAGCGTTGCTCTCCGATCAGACCGAGTTGTAGTTGTAAGGCCGTGCTTCCGAAAGGAATGCGCGGCTTATTTTCGTTTATGAA  >  NC_000913/3210640‑3211087
                                                                                                                                                                                                                              |                                                                                                                                                                                                                                     
ctcACCCTTATAAAAGTCCCTTTCAAAAAAGGCCGCGGTGCTTTACAAAGCAGCAGCAATTGCAGTAAAATTCCGCACCATTTTGAAATAAGCTGGCGTTGATGCCAGCGACAAACCGAATTAATCAAAGGTGAGAGGCACATGCCGGTAATTAAAGTACGTGAAAACAATCCGTTAGACGTAACTCTGAGTCGCTTCAAGCGTTCCTGCGAACAAACAGGGGTTCTg                                                                                                                                                                                                                                >  2:61093/1‑228 (MQ=255)
 tcACTCTTATAAAAGTCCCTTTCAAAAAAGGCCGCGGTGCTTTACAAAGCAGCAGCAATTGCAGTAAAATTCCGCACCATTTTGAAATAAGCTGGCGTTGATGCCAGCGGCAAACCGAATTAATCAAAGGTAAACGGCCCATGCCGGTAATTAAAACACCTGAAACCAAACCGTTCGACCAAGCCCTGCGGTCCTTCAAGCCTGCCCGCGAAAAAACCGGGGtttgg                                                                                                                                                                                                                                 >  2:694766/1‑224 (MQ=255)
  cACCCTTATAAAAGTCCCTTTCAAAAAAGGCCGCGGTGCTTTACAAAGCAGCAGCAATTGCAGTAAAATTCCGCACCATTTTGAAATAAGCTGGCGTTGATGCCAGCGGCAAACCGAATTAATCAAAGGTTAGAGTCACATGCCGGTAATTAAAGTACGTGAAAACAAGCCGTTAGACGTAACTCTGCGTAGCTTCAAGCTTTCCGCGGCAAAAGAAGGTGTTTTgcg                                                                                                                                                                                                                              >  1:642754/1‑226 (MQ=255)
   aCCCTTATAAAAGTCCCTTTCAAAAAAGGCCGCGGTGCTTTACAAAGCAGCAGCAATTGCAGTAAAATTCCGCACCATTTTGAAATAAGCTGGCGTTGCTGCCAGCGGCAAACCGAATTAATCAAAGGTGAGAGGCAAATGCCGGGAATTAAGGTACGAGAAACCGAGCCGTTCGACGCAGCTCGCCGGCGCTTCAGACGTTCCTTCGTAAAACCGGGTGTTTTGCCg                                                                                                                                                                                                                             >  1:391539/1‑228 (MQ=255)
   aCCCTTATAAAAGTCCCTTTCAAAAAAGGCCGCGGTGCTTTACAAAGCAGCAGCAATTGCAGTAAAATTCCGCACCATTTTGAAATAAGCTGGCGTTGATGCCAGCGGCAAACCGAATTAATCAAAGGTGAGAGGCACATGCCGGTAATTAAAGTACGTGAAAACGAGCCGTTCGACGTAGCTCTGCGTCGCTTCAAGCGTTCCTGCGAAAAAGCAGGTGTTCTGGc                                                                                                                                                                                                                              >  1:422301/1‑227 (MQ=255)
   aCCCTTATAAAAGTCCCTTTCAAAAAAGGCCGCGGTGCTTTACAAAGCAGCAGCAATTGCAGTAAAATTCCGCACCATTTTGAAATAAGCTGGCGTTGATGCCAGCGGCAAACCGAATTAATCAAAGGTGAGAGGCACATGCCGGTAATTAAAGTACGTGAAAACAAGCCGTTCGACGTAGCTCTGCGTCGCTTCAAGCGTTCCTGCGAAAAAGCAGGGGTTCTgggg                                                                                                                                                                                                                             >  2:489255/1‑226 (MQ=255)
        tataAAAGTCCCTTTCAAAAAAGGCCGCGGTGCTTTACAAAGCAGCAGCAATTGCAGTAAAATTCCGCACCATTTTGAAATAAGCTGGCGTTGATGCCAGCGGCAAACCGAATTAATCAAAGGTGAGAGGCACATGCCGGTAATTAAAGTACGTGAAAACGAGCCGTTCGACGTAGCTCTGCGTCGCTTCAAGCGTTCCTGCGAAAAAACAGGGGTTATGGCGGCAgg                                                                                                                                                                                                                        >  2:375476/1‑227 (MQ=255)
        tataAAAGTCCCTTTCAAAAAAGGCCGCGGTGCTTTACAAAGCAGCAGCAATTGCAGTAAAATTCCGCACCATTTTGAAATAAGCTGGCGTTGATGCCAGCGGCAAACCGAATTAATCAAAGGTGAGAGGCACATGCCGGTAATTAAAGTACGTGAAAACGAGCCGTTCGACGAAGCTCTGCGTCGCTTCAAGCGTTCCTGCGAAAAAGCAGGTGTGCTGGCGGAAGt                                                                                                                                                                                                                        >  1:427007/1‑228 (MQ=255)
        tataAAAGTCCCTTTCAAAAAAGGCCGCGGTGCTTTACAAAGCAGCAGCAATTGCAGTAAAATTCCGCACCATTTTGAAATAAGCTGGCGTTGATGCCAGCGGCAAACCGAATTAATCAAAGGTGAGAGGCACATGCCGGTAATTAAAGTACGTGAAAACGAACCGTTCGACTTAGCTCTGCGTCCCTTCAACCGTTCCTGCGAAAAACAAGGGGTTCTGGCGGAAGt                                                                                                                                                                                                                        >  2:653326/1‑228 (MQ=255)
        tataAAAGTCCCTTTCAAAAAAGGCCGCGGTGCTTTACAAAGCAGCAGCAATTGCAGTAAAATTCCGCACCATTTTGAAATAAGCTGGCGTTGATGCCAGCGGCAAACCGAATTAATCAAAGGTGAGAGGCACATGCCGGTAATTAAAGTACGGGAAAACAAGCCGTTCAACGAAGCCCTGGGGCGCTTCAGGCGTTCCTGCCAAAAAGCAGGGGTTCGGGGGGAAGt                                                                                                                                                                                                                        >  2:343395/1‑228 (MQ=255)
         ataAAAGTCCCTTTCAAAAAAGGCCGCGGTGCTTTACAAAGCAGCAGCAATTGCAGTAAAATTCCGCACCATTTTGAAATAAGCTGGCGTTGATGCCAGCGGCAAACCGAATTAATCAAAGGTAAGAGGCACATGCCGGTAATTAAAGTACGCGAAAACGAGCCGTTCGACGTAGCTCTGCATCGCTTCAAGCGTTTCTGCTGAAAAGAAGGTGTTCTGCCGCAAGtt                                                                                                                                                                                                                       >  1:234339/1‑228 (MQ=255)
          taAAAGTCCCTTTCAAAAAAGGCCGCGGTGCTTTACAAAGCAGCAGCAATTGCAGTAAAATTCCGCACCATTTTGAAATAAGCTGGCGTTGATGCCAGCGGCAAACCGAATTAATCAAAGGTGAGAGGCCCATGCCGGTAATTAAAGTACGTGAAAACGAGCCGTTCGCCGTAGCTCTGCGTCGCTTCAAGCGTTCATCCGAAAAAGCAGGGGTTCTGCCGGAAGTtc                                                                                                                                                                                                                      >  1:263820/1‑228 (MQ=255)
          taAAAGTCCCTTTCAAAAAAGGCCGCGGTGCTTTACAAAGCAGCAGCAATTGCAGTAAAATTCCGCACCATTTTGAAATAAGCTGGCGTTGATGCCAGCGGCAAACCGAATTAATCAAAGGTGAGAGGCACATGCCGGTAATTAAAGTACGTGAAAACGAGCCGTTCGACGTAGCTCTGCGTCGCTTCAAGCGTTCCTGCGAAAAAGCAGGGGTTCTGGCGGAAGCtc                                                                                                                                                                                                                      >  1:385472/1‑228 (MQ=255)
           aaaaATCCATTTAAAAAAAGGCCGCGGGGTTTTAAAAAGCGGAGCAAATGGCAAGAAAAGTCCGGAACAATTGGAAAAAAGCTGGGGTTGATGCAAGCGGAAAACCGAATAAATCAAAGGTGAGAGGCACATGCCGGTAATTAAAGTACGTGAAAACGAGCCGTTCGACGTAGCTCTGCGTCGCTTCAAGCGTTCCTGCGAAAAAGCAGGTGTTCTGGCGGAAGTtcg                                                                                                                                                                                                                     <  2:2848/228‑1 (MQ=255)
             aaGTCCCTTACAAAAACGGCCGCGGTGCTTTACAAAGCAGCAGCAATTGACGTAAAATTCCGCCCCATTTTGAAATAAGCTGGCGTTGATGCCAGCGGCAAACCGAATTAATCAAAGGTGAGAGGCACATGCCGGTAATTAAAGTACGTGAAAACGAGCCGTTCGACGTAGCTCTGCGTCGCTTCAAGCGTTCCTGCGAAAAAGCAGGTGTTCTGGCGGAAGTtcgtc                                                                                                                                                                                                                   <  1:666349/228‑1 (MQ=255)
               gTCCCTTTCAAAAAAGGCCGCGGTGCTTTACAAAGCAGCAGCAATTGCAGTAAAATTCCGCACCATTTTGAATTAACTTGGCGTTGATACCAGCGGCAAACCGAATTAATCAAAGGCGAGAGGCACATGCTGGTAATTAAAGTAAGTGATAACAAGCGGTTCGAAGTAGCTCTTCTTCGCTTCAAACGTTCCTGCGAACAATCAAGAGATCTGAGGGAAGttgccg                                                                                                                                                                                                                   >  1:99570/1‑222 (MQ=255)
               gTCCCTTTCAAAAAAGGCCGCGGTGCTTTACAAAGCAGCAGCAATTGCAGTAAAATTCCGCACCATTTTGAAATAAGCTGGCGTTGATGCCAGCGGCAAACCGAATTAATCAAAGGTGAGAGGCACATGCCGGTAATTAAAGTACGTGAAAACGAGCCGTTCGACGTAGCTCTGCGTCGCTTCAACCGTTCCTGCGAAAAAGCAGTGGTTCTGGCGGAAGTtcctcgt                                                                                                                                                                                                                 >  1:358697/1‑228 (MQ=255)
               gTCCCTTTCAAAAAAGGCCGCGGTGCTTTACAAAGCAGCAGCAATTGCAGTAAAATTCCGCACCATTTTGAAATAAGCTGGCGTTGATGCCAGCGGCAAACCGAATTAATCAAAGGTGAGAGGCACATGCCGGCAATTAAAGAACGTGAAAACGAGCCGTTCAACGTAGCTCTGCGTCCCTTTAAGCGTTACCGCGAAAAAGCAGGGTTTCTGGCGGAAATtcgtcgc                                                                                                                                                                                                                 >  2:662862/1‑227 (MQ=255)
               gTCCCTTTCAAAAAAGGCCGCGGTGCTTTACAAAGCAGCAGCAATTGCAGTAAAATTCCGCACCATTTTGAAATAAGCTGGCGTTGATGCCAGCGGCAAACCGAATTAATCAAAGGTAAGAGGCACATGCCGGCAATTAAAGTACGTGAAAACGAGCCGTTCGACGTAGCTCTGCGTCGCTTCAAGCGTTCCTGCGAAAAAGACGGTGTTCTGGCGGAAGGttgtctt                                                                                                                                                                                                                 >  2:496515/1‑226 (MQ=255)
                  tctgtCAAAAAAGGCCGCGGTGCTTTACAAAGCGGCAGCAATGGCAGTAAAATACCGAAACATTTTGAAATAAGCTGGCGTTGATTCCAGCGGCAAACCGAATTAGTCAAGGGTGAGGGGCACATGCCGGTAATTAAAGTACGTGAAAACGAGCCGTTCGACGTAGCTCTGCGTCGCTTCAAGCGTTCCTGCGAAAAAGCAGGTGTTCTGGCGGAAGTTcgtcgtcgt                                                                                                                                                                                                              <  2:639113/224‑1 (MQ=255)
                  ccTTTCAAAAAAGGCCGCGGTGCTTTACAAAGCAGCAGCAATTGCAGTAAAATTCCGCACCATTTTGAAATAAGCTGGCGTTGATGCCAGCGGCAAACCGAATTAATCAAAGGTTAGAGGCCCATGCCGGTAATTAAAGTACGGTAAAACGAGCCGGTCAACGTAACTCTACGTCGCTTCAATCGTTCCTCCGCAAAAGCAGGGGTCTAGGCGGAAGttcttctttgt                                                                                                                                                                                                              >  2:412608/1‑223 (MQ=255)
                    tttCAAAAAAGGCCGCGGTGCTTTACAAAGCAGCAGCAATTGCAGTAAAATTCCGCACCATTTTGAAATAAGCTGGCGTTGATGCCAGCGGCAAACCGAATTAATCAAAGGTGAGAGGCACATGCCGGTAATTAAAGTACGTGAAAACGACCCGTTCGACGTAGCTCTGCGTCGCTTCAAGCGTTCCTGCGAAAAAGCAGGGTTTCTGCCGGAAGTtcgtcgcccggg                                                                                                                                                                                                            >  2:147943/1‑222 (MQ=255)
                     ttaaaaaatGGCCCAGGTGCATTACAAGGAAGCAGCAAGAGCGGTAAAATCCCACACTTTTTAGAAATAAGCGAGCGTGGATGCGACCGGCAAACCGAATAAGTCAAAGGTGGAAGGGACATGCCGGAAATTAAAGTACGTGAAAACAAGCCGTACGACGTAGCTCTGGGTCGCTTCAATCGTTCCTGCGAAAAAGAAGGTGTTCTGGCGGAAGTTCGTCGTCGTGAg                                                                                                                                                                                                           <  2:99570/225‑1 (MQ=255)
                     ttCAAAAAAGGCCGCGGTGCTTTACAAAGCAGCAGCAATTGCAGTAAAATTCCGCACCATTTTGAAATAAGCTGGCGTTGATGCCAGCGGCAAACCGAATTAATCAAAGGAGAGAGGCACATGCCGGTAATTAAAGTACGTGAAAACTAGCCGTTCGACGTAGCTCTGCGTCGCTTCAAGCGTTACTGCGCAAAAGCAAGGGGTCTGGCGGAAGTTCGTCGTCGTggg                                                                                                                                                                                                           >  1:672154/1‑226 (MQ=255)
                                ccGGGGTGGTTTACAAAGCAGCAGCAAGTGCAGAAAAATTCCGCCCAATTTTGAAATAAGCTGGCGGTTATGCCAGCGGCAAACCAAATTAATCAAAGGTGAGAGGCACATGCCGGTAATTAAAGTACGTGAAACCGAGCCGTTCGACGTAGCTCTGCGTCGCTTCAAGCGTTCCTGCGAAAAAGCAGGTGTTCTGGCGGAAGTTCGTCGTCGTGAGTTCTATGaaaa                                                                                                                                                                                               <  1:303338/228‑1 (MQ=255)
                                   cgGTGCTTTACAAAGCAGCAGCAATTGCAGTAAAATTCCGCACCATTTTGAAATAAGCTGGCGTTGCTCCCAGCGACAAACCGAATTAATCACAGGTGAGAGGCACATGCCGGTAATTAAAGTACGTGAAAACGAACCGTTCGCCGGAGCGCTGCGTCGATTCACGAGTTCCTCCGAAAAATCCGGTGTGCTGGCGGAGGATCGTGGTGGGAGTTTCGATGAAAAAcc                                                                                                                                                                                            >  2:419631/1‑228 (MQ=255)
                                        ttttAAACAGCAGCAGCAATTGCAGTAAAAACCCGCAACATTTAGAAATAAGCGGGGTTTTATGCGAAGGACAAACGCAATAAATCAAAGGGGGGAGGCTACTGCCGGTAATTAAAGTACGTGAAAACGAGCCGTTCGACGTAGCTCTGCGTCGCTTCAAGCGTTCCTGCGAAAAAGCAGGTGTTCTGGCGGAAGTTCGTCGTCGTGAGTTCTATGAAAAACCGACTa                                                                                                                                                                                       <  2:480795/227‑1 (MQ=255)
                                        cTTTACAAAGCAGCAGCAATTGCAGTAAAATTCCGCACCATTTTGAAATAAGCTGGCGTTGATGCCAGCGGCAAACCGAATTAATCAAAGGTGAGAGGCACATGCCGGTAATTAAAGTACGTGAAAACGAGCCGTTCGACGTAGCTCTGCGTCGCTTCAAGCGTTCCTGCGAAAAAGCAGGGGTTCTGGCGGAAGTTCGTCGTCGTGAGTTCTATGAAAAAccactta                                                                                                                                                                                       >  1:243325/1‑223 (MQ=255)
                                                agcagcagcAATTGCAGTAAAATTCCGCACCATTTTGAAATAAGCTGGCGTTGATGCCAGCGGCAAACCGAATTAATCAAAGGTGAGAGGCACATGCCGGTAATTAAAGTACGTGAAAACAAGCCGTTCGACGTAGCTCTGCGTCGCTTCAAGCGTTCCTGCGAAAAAGCAGGTGTTCTGGCGGAAGTTCGTCGTCGTGAGTTCTATGAAAAACCGACTACCGCACGt                                                                                                                                                                               >  1:508117/1‑228 (MQ=255)
                                                  cagcagcaATTGCAGTAAAATTCCGCACCATTTTGAAATAAGCTGGCGTTGATGCCAGCGGCAAACCGAATTAATCAAAGGTGAGAGGCACATGCCGGTAATTAAAGTACGTGAAAACGAGCCGTTCGACGTAGCTCTGCGTCGCTTCAAGCGTTCCTGCGAAAAAGCAGGTGTTCTGACGGAAGTTCGTCGCCGTGTGTTCTTTGAAAAACGGACCACCGAAccaa                                                                                                                                                                              >  1:23126/1‑224 (MQ=255)
                                                  cagcagcaATTGCAGTAAAATTCCGCACCATTTTGAAATAAGCTGGCGTTGATGCCAGCGGCAAACCGAATTAATCAAAGGTGAGAGGCACATGCCGGTAATTAAAGTACGTGAAAACGAGCCGTTCGACGTAGCTCTGCGTCGCTTCAAGCGTTACTGCGAAAAAGCAGGAGTTCTAGCGGAACTTCGTCGTCGGGAGTTCTATGAAAAACCGCCTACCGAACATaa                                                                                                                                                                             >  1:363051/1‑228 (MQ=255)
                                                          aaTTGCAGTAAAATTCCGCACCATTTTGAAAAAAGCGGGCGTTGATGCCAGCGGCAAACCGAATTAATCAAAGGTGAGAGGCACATGCCGTTAATTAAAGTACGTGAAAACGAGCCGTTCGACGTAGCTCTGCGTCGCTTCAAGCGTTCCTGCGAAAAAGCAGGTGTTCTGGCGGAAGTTCGTCGTCGTGAGTTCTATGAAAAACCGACTACCGAACGTAAGCGCGCt                                                                                                                                                                      <  1:372063/228‑1 (MQ=255)
                                                           atggCAATAAAATTCCGAACAATTTTGAAATAAGCGGGCGTTGACGCCAACGGCAAACCGAATTAATCAAAGGTGAGAGGCACATGCCGGTAATTAAAGTACGTGAAAACGAGCCGTTCGACGTAGCTCTGCGTCGCTTCAAGCGTTCCTGCGAAAAAGCAGGTGTTCTGGCGGAAGTTCGTCGTCGTGAGTTCTATGAAAAACCGACTACCGAACGTAAGCGCGCTa                                                                                                                                                                     <  1:652716/225‑1 (MQ=255)
                                                            ttGCAGTAAAATTCCGCACCATTTTGAAATAAGCGTGCGTTGATGCCAGCGGCAAACCGAATTAATCAAAGGTGAGAGGCACATGCCGGTAATTAAAGTACGTGAAAACGAGCCGTTCGACGTAGCTCTGCGTCGCTTCAAGCGTTCCTGCGAAAAAGCAGGTGTTCTGGCGGAAGTTCGTCGTCGTGAGTTCTATGAAAAACCGACTACCGAACGTAAGCGCGCTaa                                                                                                                                                                    <  1:148238/228‑1 (MQ=255)
                                                                     aaTTCCGCACCATTTTGAAATAAGCTGCCGTTGATGCCAGCGGCAAACCGAATTAATCAAAGGTGAGAGGCACATGCCGGTAATTAAAGTACGTGAAAACGAGCCGTGCGACGTAGCTCTGCGTCGCTTCAAGCGTTCCTGCGAAAAAGCAGGTGTTCTGGCGGAAGTTCGTCGTCGTGAGTTCTATGAAAAACCGACTACCGAACGTAAGCGCGCTAAAGCTTCTg                                                                                                                                                            <  1:296387/227‑1 (MQ=255)
                                                                            cACCATTTTGAAATAAGCTGGCGTTGATGCCAGCGGCAAACCGAATTAATCAAAGGTGAGAGGCACATGCCGGTAATTAAAGTACGTGAAAACGAGCCGTTCGACGTAGCTCTGCGTCGCTTCAAGCGTTCCTGCGAAAAAGCAGGTGTTCTGGCGGAAGTCCGTCGTCGTTAGTTCTATGAAAAACCGACTACCGAACGAAAGCGCGCTAAAGCTTCTgacgggaac                                                                                                                                                    >  1:478544/1‑220 (MQ=255)
                                                                               aaTTTGGAAATAAGGTGGCTGGGAAGCCAGCGGAAAACCGAAATAAACAAGGGGGAGGGGCACATGCCGGTAATTAAAGTACGTGAAAACGAGCCGTTCGACGTAGCTCTGCGTCGCTTCAAGCGTTCCTGCGAAAAAGCAGGTGTTCTGGCGGAAGTTCGTCGTCGTGAGTTCTATGAAAAACCGACTACCGAACGTAAGCGCGCTAAAGCTTCTGCAGTGAAACGt                                                                                                                                                 <  2:90696/227‑1 (MQ=255)
                                                                                 ttttGAAATAAGCTGGCGTTGATGCCAGCGGCAAACCGAATTAATCAAAGGTGAGAGGCACATGCCGGTAATTAAAGGACGTGAAAACGAGCCGTTCGACGTAGCTCTGCGTCGCTTCAAGCGTTCCTGCGAAAAAGCAGGTGTTCTGGCGGAAGTTCGTCGTCGTGAGTTCTATGAAAAACCGACTACCGAACGTAAGCGCGCGAAAGCTTCCGCAGGCAAACGTCa                                                                                                                                               >  1:379443/1‑228 (MQ=255)
                                                                                           aGCTGGCGTTGATGCCAGCGGCAAACCGAATTAATCAAAGGTGAGAGGCACATGCCGGTAATTAAAGTACGTGAAAACGAGCCGTTCGACGTAGCTCTGCGTCGCTTCAAGCGTTCCTGCGAAAAAGCAGGTGTTCTGGCGGAAGTTCGTCGTCGTGAGTTCTATGAAAAACCGACTACCGAACGTaa                                                                                                                                                                             >  1:5268/1‑188 (MQ=255)
                                                                                                 ggtgagtGCCAGCGGAAAACACGATAAATCAAGGGGGAGGGGCACTGGCCGGTAATTAAAGTAGGTGAAAACGAGCCGTTCGACGTAGCTCTGCGTCGCTTCAAGCGTTCCTGCGAAAAAGCAGGGGTTCTGGCGGAAGTTCGTCGTCGTGAGTTCTATGAAAAACCGACTACCGAACGTAAGCGCGCTAAAGCTTCTGCAGTGAAACGTCACGCGAAGAAACTGGCt                                                                                                                               <  2:380262/222‑1 (MQ=255)
                                                                                                                       aaTTAATCAAAGGTGAGAGGCACATGCCGGTAATTAAAGTACGTGAAAACGAGCCGTTCGACGTCGCTCTGCGTCGCTTCAAGCGTTCCTGCGAAAAAGCAGGGGTTCTGGCGGACGTTCGTCGTCGTGAGTTCTATGAAAAACCGACAACCGAACGAAAGCGCGCTAAAGCTTCTTAAGTGAAACTACACGCGAAGAAACAGGCTCGCCGAAAACACAGCCGCAccc                                                                                                        >  2:618530/1‑226 (MQ=255)
                                                                                                                       aaTTAATCAAAGGTGAGAGGCACATGCCGGTAATTAAAGTACGTGAAAACGAGCCGTTCGACGTAGCTCTGCGTCGCTTCAAGCGTTCCTGCGAAAAAGCAGGTGTTCTGGCGGAAGTTCGTCGTCGTGAGTTCTATGAAAAACCGACTACCGAACGTAAg                                                                                                                                                                            >  1:570531/1‑161 (MQ=255)
                                                                                                                        atgaaTCAAAGGTGAGAGGCACATGCCGGTAATTAAAGTACGTGAAAACGAGCCGTTCGACGTAGCTCTGCGTCGCTTCAAGCGTTCCTGCGAAAAAGCAGGTGTTCTGGCGGAAGTTCGTCGTCGTGAGTTCTATGAAAAACCGACTACCGAACGTAAgc                                                                                                                                                                           <  2:570531/158‑1 (MQ=255)
                                                                                                                         ttaacaaaaGGTGGAGGGCACATGCCGGTAATTAAAGTATGTGAAAACGAGCCGTTCGACGTAGCTCTGCGTAGCTTCAAGCGTTCCTGCGAAAAAGCAGGTGTTCTGGCGGAAGTTCGTCGTCGTGAGTTCTATGAAAAACCGACTACCGAACGTAAGCGCGCTAAAGCTTCTGCAGTGAAACGTCACGCGAAGAAACTGGCTCGCGAAAACGCACGCCGCACTCGt                                                                                                      <  2:385472/222‑1 (MQ=255)
                                                                                                                           aaTCAAGGGGGAGAGCGACATGCGGGTAATTAAAGTAGTTGAAAAGGAGCGTGTGGACGTAGGTCTGGGGCGCTTCAAGGTTTCCTGCGAAAAAGCAGGTGTTCTGGCGGAAGTTCGTCGTCGTGAGTTCTATGAAAAACCGACTACCGAACGTAAGCGCGCTAAAGCTTCTGCAGTGAAACGTCACGCGAAGAAACTGGCTCGCGAAAACGCACGCCGCACTCGTc                                                                                                     <  2:285829/227‑1 (MQ=255)
                                                                                                                                   ggaagagGCACATGCCGGTAATTAAAGTACGTGAAAACGAGCCGTTCGACGTAGCTCTGCGTCGCTTCAAGCGTTCCTGCGAAAAAGCAGGTGTTCTGGCGGAAGTTCGTCGTCGTGAGTTCTATGAAAAACCGACTACCGAACGTAAGCGCGCTAAAGCTTCTGCAGTGAAACGTCACGCGAAGAAACTGGCTCGCGAAAACGCACGCCGCACTCGTCTGTACTaa                                                                                             <  1:105912/224‑1 (MQ=255)
                                                                                                                                                      aaTTAAAGGTCGGTAAAACGGGCCGTTCGGCGTAGCCCTGCGTTCCTGCAAGCGTTCCTTCGAAAAAGCAGGGGTTCTGGCGGGAGTTCGTCGTCGTGAGTTCTATGAAAAACCGACTACCGAACGTAAGCGCGCTAAAGCTTCTGCAGTGAAACGTCACGCGAAGAAACTGGCTCGCGAAAACGCACGCCGCACTCGTCTGTACTAATTCCCCGAGAGCGTTGctc                                                                          <  2:430552/227‑1 (MQ=255)
                                                                                                                                                            aGTACGTGAAAACGAGCCGTTCGACGTAGCTCTGCGTCGCTTCAAGCGTTCCTGCGAAAAAGCAGGTGTTCTGGCGGAAGTTCGTCGTCGTGAGTTCTATGAAAAACCGACTACCGAACGTAAGCGCGCTAAAGCTTCTGCAGTGAAACGTCACGCGAAGAAACTGGCTCGCGAAAACGCACGCCGCACTCGTCTGTACTAATTCCCCCAGAGCGTTGCTCCCCgcgc                                                                   >  1:514340/1‑225 (MQ=255)
                                                                                                                                                                cgggAAAACGAGCCGTTCGACGTAGCTTTGCGTCCCTTCAAGCGTTCCTGCGAAAAAGCAGGTGTTCTGGCGGAAGTTCGTCGTCGTGAGTTCTATGAAAAACCGACTACCGAACGTAAGCGCGCTAAAGCTTCTGCAGTGAAACGTCACGCGAAGAAACTGGCTCGCGAAAACGCACGCCGCACTCGTCTGTACTAATTCCCCGAGAGCGTTGCTCTCCGATCAGa                                                                <  1:375476/224‑1 (MQ=255)
                                                                                                                                                                cGTGAAAACGAGCCGTTCGACGTAGCTCTGCGTCGCTTCAAGCGTTCCTGCGAAAAAGCAGGTGTTCTGGCGGAAGTTCGTCGTCGTGAGTTCTATGAAAAACCGACTACCGAACGTAAGCGCGCTAAAGCTTCTGCAGTGAAACGTCACGCGACGAAACTGGCTCGCGAAAACGCACGCCCCCCTCGCCTGCCGAAATTACCCCAGAACGTTGCCCTCCGATCaacc                                                               >  2:514247/1‑225 (MQ=255)
                                                                                                                                                                cGTGAAAACGAGCCGTTCGACGTAGCTCTGCGTCGCTTCAAGCGTTCCTGCGAAAAAGCAGGTGTTCTGGCGGAAGTTCGTCGTCGTGAGTTCTATGAAAAACCGACTACCGAACGTAAGCGCGCTAAAGCTTCTGCAGTGAAACGTCACGCGAAGAACCTGGCTCGCAAAAACCCACGCCGCACCCGTCCGTCCTAATTCCCCCACAGCGGTACCCTCCgcacacac                                                               >  2:255865/1‑221 (MQ=255)
                                                                                                                                                                cGTGAAAACGAGCCGTTCGACGTAGCTCTGCGTCGCTTCAAGCGTTCCTGCGAAAAAGCAGGTGTTCTGGCGGAAGTTCGTCGTCGTGAGTTCTATGAAAAACCGACTACCGAACGTAAGCGCGCTAAAGCTTCTGCAGTGAAACGTCACGCGAAGAAACTGGCTCGCGAAAACGCACGCCGCACTCGTCTGTACTAATTCCCCGAGAGCGTTGCTCTCCGATCAGAc                                                               >  1:234487/1‑228 (MQ=255)
                                                                                                                                                                   gAAAACGAGCCGTTCGACGTAGCTCTGCGTCGCTTCAAGCGTTCCTGCGAAAAAGCAGGTGTTCTGGCGGAAGTTCGTCGTCGTGAGTTCTATGAAAAACCGACTGCCGAACGTAAGCGCGCTAAAGCTTCTGCAGTGAAACGTCACGCGAAGAAACTGGCTCGCGAAAACGCAAGCCGCACTCGTCTGCAATAATTCCCCGAGAGGGTTGCTCTCCGATCAGACCgg                                                            >  1:37881/1‑227 (MQ=255)
                                                                                                                                                                   gAAAACGAGCCGTTCGACGTAGCTCTGCGTCGCTTCAAGCGTTCCTGCGAAAAAGCAGGTGTTCTGGCGGAAGTTCGTCGTCGTGAGTTCTATGAAAAACCGACTACCGAACGTAAGCGCGCTAAAGCt                                                                                                                                                                >  2:135462/1‑129 (MQ=255)
                                                                                                                                                                    aaaaCGAGCCGTTCGACGTAGCTCTGCGTCGCTTCAAGCGTTCCTGCGAAAAAGCAGGTGTTCTGGCGGAAGTTCGTCGTCGTGAGTTCTATGAAAAACCGACTACCGAACGTAAGCGCGCTAAAGCtt                                                                                                                                                               <  1:135462/129‑1 (MQ=255)
                                                                                                                                                                          atccGTTCGACGTAGCGCTGCTTCGCTGCAAGCGTTCCTGCGAAAAAGCAGGTGTTATGGCGGAAGTTCGTCGTCGTGATTTCTATGAAAAACCGACTACCGAACGTAAGCGCGCTAAAGCTTCTGCAGTGAAACGTCACGCGAAGAAACTGGCTCGCGAAAACGCACGCCGCACTCGTCTGTACTAATTCCCCGAGAGCGTTGCTCTCCGATCAGACCGagttgt                                                       <  1:412608/224‑1 (MQ=255)
                                                                                                                                                                              gTTTGACTTTGCTCTGCGCTGCTTCAAGCGTTCCTGCAAAAAAGCAGGTGTTCTGGCGTAAGTTCGTCGTCGTGAGTTCTATGAAAACCCGACTACCGAACGTAAGCGCGCTAAAGCTTCTGCAGTGAAACGTCACGCGAAGAAACTGGCTCGCGAAAACGCACGCCGCACTCGTCTGTACTAATTCCCCGAGAGCGTTGCTCTCCGATCAGACCGagttgtagttg                                                  <  2:150760/227‑1 (MQ=255)
                                                                                                                                                                                        ggttggCGTCCCTTCAAGCGTTCCTGCGAAACAGCAGGTGTTCTGGCGGAAGTCCGTCGTCTTGAGTTCGAAGAAAACCCGACTACCGAACGTAAGCGCGCTAAAGCTTCTGCAGTGAAACGTCACGCGAAGAAACTGGCTCGCGAAAACGCACGCCGCACTCGTCTGTACTAATTCCCCGAGAGCGTTGCTCTCCGATCAGACCGAGTTGTAGTTGTAAGGCCGTg                                        <  2:598827/222‑1 (MQ=255)
                                                                                                                                                                                         ctctGCGTCGCTTCAAGCGTTCCTGCGAAAAAGCAGGTGTTCTGGCGGAAGTTCGTCGTCGTGAGTTCTATGAAAAACCGACTACCGAACGTAAGCGCGCTAAAGCTTCTGCAGTGAAACGTCACGCGAAGAAACTGGCTCGCGAAAACGCACGCCGCACTCGTCTGTACTAATTCCCCGAGAGCGTTGCTCTCCGGTCAGACCGAGTTGTAGTTGTAAGGcctgggt                                      >  1:408866/1‑223 (MQ=255)
                                                                                                                                                                                         ctctGCGTCGCTTCAAGCGTTCCTGCGAAAAAGCAGGTGTTCTGGCGGAAGTTCGTCGTCGTGAGTTCTATGAAAAACCGACTACCGAACGTAAGCGCGCTAAAGCTTCTGCAGTGAAACGCCACGCGAAGAAACTGGCTCGCGAAAACGCACGCCGCAATCGTCTGTACTAATTCCCCGACAGCGTTGCTCTCCGATCAGACCGCGTTGTAGTTGTAAGGcctggct                                      >  1:79138/1‑223 (MQ=255)
                                                                                                                                                                                         ctctGCGTCGCTTCAAGCGTTCCTGCGAAAAAGCAGGGGTTCTGGCGGAAGTTCGTCGTCGTGAGTTCTATGAAAAACCGACTACCGAACGTAAGCGCGCTAAAGCTTCTGCAGTGAAACGTCACGCGAAGAAACAGGCTCGCGAAAACGCACGCCGCCCGCGGCTGTGCTAATTCCCCGAGAGCGGTTCTCCCCGGCCAAACCCAGTTTTAGTTGTAAGGCCGGGcg                                      >  2:419657/1‑227 (MQ=255)
                                                                                                                                                                                             gCGTCGCTTGAAGCGTTCGTGAGAAAAAGCAGGTGTTCTGGCGGAAGTTCGTCGTCGTGAGTTCTAAGAAAAACCGACTACCGAACGTAAGCGCGCTAAAGCTTCTGCAGTGAAACGTCACGCGAAGAAACTGGCTCGCGAAAACGCACGCCGCACTCGTCTGTACTAATTCCCCGAGAGCGTTGCTCTCCGATCAGACCGAGTTGTAGTTGTAAGGCCGTg                                        <  1:264454/222‑1 (MQ=255)
                                                                                                                                                                                              ggTCGTTTCAAGCGTTCCTGCAAAAAAGCAGGTTTTCCGGCGGAAGGTCGGCGTCGTGAGTTCTAGGAAAAACCGACTACCGAACGTAAGCGCGCTAAAGCTTCTGCAGTGAAACGTCACGCGAAGAAACTGGCTCGCGAAAACGCACGCCGCACTCGTCTGTACTAATTCCCCGAGAGCGTTGCTCTCCGATCAGACCGAGTTGTAGTTGTAAGGCCGTGCTTCCGa                                 <  2:408866/227‑1 (MQ=255)
                                                                                                                                                                                               gTCGCTTCAAGCGTTCCTGCGAAAGAGCAGGGATGCTGGCGGAAGTTCGTCGTCGTAAGTTCTAGGAAAAACCGACTACCGAACGTAAGCGCGCTAAAGCTTCTGCAGTAAAACGTCACGCGAAGAAACTGGCTCGCGAAAACGCACGCCGCACTCGTCTGTACTAATTCCCCGAGAGCGTTGCTCTCCGATCAGACCGAGTTGTAGTTGTAAGGCCGTGCTTCCGa                                 <  1:248100/227‑1 (MQ=255)
                                                                                                                                                                                               gTCGCTTCAAGCGTTCCTGCGAAAAAGCAGGTGTTCTGGCGGAAGTTCGTCGTCGTGAGTTCTATGAAAAACCGACTACCGAACGTAAGCGCGCTAAAGCTTCTGCAGTGAAACGTCACGCGAAGAAACTGGCTCGCGAAAACGCACGCCGCACTCGTCTGTACCAATTCCCCCAGAGCGTGGCTCCCCGACCCGCCCCGGTTGTAGTTGTAAGGCCGAGCTTCCGaa                                >  2:41038/1‑228 (MQ=255)
                                                                                                                                                                                                   cTTCAAGCGTTCCTGCGAAAAAGCAGGTGTTCTGGCGGAAGTTCGTCGTCGTGAGTTCTATGAAAAACCGACTACCGAACGTAAGCGCGCTAAAGCTTCTGCAGTGAAACGTCACGCGAAGAAACTGGCTCGCGAAAACGCACGCCGCACTCGTCTGTACTAATTCCCCGAGAGCGTTGCTCTCCGATCAGACCGAGTTGTGGTTGTAAGGCCGTGCTTCCGAAAGGa                            >  1:434442/1‑228 (MQ=255)
                                                                                                                                                                                                        atggTTCCAGCGAAAAAGCGGTGTTTTCGGCGGAAGTTCGGCGGAGTGAGTTCTATGAAAAACCGACTACCGAACTTAAGCCCGCTAAAGCTTCTGCAGTGAAACGGCACGCGAAGAAACTGGCTCGCGAAAACGCACGCCGCACTCGTCTGTACTAATTCCCCGAGAGCGTTGCTCTCCGATCAGACCGAGTTGTAGTTGTAAGGCCGTGCTTCCGAAAGGAATgcg                       <  2:400228/225‑1 (MQ=255)
                                                                                                                                                                                                                 gCGAAAAAGCAGGTGTTCTGGCGGAAGTTCGTCGTCGTGAGTTCTATGAAAAACCGACTACCGAACGTAAGCGCGCTAAAGCTTCTGCAGTGAAACGTCACGCGAAGAAACTGGCTCGCGAAAACGCACGCCGCACTCGTCTGTACTCATTCCCCGAGAGCGTTCCTCTCCGATCGGACCGGGTTGTAGTTGTAAGGCCGTGCTCCCGAAAGGAACTCGCGGCTTAtt              >  1:558931/1‑228 (MQ=255)
                                                                                                                                                                                                                             gtgtTCTGGGGGAAGGTCGTCGTGGTGAGTTCTATGAAAAACCGACTACCGAACGTAAGCGCGCTAAAGCTTCTGCAGGGAAACGTCACGCGAAGAAACTGGCTCGCGAAAACGCACGCCGCACTCGTCTGTACTAATTCCCCGAGAGCGTTGCTCTCCGATCAGACCGAGTTGTAGTTGTAAGGCCGTGCTTCCGAAAGGAATGCGCGGCTTATTTTCGTTTATGaa  <  2:422301/228‑1 (MQ=255)
                                                                                                                                                                                                                              |                                                                                                                                                                                                                                     
CTCACCCTTATAAAAGTCCCTTTCAAAAAAGGCCGCGGTGCTTTACAAAGCAGCAGCAATTGCAGTAAAATTCCGCACCATTTTGAAATAAGCTGGCGTTGATGCCAGCGGCAAACCGAATTAATCAAAGGTGAGAGGCACATGCCGGTAATTAAAGTACGTGAAAACGAGCCGTTCGACGTAGCTCTGCGTCGCTTCAAGCGTTCCTGCGAAAAAGCAGGTGTTCTGGCGGAAGTTCGTCGTCGTGAGTTCTATGAAAAACCGACTACCGAACGTAAGCGCGCTAAAGCTTCTGCAGTGAAACGTCACGCGAAGAAACTGGCTCGCGAAAACGCACGCCGCACTCGTCTGTACTAATTCCCCGAGAGCGTTGCTCTCCGATCAGACCGAGTTGTAGTTGTAAGGCCGTGCTTCCGAAAGGAATGCGCGGCTTATTTTCGTTTATGAA  >  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: