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

GGGCTTCACTACTCGTGCTAATTCTTCTGCTTTACACGGCGCGTAAATACGTATTATGGCGTCCATACTGGTATCGGAAAACGGCAAACGGTGGCTGGAAG‑CGACACAAAAAGTGACCTGCGGATAGCGTTTCGCCGCCGCTTTTATCGCTACCTTCGAAACATCCAG‑ACCAAACGTGGTGATTTCGGGCAACGCATCGGCAAATGCGTGTGTGTAATACCCTTCACCACAGCCAATATCCAGCACCGCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTGGGCGACAATTGCATCACGCAGCGGCTGATAATGTCCGGCATCTAAGAATGCGCGGCGTGCTTGCATCATTTCCGCGCTGTCGCCCGGATCACGAGACCGTTTATGCTGAACGGGCAGCAGATTGACATACCCTTCTTTCGCCATATCAAACTGATGTCGCTGGGGACAGATATAGCTGTTTTTTTCACGCGAAAGAGGCTGATG  >  NC_000913/1906551‑1907036
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                
gggCTTCACTACTCGTGCTAATTCTTCTGCTTTACACGGCGCGTAAATACGTATTATGGCGTCCATACTGGTATCGGAAAACGGCAAACGGTGGCTGGAAG‑CGACACAAAAAGTGACCTGCGGATAGCGTTTCGCCGCCGCTTTTATCGCTACCTTCGAAACATCCAG‑ACCAAACGTGGTGATTTCGGGCAACGCATCGGCAAATGCGTGTGTGTAATACCCTTCACCACAGCCAATATCCAGCAtcgccg                                                                                                                                                                                                                                             >  1:552721/1‑251 (MQ=255)
 ggattCACTACTCGTGCTAATTCTTCTGCGTTACACGGCGCGTAAATACGTATTATGGCGTCCATACTGGTACCGGAAAACGGCAAACGGTGGCTGGAAG‑CGACACAAAAAGTGACCTGCGGATAGCGTTTCGCCGCCGCTTTTATCGCTACCTTCGAAACATCCAG‑ACCAAACGTGGTGATTTCGGGCAACGCATCGGCAAATGCGTGTGTGTAATACCCTTCACCACAGCCAATATCCAGCATCGCCGt                                                                                                                                                                                                                                            <  1:432966/248‑1 (MQ=255)
      cACTACTCGTGCTAATTCTTCTGCTTTACACGGCGCGTAAATACGTATTATGGCGTCCATACTGGTATCGGAAAACGGCAAACGGTGGCTGGAAG‑CGACACAAAAAGTGACCTGCGGATAGCGTTTCGCCGCCGCTTTTATCGCTACCTTCGAAACATCCAG‑ACCAAACGTGGTGATTTCGGGCAACGCATCGGCAAATGCGTGTGTGTAATACCCTTCACCACAGCCAATATCCAGCATCGCCGTGGCCt                                                                                                                                                                                                                                       >  2:463061/1‑251 (MQ=255)
       cctactCGTGCTAATTCTTCTGCTTTACACGGCGCGTAAATACGTATTATGGCGTCCATACTGGTATCGGAAAACGGCAAACGGTGGCTGGAAG‑CGACACAAAAAGTGACCTGCGGATAGCGTTTCGCCGCCGCTTTTATCGCTACCTTCGAAACATCCAG‑ACCAAACGTGGTGATTTCGGGCAACGCATCGGCAAATGCGTGTGTGTAATACCCTTCACCACAGCCAATATCCAGCATCGCCGTGGCCtt                                                                                                                                                                                                                                      <  2:486806/250‑1 (MQ=255)
             cGTGCTAATTCTTCTGCTTTACACGGCGCGTAAATACGTATTATGGCGTCCATACTGGTATCGGAAAACGGCAAACGGTGGCGGGAAG‑CGACACAAAAAGTGACCTGCGGATAGCGTTTCGCCGCCGCTTTTATCGCTACCTTCGAAACATCCAG‑ACCAAACGTGGTGATTTCGGGCAACGCATCGGCAAATGCGTGTGTGTAATACCCTTCACCACAGCCAATATCCAGCATCGCCGTGGCCTTatcatc                                                                                                                                                                                                                                <  1:585635/251‑1 (MQ=255)
                         tctGCTTTACACGGCGCGTAAATACGTATTATGGCGTCCATACTGGTATCGGAAAACGGCAAACGGTGGCTGGAAG‑CGACACAAAAAGTGACCTGCGGATAGCGTTTCGCCGCCGCTTTTATCGCTACCTTCGAAACATCCAG‑ACCAAACGTGGTGATTTCGGGCAACGCATCGGCAAATGCGTGTGTGTAATACCCTTCACCACAGCCAATATCCAGCATCGCCGTGGCCTTATCATCAAGCCGTTccc                                                                                                                                                                                                                     <  1:481125/250‑1 (MQ=255)
                         tctGCTTTACACGGCGCGTAAATACGTATTATGGCGTCCATACTGGTATCGGAAAACGGCAAACGGTGGCTGGAAG‑CGACACAAAAAGTGACCTGCGGATAGCGTTTCGCCGCCGCTTTTATCGCTACCTTCGAAACATCCAG‑ACCAAACGTGGTGATTTCGGGCAACGCATCGGCAAATGCGTGTGTGTAATACCCTTCACCACAGCCAATATCCAGCATCGCCGTGGCCTTATCATCAAGCCGTTc                                                                                                                                                                                                                       <  1:421754/248‑1 (MQ=255)
                                      gcgcGTAAATACGTATTATGGCGTCCATACTGGTATCGGAAAACGGCAAACGTTGGCTGGAAG‑CGACACAAAAAGTGACCGGCGGATAGCGTTTCGCCGCCGCTTTTATCGCTACCTTCGAAACATCCAG‑ACCAAACGTGGTGATTTCGGGCAACGCATCGGCAAATGCGTGTGTGTAATACCCTTCACCACAGCCAATATCCAGCATCGCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTGGGCGACa                                                                                                                                                                                                       <  2:507022/251‑1 (MQ=255)
                                          gTAAATACGTATTATGGCGTCCATACTGGTATCGGAAAACGGCAAACGGTGGCTGGAAG‑CGACACAAAAAGTGACCTGCGGATAGCGTTTCGCCGCCGCTTTTATCGCTACCTTCGAAACATCCAG‑ACCAAACGTGGTGATTTCGGGCAACGCATCGGCAAATGCGTGTGTGTAATACCCTTCACCACAGCCAATATCCAGCATCGCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTGGGCGACAATTg                                                                                                                                                                                                   >  2:411431/1‑251 (MQ=255)
                                              aTACGTATTATGGCGTCCATACTGGTATCGGAAAACGGCAAACGGTGGCTGGAAG‑CGACCCAAAACGCGACAAGCGGATAGCGTTTCGCCGCCGCTTTTATCGCTACCTTCGCAACATCCAG‑ACCAAACGTGGTGATTTCGGGCAACGCATCGGCACATGCGTGTGTGTAATACCCTTCACCACAGCCAATATCCAGCATCGCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTGGGCGACAATTGCATc                                                                                                                                                                                               >  1:358815/1‑251 (MQ=255)
                                              aTACGTATTATGGCGTCCATACTGGTATCGGAAAACGGCAAACGGTGGCTGGAAG‑CGACACAAAAAGTGACCTGCGGATAGCGTTTCGCCGCCGCTTTTATCGCTACCTTCGAAACATCCAG‑ACCAAACGTGGTGATTTCGGGCAACGCATCGGCAAATGCGTGTGTGTAATACCCTTCACCACAGCCAATATCCAGCATCGCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTGGGCGACAATTGCATc                                                                                                                                                                                               >  2:290784/1‑251 (MQ=255)
                                                      tatGGCGTCCATAGTGGTATCGGAGAACGGCGAACGTTGG‑TGGGAGTCGCCACAAAAAGTGCCCGGCGGATAGCGTTGCGCCGCCGCTTTTATCGCTACCTTCGAAACATCCAG‑ACCTACCGTGTTGATTTCGGGCAACGCATCGGCAAATGCGTTTGTTTATTACCCTTCGCCACATCCAATATCCGGTATCGCCGTGGCCTTCTCTTTAAGCCGGTCTCACAGTTTTGCTACAGTTGCATCACGCAGc                                                                                                                                                                                        <  2:273039/250‑1 (MQ=255)
                                                       atGGCGTCCATACTGGCATCGGAAAACGGCAAACGGTGGCTGGAAG‑CGACACACAAAGTGACCTGCGGATAGCGTTTCGCCGCCGCTTTTATCGCTACCTTCGAAACATCCGG‑ACCAAACGTGGTGATTTCGGGCAACGCATCGGCAAATGCGTGTGTGTAATACCCTTAACCACAGCCAATATCCAGCATCGCCGTGGCCTTATCATCAAGCCGTTCCCTCCGTTGGGCGACAATTGCATCACGCAGCgg                                                                                                                                                                                      >  2:565021/1‑251 (MQ=255)
                                                               cATACTGGTATCGGAAAACGGCAAACGGTGGCTGGAAG‑CGACACAAAAAGTGACCTGCGGATAGCGTTTCGCCGCCGCTTTTATCGCTACCTTCGAAACATCCAG‑ACCAAACGTGGTGATTTCGGGCAACGCATCGGCAAATGCGTGTGTGTAATACCCTTCACCACAGCCAATATCCAGCATCGCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTGGGCGACAATTGCATCACGCAGCGGCTGATAAt                                                                                                                                                                              <  1:290784/251‑1 (MQ=255)
                                                                  aCTGGTATCGGAAAACGGCAAACGGTGGCTGGAAG‑CGACACAAAAAGTGACCTGCGGATAGCGTTTCGCCTCCGCTTTTATCGCTACCTTCGAAACATCCAG‑ACCAAACGTGGTGATTTCGGGCAACGCATCGGCAAATGCGTGTGTGTAATACCCTTCACCACAGCCAATATCCAGCATCGCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTGGGCGACAATTGCATCACGCAGCGGCTGATAATGt                                                                                                                                                                            <  1:31339/250‑1 (MQ=255)
                                                                     ggTATCGGAAAACGGCAAACGGTGGCTGGAAG‑CGACACAAAAAGTGACCTGCGGATAGCGTTTCGCCGCCGCTTTTATCGCTACCTTCGAAACATCCAG‑ACCAAACGTGGTGATTTCGGGCAACGCATCGGCAAATGCGTGTGTGTAATACCCTTCACCACAGCCAATATCCAGCATCGCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTGGGCGACAATTGCATCACGCAGCGGCTGATAATGTCCgg                                                                                                                                                                        <  1:204889/251‑1 (MQ=255)
                                                                      gTATCGGAAAACGGCAAACGGTGGCTGGAAG‑CGACACAAAAAGTGACCTGCGGATAGCGTTTCGCCGCCGCTTTTATCGCTACCTTCGAAACATCCAG‑ACCAAACGTGGTGATTTCGGGCAACGCATCGGCAAATGCGTGTGTGTAATACCCTTCACCACAGCCAATATCCAGCATCGCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTGGGCGACAATTGCATCACGCa                                                                                                                                                                                          >  2:148965/1‑232 (MQ=255)
                                                                      gTATCGGAAAACGGCAAACGGTGGCTGGAAG‑CGACACAAAAAGTGACCTGCGGATAGCGTTTCGCCGCCGCTTTTATCGCTACCTTCGAAACATCCAG‑ACCAAACGTGGTGATTTCGGGCAACGCATCGGCAAATGCGTGTGTGTAATACCCTTCACCACAGCCAATATCCAGCATCGCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTGGGCGACAATTGCATCACGCa                                                                                                                                                                                          <  1:148965/232‑1 (MQ=255)
                                                                           ggAAAACGGCAAACGGTGGCGGGAAG‑CGACACAAAAAGTGACCTGCGGATAGCGTTTCGCCGCCGCTTTTATCGCTACCTTCGAAACATCCAG‑ACCAAACGTGGTGATTTCGGGCAACGCATCGGCAAATGCGTGTGTGTAATACCCTTCACCACAGCCAATATCCAGCATCGCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTGGGCGACAATTGCATCACGCAGCGGCTGATAATGTCCGGCATCTa                                                                                                                                                                  <  1:483619/251‑1 (MQ=255)
                                                                                   gCAAACGGTGGCTGGAAG‑CGACACAAAAAGTGACCTGCGGATAGCGTTTCGCCGCCGCTTTTATCGCTACCTTCGAAACATCCAG‑ACCAAACGTGGTGATTTCGGGCAACGCATCGGCAAATGCGTGTGTGTAATACCCTTCACCACAGCCAATATCCAGCATCGCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTGGGCGACAATTGCATCACGCAGCGGCTGATAATGTCCGGCATCTAAGAATgcg                                                                                                                                                          <  1:463061/251‑1 (MQ=255)
                                                                                        cGGTGGCTGGAAG‑CGACACAAAAAGTGACCTGCGGATAGCGTTTCGCCGCCGCTTTTATCGCTACCTTCGAAACATCCAG‑ACCAAACGTGGTGATTTCGGGCAACGCATCGGCAAATGCGTGTGTGTAATACCCTTCTCCACAGCCAATATCCAGCATCGCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTGGGCGACAATTGCATCACGCAGCGGCTGATAATGTCCGGCATCTAAGAATGCgcggcg                                                                                                                                                     <  1:257943/251‑1 (MQ=255)
                                                                                             gCTGGACG‑CGACACAAAAAGTGACCTGCGGATAGCGTTTCGCCGCCGCTTTTATCGCTACCTTCGAAACATCCAG‑ACCAAACGTGGTGATTTCGGGCAACGCATCGGCAAATGCGTGTGTGTAATACCCTTCACCACAGCCAATATCCAGCATCGCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTGGGCGACAATTGCATCACGCAGCGGCTGATAATGTCCGGCATCTAAGAATGCGCGGCGTGCtt                                                                                                                                                >  1:107371/1‑251 (MQ=255)
                                                                                             gCTGGAAG‑AGAAACAAAAAGTGACCTGGGGATCGCGTTTCGCCGCCGCTTTTATCGCTACCTTCGAAACATCCAG‑ACCAAACGTGGTGATTTCGGGCAACGCATCGGCAAATGCGTGTGTGTAATACCCTTCACCACAGCCAATATCCAGCATCGCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTGGGCGACAATTGCATCACGCAGCGGCTGATAATGTCCGGCATCTAAGAATGCGCGGCGTGCtt                                                                                                                                                <  2:452350/251‑1 (MQ=255)
                                                                                                ggAAG‑CGACACAAAAAGTAAACTGCGGATAGCGTTTCGCCGCAGCTTTTATCGCTACCTTCGAAACATCCAG‑ACCAAACCAGGTGATTTCGGGCAACGCATCGGCAAATGCGTGTGTGTAATACCCTTCACCACAGCCAATATAAAGCATCGCCGTGGCCTAATCATCAAGCCGTGCCCTAAGTTGAGAGACAATTGCATCAAGCAGCGGCCGATAATGTCCGGCATCTAAGAATGCGCGGCGTGCTTGc                                                                                                                                              >  2:218745/1‑250 (MQ=255)
                                                                                                                   gACCTGCGGATAGCGTTTCGCCGCCGCTTTTATCGCTACCTTCGAAACATCCAG‑ACCAAACGTGGTGATTTCGGGCAACGCATCGGCAAATGCGTGTGTGTAATACCCTTCACCACAGCCAATATCCAGCATCGCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTGGGCGACAATTGCATCACGCAGCGGCTGATAATGTCCGGCATCTAAGAATGCGCGGCGTGCTTGCATCATTTCCGCGCTGTCGc                                                                                                                           <  1:555920/251‑1 (MQ=255)
                                                                                                                          ggATAGCGTTTAGCCGCCGCTTTTATCGCTACCTTCGAAACATCCAG‑ACCAAACGTGGTGATTTCGGGCAACGCATCGGCAAATGCGTGTGTGTAATACCCTTCACCACAGCCAATATCCAGCATCGCCGTGGCCTTATCATCAAGCAGTTCCCTCAGTTGGGCGACAATTGCATCACGCAGCGGCTGATAATGTCCGCCATCTAAGAATGCGCGGCGTGCTTGCATCATTTCCGCGCTGTCGCCCGGATc                                                                                                                    >  2:68502/1‑251 (MQ=255)
                                                                                                                           gATAGCGTTTCGCCGCCGCTTTTATCGCTACCTTCGAAACATCCAG‑ACCAAACGTGGTGATTTCGGGCAACGCATCGGCAAATGCGTGTGTGTAATACCCTTCACCACAGCCAATATCCAGCATCGCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTTGGCGACAATTGCATCACGCAGCGGCTGATAATGTCCGGCATCTAAGAATGCGCGGCGTGCTTGCATCATTTCCGCGCTGTCGCCCGGATCa                                                                                                                   <  2:333537/251‑1 (MQ=255)
                                                                                                                            aTAGCGTTTCGCCGCCGCTTTTATCGCTACCTTCGAAACATCCAG‑ACCAAACGTGGTGATTTCGGGCAACGCATCGGCAAATGCGTGTGTGTAATACCCTTCACCACAGCCAATATCCAGCATCGCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTGGGCGACAATTGCATCACGCAGCGGCTGATAATGTCCGGCATCTAAGAATGCGCGGCGTGCTTGCATCATTTCCGCGCTGTCGCCCGGATCAc                                                                                                                  <  1:157064/251‑1 (MQ=255)
                                                                                                                              aGCGTTTCGCCGCCGCTTTTATCGCTACCTTCGAACCATCCAG‑ACCAAACGTGGTGATTTCGGGCAACGCATCGGCAAATGCGTGTGTGTAATCCCCTTCACCACAGCCAATATCCAGCATCGCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTGGGCGACAATTGCATCACGCAGCGGCTGATAATGTCCGGCATCTAAGAATGCGCGGCGTGCTTGCATCATTTCCGCGCTGTCGCCCGGATCACga                                                                                                                <  2:477610/251‑1 (MQ=255)
                                                                                                                                cGTTTAGCCGCCGATTTTATCGCTACCTTCGAAACATCCAG‑ACCAAACGTGGTGATTTCGGGCAACGCATCGGCAAATGCGTGTGTGTAATACCCTTCACCACAGCAAATATCCAGCATCGCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTGGGCGACAATTGCATCACGCAGCGGCTGATAATGTCCGGCATCTAAGAATGCGCGGCGTGCTTGCATCATTTCCGTGCTGTCGCCCGGATCACg                                                                                                                 >  2:389938/1‑248 (MQ=255)
                                                                                                                                   ttCGCCGCCGCTTTTATCGCTACCTTCGAAACATCCAG‑ACCAAACGTGGTGATTTCGGGCAACGCATCGGCAAATGCGTGTGTGTAATACCCTTCACCACAGCCAATATCCAGCATCGCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTGGGCGACAATTGCATCACGCAGCGGCTGATAATGTCCGGCATCTAAGAATGCGCGGCGTGCTTGCATCATTTCCGCGCTGTCGCCCGGATCACGAGAcc                                                                                                            <  1:141919/250‑1 (MQ=255)
                                                                                                                                    tCGCCGCCGCTTTTATCGCTACCTTCGAAACATCCAG‑ACCAAACGTGGTGATTTCGGGCAACGCATCGGCAAATGCGTGTGTGTAATACCCTTCACCACAGCCAATATCCAGCATCGCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTGGGCGACAATTGCATCACGCAGCGGCTGATAATGTCCGGCATCTAAGAATGCGCGGCGTGCTTGCATCATTTCCGCGCTGTCGCCCGGATCACGAGACCg                                                                                                           <  1:105983/250‑1 (MQ=255)
                                                                                                                                             cTTTTATCGCTACCTTCGAAACATCCAG‑ACCAAACGTGGTGATTTCGGGCAACGCATCGGCAAATGCGTGTGTGTAATACCCTTCACCACAGCCAATATCCAGCATCGCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTGGGCGACAATTGCATCACGCAGCGGCTGATAATGTCCGGCATCTAAGAATGCGCGGCGTGCTTGCATCATTTCCGCGCTGTCGCCCGGATCACGAGACCGTTTATGCTg                                                                                                  >  1:89342/1‑250 (MQ=255)
                                                                                                                                                  aTCGCTACCTTCGAAACATCCAG‑ACCAAACGTGGTGATTTCGGGCAACGCATCGGCAAATGCGTGTGTGTAATACCCTTCACCACAGCCAATATCCAGCATCGCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTGGGCGACAATTGCATCACGCAGCGGCTGATAATGTCCGGCATCTAAGAATGCGCGGCGTGCTTGCATCATTTCCGCGCTGTCGCCCGGATCACGAGACCGTTTATGCTGAACggg                                                                                            >  2:11841/1‑251 (MQ=255)
                                                                                                                                                      cTACCTTCGAAACATCCAG‑ACCAAACGTGGTGATTTCGGGCAACGCATCGGCAAATGCGTGTGTGTAATACCCTTCACCACAGCCAATATCCAGCATCGCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTGGGCGACAATTGCATCACGCAGCGGCTGATAATGTCCGGCATCTAAGAATGCGCGGCGTGCTTGCATCATTTCCGCGCTGTCGCCCGGATCACGAGACCGTTTATGCTGAACGGgcagc                                                                                        >  1:474627/1‑251 (MQ=255)
                                                                                                                                                      cTACCTTCGAAACATACAGCACCAA‑CGTAGTGATTTCGGGCAACGCATCGGCAAATGCGTGTGTGTAATACCCTTCACCACAGCCAATATCCAGCATCGCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTGGGCGACAATTGCATCACGCAGCGGCTGATAATGTCCGGCATCTAAGAATGCGCGGCGTGCTTGCATCATTTACGCGCTGTCGCCCGGATCACGAGACCGTTTCTGCTGAACGGgcagc                                                                                        >  2:357967/1‑251 (MQ=255)
                                                                                                                                                               aaaCATCCAG‑ACCAAACGTGGTGATTTCGGGCAACGCATCGGCAAATGCGTGTGTGTAATACCCTTCACCACAGCCAATATCCAGCATCGCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTGGGCGACAATTGCATCACGCAGCGGCTGATAATGTCCGGCATCTAAGAATGCGCGGCGTGCTTGCATCATTTCCGCGCTGTCGCCCGGATCACGAGACCGTTTATGCTGAACGGGCAGCAGATTGAc                                                                                >  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)
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                
GGGCTTCACTACTCGTGCTAATTCTTCTGCTTTACACGGCGCGTAAATACGTATTATGGCGTCCATACTGGTATCGGAAAACGGCAAACGGTGGCTGGAAG‑CGACACAAAAAGTGACCTGCGGATAGCGTTTCGCCGCCGCTTTTATCGCTACCTTCGAAACATCCAG‑ACCAAACGTGGTGATTTCGGGCAACGCATCGGCAAATGCGTGTGTGTAATACCCTTCACCACAGCCAATATCCAGCACCGCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTGGGCGACAATTGCATCACGCAGCGGCTGATAATGTCCGGCATCTAAGAATGCGCGGCGTGCTTGCATCATTTCCGCGCTGTCGCCCGGATCACGAGACCGTTTATGCTGAACGGGCAGCAGATTGACATACCCTTCTTTCGCCATATCAAACTGATGTCGCTGGGGACAGATATAGCTGTTTTTTTCACGCGAAAGAGGCTGATG  >  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: ‑