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,7960CT98.0% 151.3 / NA 49V89M (GTG→ATG) rlmA23S rRNA m(1)G745 methyltransferase
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base T (20/28);  total (21/28)
Fisher's exact test for biased strand distribution p-value = 4.29e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.59e-01

GGGCTTCACTACTCGTGCTAATTCTTCTGCTTTACACGGCGCGTAAATACGTATTATGGCGTCCATACTGGTATCGGAAAACGGCAAACGGTGGCTGGAAGCGACACAAAAAGTGACCTGCGGATAGCGTTTCGCCGCCGCTTTTATCGCTACCTTCGAAACATCCAGACCAAACGTGGTGATTTCGGGCAACGCATCGGCAAATGCGTGTGTGTAATACCCTTCACCACAGCCAATATCCAGCACCGCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTGGGCGACAATTGCATCACGCAGCGGCTGATAATGTCCGGCATCTAAGAATGCGCGGCGTGCTTGCATCATTTCCGCGCTGTCGCCCGGATCACGAGACCGTTTATGCTGAACGGGCAGCAGATTGACATACCCTTCTTTCGCCATATCAAACTGATGTCGCTGGGGACAGATATAGCTGTTTTTTTCACGCGAAAGAGGCTGATG  >  NC_000913/1906551‑1907036
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                
gggCTTCACTACTCGTGCTAATTCTTCTGCTTTACACGGCGCGTAAATACGTATTATGGCGTCCATACTGGTATCGGAAAACGGCAAACGGTGGCTGGAAGCGACACAAAAAGTGACCTGCGGATAGCGTTTCGCCGCCGCTTTTATCGCTACCTTCGAAACATCCAGACCAAACGTGGTGATTTCGGGCAACGCATCGGCAAATGCGTGTGTGTAATACCCTTCACCACAGCCAATATCCAGCAtcgccg                                                                                                                                                                                                                                             >  1:552721/1‑251 (MQ=255)
 ggattCACTACTCGTGCTAATTCTTCTGCGTTACACGGCGCGTAAATACGTATTATGGCGTCCATACTGGTACCGGAAAACGGCAAACGGTGGCTGGAAGCGACACAAAAAGTGACCTGCGGATAGCGTTTCGCCGCCGCTTTTATCGCTACCTTCGAAACATCCAGACCAAACGTGGTGATTTCGGGCAACGCATCGGCAAATGCGTGTGTGTAATACCCTTCACCACAGCCAATATCCAGCATCGCCGt                                                                                                                                                                                                                                            <  1:432966/248‑1 (MQ=255)
      cACTACTCGTGCTAATTCTTCTGCTTTACACGGCGCGTAAATACGTATTATGGCGTCCATACTGGTATCGGAAAACGGCAAACGGTGGCTGGAAGCGACACAAAAAGTGACCTGCGGATAGCGTTTCGCCGCCGCTTTTATCGCTACCTTCGAAACATCCAGACCAAACGTGGTGATTTCGGGCAACGCATCGGCAAATGCGTGTGTGTAATACCCTTCACCACAGCCAATATCCAGCATCGCCGTGGCCt                                                                                                                                                                                                                                       >  2:463061/1‑251 (MQ=255)
       cctactCGTGCTAATTCTTCTGCTTTACACGGCGCGTAAATACGTATTATGGCGTCCATACTGGTATCGGAAAACGGCAAACGGTGGCTGGAAGCGACACAAAAAGTGACCTGCGGATAGCGTTTCGCCGCCGCTTTTATCGCTACCTTCGAAACATCCAGACCAAACGTGGTGATTTCGGGCAACGCATCGGCAAATGCGTGTGTGTAATACCCTTCACCACAGCCAATATCCAGCATCGCCGTGGCCtt                                                                                                                                                                                                                                      <  2:486806/250‑1 (MQ=255)
             cGTGCTAATTCTTCTGCTTTACACGGCGCGTAAATACGTATTATGGCGTCCATACTGGTATCGGAAAACGGCAAACGGTGGCGGGAAGCGACACAAAAAGTGACCTGCGGATAGCGTTTCGCCGCCGCTTTTATCGCTACCTTCGAAACATCCAGACCAAACGTGGTGATTTCGGGCAACGCATCGGCAAATGCGTGTGTGTAATACCCTTCACCACAGCCAATATCCAGCATCGCCGTGGCCTTatcatc                                                                                                                                                                                                                                <  1:585635/251‑1 (MQ=255)
                         tctGCTTTACACGGCGCGTAAATACGTATTATGGCGTCCATACTGGTATCGGAAAACGGCAAACGGTGGCTGGAAGCGACACAAAAAGTGACCTGCGGATAGCGTTTCGCCGCCGCTTTTATCGCTACCTTCGAAACATCCAGACCAAACGTGGTGATTTCGGGCAACGCATCGGCAAATGCGTGTGTGTAATACCCTTCACCACAGCCAATATCCAGCATCGCCGTGGCCTTATCATCAAGCCGTTccc                                                                                                                                                                                                                     <  1:481125/250‑1 (MQ=255)
                         tctGCTTTACACGGCGCGTAAATACGTATTATGGCGTCCATACTGGTATCGGAAAACGGCAAACGGTGGCTGGAAGCGACACAAAAAGTGACCTGCGGATAGCGTTTCGCCGCCGCTTTTATCGCTACCTTCGAAACATCCAGACCAAACGTGGTGATTTCGGGCAACGCATCGGCAAATGCGTGTGTGTAATACCCTTCACCACAGCCAATATCCAGCATCGCCGTGGCCTTATCATCAAGCCGTTc                                                                                                                                                                                                                       <  1:421754/248‑1 (MQ=255)
                                      gcgcGTAAATACGTATTATGGCGTCCATACTGGTATCGGAAAACGGCAAACGTTGGCTGGAAGCGACACAAAAAGTGACCGGCGGATAGCGTTTCGCCGCCGCTTTTATCGCTACCTTCGAAACATCCAGACCAAACGTGGTGATTTCGGGCAACGCATCGGCAAATGCGTGTGTGTAATACCCTTCACCACAGCCAATATCCAGCATCGCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTGGGCGACa                                                                                                                                                                                                       <  2:507022/251‑1 (MQ=255)
                                          gTAAATACGTATTATGGCGTCCATACTGGTATCGGAAAACGGCAAACGGTGGCTGGAAGCGACACAAAAAGTGACCTGCGGATAGCGTTTCGCCGCCGCTTTTATCGCTACCTTCGAAACATCCAGACCAAACGTGGTGATTTCGGGCAACGCATCGGCAAATGCGTGTGTGTAATACCCTTCACCACAGCCAATATCCAGCATCGCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTGGGCGACAATTg                                                                                                                                                                                                   >  2:411431/1‑251 (MQ=255)
                                              aTACGTATTATGGCGTCCATACTGGTATCGGAAAACGGCAAACGGTGGCTGGAAGCGACCCAAAACGCGACAAGCGGATAGCGTTTCGCCGCCGCTTTTATCGCTACCTTCGCAACATCCAGACCAAACGTGGTGATTTCGGGCAACGCATCGGCACATGCGTGTGTGTAATACCCTTCACCACAGCCAATATCCAGCATCGCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTGGGCGACAATTGCATc                                                                                                                                                                                               >  1:358815/1‑251 (MQ=255)
                                              aTACGTATTATGGCGTCCATACTGGTATCGGAAAACGGCAAACGGTGGCTGGAAGCGACACAAAAAGTGACCTGCGGATAGCGTTTCGCCGCCGCTTTTATCGCTACCTTCGAAACATCCAGACCAAACGTGGTGATTTCGGGCAACGCATCGGCAAATGCGTGTGTGTAATACCCTTCACCACAGCCAATATCCAGCATCGCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTGGGCGACAATTGCATc                                                                                                                                                                                               >  2:290784/1‑251 (MQ=255)
                                                      tatGGCGTCCATAGTGGTATCGGAGAACGGCGAACGTTGGTGGGAGTCGCCACAAAAAGTGCCCGGCGGATAGCGTTGCGCCGCCGCTTTTATCGCTACCTTCGAAACATCCAGACCTACCGTGTTGATTTCGGGCAACGCATCGGCAAATGCGTTTGTTTATTACCCTTCGCCACATCCAATATCCGGTATCGCCGTGGCCTTCTCTTTAAGCCGGTCTCACAGTTTTGCTACAGTTGCATCACGCAGc                                                                                                                                                                                        <  2:273039/250‑1 (MQ=255)
                                                       atGGCGTCCATACTGGCATCGGAAAACGGCAAACGGTGGCTGGAAGCGACACACAAAGTGACCTGCGGATAGCGTTTCGCCGCCGCTTTTATCGCTACCTTCGAAACATCCGGACCAAACGTGGTGATTTCGGGCAACGCATCGGCAAATGCGTGTGTGTAATACCCTTAACCACAGCCAATATCCAGCATCGCCGTGGCCTTATCATCAAGCCGTTCCCTCCGTTGGGCGACAATTGCATCACGCAGCgg                                                                                                                                                                                      >  2:565021/1‑251 (MQ=255)
                                                               cATACTGGTATCGGAAAACGGCAAACGGTGGCTGGAAGCGACACAAAAAGTGACCTGCGGATAGCGTTTCGCCGCCGCTTTTATCGCTACCTTCGAAACATCCAGACCAAACGTGGTGATTTCGGGCAACGCATCGGCAAATGCGTGTGTGTAATACCCTTCACCACAGCCAATATCCAGCATCGCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTGGGCGACAATTGCATCACGCAGCGGCTGATAAt                                                                                                                                                                              <  1:290784/251‑1 (MQ=255)
                                                                  aCTGGTATCGGAAAACGGCAAACGGTGGCTGGAAGCGACACAAAAAGTGACCTGCGGATAGCGTTTCGCCTCCGCTTTTATCGCTACCTTCGAAACATCCAGACCAAACGTGGTGATTTCGGGCAACGCATCGGCAAATGCGTGTGTGTAATACCCTTCACCACAGCCAATATCCAGCATCGCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTGGGCGACAATTGCATCACGCAGCGGCTGATAATGt                                                                                                                                                                            <  1:31339/250‑1 (MQ=255)
                                                                     ggTATCGGAAAACGGCAAACGGTGGCTGGAAGCGACACAAAAAGTGACCTGCGGATAGCGTTTCGCCGCCGCTTTTATCGCTACCTTCGAAACATCCAGACCAAACGTGGTGATTTCGGGCAACGCATCGGCAAATGCGTGTGTGTAATACCCTTCACCACAGCCAATATCCAGCATCGCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTGGGCGACAATTGCATCACGCAGCGGCTGATAATGTCCgg                                                                                                                                                                        <  1:204889/251‑1 (MQ=255)
                                                                      gTATCGGAAAACGGCAAACGGTGGCTGGAAGCGACACAAAAAGTGACCTGCGGATAGCGTTTCGCCGCCGCTTTTATCGCTACCTTCGAAACATCCAGACCAAACGTGGTGATTTCGGGCAACGCATCGGCAAATGCGTGTGTGTAATACCCTTCACCACAGCCAATATCCAGCATCGCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTGGGCGACAATTGCATCACGCa                                                                                                                                                                                          >  2:148965/1‑232 (MQ=255)
                                                                      gTATCGGAAAACGGCAAACGGTGGCTGGAAGCGACACAAAAAGTGACCTGCGGATAGCGTTTCGCCGCCGCTTTTATCGCTACCTTCGAAACATCCAGACCAAACGTGGTGATTTCGGGCAACGCATCGGCAAATGCGTGTGTGTAATACCCTTCACCACAGCCAATATCCAGCATCGCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTGGGCGACAATTGCATCACGCa                                                                                                                                                                                          <  1:148965/232‑1 (MQ=255)
                                                                           ggAAAACGGCAAACGGTGGCGGGAAGCGACACAAAAAGTGACCTGCGGATAGCGTTTCGCCGCCGCTTTTATCGCTACCTTCGAAACATCCAGACCAAACGTGGTGATTTCGGGCAACGCATCGGCAAATGCGTGTGTGTAATACCCTTCACCACAGCCAATATCCAGCATCGCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTGGGCGACAATTGCATCACGCAGCGGCTGATAATGTCCGGCATCTa                                                                                                                                                                  <  1:483619/251‑1 (MQ=255)
                                                                                   gCAAACGGTGGCTGGAAGCGACACAAAAAGTGACCTGCGGATAGCGTTTCGCCGCCGCTTTTATCGCTACCTTCGAAACATCCAGACCAAACGTGGTGATTTCGGGCAACGCATCGGCAAATGCGTGTGTGTAATACCCTTCACCACAGCCAATATCCAGCATCGCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTGGGCGACAATTGCATCACGCAGCGGCTGATAATGTCCGGCATCTAAGAATgcg                                                                                                                                                          <  1:463061/251‑1 (MQ=255)
                                                                                        cGGTGGCTGGAAGCGACACAAAAAGTGACCTGCGGATAGCGTTTCGCCGCCGCTTTTATCGCTACCTTCGAAACATCCAGACCAAACGTGGTGATTTCGGGCAACGCATCGGCAAATGCGTGTGTGTAATACCCTTCTCCACAGCCAATATCCAGCATCGCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTGGGCGACAATTGCATCACGCAGCGGCTGATAATGTCCGGCATCTAAGAATGCgcggcg                                                                                                                                                     <  1:257943/251‑1 (MQ=255)
                                                                                             gCTGGACGCGACACAAAAAGTGACCTGCGGATAGCGTTTCGCCGCCGCTTTTATCGCTACCTTCGAAACATCCAGACCAAACGTGGTGATTTCGGGCAACGCATCGGCAAATGCGTGTGTGTAATACCCTTCACCACAGCCAATATCCAGCATCGCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTGGGCGACAATTGCATCACGCAGCGGCTGATAATGTCCGGCATCTAAGAATGCGCGGCGTGCtt                                                                                                                                                >  1:107371/1‑251 (MQ=255)
                                                                                             gCTGGAAGAGAAACAAAAAGTGACCTGGGGATCGCGTTTCGCCGCCGCTTTTATCGCTACCTTCGAAACATCCAGACCAAACGTGGTGATTTCGGGCAACGCATCGGCAAATGCGTGTGTGTAATACCCTTCACCACAGCCAATATCCAGCATCGCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTGGGCGACAATTGCATCACGCAGCGGCTGATAATGTCCGGCATCTAAGAATGCGCGGCGTGCtt                                                                                                                                                <  2:452350/251‑1 (MQ=255)
                                                                                                ggAAGCGACACAAAAAGTAAACTGCGGATAGCGTTTCGCCGCAGCTTTTATCGCTACCTTCGAAACATCCAGACCAAACCAGGTGATTTCGGGCAACGCATCGGCAAATGCGTGTGTGTAATACCCTTCACCACAGCCAATATAAAGCATCGCCGTGGCCTAATCATCAAGCCGTGCCCTAAGTTGAGAGACAATTGCATCAAGCAGCGGCCGATAATGTCCGGCATCTAAGAATGCGCGGCGTGCTTGc                                                                                                                                              >  2:218745/1‑250 (MQ=255)
                                                                                                                   gACCTGCGGATAGCGTTTCGCCGCCGCTTTTATCGCTACCTTCGAAACATCCAGACCAAACGTGGTGATTTCGGGCAACGCATCGGCAAATGCGTGTGTGTAATACCCTTCACCACAGCCAATATCCAGCATCGCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTGGGCGACAATTGCATCACGCAGCGGCTGATAATGTCCGGCATCTAAGAATGCGCGGCGTGCTTGCATCATTTCCGCGCTGTCGc                                                                                                                           <  1:555920/251‑1 (MQ=255)
                                                                                                                          ggATAGCGTTTAGCCGCCGCTTTTATCGCTACCTTCGAAACATCCAGACCAAACGTGGTGATTTCGGGCAACGCATCGGCAAATGCGTGTGTGTAATACCCTTCACCACAGCCAATATCCAGCATCGCCGTGGCCTTATCATCAAGCAGTTCCCTCAGTTGGGCGACAATTGCATCACGCAGCGGCTGATAATGTCCGCCATCTAAGAATGCGCGGCGTGCTTGCATCATTTCCGCGCTGTCGCCCGGATc                                                                                                                    >  2:68502/1‑251 (MQ=255)
                                                                                                                           gATAGCGTTTCGCCGCCGCTTTTATCGCTACCTTCGAAACATCCAGACCAAACGTGGTGATTTCGGGCAACGCATCGGCAAATGCGTGTGTGTAATACCCTTCACCACAGCCAATATCCAGCATCGCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTTGGCGACAATTGCATCACGCAGCGGCTGATAATGTCCGGCATCTAAGAATGCGCGGCGTGCTTGCATCATTTCCGCGCTGTCGCCCGGATCa                                                                                                                   <  2:333537/251‑1 (MQ=255)
                                                                                                                            aTAGCGTTTCGCCGCCGCTTTTATCGCTACCTTCGAAACATCCAGACCAAACGTGGTGATTTCGGGCAACGCATCGGCAAATGCGTGTGTGTAATACCCTTCACCACAGCCAATATCCAGCATCGCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTGGGCGACAATTGCATCACGCAGCGGCTGATAATGTCCGGCATCTAAGAATGCGCGGCGTGCTTGCATCATTTCCGCGCTGTCGCCCGGATCAc                                                                                                                  <  1:157064/251‑1 (MQ=255)
                                                                                                                              aGCGTTTCGCCGCCGCTTTTATCGCTACCTTCGAACCATCCAGACCAAACGTGGTGATTTCGGGCAACGCATCGGCAAATGCGTGTGTGTAATCCCCTTCACCACAGCCAATATCCAGCATCGCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTGGGCGACAATTGCATCACGCAGCGGCTGATAATGTCCGGCATCTAAGAATGCGCGGCGTGCTTGCATCATTTCCGCGCTGTCGCCCGGATCACga                                                                                                                <  2:477610/251‑1 (MQ=255)
                                                                                                                                cGTTTAGCCGCCGATTTTATCGCTACCTTCGAAACATCCAGACCAAACGTGGTGATTTCGGGCAACGCATCGGCAAATGCGTGTGTGTAATACCCTTCACCACAGCAAATATCCAGCATCGCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTGGGCGACAATTGCATCACGCAGCGGCTGATAATGTCCGGCATCTAAGAATGCGCGGCGTGCTTGCATCATTTCCGTGCTGTCGCCCGGATCACg                                                                                                                 >  2:389938/1‑248 (MQ=255)
                                                                                                                                   ttCGCCGCCGCTTTTATCGCTACCTTCGAAACATCCAGACCAAACGTGGTGATTTCGGGCAACGCATCGGCAAATGCGTGTGTGTAATACCCTTCACCACAGCCAATATCCAGCATCGCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTGGGCGACAATTGCATCACGCAGCGGCTGATAATGTCCGGCATCTAAGAATGCGCGGCGTGCTTGCATCATTTCCGCGCTGTCGCCCGGATCACGAGAcc                                                                                                            <  1:141919/250‑1 (MQ=255)
                                                                                                                                    tCGCCGCCGCTTTTATCGCTACCTTCGAAACATCCAGACCAAACGTGGTGATTTCGGGCAACGCATCGGCAAATGCGTGTGTGTAATACCCTTCACCACAGCCAATATCCAGCATCGCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTGGGCGACAATTGCATCACGCAGCGGCTGATAATGTCCGGCATCTAAGAATGCGCGGCGTGCTTGCATCATTTCCGCGCTGTCGCCCGGATCACGAGACCg                                                                                                           <  1:105983/250‑1 (MQ=255)
                                                                                                                                             cTTTTATCGCTACCTTCGAAACATCCAGACCAAACGTGGTGATTTCGGGCAACGCATCGGCAAATGCGTGTGTGTAATACCCTTCACCACAGCCAATATCCAGCATCGCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTGGGCGACAATTGCATCACGCAGCGGCTGATAATGTCCGGCATCTAAGAATGCGCGGCGTGCTTGCATCATTTCCGCGCTGTCGCCCGGATCACGAGACCGTTTATGCTg                                                                                                  >  1:89342/1‑250 (MQ=255)
                                                                                                                                                  aTCGCTACCTTCGAAACATCCAGACCAAACGTGGTGATTTCGGGCAACGCATCGGCAAATGCGTGTGTGTAATACCCTTCACCACAGCCAATATCCAGCATCGCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTGGGCGACAATTGCATCACGCAGCGGCTGATAATGTCCGGCATCTAAGAATGCGCGGCGTGCTTGCATCATTTCCGCGCTGTCGCCCGGATCACGAGACCGTTTATGCTGAACggg                                                                                            >  2:11841/1‑251 (MQ=255)
                                                                                                                                                      cTACCTTCGAAACATCCAGACCAAACGTGGTGATTTCGGGCAACGCATCGGCAAATGCGTGTGTGTAATACCCTTCACCACAGCCAATATCCAGCATCGCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTGGGCGACAATTGCATCACGCAGCGGCTGATAATGTCCGGCATCTAAGAATGCGCGGCGTGCTTGCATCATTTCCGCGCTGTCGCCCGGATCACGAGACCGTTTATGCTGAACGGgcagc                                                                                        >  1:474627/1‑251 (MQ=255)
                                                                                                                                                      cTACCTTCGAAACATACAGCACCAACGTAGTGATTTCGGGCAACGCATCGGCAAATGCGTGTGTGTAATACCCTTCACCACAGCCAATATCCAGCATCGCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTGGGCGACAATTGCATCACGCAGCGGCTGATAATGTCCGGCATCTAAGAATGCGCGGCGTGCTTGCATCATTTACGCGCTGTCGCCCGGATCACGAGACCGTTTCTGCTGAACGGgcagc                                                                                        >  2:357967/1‑251 (MQ=255)
                                                                                                                                                               aaaCATCCAGACCAAACGTGGTGATTTCGGGCAACGCATCGGCAAATGCGTGTGTGTAATACCCTTCACCACAGCCAATATCCAGCATCGCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTGGGCGACAATTGCATCACGCAGCGGCTGATAATGTCCGGCATCTAAGAATGCGCGGCGTGCTTGCATCATTTCCGCGCTGTCGCCCGGATCACGAGACCGTTTATGCTGAACGGGCAGCAGATTGAc                                                                                >  1:439958/1‑250 (MQ=255)
                                                                                                                                                                          aCCAAACGTGGTGATTGCGGGCAACGCATCGGCAAATGCGTGTGTGTAATCCCCTTCACCACAGCCAATAGCCAGCATCGCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTGGGCGACAAGTGCATCACGCAGCGGCTGATAATGTCCGGCATCTAAGAATGCGCGGCGTGCTTGCATCATTTCCGCGCTGTCGCCCGGATCACGAGACCGTTTATGCTGAACGGGCAGCAGATTGACATACCcttctt                                                                     <  1:68502/251‑1 (MQ=255)
                                                                                                                                                                            cAAACGTGGTGATTTCGGGCAACGCATCGGCAAATGCGTGGGTGTAATACCCTTCACCACAGCCAATATCCAGCATCGCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTGGGCGACAATTGCATCACGCAGCGGCTGATAATGTCCGGCATCTAAGAATGCGCGGCGTGCTTGCATCATTTCCGCGCTGTCGCCCGGAGCACGAGACCGTTTATGCTGAACGGGCAGCAGATTGACATACCCGTCTTTc                                                                   <  1:391972/251‑1 (MQ=255)
                                                                                                                                                                                 gtggtgATTTCGGGCAACGCATCGGCAAATGCGTGTGTGTAATACCCTTCACCACAGCCAATATCCAGCATCGCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTGGGCGACAATTGCATCACGCAGCGGCTGATAATGTCCGGCATCTAAGAATGCGCGGCGTGCTTGCATCATTTCCGCGCTGTCGCCCGGATCACGAGACCGTTTATGCTGAACGGGCAGCAGATTGACATACCCTTCTTTCGCCat                                                              <  1:197152/251‑1 (MQ=255)
                                                                                                                                                                                           cGGGCAACGCATAGGCAAATGCGGGTGAGTAATACCCTTCACCACAGCAAATATCCAGCATCGCAGTGGCCTTATCATCAAGCCGTTCCCTCAGTTGGGCGACAATTGCATCACGCAGCGGCTGATAATGTCCGGCATCTAAGAATGCGCGGCGTGCTTGCATCATTTCCGCGCTGTCGCCCGGATCACGAGACCGTTTATGCTGAACGGGCAGCAGATTGACATCCCCTTCTTTCGCCATATCAAACTgc                                                    >  2:306529/1‑250 (MQ=255)
                                                                                                                                                                                            gggCAACGCATCGGCAAATGCGTGTGTGTAATCCCCTTCACCACAGCCAATATCCAGCATCGCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTGGGCGACAATTGCATCACGCAGCGGCTGATAATGTCCGGCATCTAAGAATGCGCGGCGTGCTTGCATCATTTCCGCGCTGTCGCCCGGATCACGAGACCGTTTATGCTGAACGGGCAGCAGATTGACATACCCTTCTTTCGCCATATCAAACTGAt                                                   <  1:33987/251‑1 (MQ=255)
                                                                                                                                                                                                             aTGCGTGTGTGTAATCCCCTTCACCACAGCCAATATCCAGCATCGCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTGGGCGACAATTGCATCACGCAGCGGCTGATAATGTCCGGCATCTAAGAATGCGCGGCGTGCTTGCATCATTTCCGCGCTGTCGCCCGGATCACGAGACCGTTTATGCTGAACGGGCAGCAGATTGACATACCCTTCTTTCGCCATATCAAACTGATGTCGCTGGGGACAGata                                  <  2:89342/251‑1 (MQ=255)
                                                                                                                                                                                                                 gtgtgtgtAATACCCTTCACCACAGCCAATATCCAGCATCGCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTGGGCGACAATTGCATCACGCAGCGGCTGATAATGTCCGGCATCTAAGAATGCGCGGCGTGCTTGCATCATTTCCGCGCTGTCGCCCGGATCACGAGACCGTTTATGCTGAACGGGCAGCAGATTGACATACCCTTCTTTCGCCATATCAAACTGATGTCGCTGGGGACAGATATAGc                              >  1:171427/1‑251 (MQ=255)
                                                                                                                                                                                                                            aCCCTTCACCACAGACAATATCCAGCATCGCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTAGGCGACAATTGCATCACGCAGCGGCTGATACTGTCCGGCATCTAAGAATGCGCGGCGTGCTTGCATCATTTCCGCGCTGTCGCACGGATCACGAGACCGTTTATGCTGAACGGGCAGCAGATTGACATACCCTTATTTCGCCATATCAAACTGATGTCGCTGGGGAAAGAGATAGCTGTTTTTTTcc                   >  2:76813/1‑250 (MQ=255)
                                                                                                                                                                                                                                   accacAGCCAATATCCAGCATCGCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTGGGCGACAATTGCATCACGCAGCGGCTGATAATGTCCGGCATCTAAGAATGCGCGGCGTGCTTGCATCATTTCCGCGCTGTCGCCCGGATCACGAGACCGTTTATGCTGAACGGGCAGCAGATTGACATACCCTTCTTTCGCCATATCAAACTGATGTCGCTGGGGACAGATATAGCTGTTTTTTTCACGCGaaa            <  1:411431/251‑1 (MQ=255)
                                                                                                                                                                                                                                          ccAATATCCAGCATCGCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTGGGCGACAATTGCATCACGCAGCGGCTGATAATGTCCGGCATCTAAGAATGCGCGGCGTGCTTGCATCATTTCCGCGCTGTCGCCCGGATCACGAGACCGTTTATGCTGAACGGGCAGCAGATTGACATACCCTTCTTTCGCCATATCAAACTGCTGTCGCTGGGGACAGATATAGCTGTTTTTTTCACGCGAAAGAggatg     >  2:128673/1‑248 (MQ=255)
                                                                                                                                                                                                                                          ccAATATCCAGCATCGCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTGGGCGACAATTGCATCACGCAGCGGCTGATAATGTCCGGCATCTAAGAATGCGCGGCGTGCTTGCATCATTTCCGCGCTGTCGCCCGGATCACGAGACCGTTTATGCTGAACGGGCAGCAGATTGACATACCCTTCTTTCGCCATATCAAACTGATGTCGCTGGGGACAGATATAGCTGTTTTTTTCACGCGAAAGAGGCTg     >  1:537123/1‑251 (MQ=255)
                                                                                                                                                                                                                                             atatCCAGCATCGCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTGGGCGACAATTGCATCACGCAGCGGCTGATAATGTCCGGCATCTAAGAATGCGCGGCGTGCTTGCATCATTTCCGCGCTGTCGCCCGGATCACGAGACCGTTTATGCTGAACGGGCAGCAGATTGACATACCCTTCTTTCGCCATATCAAACTGATGTCGCTGGGGACAGATATAGCTGTTTTTTTCACGCGAAAGAGGCTGATg  <  1:565021/251‑1 (MQ=255)
                                                                                                                                                                                                                                             atatCCAACACAGCCGTGGCATTATCATAACGTCGTTCCATCAGTTGGGCGAAACTTGCATCACGCAGCGGCTGATAATGTCCGGCATCTAAGAATGCGCGGAGTGCTTGCATCATTTCCGCGCTGTCGCCCGGATCACGAGACCGTTTATGCTGAACGGGCAGCAGATTGACATACCCTTCTTTCGCCATATCAAACTGATGTCGCTGGGGACGGATATAGCTGTTTTTTTCACGCGAAAGAGGCTGATg  >  2:364198/1‑251 (MQ=255)
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                
GGGCTTCACTACTCGTGCTAATTCTTCTGCTTTACACGGCGCGTAAATACGTATTATGGCGTCCATACTGGTATCGGAAAACGGCAAACGGTGGCTGGAAGCGACACAAAAAGTGACCTGCGGATAGCGTTTCGCCGCCGCTTTTATCGCTACCTTCGAAACATCCAGACCAAACGTGGTGATTTCGGGCAACGCATCGGCAAATGCGTGTGTGTAATACCCTTCACCACAGCCAATATCCAGCACCGCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTGGGCGACAATTGCATCACGCAGCGGCTGATAATGTCCGGCATCTAAGAATGCGCGGCGTGCTTGCATCATTTCCGCGCTGTCGCCCGGATCACGAGACCGTTTATGCTGAACGGGCAGCAGATTGACATACCCTTCTTTCGCCATATCAAACTGATGTCGCTGGGGACAGATATAGCTGTTTTTTTCACGCGAAAGAGGCTGATG  >  NC_000913/1906551‑1907036

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: