Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,906,796 C→T V89M (GTG→ATG)  rlmA ← 23S rRNA m(1)G745 methyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,906,7960CT97.3% 216.0 / NA 73V89M (GTG→ATG) rlmA23S rRNA m(1)G745 methyltransferase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (27/44);  minor base G (0/2);  total (27/46)
Fisher's exact test for biased strand distribution p-value = 5.27e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCACTACTCGTGCTAATTCTTCTGCTTTACACGGCGCGTAAATACGTATTATGGCGTCCATACTGGTATCGGAAAACGGCAAACGGTGGCTGGAAGCGACACAAAAAGTGACCTGCGGATAGCGTTTCGCCGCCGCTTTTATCGCTACCTTCGAAACATCCAGACCAAACGTGGTGATTTCGGGCAACGCATCGGCAAATGCGTGTGTGTAATACCCTTCACCACAGCCAATATCCAGCACCGCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTGGGCGACAATTGCATCACGCAGCGGCTGATAATGTCCGGCATCTAAGAATGCGCGGCGTGCTTGCATCATTTCCGCGCTGTCGCCCGGATCACGAGACCGTTTATGCTGAACGGGCAGCAGATTGACATACCCTTCTTTCGCCATATCAAACTGATGTCGCTGGGGACAGATATAGCTGTTTTTTTCACGCGAAAGAGGCTGATGGCAAAG  >  NC_000913/1906556‑1907042
                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                      
tCACTACTCGTGCTAATTCTTCTGCTTTACACGGCGCGTAAATACGTATTATGGCGTCCATACTGGTATCGGAAAACGGCAAACGGTGGCTGGAAGCGACACAAAAAGTGACCTGCGGATAGCGTTTCGCCGCCGCTTTTATCGCTACCTTCGAAACATCCAGACCAAACGTGGTGATTTCGGGCAACGCATCGGCAAATGCGTGTGTGTAATACCCTTCACCACAGCCAATATCCAGCATCGCCGTGGcc                                                                                                                                                                                                                                              >  2:209672/1‑251 (MQ=255)
tCACTACTCGTGCTAATTCTTCTGCTTTACACGGCGCGTAAATACGTATTATGGCGTCCATACTGGTATCGGAAAACGGCAAACGGTGGCTGGAAGCGACACAAAAAGTGACCTGCGGATAGCGTTTCGCCGCCGCTTTTATCGCTACCTTCGAAACATCCAGACCAAACGTGGTGATTTCGGGCAACGCATCGGCAAATGCGTGTGTGTAATACCCGTCACCACAGCCAATATCCAGCATCGCCGTgg                                                                                                                                                                                                                                                <  1:480875/249‑1 (MQ=255)
     actCGTGCTAATTCTTCTGCTTTACACGGCGCGTAAATACGTATTATGGCGTCCATACTGGTATCGGAAAACGGCAAACGGTGGCTGGAAGCGACACAAAAAGTGACCTGCGGATAGCGTTTCGCCGCCGCTTTTATCGCTACCTTCGAAACATCCAGACCAAACGTGGTGATTTCGGGCAACGCATCGGCAAATGCGTGTGTGTAATACCCTTCACCACAGCCAATATCCAGCATCGCCGTGGCCTTatc                                                                                                                                                                                                                                         <  2:98500/251‑1 (MQ=255)
     actCGTGCTAATTCTTCTGCTTTACACGGCGCGTAAATACGTATTATGGCGTCCATACTGGTATCGGAAAACGGCAAACGGTGGCTGGAAGCGACACAAAAAGTGACCTGCGGATAGCGTTTCGCCGCCGCTTTTATCGCTACCTTCGAAACATCCAGACCAAACGTGGTGATTTCGGGCAACGCATCGGCAAATGCGTGTGTGTAATACCCTTCACCACAGCCAATATCCAGCATCGCCGTGGCCTTat                                                                                                                                                                                                                                          <  1:180201/250‑1 (MQ=255)
              aaTTCTTCTGCTTTACACGGCGCGTAAATACGTATTATGGCGTCCATACTGGTAGCGGAAAACGGCAAACGGTGGCTGGAAGCGACACAAAAAGTGACCTGCGGATAGCGTTTCGCCGCCGCTTTTATCGCTACCTTCGAAACATCCAGACCAAACGTGGTGATTTCGGGCAACTCATCGGCAAATGCGTGTGTGTAATACCCTTCACCACAGCCAATATCCAGCATCGCCGTGGCCTTATCATCAAGcc                                                                                                                                                                                                                                 <  2:526494/250‑1 (MQ=255)
                    tctGCTTTACACGGCGCGTAAATACGTATTATGGCGTCCATACTGGTATCGGAAAACGGCAAACGGTGGCTGGAAGCGACACAAAAAGTGACCTGCGGATAGCGTTTCGCCGCCGCTTTTATCGCTACCTTCGAAACATCCAGACCAAACGTGGTGATTTCGGGCAACGCATCGGCAAATGCGTGTGTGTAATACCCTTCACCACAGCCAATATCCAGCATCGCCGTGGCCTTATCATCAAGCCGTTcc                                                                                                                                                                                                                            <  1:197406/249‑1 (MQ=255)
                    tctGCTTTACACGGCGCGTAAATACGTATTATGGCGTCCATACTGGTATCGGAAAACGGCAAACGGTGGCTGGAAGCGACACAAAAAGTGACCTGCGGATAGCGTTTCGCCGCCGCTTTTATCGCTACCTTCGAAACATCCAGACCAAACGTGGTGATTTCGGGCAACGCATCGGCAAATGCGTGTGTGTAATACCCTTCACCACAGCCAATATCCAGCATCGCCGTGGCCTTATCATCAAGCCGTTCCCt                                                                                                                                                                                                                          >  1:523118/1‑251 (MQ=255)
                     ctGCTTTACCCGGCGCGTAAATACGTATTATGGCGTCCATACTGGTATCGGAACACGGCAAACGGTGGCTGGAAGCGACACAAAAAGTGACCTGCGGATAGCGTTTCGCCGCCGCTTTTATCGCTACCTTCGAAACATCCAGACCCAACGTGGTGATTTCGGGCAACGCATCGGCAAATGCGTGTGTGTAATACCCTTCACCACAGCCAATATCCAGCATCGCCGTGGCCTTATCATCAAGCCGTTCCCTc                                                                                                                                                                                                                         >  2:492360/1‑251 (MQ=255)
                       gCTTTACACGGCGCGTAAATACGTATTATGGCGTCCATACTGGTATCGGAAAACGGCAAACGGTGGCTGGAAGCGACACAAAAAGTGACCTGCGGATAGCGTTTCGCCGCCGCTTTTATCGCTACCTTCGAAACATCCAGACCAAACGTGGTGATTTCGGGCAACGCATCGGCAAATGCGTGTGTGTAATACCCTTCACCACAGCCAATATCCAGCATCGCCGTGGCCTTATCATCAAGCCGTTCCCTc                                                                                                                                                                                                                         <  1:15456/249‑1 (MQ=255)
                           tACACGGCGCGTAAATACGTATTATGGCGTCCATACTGGTATCGGAAAACGGCAAACGGTGGCTGGAAGCGACACAAAAAGTGACCTGCGGATAGCGTTTCGCCGCCGCTTTTATCGCTACCTTCGAAACATCCAGACCAAACGTGGTGATTTCGGGCAACGCATCGGCAAATGCGTGTGTGTAATACCCTTCACCACAGCCAATATCCAGCATCGCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTg                                                                                                                                                                                                                    <  1:94455/250‑1 (MQ=255)
                            acacGGCGCGTAAATACGTATTATGGCGTCCATACTGGTATCGGAAACCGGCAAACGGTGGCGGGAAGCGACACAAAAAGTGACCTGCGGATAGCGTTTCGCCGCCGCTTTTATCGCTACCTTCGAAACATCCAGACCAAACGTGGTGATTTCGGGCAACGCATCGGCAAATGCGTGTGTGTAATACCCTTCACCACAGCCAATATCCAGCATCGCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTgg                                                                                                                                                                                                                   <  2:33654/250‑1 (MQ=255)
                                   gcttAAATACGTATTATGGCGTCCATACTGGTATCGGAAAACGGCAAACGGTGGCTGGAAGCGACACAAAAAGTGACCTGCGGATAGCGTTTCGCCGCCGCTTTTATCGCTACCTTCGAAACATCCAGACCAAACGTGGTGATTTCGGGCAACGCATCGGCAAATGCGTGTGTGTAATACCCTTCACCACAGCCAATATCCAGCATCGCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTGGGCGACAAt                                                                                                                                                                                                           <  1:170950/248‑1 (MQ=255)
                                     ggaaaTACGTATTATGGCGTCCATACTGGTAGCGGAAAACGGCAAACGGTGGCTGGAAGCGACACAAAAAGGGACCTGCGGATAGCGTTTCGCCGCCGCTTTTATCGCTACCTGCGAAACAGCCAGTCCAAACGTGGTGATTTCGGGCAACGCAGCGGCAAATGCGTGTGTGTAATACCCTTCACCACAGCCAATATCCAGCAGCGCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTGGGCGACAATTg                                                                                                                                                                                                         <  1:571253/249‑1 (MQ=255)
                                     gTAAATACGTATTATGGCGACCATACTGGTCTCGGAAACCGACAAACGGACGATGGAAGCGACACACAAAGTGACCTGCGGATAGCGTTTCGCCGCCGCTTTTATCGCTACCTTAGAAACATACAGACCAAACGTGGTGATTTCGGGCAACGCATCGGCAAATGCGTGTGTGTAATAACCTTCACCAAAGACAATATCCAGCATCGCCGTGGCCTTATCATCAAGCCGTTCCCTAAGTTGGGCGACAATTg                                                                                                                                                                                                         >  2:75408/1‑251 (MQ=255)
                                       aaaTACGTATTATGGCGTCCATACGGGTATCGGAAAACGGCAAACGGTGGCTGGAAGCGACACAAAAAGGGACCTGCGGATAGCGTTTCGCCGCCGCTTTTATCGCTACCTTCGAAACATCCAGACCAAACGTGGTGATTTCGGGCAACGCATCGGCAAATGCGTGTGTGTAATACCCTTCACCACAGCCAATATCCAGCATCGCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTGGGCGTCAATTGc                                                                                                                                                                                                        <  2:286647/250‑1 (MQ=255)
                                        aaTACGTATTATGGCGTCCATACTGGTATCGGAAAACGGCAAACGGTGGCTGGAAGCGACACAAAAAGTGACCTGCGGATAGCGTTTCGCCGCCGCTTTTATCGCTACCTTCGAAACATCCAGACCAAACGTGGTGATTTCGGGCAACGCATCGGCAAATGCGTGTGTGTAATACCCTTCACCACAGCCAATATCCAGCATCGCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTGGGCGACAATTGCAt                                                                                                                                                                                                      >  1:464768/1‑251 (MQ=255)
                                         aTACGGATTGTGGCGTCCATATTGGTAGCGGAAAACGGCAAACGGTGGCTGGAAGCGACACAAAAAGGGACCTGCGGATAGCGTTTCGCCGCCGCTTTTATCGCTACCTTCGAAACACCCATACAAAACGTGGTGATTGCGGGCAACGCATCGGCAAATGCGTGTGTGTAATACCCTTCACCACAGCCAATATCCAGCATCGCCGTTGACTTATCATCAAGCCGTTCCCTCAGTTGGGCGACAATTGCATc                                                                                                                                                                                                     <  1:348364/251‑1 (MQ=255)
                                          agcGTATTATGGCGTCCATACGGGTATCGGAAAACGGCAAACGGTGGCTGGAAGCGACACAAAAAGTGACCTGCGGATAGCGTTTCGCCGCCGCTTTTATCGCTACCTTCGAAACATCCAGACCAAACGTGGTGATTTCGGGCAACGCATCGGCAAATGCGTGTGTGTAATACCCTTCACCACAGCCAATATCCAGCATCGCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTGGGCGACAATTGCATCa                                                                                                                                                                                                    <  2:90920/249‑1 (MQ=255)
                                             gTATTATGGCGTCCATACTGGTATCGGAAAACGGCAAACGGTGGCTGGAAGCGACACAAAAAGTGACCTGCGGATAGCGTTTCGCCGCCGCTTTTATCGCTACCTTCGAAACATCCAGACCAAACGTGGTGATTTCGGGCAACGCATCGGCAAATGCGTGTGTGTAATACCCTTCACCACAGCCAATATCCAGCATCGCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTGGGCGACAATTGCATCACGc                                                                                                                                                                                                 <  1:267911/251‑1 (MQ=255)
                                               attatGGCGTCCATACTGGTATCGGAAAACGGCAAACGGTGGCTGGAAGCGACACAAAAAGTGACCTGCGGATAGCGTTTCGCCGCCGCTTTTATCGCTACCTTCGAAACATCCAGACCAAACGTGGTGATTTCGGGCAACGCATCGGCAAATGCGTGTGTGTAATACCCTTCACCACAGCCAATATCCAGCATCGCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTGGGCGACAATTGCATCACGc                                                                                                                                                                                                 >  2:132667/1‑249 (MQ=255)
                                                 tatGGCGTCCATACTGGTATCGGAAAACGGCAACCGTTGGCTGGAAGCGACACAAAAAGTGACCTGCGGATAGCTTTTCGCCGCCGCTTTTATTGCTACCTTCGAAACATCCAGACCAAACGTGGTGATTGCGGGCAACGCATCGTCAAATGCGTGTGTGTAATACCCTTCACCACAGCCAATATCCAGCATCGCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTGGGCGACAATTGCATCACGCAGc                                                                                                                                                                                              <  2:514599/250‑1 (MQ=255)
                                                    ggCGTCCATACTGGTATCGGAAAACGGCAAACGGTGGCTGGAAGCGACACAAAAAGTGACCTGCGGATAGCGTTTCGCCGCCGCTTTTATCGCTACCTTCGAAACATCCAGACCAAACGTGGTGATTTCGGGCAACGCATCGGCAAATGCGTGTGTGTAATACCCTTCACCACAGCCAATATCCAGCATCGCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTGGGCGACAATTGCATCACGCAGCGGCt                                                                                                                                                                                          <  2:213515/251‑1 (MQ=255)
                                                    ggCGTCCATACTGGTATCGGAAAACGGCAAACGGTGGCTGGAAGCGACACAAAAAGTGACCTGCGGATAGCGTTTCGCCGCCGCTTTTATCGCTACCTTCGAAACATCCAGACCAAACGTGGTGATTTCGGGCAACGCATCGGCAAATGCGTGTGTGTAATACCCTTCACCACAGCCAATATCCAGCATCGCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTGGGCGACAATTGCATCACGCAGCGGCt                                                                                                                                                                                          <  1:265153/251‑1 (MQ=255)
                                                             aCTGGTATCGGAAAACGGCAAACGGTGGCTGGAAGCGACACAAAAAGTGACCTGCGGATAGCGTTTCGCCGCCGCTTTTATCGCTACCTTCGAAACATCCAGACCAAACGTGGTGATTTCGGGCAACGCATCGGCAAATGCGTGTGTGTAATACCCTTCACCACAGCCAATATCCAGCATCGCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTGGGCGACAATTGCATCACGCAGCGGCTGATAATGTc                                                                                                                                                                                 >  1:233357/1‑251 (MQ=255)
                                                                          accggaaaaCGGTGGCTGGAGGCGACACAAAAAGTGACCTGCGGATAGCGTTTCGCCGCCGCTTTTATCGCTACCTTCGAAACATCCAGACCAAACGTGGTGATTTCGGGCAACGCATCGGCAAATGCGTGTGTGTAATACCCTTCACCACAGCCAATATCCAGCATCGCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTGGGCGACAATTGCATCACGCAGCGGCTGATAATGTCCGGCAGCTAAGaa                                                                                                                                                                    <  2:142334/245‑1 (MQ=255)
                                                                                  aCGGTGGCTGGAAGCGACACAAAAAGTGACCTGCGGATAGCGTTTCGCCGCCGCTTTTATCGCTACCTTCGAAACATCCAGACCAAACGTGGTGATTTCGGGCAACGCATCGGCAAATGCGTGTGTGTAATACCCTTCACCACAGCCAATATCCAGCATCGCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTGGGCGACAATTGCATCACGCAGCGGCTGATAATGTCCGGCATCTAAGAATGCgcggc                                                                                                                                                            <  1:167093/251‑1 (MQ=255)
                                                                                   cGGTGGCTGGAAGCGACACAAAAAGTGACCTGCGGATAGCGTTTCGCCGCCGATTTTATCGCTACCTTCGAAACATCCAGACCAAACGTGGTGATTTCGGGCAACGCATCGGCAAATGCGTGTGTGTAATACCCTTCACCACAGCCAATATCCAGAATCGCAGTGGCCTTATCATCAAGCCGTTCCCTCAGTTGGGCGACAATTGCATCACGCAGCGGCTGATAATGTCCGGCATCTAAGAATGCgcggcg                                                                                                                                                           >  1:335882/1‑251 (MQ=255)
                                                                                          tGGAAGCGACACAAAAAGAGACCTGCGGATAGCGTTTCGCCGCCGCTTTTATCGCTACCTTCGAAACATCCAGACCAAACGTGGAGATTTCGGGCAACGCATCGGCAAATGCGTGTGTGTAATACCCTTCACCACAGCCAATATCCAGCATCGCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTGGGCGACAATTGCATCACGCAGCGGCTGATAATGTCCGGCATCTAAGAATGCGCGGCGTTCTTGc                                                                                                                                                    >  2:563403/1‑251 (MQ=255)
                                                                                                cGACACAAAAAGTGACCTGCGGATAGCGTTTCGCCGCCGCTTTTATCGCTACCTTCGAAACATCCAGACCAAACGTGGTGATTTCGGGCAACGCATCGGCAAATGCGTGTGTGTAATACCCTTCACCACAGCCAATATCCAGCATCGCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTGGGCGACAATTGCATCACGCAGCGGCTGATAATGTCCGGCATCTAAGAATGCGCGGCGTGCTGGCATCAtt                                                                                                                                              <  1:431007/251‑1 (MQ=255)
                                                                                                        aaaGTGACCTGCGGATAGCGTTTCGCCGCCGCTTTTATCGCTACCTTCGAAACATCCAGACCAAACGTGGTGATTTCGGGCAACGCATCGGCAAATGCGTGTGTGTAATACCCTTCACCACAGCCAATATCCAGCATCGCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTGGGCGACAATTGCATCACGCAGCGGCTGATAATGTCCGGCATCTAAGAATGCGCGGCGTGCTTGCATCATTTCCGCGCt                                                                                                                                      <  2:459602/251‑1 (MQ=255)
                                                                                                         aaGTGACCGGCGGATTGCGTTTCGCCGCCGCTTTTATCGCTACCGTCGAAACATCCAGACCAAACGTGGTGATTTCGGGCAACGCATCGGCAAATGCGTGTGTGTAATACCCTTCACCACAGCCAATATCCAGCATCGCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTGGGCGACAATTGCATCACGCAGCGGCTGATAATGTCCGGCATCTAAGAATGCGCGGCGTGCTTGCATCATTTCCGCGCTg                                                                                                                                     <  2:315715/251‑1 (MQ=255)
                                                                                                           gTGACCTGCGGATAGCGTTTCGCCGCCGCTTTTATCGCTACCTTCGAAACATCCAGACCAAACGTGGTGATTTCGGGCAACGCATCGGCAAATGCGTGTGTGTAATACCCTTCACCACAGCCAATATCCAGCATCGCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTGGGCGACAATTGCATCACGCAGCGGCTGATAATGTCCGGCATCTAAGAATGCGCGGCGTGCTTGCATCATTTCCGCGCTGTc                                                                                                                                   >  2:183456/1‑251 (MQ=255)
                                                                                                             gACCTGCGGATAGCGTTTCGCCGCCGCTTTTATCGCTACCTTCGAAACATCCAGACCAAACGTGGTGATTTCGGGCAACGCATCGGCAAATGCGTGTGTGTAATACCCTTCCCCACAGCCAATATCCAGCATCGCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTGGGCGACAATTGCATCACGCAGCGGCTGATAATGTCCGGCATCTAAGAATGCGCGGCGTGCTTGCATCATTTCCGCGCTGTCGc                                                                                                                                 >  2:381083/1‑251 (MQ=255)
                                                                                                              aCCTGCGGATAGCGTTTCGCCGCCGCTTTTATCGCTACCTTCGAAACATCCAGACCAAACGTGGTGATTTCGGGCAACGCATCGGCAAATGCGTGTGTGTAATACCCTTCACCACAGCCAATATCCAGCATCGCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTGGGCGACAATTGCATCACGCAGCGGCTGATAATGTCCGGCATCTAAGAATGCGCGGCGTGCTTGCATCATTTCCGCGCTGTCGcc                                                                                                                                >  1:337917/1‑251 (MQ=255)
                                                                                                               ccTGCGGAGAGCGTTTCGCCGCCGCTTTTATCGCGACCTTCGAAACAACCAGACCAAACGTGGTGATTTCGGGCAACGCATCGGCAAATGCGTGTGTGTAATACCCTTCACCACCGCCAATATCCAGCATCGCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTGGGCGACAATTGCATCACGCAGCGGCTGATAATGTCCGGCATCTAAGAATGCGCGGCGTGCTTGCATCATTTCCGCGCTGTCGccc                                                                                                                               <  2:314690/251‑1 (MQ=255)
                                                                                                                   cGGATAGCTTTTCGCCGCCGCTTTTATCGCTACCTGCGAAACAGCCAGACCAAACGGGGTGATTTCGGGCAACGCATCGGCAAAGGCGTGTGTGTAATTCCCTTCACCACAGCCAATATCCAGCATCGCCGTGGCCTTATCATCAAGCCTTTCCCTCAGTTGGGCGACAATTGCATCACGCAGCGGCTGATAATGTCCGGCATCTAAGAATGCGCGGCGTGCTTGCATCATTTCCGCGCTGTCGCCAGGAt                                                                                                                           <  2:554909/251‑1 (MQ=255)
                                                                                                                   cGGATAGCGTTTCGCCGCCGCTTTTATCGCTACCTTCGAAACATCCAGACCAAACGTGGTGATTTCGGGCAACGCATCGGCAAATGCGTGTGTGTAATACCCTTCACCACAGCCAATATCCAGCATCGCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTGGGCGACAATTGCATCACGCAGCGGCTGATAATGTCCGGCATCTAAGAATGCGCGGCGTGCTTGCATCATTTCCGCGCTGTCGCCCGGAt                                                                                                                           >  1:420986/1‑251 (MQ=255)
                                                                                                                       tAGCGTTTCGCCGCCGCTTTTATCGCTACCTTCGAAACATCCAGACCAAACGTGGTGATTTCGGGCAACGCATCGGCAAATGCGTGTGTGTAATACCCTTCACCACAGCCAATATCCAGCATCGCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTGGGCGACAATTGCATCACGCAGCGGCTGATAATGTCCGGCATCTAAGAATGCGCGGCGTGCTTGCATCATTTCCGCGCTGTCGCCCGGATCACg                                                                                                                       >  1:136250/1‑251 (MQ=255)
                                                                                                                          ggttgcgccgccgcTTTTATCGCTGCCTTCGAAACATCCAGACCAAACGTGGTGATTTCGGGCAACGCATCGGCAAATGCGTGTGTGTAATACCCTTCACCTCAGCCAATATCCAGCAGCTCCGTGGCCTTATCAGCAAGCCGTTCCCTCAGGTGGGCGACAATTGCGTCACGCAGCGGCTGAGGATGTCCGGCGTCTAAGAATGCGCGGCGGGCTTGCGTCAGTTCCGCTCTGTCGCCCGGATCACgaga                                                                                                                    <  2:368280/246‑1 (MQ=255)
                                                                                                                            tttCGCCGCCGCTTTTATCGCTACCTTCGAAACATCCAGACCAAACGTGGTGATTTCGGGCAACGCATCGGCAAATGCGTGTGTGTAATACCCTTCACCACAGCCAATATCCAGCATCGCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTGGGCGACAATTGCATCACGCAGCGGCTGATAATGTCCGGCATCTAAGAATGCGCGGCGTGCTTGCATCATTTCCGCGCTGTCGCCCGGATCACGAGAcc                                                                                                                  >  1:279642/1‑251 (MQ=255)
                                                                                                                            gttCGCCGCCGCTTTTATCGCTACCTTCGAAACATCCAGACCAAACGTGGTGATTTCGGGCAACGCATCGGCAAATGCGTGTGTGTAATACCCTTCACCACAGCCAATATCCAGCATCGCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTGGGCGACAATTGCATCACGCAGCGGCTGATAATGTCCGGCATCTAAGAATGCGCGGCGTGCTTGCATCATTTCCGCGCTGTCGCCCGGATCACGAGAcc                                                                                                                  <  2:233357/250‑1 (MQ=255)
                                                                                                                               cgccgccgcTTTTATCGCTACCTTCGAAACATCCAGACCAAACGTGTTGATTTCGGGCAACGCATCGGCAAATGCGTGTGTGTAATACCCTTCACCACAGCCAATATCCAGCATCGCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTGGGCGACAATTGCATCACGCAGCGGCTGATAATGTCCGGCATCTAAGAATGCGCGGCGTGCTTGCATCATTTCCGCGCTGTCGCCCGGATCACGAGACCGtt                                                                                                               <  2:186177/251‑1 (MQ=255)
                                                                                                                               cgccgccgcTTTTATCGCTACCTTCGAAACATCCAGACCAAACGTGGTGATTTCGGGCAACGCATCGGCAAATGCGTGTGTGTAATACCCTTCACCACAGCCAATATCCAGCATCGCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTGGGCGACAATTGCATCACGCAGCGGCTGATAATGTCCGGCATCTAAGAATGCGCGGCGTGCTTGCATCATTTCCGCGCTGTCGCCCGGATCACGAGACCGt                                                                                                                <  1:504507/250‑1 (MQ=255)
                                                                                                                                       cTTTTCTCGCTACCTTAGAAACATCCAGACCAAACGTGGTGATTTCGGGCAACGCATCGGCAAATGCGTGTGTGTAATACCCTTCACCACAGCCAATATCCAGCATCGCCGTGGCCTTATCATCAAGCCGTTCCCTCCGTTGGGCGACAATTGCATCACGCAGCGGCTGATAATGTCCGGCATCTAAGAATGCGCGGCGTGCTTGCATCATTTCCGCGCTGTCGCCCGGATCACGAGACCGTTTATGCTg                                                                                                        >  2:29963/1‑250 (MQ=255)
                                                                                                                                              cgcgaCCTTCGAAACATCCAGACCAAACGTGGTGATTTCGGGCAACGCATCGGCAAATGCGTGTGTGTAATACCCTTCACCACAGCCAATATCCAGCATCGCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTGGGCGACAATTGCATCACGCAGCGGCTGATAATGTCCGGCATCTAAGAATGCGCGGCGTGCGTGCATCATGGCCGCGCTGTCGCCCGGATCACGAGACCGTTTATGGTGAACGGgca                                                                                                <  2:245346/247‑1 (MQ=255)
                                                                                                                                              cGCTACCTTAGAAACATCCAGACCAAACGTGGTGATTTCGGGCAACTCATCGGCAAATGCGTGTGTGTAATACCCTTCACCACAGCCAATATCCAGCATCGCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTGGGCGACAATTGCATCACGCAGCGGCTGATAATGTCCGGCATCTAAGAATGCGCGGCGTGCTTGCATCATTGCCGCGCTGTCGCCCGGATCACGAGATCGTTTTTGCTGAACGGgca                                                                                                <  2:410005/251‑1 (MQ=255)
                                                                                                                                               gCTACATTCGAAACATCCAGACCAAACGTGGTGCTTTCGGGCAACGCATCGGCAAATGCGTGTGTGTAATACCCTTCACCACAGCCAATATCCAGCATCGCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTGGGCGACAATTGCATCACGCAGCGGCTGATAATGTCCGGCATCTAAGAATGCGCGGCGTGCTTGCATCATTTCCGCGCTGTCGCCCGGATCACGAGACCGTTTATGCTGAACGGgc                                                                                                 >  2:572912/1‑249 (MQ=255)
                                                                                                                                                      tCGAAACATCCAGTCCATACGTGTTGATTTCGGGCAACTGATCGGCAGATTCGTGTGTGTAATACCCTCCACCACAGCCAATATCCAGCATCGCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTGGGCGACAATTGCATCACGCAGCGGCTGATAATGTCCGGCAGCTAAGAATGCGCGGCGTGCTTGCATCATTTCCGCGCTGTCTCCCGGATCACGAGACCGTTTATGCTGAACGGGCAGCAGAt                                                                                          <  2:20453/249‑1 (MQ=255)
                                                                                                                                                             aTCCAGACCAAACGTGGTGATTTCGGGCAACTCATCGGCAAATGCGTGTGTGTAATACCCTTCACCACAGCCAATATCCAGCATCGCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTGGGCGACAATTGCATCACGCAGCGGCTGATAATGTCCGGCATCTAAGAATGCGCGGCGTGCTTGCATCATTTCCGCGCTGTCGCCCGGATCACGAGACCGTTTCTGCTGAACGGGCAGCAGATTGACATAcc                                                                                 >  2:66102/1‑251 (MQ=255)
                                                                                                                                                                 aGACCAAACGTGGTGATTTCGGGCAACGCATCGGCAAATGCGTGTGTGTAATACCCTTCACCACAGCCAATATCCAGCATCGCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTGGGCGACAATTGCATCACGCAGCGGCTGATAATGTCCGGCATCTAAGAATGCGCGGCGTGCTTGCATCATTTCCGCGCTGTCGCCCGGATCACGAGACCGTTTATGCTGAACGGGCAGCAGATTGACATACCcttc                                                                             <  1:278050/251‑1 (MQ=255)
                                                                                                                                                                    ccAAACGTGGTGATTTCGGGCAACGCATCGGCAAATGCGTGTGTGTAATACCCTTCACCACAGCCAATATCCAGCATCGCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTGGGCGACAATTGCATCACGCAGCGGCTGATAATGTCCGGCATCTAAGAATGCGCGGCGTGCTTGCATCATTTCCGCGCTGTCGCCCGGATCACGAGACCGTTTATGCTGAACGGGCAGCAGATTGACATACCCTTCttt                                                                          <  2:494434/251‑1 (MQ=255)
                                                                                                                                                                     cAAACGGGGTGATTTCGGGCAACGCATCGGCAAATGCGGGTGTGTAATACCCTTCACCACAGCCAATATCCAGCATCGCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTGGGCGACAATTGCATCACGCAGCGGCTGATAATGTCCGGCATCTAAGAATGCGCGGCGTGCTTGCATCATTTCCGCGCTGTCGCCCGGATCACGAGACCGTTTATGCTGAACGGGCAGCAGATTGACATACCCTTCTTTc                                                                         <  2:341076/251‑1 (MQ=255)
                                                                                                                                                                               gATTTCGGGCAACGCATCGGCAAATGCGTGTGTGTAATACCCTTCACCACAGCCAATATCCAGCATCGCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTGGGCGACAATTGCATCACGCAGCGGCTGATAATGTCCGGCATCTAAGAATGCGCGGCGTGCTTGCATCATTTCCGCGCTGTCGCCCGGATCACGAGACCGTTTATGCTGAACGGGCAGCAGATTGACATACCCTTCTTTCGCCATATCa                                                                >  2:270048/1‑250 (MQ=255)
                                                                                                                                                                                   tCGGGCAACGCATCGGCAAATGCGTGTGTGTAATCCCCGTCACCACAGCCAATATCCAGCATCGCCGTGGCCGTATCATCAAGCCGTTCCCTCAGTTGGGCGACAATTGCATCACGCAGCGGCTGATAATGTCCGGCATCTAAGAATGCGCGGCGTGCTTGCATCATTTCCGCGCTGTCGCCCGGATCACGAGACCGTTTATGCTGAACGGGCAGCAGATTGACATACCCTTCTTTCGCCATATCAAACt                                                            <  1:565267/250‑1 (MQ=255)
                                                                                                                                                                                    cGGGCAACGCATCGGCAAATGCGTGTGTGTAATACCCTTCACCACAGCCAATATCCAGCATCGCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTGGGCGACAATTGCATCACGCAGCGGCTGATAATGTCCGGCATCTAAGAATGCGCGGCGTGCTTGCATCATTTCCGCGCTGTCGCCCGGATCACGAGACCGTTTATGCTGAACGGGCAGCAGATTGACATACCCTTCTTTCGCCATATCAAACTGa                                                          <  1:183456/251‑1 (MQ=255)
                                                                                                                                                                                             cATCGGCAAATGCGTGTGTGTAATACCCTTCACCACAGCCAATATCCAGCATCGCCGTTGCCTTATCATCAAGCCGTTCCCTCAGTTGGGCGACAATTGCATCACGCAGCGGCTGATAATGTCCGGCATCTAAGAATGCGCGGCGTGCTTGCATCATTTCCGCGCTGTCGCCCGGATCACGAGACCGTTTATGCTGAACGGGCAGCAGATTGACATACCCTTCTTTCGCCATATCAAACTGATGTCGCTgg                                                 <  1:282296/251‑1 (MQ=255)
                                                                                                                                                                                                    aaaTGCGTGTGTGTAATACCCTTCACCACAGCCAATATCCAGCATCGCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTGGGCGACAATTGCATCACGCAGCGGCTGATAATGTCCGGCATCTAAGAATGCGCGGCGTGCTTGCATCATTTCCGCGCTGTCGCCCGGATCACGAGACCGTTTATGCTGAACGGGCAGCAGATTGACATACCCTTCTTTCGCCATATCAAACTGATGTCGCTGGGGACAGa                                          <  1:92761/251‑1 (MQ=255)
                                                                                                                                                                                                     aaTGCGTGTGTGTAATACCCTTCACCACAGCCAATATCCAGCATCGCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTGGGCGACAATTGCATCACGCAGCGGCTGATAATGTCCGGCATCTAAGAATGCGCGGCGTGCTTGCATCATTTCCGCGCTGTCGCCCGGATCACGAGACCGTTTATGCTGAACGGGCAGCAGATTGACATACCCTTCTTTCGCCATATCAAACTGATGTCGCTGGGGACAGat                                         <  1:233621/251‑1 (MQ=255)
                                                                                                                                                                                                       tGCGTGTGTGTAATACCCTTCACCACAGCCAATATCCAGCATCGCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTGGGCGACAATTGCATCACGCAGCGGCTGATAATGTCCGGCATCTAAGAATGCGCGGCGTGCTTGCATCATTTCCGCGCTGTCGCCCGGATCACGAGACCGTTTATGCTGAACGGGCAGCAGATTGACATACCCTTCTTTCGCCATATCAAACTGATGTCGCTGGGGACAGatat                                       >  2:100357/1‑251 (MQ=255)
                                                                                                                                                                                                        gCGTGTGTGTATTCCACTTCACCAAAGCCTATATCCAGCATCGCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTGGGCGACAATTGCATCACGCAGCGGCTGATAATGTCCGGCATCTAAGAATGCGCGGCGTGCTTGCATCATTTCCGCGCTGTCGCCCGGATCACGAGACCGTTTATGCTGAACGGGCAGCAGATTGACATACCCTTCTTTCGCCATATCAAACCACTGTCGCTggggcacgatat                                       >  2:300149/1‑242 (MQ=255)
                                                                                                                                                                                                         cGTGTGTGTAATACCCTTCACCACAGCCAATATCCAGCATCGCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTGGGCGACAATTGCATCACGCAGCGGCTGATAATGTCCGGCATCTAAGAATGCGCGGCGTGCTTGCATCATTTCCGCGCTGTCGCCCGGATCACGAGACCGTTTATGCTGAACGGGCAGCAGATTGACATACCCTTCTTTCGCCATATCAAACTGATGTCGCTGGGGACAGatat                                       >  2:400184/1‑249 (MQ=255)
                                                                                                                                                                                                           tgtgtgtAATACCCTTCACCACAGCCAATATCCAGCATCGCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTGGGCGACAATTGCATCACGCAGCGGCTGATAATGTCCGGCATCTAAGAATGCGCGGCGTGCTTGCATCATTTCCGCGCTGTCGCCCGGATCACGAGACCGTTTATGCTGAACGGGCAGCAGATTGACATACCCTTCTTTCGCCATATCAAACTGATGTCGCTGGGGACAGATATAGCt                                   >  2:413196/1‑251 (MQ=255)
                                                                                                                                                                                                            gtgtgtAATACCCTTCACCACAGCCAATATCCAGCATCGCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTGGGCGACAATTGCATCACGCAGCGGCTGATAATGTCCGGCATCTAAGAATGCGCGGCGTGCTTGCATCATTTCCGCGCTGTCGCCCGGATCACGAGACCGTTTATGCTGAACGGGCAGCAGATTGACATACCCTTCTTTCGCCATATCAAACTGATGTCGCTGGGGACAGATATAGCt                                   <  1:336173/250‑1 (MQ=255)
                                                                                                                                                                                                                         ttCACCACAGCCAATATCCAGCATCGCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTGGGCGACAATTGCATCACGCAGCGGCTGATAATGTCCGGCATCTAAGAATGCGCGGCGTGCTTGCATCATTTCCGCGCTGTCGCCCGGATCACGAGACCGTTTATGCTGAACGGGCAGCAGATTGACATACCCTTCTTTCGCCATATCAAACTGATGTCGCTGGGGACAGATATAGCTGTTTTTTTCAcgcg                     >  1:149378/1‑251 (MQ=255)
                                                                                                                                                                                                                            accacAGCCAATATCCAGCATCGCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTGGGCGACAATTGCATCACGCAGCGGCTGATAATGTCCGGCATCTAAGAATGCGCGGCGTGCTTGCATCATTTCCGCGCTGTCGCCCGGATCACGAGACCGTTTATGCTGAACGGGCAGCAGATTGACATCCCCTTCTTTCGCCATATCAAACTGATGTCGCTGGGGACAGATATAGCTGTTTTTTTCACGCGaaa                  <  2:464768/251‑1 (MQ=255)
                                                                                                                                                                                                                                  gcccatatCCAGCATCGCCGTGGCCGTATCATCAAGCCGTTCCCTCAGTTGGGCTACAATTGCATCACGCAGCGGCTGATAATGTCCGGCATCTAAGAATGCGCGGCGTGCTTGCATCATTTCCGCGCTGTCGCCCGGATCACGAGACCGTTTATGCTGAACGGGCAGCAGATTGACATACCCTTCTTTCTCCATATCAAACTGATGTCGCTGGGGACAGATATAGCTTTTTTTTTCACTCCAAAGAGGCt            <  2:335882/247‑1 (MQ=255)
                                                                                                                                                                                                                                  gCCAATATCCAGCATCGCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTGGGCGACAATTGCATCACGCAGCGGCTGATAATGTCCGGCATCTAAGAATGCGCGGCGTGCTTGCATCATTTCCGCGCTGTCGCCCGGATCACGAGACCGTTTATGCTGAACGGGCAGCAGATTGACATACCCTTCTTTCGCCATATCAAACTGATGTCGCTGGGGACAg                                           <  2:396988/220‑1 (MQ=255)
                                                                                                                                                                                                                                  gCCAATATCCAGCATCGCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTGGGCGACAATTGCATCACGCAGCGGCTGATAATGTCCGGCATCTAAGAATGCGCGGCGTGCTTGCATCATTTCCGCGCTGTCGCCCGGATCACGAGACCGTTTATGCTGAACGGGCAGCAGATTGACATACCCTTCTTTCGCCATATCAAACTGATGTCGCTGGGGACAg                                           >  1:396988/1‑220 (MQ=255)
                                                                                                                                                                                                                                     aaTATCCAGCATCGCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTGGGCGACAATTGCATCACGCAGCGGCTGATAATGTCCGGCATCTAAGAATGCGCGGCGTGCTTGCATCATTTCCGCGCTGTCGCCCGGATCACGAGACCGTTTATGCTGAACGGGCAGCAGATTGACATACCCTTCTTTCGCCATATCAAACTGATGTCGCTGGGGACAGATATAGCTGTTTTTTTCACGCGAAAGAGGCTGAt         >  1:404303/1‑251 (MQ=255)
                                                                                                                                                                                                                                      atatCCAGCATCGCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTGGGCGACAATTGCATCACGCAGCGGCTGATAATGTCCGGCATCTAAGAATGCGCGGCGTGCTTGCATCATTTCCGCGCTGTCGCCCGGATCACGAGACCGTTTATGCTGAACGGGCAGCAGATTGACATACCCTTCTTTCGCCATATCAAACTGATGTCGCTGGGGACAGATATAGCTGTTTTTTTCACGCGAAAGAGGCTGATg        <  1:538701/251‑1 (MQ=255)
                                                                                                                                                                                                                                        atCCAGCATCGCCGTGGCCTTATCATCAAGCCGTTCCCTCAGGTGGGCGACAATTGCATCACGCAGCGGCTGATAATGTCCGGCATCTAAGAATGCGCGGCGTGCTTGCATCATTTCCGCGCTGTCGCCCGGATCACGAGACCGTTTATGCTGAACGGGCAGCAGATTGACATACCCTTCTTTCGCCATATCAAACTGATGGCGCTGGGGACAGATATAGCTGTTTTTTTCACGCGAAAGAGGCTGATGGc      <  2:432229/251‑1 (MQ=255)
                                                                                                                                                                                                                                          ccAGCATCGCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTGGGCGACAATTGCATCACGCAGCGGCTGATAATGTCCGGCATCTAAGAATGCGCGGCGTGCTTGCATCATTTCCGCGCTGTCGCCCGGATCACGAGACCGTTTATGCTGAACGGGCAGCAGATTGACATACCCTTCTTTCGCCATATCAAACTGATGTCGCTGGGGACAGATATAGCTGTTTTTTTCACGCGAAAGAGGCTGATGGCa     >  1:70550/1‑250 (MQ=255)
                                                                                                                                                                                                                                            aGCATCGCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTGGGCGACAATTGCATCACGCAGCGGCTGATAATGTCCGGCATCTAAGAATGCGCGGCGTGCTTGCATCATTTCCGCGCTGTCGCCCGGATCACGAGACCGTTTATGCTGAACGGGCAGCAGATTGACATACCCTTCTTTCGCCATATCAAACTGATGTCGCTGGGGACAGATATAGCTTTTTTTTTCACTCGAAAGAGGCTGATGGCAAAg  <  1:464573/251‑1 (MQ=255)
                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                      
TCACTACTCGTGCTAATTCTTCTGCTTTACACGGCGCGTAAATACGTATTATGGCGTCCATACTGGTATCGGAAAACGGCAAACGGTGGCTGGAAGCGACACAAAAAGTGACCTGCGGATAGCGTTTCGCCGCCGCTTTTATCGCTACCTTCGAAACATCCAGACCAAACGTGGTGATTTCGGGCAACGCATCGGCAAATGCGTGTGTGTAATACCCTTCACCACAGCCAATATCCAGCACCGCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTGGGCGACAATTGCATCACGCAGCGGCTGATAATGTCCGGCATCTAAGAATGCGCGGCGTGCTTGCATCATTTCCGCGCTGTCGCCCGGATCACGAGACCGTTTATGCTGAACGGGCAGCAGATTGACATACCCTTCTTTCGCCATATCAAACTGATGTCGCTGGGGACAGATATAGCTGTTTTTTTCACGCGAAAGAGGCTGATGGCAAAG  >  NC_000913/1906556‑1907042

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