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

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