Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 90,720 C→A 7.8% R209R (CGC→CGA rsmH → 16S rRNA m(4)C1402 methyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091390,7200CA 7.8% 221.0 / 3.3 65R209R (CGC→CGArsmH16S rRNA m(4)C1402 methyltransferase
Reads supporting (aligned to +/- strand):  ref base C (31/28);  new base A (2/4);  total (33/32)
Fisher's exact test for biased strand distribution p-value = 4.27e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GATGGTCCGCTGGACATGCGTATGGACCCAACCCGTGGGCAGTCAGCCGCTGAATGGCTACAAACCGCAGAAGAAGCCGATATCGCCTGGGTATTGAAAACCTATGGTGAAGAGCGTTTTGCCAAACGCATTGCCCGCGCCATTGTCGAGCGTAACCGCGAACAGCCGATGACCCGCACCAAAGAACTGGCGGAAGTCGTGGCTGCTGCAACGCCGGTGAAAGATAAGTTTAAACATCCCGCGACCCGTACCTTCCAGGCGGTGCGCATTTGGGTAAACAGTGAACTGGAGGAGATAGAGCAGGCGCTAAAAAGCTCGCTCAACGTGCTGGCCCCGGGTGGGCGGCTTTCGATCATCAGCTTCCACTCGCTGGAAGACCGTATTGTGAAACGTTTTATGCGTGAAAACAGCCGCGGTCCGCAAGTTCCGGCAGGGTTACCGATGACTGAAGAGCAGCTCAAAAAACTGGGTGGCCGTCAGCTGCGAGCACTAGGCAAGTTAATGCC  >  NC_000913/90454‑90959
                                                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                               
gATGGTCCGCTGGACATGCGTATGGACCCAACCCGTGGGCAGTCAGCCGCTGAATGGCTACAAACCGCAGAAGAAGCCGATATCGCCTGGGTATTGAAAACCTATGGTGAAGAGCGTTTTGCAAAACGCATTGCCCGCGCCATTGTCGAGCGTAACCGCGAACAGCCGATGACCCGCACCAAAGAACTGGCGGAAGTCGTGGCTGCTGCAACGCCGGTGAAAGATAAGTTTAAACATCCCGCGACCCGTACCTTCCAGGCGGTGCGCATTTg                                                                                                                                                                                                                                            <  1:129564/272‑1 (MQ=255)
             agaagggTAGGGACCCAACCCGTGGGCAGACACCCGCTGAAGGGCTACAAACCGCAGAAGAAGCCGATATCGCCTGGGTATTGAAAACCTATGGTGAAGAGCGTTTTGCCAAACGCATTGCCCGCGCCATTGTCGAGCGTAACCGCGAACAGCCGATGACCCGCACCAAAGAACTGGCGGAAGTCGTGGCTGCTGCAACGCCGGTGAAAGATAAGTTTAAACATCCCGCGACCCGTACCTTCCAGGCGGTGCGCATTTGGGTAAACAGt                                                                                                                                                                                                                                  <  2:65749/263‑1 (MQ=255)
               aacagttgtgACCAATCCCGTGGGCATTCATCCGTTGAATGGTCTCAAACCGCAGAAGAAGCCGATATCGCCTGGGTATTGAAAACCTATGGTGAAGAGCGTTTTGCCAAACGCATTGCCCGCGCCATTGTCGAGCGTAACCGCGAACAGCCGATGACCCGCACCAAAGAACTGGCGGAAGTCGTGGCTGCTGCAACGCCGGTGAAAGATAAGTTTAAACATCCCGCGACCCGTACCTTCCAGGCGGTGCGCATTTGGGTAAACAGTGAACt                                                                                                                                                                                                                             <  2:129292/263‑1 (MQ=255)
                                ccGTGGGCAGTCAGCCGCTGAATGGCTACAAACCGCAGAAGAAGCCGATATCGCCTGGGTATTGAAAACCTATGGTGAAGAGCGTTTTGCCAAACGCATTGCCCGCGCCATTGTCGAGCGTAACCGCGAACAGCCGATGACCCGCACCAAAGAACTGGCGGAAGTCGTGGCTGCTGCAACGCCGGTGAAAGATAAGTTTAAACATCCCGCGACCCGTACCTTCCAGGCGGTGCGCATTTTGGGAAAAAGGGAACAGGGGGGGAAAAaggagg                                                                                                                                                                                                            >  2:470836/1‑272 (MQ=255)
                                 cGTGGGCAGTCAGCCGCTGAATGGCTACAAACCGCAGAAGAAGCCGATATCGCCTGGGTATTGAAAACCTATGGTGAAGAGCGTTTTGCCAAACGCATTGCCCGCGCCATTGTCGAGCGTAACCGCGAACAGCCGATGACCCGCACCAAAGAACTGGCGGAAGTCGTGGCTGCTGCAACGCCGGTGAAAGATAAGTTTAAACATCCCGCGACCCGTACCTTCCAGGCGGTGCGCATTTGGGTAAACAg                                                                                                                                                                                                                                   >  1:65749/1‑248 (MQ=255)
                                 cGTGGGCAGTCAGCCGCTGAATGGCTACAAACCGCAGAAGAAGCCGATATCGCCTGGGTATTGAAAACCTATGGTGAAGAGCGTTTTGCCAAACGCATTGCCCGCGCCATTGTCGAGCGTAACCGCGAACAGCCGATGACCCGCACCAAAGAACTGGCGGAAGTCGTGGCTGCTGCAACGCCGGTGAAAGATAAGTTTAAACATCCCGCGACCCGTACCTTCCAGGCGGTGCGCATTTGGGTAAACAGTGAAc                                                                                                                                                                                                                              >  1:129292/1‑253 (MQ=255)
                                 cGTGGGCAGTCAGCCGCTGAATGGCTACAAACCGCAGAAGAAGCCGATATCGCCTGGGTATTGAAAACCTATGGTGAAGAGCGTTTTGCCAAACGCATTGCCCGCGCCATTGTCGAGCGTAACCGCGAACAGCCGATGACCCGCACCAAAGAACTGGCGGAAGTCGTGGCTGCTGCAACGCCGGTGAAAGATAAGTTTAAACATCCCGCGACCCGTACCTTCCAGGCGGTGCGCATTTGGGTAAACAGTGAACTGGAGGAGATAGAGCAGgc                                                                                                                                                                                                           >  1:619276/1‑272 (MQ=255)
                                                                      aagaagCCGATATCGCCTGGGTATTGAAAACCTATGGTGAAGAGCGTTTTGCCAAACGCATTGCCCGCGCCATTGTCGAGCGTAACCGCGAACAGCCGATGACCCGCACCAAAGAACTGGCGGAAGTCGTGGCTGCTGCAACGCCGGTGAAAGATAAGTTTAAACATCCCGCGACCCGTACCTTCCAGGCGGTGCGCAGTTGGGTAAACAGTGAACTGGAGGAGATAGAGAAGGCGCTAAAAAGCTCGCCCAACGTGATGGCCCCGGGTgg                                                                                                                                                                       >  1:516781/1‑271 (MQ=255)
                                                                              gATATCGCCTGGGTATTGAAAACCTATGGTGAAGAGCGTTTTGCCAAACGCATTGCCCGCGCCATTGTCGAGCGTAACCGCGAACAGCCGATGACCCGCACCAAAGAACTGGCGGAAGTCGTGGCTGCTGCAACGCCGGTGAAAGATAAGTTTAAACATCCCGCGACCCGTACCTTCCAGGCGGTGCGCATTTGGGTAAACAGTGAACTGGAGGAGATAGAGCAGGCGCTAAAAAGCTCGCTCAACGTGCTGGCCCCGGGTGggcggcatca                                                                                                                                                              >  1:307682/1‑268 (MQ=255)
                                                                               atatCGCCTGGGTAGTGAAAACCTATGGTGAAGAGCGTTTTGCCAAACGCATTGCCCGCGCCATTGTCGAGCGTAACCGCGAACAGCCGATGACCCGCACCAAAGAACTGGCGGAAGTCGTGGCTGCTGCAACGCCGGTGAAAGATAAGTTTAAACATCCCGCGACCCGTACCTTCCAGGCGGTGCGCATTTGGGTAAACAGTGAACTGGAGGAGATAGAGCAGGCGCTAAAAAGCTCGCTCAACGTGCTGGCCCCGGGTGGGCGGCTTTCg                                                                                                                                                             <  1:419324/272‑1 (MQ=255)
                                                                                  cgcccTGTGTTTTGGAAACCTTTTGGTGACATGCGTTTTGCCAAACGCTTGGCCCGCGCCATTGTCGAGCTTAACCGCGAACAGCCGATGACCCGAACCAAAGAACTGGCGGAAGTCGTGGCTGCTGCAACGCCGGTGAAAGATAAGTTTAAACATCCCGCGACCCGTACCTTCCAGGCGGTGCGCATTTGGGTAAACAGTGAACTGGAGGAGATAGAGCAGGCGCTAAAAAGCTCGCTCAACGTGCTGGCCCCGGGTGGGCGGCTTTCGat                                                                                                                                                           <  2:917256/269‑1 (MQ=255)
                                                                                                 aaaCCTATGGTGAAGAGCGTTTTGCCAAACGCATTGCCCGCGCCATTGTCGAGCGTAACCGCGAACAGCCGATGACCCGCACCAAAGAACTGGCGGAAGTCGTGGCTGCTGCAACGCCGGTGAAAGATAAGTTTAAACATCCCGCGACCCGTACCTTCCAGGCGGTGCGCATTTGGGTAAACAGTGAACTGGAGGAGATAGAGCAGGCGCTAAAAAGCTCGCCCaa                                                                                                                                                                                         >  1:724392/1‑226 (MQ=255)
                                                                                                 aaaCCTATGGTGAAGAGCGTTTTGCCAAACGCATTGCCCGCGCCATTGTCGAGCGTAACCGCGAACAGCCGATGACCCGCACCAAAGAACTGGCGGAAGTCGTGGCTGCTGCAACGCCGGTGAAAGATAAGGTTAAACATCCCGCGACCCGTACCTTCCAGGCGGTGCGCATTTGGGTAAACAGTGAACTGGAGGAGATAGAGCAGGCGCTAAAAAGCTCGCTCAACGTGCTGGCCCCGGGTGGGCGGCTTTCGAACAACAGCTTCAACTCg                                                                                                                                           >  1:691943/1‑272 (MQ=255)
                                                                                                          ggTGAAGAGCGTTTTGCCAAACGCATTGCCCGCGCCATTGTCGAGCGTAACCGCGAACAGCCGATGACCCGCACCAAAGAACTGGCGGAAGTCGTGGCTGCTGCAACGCCGGTGAAAGATAAGTTTAAACATCCCGCGACCCGTACCTTCCAGGCGGTGCGCATTTGGGTAAACAGTGAACTGGAGGAGATAGAGCAGGCGCTAAAAAGCTCGCTCAACGTGCt                                                                                                                                                                                   <  1:674177/224‑1 (MQ=255)
                                                                                                           gTGAAGAGCGTTTTGCCAAACGCATTGCCCGCGCCATTGTCGAGCGTAACCGCGAACAGCCGATGACCCGCACCAAAGAACTGGCGGAAGTCGTGGCTGCTGCAACGCCGGTGAAAGATAAGTTTAAACATCCCGCGACCCGTACCTTCCAGGCGGTGCGCATTTGGGTAAACAGTGAACTGGAGGAGATAGAGCAGGCGCTAAAAAgcccg                                                                                                                                                                                              >  2:392632/1‑212 (MQ=255)
                                                                                                            tGAAGAGCGTTTTGCCAAACGCATTGCCCGCGCCATTGTCGAGCGTAACCGCGAACAGCCGATGACCCGCACCAAAGAACTGGCGGAAGTCGTGGCTGCTGCAACGCCGGTGAAAGATAAGTTTAAACATCCCGCGACCCGTACCTTCCAGGCGGTGCGCATTTGGGTAAACAGTGAACTGGAGGAGATAGAGCAGGCGCTAAAAAgctcgc                                                                                                                                                                                             <  1:392633/212‑1 (MQ=255)
                                                                                                                            aaccGCATTGCCCGCGCAATGGTCGAGCGTAACCGCGAACAGCCGATGACCCGCACCAAAGAACTGGCGGAAGTCGTGGCTGCTGAAACGCCGGTGAAAGATAAGTTTAACCATCCCGCGACCCGTACCTTCCAGGAGTTGCGAATTTGGGAAAACAGTGAAGTGGAGGAGATAGAGCAGGCGCTAAAAAGCTCGCTCAACGTGCTGGCCCGGGGTGGGCGGCTTTCGATCATCAGCTTCCACTCGCTGGAAGACCGTATTGTGAAACGtt                                                                                                                  <  1:470837/268‑1 (MQ=255)
                                                                                                                                       ccGCGCCATTGTCGAGCGTAACCGCGAACAGCCGATGACCCGCACCAAAGAACTGGCGGAAGTCGTGGCTGCTGCAACGCCGGTGAAAGATAAGTTTAAACATCCCGCGACCCGTACCTTCCAGGCGGTGCGCATTTGGGTAAACAGTGAACTGGAGGAGATAGAGCAGGCGCTAAAAAGCTCGCTCAACg                                                                                                                                                                                       >  1:388000/1‑191 (MQ=255)
                                                                                                                                        cgcgCCATTGTCGAGCGTAACCGCGAACAGCCGATGACCCGCACCAAAGAACTGGCGGAAGTCGTGGCTGCTGCAACGCCGGTGAAAGATAAGTTTAAACATCCCGCGACCCGTACCTTCCAGGCGGTGCGCATTTGGGTAAACAGTGAACTGGAGGAGATAGAGCAGGCGCTAAAAAGCTCGCTCAACGt                                                                                                                                                                                      <  2:387999/191‑1 (MQ=255)
                                                                                                                                                        gTAACCGCGAACAGCCGATGACCCGCACCAAAGAACTGGCGGAAGTCGTGGCTGCTGCAACGCCGGTGAAAGATAAGTTTAAACATCCCGCGACCCGTACCTTCCAGGCGGTGCGCATTTGGGTAAACAGTGAACTGGAGGAGATAGAGCAGGCGCTAAAAAGCTCGCTCAACGTGCTGGCCCCggg                                                                                                                                                                          >  1:431703/1‑187 (MQ=255)
                                                                                                                                                        gTAACCGCGAACAGCCGATGACCCGCACCAAAGAACTGGCGGAAGTCGTGGCTGCTGCAACGCCGGTGAAAGATAAGTTTAAACATCCCGCGACCCGTACCTTCCAGGCGGTGCGCATTTGGGTAAACAGTGAACTGGAGGAGATAGAGCAGGCGCTAAAAAGCTCGCTCAACGTGCTGGCCCCGGGTGGGCGGCTTTCGATCATCAGCTTCCACTCGCTGGAAGACCGTATTGTGAAACGTTTTATGCGTGAAAACAGCCGCGGTCCGCaa                                                                                     >  1:229415/1‑272 (MQ=255)
                                                                                                                                                         tAACCGCGAACAGCCGATGACCCGCACCAAAGAACTGGCGGAAGTCGTGGCTGCTGCAACGCCGGTGAAAGATAAGTTTAAACATCCCGCGACCCGTACCTTCCAGGCGGTGCGCATTTGGGTAAACAGTGAACTGGAGGAGATAGAGCAGGCGCTAAAAAGCTCGCTCAACGTGCTGGCCCCGGGt                                                                                                                                                                         <  2:431702/187‑1 (MQ=255)
                                                                                                                                                            ccGCGAACAGCCGATGACCCGCACCAAAGAACTGGCGGAAGTCGTGGCTGCTGCAACGCCGGTGAAAGATAAGTTTAAACATCCCGCGACCCGTACCTTCCAGGCGGTGCGCATTTGGGTAAACAGTGAACTGGAGGAGATAGAGCAGGCGCTAAAAAGCTCGCTCAACGTGCTGGCCCCGGGTGGGCGGCTTTCGATCATCAGCTTCCACTCGCTGGAAGACCGTATTGTGAAACGTTTTATGCGTGAAAACAGCCGCGGTCCGCAAGTTc                                                                                 >  1:498722/1‑272 (MQ=255)
                                                                                                                                                             cgcgAACAGCCGATGACCCGCACCAAAGAACTGGCGGAAGTCGTGGCTGCTGCAACGCCGGTGAAAGATAAGTTTAAACATCCCGCGACCCGTACCTTCCAGGCGGTGCGCATTTGGGTAAACAGTGAACTGGAGGAGATAGAGCAGGCGCTAAAAAGCTCGCTCAACGTGCTGGCCCCGGGTGGGCGGCTTTCGATCATCAGCTTCCACTCGCTGGAAGACCGTATTGTGAAACGTTTTATGCGTGAAAACAGCcg                                                                                               <  1:611346/257‑1 (MQ=255)
                                                                                                                                                               cgAACAGCCGATGACCCGCACCAAAGAACTGGCGGAAGTCGTGGCTGCTGCAACGCCGGTGAAAGATAAGTTTAAACATCCCGCGACCCGTACCTTCCAGGCGGTGCGCATTTGGGTa                                                                                                                                                                                                                                        >  2:760623/1‑118 (MQ=255)
                                                                                                                                                                gAACAGCCGATGACCCGCACCAAAGAACTGGCGGAAGTCGTGGCTGCTGCAACGCCGGTGAAAGATAAGTTTAAACATCCCGCGACCCGTACCTTCCAGGCGGTGCGCATTTGGGTaa                                                                                                                                                                                                                                       <  1:760624/118‑1 (MQ=255)
                                                                                                                                                                                   ccAAAGAACTGGCGGAAGTCGTGGCTGCTGCAACGCCGGTGAAAGATAAGTTTAAACATCCCGCGACCCGTACCTTCCAGGCGGTGCGCATTTGGGTAAACAGTGAAc                                                                                                                                                                                                                              >  1:917490/1‑108 (MQ=255)
                                                                                                                                                                                   ccAAAGAACTGGCGGAAGTCGTGGCTGCTGCAACGCCGGTGAAAGATAAGTTTAAACATCCCGCGACCCGTACCTTCCAGGCGGTGCGCATTTGGGTAAACAGTGAACTGGAGGAGATAGAGCAGGCGCTAAAAAGCTCGCTCAACGTGCTg                                                                                                                                                                                  >  1:94064/1‑152 (MQ=255)
                                                                                                                                                                                   ccAAAGAACTGGCGGAAGTCGTGGCTGCTGCAACGCCGGTGAAAGATAAGTTTAAACATCCCGCGACCCGTACCTTCCAGGCGGTGCGCATTTGGGTAAACAGTGAACTGGAGGAGATAGAGCAGGCGCTAAAAAGCTCGCTCAACGTGCTGGCCCCGGGTGGGCGGCTTTCGATCATCAGCTTCCACTCGCTGGAAGACCGTATTGTGAAACGTTTTATGCGTGAAAACAGCCGCGGTCCGCAAGTTCCGGCAGGGTTACCGATGACTGaa                                                          >  1:597860/1‑272 (MQ=255)
                                                                                                                                                                                    cAAAGAACTGGCGGAAGTCGTGGCTGCTGCAACGCCGGTGAAAGATAAGTTTAAACATCCCGCGACCCGTACCTTCCAGGCGGTGCGCATTTGGGTAAACAGTGAACt                                                                                                                                                                                                                             <  2:917489/108‑1 (MQ=255)
                                                                                                                                                                                    cAAAGAACTGGCGGAAGTCGTGGCTGCTGCAACGCCGGTGAAAGATAAGTTTAAACATCCCGCGACCCGTACCTTCCAGGCGGTGCGCATTTGGGTAAACAGTGAACTGGAGGAGATAGAGCAGGCGCTAAAAAGCTCGCTCAACGTGCTgg                                                                                                                                                                                 <  2:94064/152‑1 (MQ=255)
                                                                                                                                                                                          aCTGGCGGAAGTCGTGGCTGCTGCAACGCCGGTGAAAGATAAGTTTAAACATCCCGCGACCCGTACCTTCCAGGCGGTGCGCATTTGGGTAAACAg                                                                                                                                                                                                                                   >  1:889480/1‑96 (MQ=255)
                                                                                                                                                                                          aCTGGCGGAAGTCGTGGCTGCTGCAACGCCGGTGAAAGATAAGTTTAAACATCCCGCGACCCGTACCTTCCAGGCGGTGCGCATTTGGGTAAACAGTGAACTGGAGGAGATAGAGCAGGCGCTAAAAAGCTCGCTCAACGTGCTGGCCCCGGGTGGGCGGCTTTCGATCATCAGCTTc                                                                                                                                                 >  2:1022671/1‑178 (MQ=255)
                                                                                                                                                                                           cTGGCGGAAGTCGTGGCTGCTGCAACGCCGGTGAAAGATAAGTTTAAACATCCCGCGACCCGTACCTTCCAGGCGGTGCGCATTTGGGTAAACAGt                                                                                                                                                                                                                                  <  2:889479/96‑1 (MQ=255)
                                                                                                                                                                                           cTGGCGGAAGTCGTGGCTGCTGCAACGCCGGTGAAAGATAAGTTTAAACATCCCGCGACCCGTACCTTCCAGGCGGTGCGCATTTGGGTAAACAGTGAACTGGAGGAGATAGAGCAGGCGCTAAAAAGCTCGCTCAACGTGCTGGCCCCGGGTGGGCGGCTTTCGATCATCAGCTTcc                                                                                                                                                <  1:1022672/178‑1 (MQ=255)
                                                                                                                                                                                              gCGGAAGTCGTGGCTGCTGCAACGCCGGTGAAAGATAAGTTTAAACATCCCGCGACCCGTACCTTCCAGGCGGTGCGCATTTGGGTAAACAGTGAACTGGAGGAGATAGAGCAGGCGCTAAAAAGCTCGCTCAACGTGCTGGCCCCggg                                                                                                                                                                          >  1:904687/1‑149 (MQ=255)
                                                                                                                                                                                               cGGAAGTCGTGGCTGCTGCAACGCCGGTGAAAGATAAGTTTAAACATCCCGCGACCCGTACCTTCCAGGCGGTGCGCATTTGGGTAAACAGTGAACTGGAGGAGATAGAGCAGGCGCTAAAAAGCTCGCTCAACGTGCTGGCCCCGGGt                                                                                                                                                                         <  2:904686/149‑1 (MQ=255)
                                                                                                                                                                                                 gAAGTCGTGGCTGCTGCAACGCCGGTGAAAGATAAGTTTAAACATCCCGCGACCCGTACCTTCCAGGCGGTGCGCATTTGGGTAAACAGTGAACTGGAGGAGATAGAGCAGGCGCTAAAAAGCTCGCTCAACGTGCTGGCCCCGGGTGGGCGGCTTTCGa                                                                                                                                                            >  2:1110519/1‑160 (MQ=255)
                                                                                                                                                                                                 gAAGTCGTGGCTGCTGCAACGCCGGTGAAAGATAAGTTTAAACATCCCGCGACCCGTACCTTCCAGGCGGTGCGCATTTGGGTAAACAGTGAACTGGAGGAGATAGAGCAGGCGCTAAAAAGCTCGCTCAACGTGCTGGCCCCGGGTGGGCGGCTTTCGATCATCAGCt                                                                                                                                                   >  1:615319/1‑169 (MQ=255)
                                                                                                                                                                                                 gAAGTCGTGGCTGCTGCAACGCCGGTGAAAGATAAGTTTAAACATCCCGCGACCCGTACCTTCCAGGCGGTGCGCATTTGGGTAAACAGTGAACTGGAGGAGATAGAGCAGGCGCGAAAAAGCTCGCTCAACGTGCTGGCCCCGGGTGGGCGGCTTTCGATCATCAGCTTCCACTCGCTGGAAGACCGTATTGTGAAACGTTTTATGCGTGAAAACAGCCGCGGTCCGCAAGTTCCGGCAGGGTTACCGATGACTGAag                                                         >  1:253418/1‑259 (MQ=255)
                                                                                                                                                                                                  aaGTCGTGGCTGCTGCAACGCCGGTGAAAGATAAGTTTAAACATCCCGCGACCCGTACCTTCCAGGCGGTGCGCATTTGGGTAAACAGTGAACTGGAGGAGATAGAGCAGGCGCTAAAAAGCTCGCTCAACGTGCTGGCCCCGGGTGGGCGGCTTTCGat                                                                                                                                                           <  1:1110520/160‑1 (MQ=255)
                                                                                                                                                                                                  aaGTCGTGGCTGCTGCAACGCCGGTGAAAGATAAGTTTAAACATCCCGCGACCCGTACCTTCCAGGCGGTGCGCATTTGGGTAAACAGTGAACTGGAGGAGATAGAGCAGGCGCTAAAAAGCTCGCTCAACGTGCTGGCCCCGGGTGGGCGGCTTTCGATCATCAGCtt                                                                                                                                                  <  2:615318/169‑1 (MQ=255)
                                                                                                                                                                                                      ccaggatgctgcAAAGCCGGTGAAGGATAAGTTTTAACATCCCGCGACCGGTCCCTTCCGAGCGGTGGGCATTTGGGTAAACAGTGAACTGGAGGAGGTAGAGCAGGCGCTAAAAAGCTCGCTCAACGTGCTGGCCCCGGGTGGGCGGCTTTCGATCATCAGCTTCCACTCGCTGGAAGACCGTGTTGTGAAACGTTTTATGCGTGAAAACAGCCGCGGTCCGCAAGTTCCGGCAGGGTTACCGATGACTGAAGAGCAGCTCAAAAAACt                                         <  2:691942/264‑1 (MQ=255)
                                                                                                                                                                                                        ttgttgctgcAACGCCGGGTAAATATAAGTTTAAACATCCCGCGACCCGTACCTTCCAGGCGGTGCGCATTTGGGTAAACAGTGAACTGGAGGAGATAGAGCAGGCGCTAAAAAGCTCGCTCAACGTGCTGGCCCCGGGTGGGCGGCTTTCGATCATCAGCTTCCACTCGCTGGAAGACCGTATTGTGAAACGTTTTATGCGTGAAAACAGCCGCGGTCCGCAAGTTCCGGCAGGGTTACCGATGACTGAAGAACAGCTCAAAAAACTgg                                       <  1:929609/266‑1 (MQ=255)
                                                                                                                                                                                                                 cAACGCCGGTGAAAGATAAGTTTAAACATCCCGCGACCCGTACCTTCCAGGCGGTGCGAATTTGGGTAAACAGTGAACTGGAGGAGATAGAGCAGGCGCTAAAAAGCTCGCTCAACGTGCTGGCCCCGGGTGGGCGGCTTTCGATCATCAGCt                                                                                                                                                   >  2:907484/1‑153 (MQ=255)
                                                                                                                                                                                                                 cAACGCCGGTGAAAGATAAGTTTAAACATCCCGCGACCCGTACCTTCCAGGCGGTGCGAATTTGGGTAAACAGTGAACTGGAGGAGATAGAGCAGGCGCTAAAAAGCTCGCTCAACGTGCTGGCCCCGGGTGGGCGGCTTTCGATCATCAGCt                                                                                                                                                   >  2:361943/1‑153 (MQ=255)
                                                                                                                                                                                                                  aaCGCCTGTGAAAGATAAGGTTAAACATCCCGCGACCCGTACCTTCCAGGCGGTGCGAATTTGGGTAAACAGTGAACTGGAGGAGATAGAGCAGGCGCTAAAAAGCTCGCTCAACGTGCTGGCCCCGGGTGGGCGGCTTTCGATCATCAGCtt                                                                                                                                                  <  1:907485/153‑1 (MQ=255)
                                                                                                                                                                                                                  aaCGCCGGTGAAAGATAAGTTTAAACATCCCGCGACCCGTACCTTCCAGGCGGTGCGAATTTGGGTAAACAGTGAACTGGAGGAGATAGAGCAGGCGCTAAAAAGCTCGCTCAACGTGCTGGCCCCGGGTGGGCGGCTTTCGATCATCAGCtt                                                                                                                                                  <  1:361944/153‑1 (MQ=255)
                                                                                                                                                                                                                      ccGGTGAAAGATAAGTTTAAACATCCCGCGACCCGTACCTTCCAGGCGGTGCGCATTTGGGTAAACAg                                                                                                                                                                                                                                   >  1:524109/1‑68 (MQ=255)
                                                                                                                                                                                                                       cGGTGAAAGATAAGTTTAAACATCCCGCGACCCGTACCTTCCAGGCGGTGCGCATTTGGGTAAACAGt                                                                                                                                                                                                                                  <  2:524108/68‑1 (MQ=255)
                                                                                                                                                                                                                            aaaGATAAGTTTAAACATCCCGCGACCCGTACCTTCCAGGCGGTGCGCATTTGGGTAAACAg                                                                                                                                                                                                                                   >  1:241486/1‑62 (MQ=255)
                                                                                                                                                                                                                            aaaGATAAGTTTAAACATCCCGCGACCCGTACCTTCCAGGCGGTGCGCATTTGGGTAAACAGTGAACTGGAGGAGATAGAGCAGGCGCTAAAAAGCTCGCTCAACGTGCTGGCCCCGGGTGGGCGGCTTTCGATCATCAGCTTCCACTCGCTGGAAGACCGTAtt                                                                                                                            >  1:562340/1‑165 (MQ=255)
                                                                                                                                                                                                                            aaaGATAAGTTTAAACATCCCGCGACCCGTACCTTCCAGGCGGTGCGCATTTGGGTAAACAGTGAACTGGAGGAGATAGAGCAGGCGCTAAAAAGCTCGCTCAACGTGCTGGCCCCGGGTGGGCGGCTTTCGATCATCAGCTTACACTCGCTGGAAGACCGTAATGTGAAACGTTTTATGCGTGAAAACAGCCGCGGTCCGCAAGTTCCGGCAGGGTTCACGATGCCTGAAGAGCAGCACAAAAAACTGGGTGGGCGTCAggccgcagccc                  >  2:660557/1‑261 (MQ=255)
                                                                                                                                                                                                                             aaGATAAGTTTAAACATCCCGCGACCCGTACCTTCCAGGCGGTGCGCATTTGGGTAAACAGt                                                                                                                                                                                                                                  <  2:241486/62‑1 (MQ=255)
                                                                                                                                                                                                                             aaGATAAGTTTAAACATCCCGCGACCCGTACCTTCCAGGCGGTGCGCATTTGGGTAAACAGTGAACTGGAGGAGATAGAGCAGGCGCTAAAAAGCTCGCTCAACGTGCTGGCCCCGGGTGGGCGGCTTTCGATCATCAGCTTCCACTCGCTGGAAGACCGTATtg                                                                                                                           <  2:562339/165‑1 (MQ=255)
                                                                                                                                                                                                                                           ctcccgccgccccGTACTTTCCGGGGGGTTCGAAGTTGGGTAAGCAGTGAACTGGGGGGGATAGAGCAGGCGCTAAAAAGCTCGCTAAACGTGCTGGCCCCGGGTGGGCGGCTTTCGATCATCAGCTTCCACTCGCTGGAAGACCGTATTGTGAAACGTTTTATGCGTGAAAACAGCCGCGGTCCGCAAGTTCCGGCAGGGTTACCGATGACTGAAGAGCAGCTCAAAAAACTGGGTGGCCGTCAGCTGCGAGCACTAGGCAAGTTAATGcc  <  2:498721/262‑1 (MQ=255)
                                                                                                                                                                                                                                            aTCCCGCGACCCGTACCTTCCAGGCGGTGCGCATTTGGGTAAACAGTGAACTGGAGGAGATAGAGCAGGCGCTAAAAAGCTCGCTCAACGTGCTGGCCCCGGGTGGGCGGCTTTCGATCATCAGCTTCCACTCGCTGGAAGACCGTATTGTGAAACGTTTTATGCGTGAAAACAGCCGCGGTCCGCAAGTTCCGGCAGGGTTAc                                                                     >  1:283730/1‑204 (MQ=255)
                                                                                                                                                                                                                                             tCCCGCGACCCGTACCTTCCAGGCGGTGCGCATTTGGGTAAACAGTGAACTGGAGGAGATAGAGCAGGCGCTAAAAAGCTCGCTCAACGTGCTGGCCCCGGGTGGGCGGCTTTCGATCATCAGCTTCCACTCGCTGGAAGACCGTATTGTGAAACGTTTTATGCGTGAAAACAGCCGCGGTCCGCAAGTTCCGGCAGGGTTAcc                                                                    <  2:283729/204‑1 (MQ=255)
                                                                                                                                                                                                                                                  cgACCCGTACCTTCCAGGCGGTGCGCATTTGGGTAAACAGTGAACTGGAGGAGATAGGGCAGGCGCTAAAAAGCTCGCTCAACGTGCTGGCCCCGGGTGGGCGGCTTTCGATCATCAGCTTCCACTCGCTGGAAGAGCGTATTGTGAAACGTTTTATGCGTGAAAACAGCCGCGGTCCGCAAGTTCCGGCAGGGTTACCGATGACTGAAGAGCAGCTCAAAAAACTGGGTGGCCGTCAGCt                          <  1:605462/241‑1 (MQ=255)
                                                                                                                                                                                                                                                  cgACCCGTACCTTCCAGGCGGTGCGCATTTGGGTAAACAGTGAACTGGAGGAGATAGAGCAGGCGCTAAAAAGCTCGCTCAACGTGCTGGCCCCGGGTGGGCGGCTTTCGATCATCAGCTTCCACTCGCTGGAAGACCGTATTGTGAAACGTTTTATGCGTGAAAACAGCCGCGGTCCGCAAGTTCCGGCAGGGTTAc                                                                     >  1:1089624/1‑198 (MQ=255)
                                                                                                                                                                                                                                                   gACCCGTACCTTCCAGGCGGTGCGCATTTGGGTAAACAGTGAACTGGAGGAGATAGAGCAGGCGCTAAAAAGCTCGCTCAACGTGCTGGCCCCGGGTGGGCGGCTTTCGATCATCAGCTTCCACTCGCTGGAAGACCGTATTGTGAAACGTTTTATGCGTGAAAACAGCCGCGGTCCGCAAGTTCCGGCAGGGTTAcc                                                                    <  2:1089623/198‑1 (MQ=255)
                                                                                                                                                                                                                                                                  ggCGGTGCGCATTTGGGTAAACAGTGAACTGGAGGAGATAGAGCAGGCGCTAAAAAGCTCGCTCAACGTGCTGGCCCCGGGTGGGCGGCTTTCGATCATCAGCTTCCACTCGCTGGAAGACCGTATTGTGAAACGTTTTATGCGTGAAAACAGc                                                                                                 >  2:608066/1‑154 (MQ=255)
                                                                                                                                                                                                                                                                   gCGGTGCGCATTTGGGTAAACAGTGAACTGGAGGAGATAGAGCAGGCGCTAAAAAGCTCGCTCAACGTGCTGGCCCCGGGTGggcgg                                                                                                                                                                   >  2:1113515/1‑87 (MQ=255)
                                                                                                                                                                                                                                                                   gCGGTGCGCATTTGGGTAAACAGTGAACTGGAGGAGATAGAGCAGGCGCTAAAAAGCTCGCTCAACGTGCTGGCCCCGGGTGGGCGGCTTTCGATCATCAGCTTCCACTCGCTGGAAGACCGTATTGTGAAACGTTTTATGCGTGAAAACAGcc                                                                                                <  1:608067/154‑1 (MQ=255)
                                                                                                                                                                                                                                                                    cGGTGCGCATTTGGGTAAACAGTGAACTGGAGGAGATAGAGCAGGCGCTAAAAAGCTCGCTCAACGTGCTGGCCCCGGGTGggcggc                                                                                                                                                                  <  1:1113516/87‑1 (MQ=255)
                                                                                                                                                                                                                                                                         cgcATTTGGGTAAACAGTGAACTGGAGGAGATAGAGCAGGCGCTAAAAAGCTCGCTCAACGTGCTGGCCCCGGGTGGGCGGCTTTCGATCATCAGCTTCCACTCGCTGGAAGACCGTATTGTGAAACGTTTTATGCGTGAAAACa                                                                                                   >  2:131252/1‑145 (MQ=255)
                                                                                                                                                                                                                                                                          gcATTTGGGTAAACAGTGAACTGGAGGAGATAGAGCAGGCGCTAAAAAGCTCGCTCAACGTGCTGGCCCCGGGTGGGCGGCTTTCGATCATCAGCTTCCACTCGCTGGAAGACCGTATTGTGAAACGTTTTATGCGTGAAAACAg                                                                                                  <  1:131252/145‑1 (MQ=255)
                                                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                               
GATGGTCCGCTGGACATGCGTATGGACCCAACCCGTGGGCAGTCAGCCGCTGAATGGCTACAAACCGCAGAAGAAGCCGATATCGCCTGGGTATTGAAAACCTATGGTGAAGAGCGTTTTGCCAAACGCATTGCCCGCGCCATTGTCGAGCGTAACCGCGAACAGCCGATGACCCGCACCAAAGAACTGGCGGAAGTCGTGGCTGCTGCAACGCCGGTGAAAGATAAGTTTAAACATCCCGCGACCCGTACCTTCCAGGCGGTGCGCATTTGGGTAAACAGTGAACTGGAGGAGATAGAGCAGGCGCTAAAAAGCTCGCTCAACGTGCTGGCCCCGGGTGGGCGGCTTTCGATCATCAGCTTCCACTCGCTGGAAGACCGTATTGTGAAACGTTTTATGCGTGAAAACAGCCGCGGTCCGCAAGTTCCGGCAGGGTTACCGATGACTGAAGAGCAGCTCAAAAAACTGGGTGGCCGTCAGCTGCGAGCACTAGGCAAGTTAATGCC  >  NC_000913/90454‑90959

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 9 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: