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,7960CT96.1% 141.4 / NA 53V89M (GTG→ATG) rlmA23S rRNA m(1)G745 methyltransferase
Reads supporting (aligned to +/- strand):  ref base C (1/1);  major base T (18/31);  minor base G (1/1);  total (20/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.61e-01

GGGCTTCACTACTCGTGCTAATTCTTCTGCTTTACACGGCGCGTAAATACGTATTATGGCGTCCATACTGGTATCGGAAAACGGCAAACGGTGGCTGGAAGCGACACAAAAAGTGACCTGCGGATAGCGTTTCGCCGCCGCTTTTATCGCTACCTTCGAAACATCCAGACCAAACGTGGTGATTTCGGGCAACGCATCGGCAAATGCGTGTGTGTAATACCCTTCACCACAGCCAATATCCAGCACCGCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTGGGCGACAATTGCATCACGCAGCGGCTGATAATGTCCGGCATCTAAGAATGCGCGGCGTGCTTGCATCATTTCCGCGCTGTCGCCCGGATCACGAGACCGTTTATGCTGAACGGGCAGCAGATTGACATACCCTTCTTTCGCCATATCAAACTGATGTCGCTGGGGACAGATATAGCTGTTTTTTTCACGCGAAAGAGGCTGATG  >  NC_000913/1906551‑1907036
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                
gggCTTCACTACTCGTGCTAATTCTTCTGCTTTACACGGCGCGTAAATACGTATTATGGCGTCCATACTGGTATCGGAAAACGGCAAACGGTGGCGGGAAGCGACACAAAAAGTGACCTGCGGATAGCGTTTCGCCGCCGCTTTTATCGCTACCTTCGAAACATCCAGACCAAACGTGGTGATTTCGGGCAACGCATCGGCAAATGCGTGTGTGTAATACCCTTCACCACAGCCAATATCCAGCAtcgccg                                                                                                                                                                                                                                             <  1:553924/251‑1 (MQ=255)
      cACTACACGTGCTAATTCTTCTGCTTTACACGGCGCGTAAATACTTATTATGGCGTCCATCCTGGTATCGGAAAACGGCAAACGGTGGCTGGAAGCGACACAAAAAGGGACCGGCGGATAGCGTTTCGCCGCCGCTTTTATCGCTACCTTCGAAACATCCAGACCAAACGTGGTGATTTCGGGCAACGCATCGGCAAATGCGTGTGTGTAATACCCTTCACCACAGCCAATATCCAGCATCGCCGTGGcc                                                                                                                                                                                                                                        <  2:213122/250‑1 (MQ=255)
       cctacgCGTGCTATTTCTGCTTCTTGACACGGCGCGTAAATTCGTATGATGGCGTCCATACTGGGATCGGAAAACGGCAAACGGTGGCTGGAAGCGACACAAAAAGGGACCTGCGGATAGCTTTTCGCCGCCGCTTTTATCGCTACCTTCGAAACATCCAGACCAAACGTGGTGATTTCGGGCAACGCATCGGCAAATGCGTGTGTGTAATACCCTTCCCCACAGCCAATATCCAGCAGCGCCGTGGCCtt                                                                                                                                                                                                                                      <  2:534154/250‑1 (MQ=255)
       actactCGTGCTAATTCTTCTGCTTTACACGGCGCGTAAATACGTATTATGGCGTCCATACTGGTATCGGAAAACGGCAAACGGTGGCTGGAAGCGACACAAAAAGTGACCTGCGGATAGCGTTTCGCCGCCGCTTTTATCGCTACCTTCGAAACATCCAGACCAAACGTGGTGATTTCGGGCAACGCATCGGCAAATGCGTGTGTGTAATACCCTTCACCACAGCCAATATCCAGCATCGCCGTGGCCtt                                                                                                                                                                                                                                      >  1:345307/1‑251 (MQ=255)
          actCGTGCTAATTCTTCTGCTTTACACGGCTCGTAAATACGTATTATGGCGTCCATACTGGTATCGGAAAACGGCAAACGGTGGCTGGAAGCGACACAAAAAGTGACCTGCGGATAGCGTTTCGCCGCCGCTTTTATCGCTACCTTCGAAACATCCAGACCAAACGTGGTGATTTCGGGCAACGCATCGGCAAATGCGTGTGTGTAATACCCTTCACCACAGCCAATATCCAGCATCGCCGTGGCCTTat                                                                                                                                                                                                                                    <  1:2688/250‑1 (MQ=255)
             cGTGCTAATTCTTCTGCTTTACACGGCGCGTAAATACGTATTATGGCGTCCATACTGGTATCGGAAAACGGCAAACGGTGGCTGGAAGCGACACAAAAAGTGACCTGCGGATAGCGTTTCGCCGCCGCTTTTATCGCTACCTTCGAAACATCCAGACCAAACGTGGTGATTTCGGGCAACGCATCGGCAAATGCGTGTGTGTAATACCCTTCACCACAGCCAATATCCAGCATCGCCGTGGCCTTatcatc                                                                                                                                                                                                                                >  1:103279/1‑251 (MQ=255)
             cGTGCTAATTCTTCTGCTTTACACGGCGCGTAAATACGTATTATGGCGTCCATACTGGTATCGGAAAACGGCAAACGGTGGCTGGAAGCGACACAAAAAGTGACCTGCGGATAGCGTTTCGCCGCCGCTTTTATCGCTACCTTCGAAACATCCAGACCAAACGTGGTGATTTCGGGCAACGCATCGGCAAATGCGTGTGTGTAATACCCTTCACCACAGCCAATATCCAGCATCGCCGTGGCCTTatcatc                                                                                                                                                                                                                                >  1:303776/1‑251 (MQ=255)
             cGTGCTAATTCTTCTGCTTTACACGGCGCGTAAATACGTATTATGGCGTCCATACTGGTATCGGAAAACGGCAAAAGGAGACTGGAAGCGACACAAAAAGTGACCTTCGGAAAGCGTTTCGCCGCCGCTTTTATCGCTACCATCGAAAAATCAAGCCCAAAAGTGGTGATTTCGGGCTACACATCGCAAAATGCGTGTGTGTATAACCCTTCCCAACAGCCAATATACTGCAGCGCAGTAGCCTTatacac                                                                                                                                                                                                                                >  2:537137/1‑247 (MQ=255)
               tGCTAATTCTTCTGCTTTACACGGCGCGTAAATACGTATTATGGCGTCCATACTGGTATCGGAAAACGGCAAACGGTGGCTGGAAGCGACACAAAAAGTGACCTGCGGATAGCGTTTCGCCGCCGCTTTTATCGCTACCTTCGAAACATCCAGACCAAACGTGGTGATTTCGGGCAACGCATCGGCAAATGCGTGTGTGTAATACCCTTCACCACAGCCAATATCCAGCATCGCCGTGGCCTTAtcatca                                                                                                                                                                                                                               >  1:101719/1‑250 (MQ=255)
                            gCTTTACACGGCGCGTAAATACGTATTATGGCGTCCATACTGGTATCGGAAAACGGCAAACGGTGGCTGGAAGCGACACAAAAAGTGACCTGCGGATAGCGTTTCGCCGCCGCTTTTATCGCTACCTTCGAAACATCCAGACCAAACGTGGTGATTTCGGGCAACGCATCGGCAAATGCGTGTGTGTAATACCCTTCACCACAGCCAATATCCAGCACCGCCGTGGCCTTATCATCAAGCCGTTCCCTc                                                                                                                                                                                                                   >  1:360033/1‑249 (MQ=255)
                            gCTTTACACGGCGCGTAAATACGTATTATGGCGCCCATACTGGTATCGGAAAACGGCAAACGGTGGCTGGAAGCGACACAAAAAGTGACCTGCGGATAGCGTTTCGCCGCCGCTTTTATCTCTACCTTCGAAACATCCAGACCAAACGTGGTGATTTCGGGCAACGCATCGGCAAATGCGTGTGTGTAATACCCTTCACCACAGCCAATATCCAGCATCGCCGTGGCCTTATCATCAAGCCGTTCCCTCAg                                                                                                                                                                                                                 <  1:74143/251‑1 (MQ=255)
                                            aaaTACGTATTATGGCGTCCATACTGGTATCGGAAAACGGCAAACGGTGGCTGGAAGCGACACAAAAAGTGACCTGCGGATAGCGTTTCGCCGCCGCTTTTATCGCTACCTTCGAAACATCCAGACCAAACGTGGTGATTTCGGGCAACGCATCGGCAAATGCGTGTGTGTAATACCCTTCACCACAGCCAATATCCAGCATCGCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTGGGCGACAATTGc                                                                                                                                                                                                  >  1:461698/1‑250 (MQ=255)
                                                    attatGGCGTCCATACTGGTATCGGAAAACGGCAAACGGTGGCTGGAAGCGACACAAAAAGTGACCTGCGGATAGCGTTTCGCCGCCGCTTTTATCGCTACCTTCGAAACATCCAGACCAAACGTGGTGATTTCGGGCAACGCATCGGCAAATGCGTGTGTGTAATACCCTTCACCACAGCCAATATCCAGCATCGCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTGGGCGACAATTGCATCACGc                                                                                                                                                                                           >  1:539824/1‑249 (MQ=255)
                                                                          cGGAAAACGGCAAACGGTGGCTGGAAGCGACACAAAAAGTGACCTGCGGATAGCGTTTCGCCGCCGCTTTTATCGCTACCTTCGAAACATCCAGACCAAACGTGGTGATTTCGGGCAACGCATCGGCAAATGCGTGTGTGTAATACCCTTCACCACAGCCAATATCCAGCATCGCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTGGGCGACAATTGCATCACGCAGCGGCTGATAATGTCCGGCATc                                                                                                                                                                    <  1:119562/250‑1 (MQ=255)
                                                                          cGGAAAAAGGCAAACGGTGGCTGGAAGCGACACAAAAAGTGACCTGCGGATAGCGTTTCGCCGCCGCTTTTATCGCTACCTTCGAAACATCCAGACCAAACGTGGTGATTTCGGGCAACGCATCGGCAAATGCGTGTGTGTAATACCCTTCACCACAGCCAATATCCAGCATCGCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTGGGCGACAATTGCATCACGCAGCGGCTGATAATGTCCGGCATc                                                                                                                                                                    <  2:103279/250‑1 (MQ=255)
                                                                                   gCAAACGGTGGCTGGAAGCGACACAAAAAGTGACCTGCGGATAGCGTTTCGCCGCCGCTTTTATCGCTACCTTCGAAACATCCAGACCAAACGTGGTGATTTCGGGCAACGCATCGGCAAATGCGTGTGTGTAATACCCTTCACCACAGCCAATATCCAGCATCGCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTGGGCGACAATTGCATCACGCAGCGGCTGATAATGTCCGGCATCTAAGAATgcg                                                                                                                                                          <  1:320645/251‑1 (MQ=255)
                                                                                      aaCGGTGGCTGGAAGCGACACAAAAAGTGACCTGCGGATAGCGTTTCGCCGCCGCTTTTATCGCTACCTTCGAAACATCCAGACCAAACGTGGTGATTTCGGGCAACGCATCGGCAAATGCGTGTGTGTAATACCCTTCACCACAGCCAATATCCAGCATCGCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTGGGCGACAATTGCATCACGCAGCGGCTGATAATGTCCGGCATCTAAGAATGCgcgg                                                                                                                                                       >  2:97996/1‑251 (MQ=255)
                                                                                         ggTGGCTGGAAGCGACACAAAAAGTGACCTGCGGATAGCGTTTCGCCGCCGCTTTTATCGCTACCTTCGAAACATCCAGACCAAACGTGGTGATTTCGGGCAACGCATCGGCAAATGCGTGTGTGTAATACCCTTCACCACAGCCAATATCCAGCATCGCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTGGGCGACAATTGCATCACGCAGCGGCTGATAATGTCCGGCATCTAAGAATGCgcggcgg                                                                                                                                                    >  2:231291/1‑250 (MQ=255)
                                                                                              cTGGAAGCGACACAAAAAGTGACCTGCGGATAGCGTTTCGCCGCCGCTTTTATCGCTACCTTCGAAACATCCAGACCAAACGTGGTGATTTCGGGCAACGCATCGGCAACTGCGTGTGTGTAATACCCTTCACCACAGCCAATATCCAGCATCGCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTGGGCGACAATTGCATCACGCAGCGGCTGATAATGTCCGGCATCTAAGAATGCGCGGCGTGCTTg                                                                                                                                               >  1:67243/1‑251 (MQ=255)
                                                                                                     cGACACAAAAAGTGACCTGCGGATAGCGTTTCGCCGCCGCTTTTATCGCTACCTTCGAAACATCCAGACCAAACGGGGTGATTTCGGGCAACGCATCGTCAAATGCGTGTGTGTAATACCCTTCACCACAGCCAATATCCAGCATCGCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTGGGCGACAATTGCATCACGCAGCGGCTGATAATGTCCGGCATCTAAGAATGCGCGGCGTGCTTGCATCAtt                                                                                                                                        <  1:245633/251‑1 (MQ=255)
                                                                                                        cacaAAAAGTGACCTGCGGATAGCGTTTCGCCGCCGCTTTTATCGCTACCTTCGAAACATCCAGACCAAACGTGGTGATTTCGGGCAACGCAGCGGCAAATGCGTGTGTGTAATACCCTTCACCACAGCCAATATCCAGCATCGCCGTGGCCGTATCATCAAGCCGTTCCCTCAGTTGGGCGACAATTGCATCACGCAGCGGCTGATAATGTCCGGCATCTAAGAATGCGCGGCGTGCTTGCATCATTTcc                                                                                                                                     <  1:267494/251‑1 (MQ=255)
                                                                                                                    cctgcgcatcgctttgccccgccgcTTTTATCGCTACCTTCGAAACATCCAGACCAAACGTGGTGATTTCGGCCAACGCAGCGGAAAATGCGTGTGTGTAATCCCCTTCACCCCAGCCAATACCCAGCATCGCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTGGGCGACAATTGCATCACGCAGCGGCTGATAATGTCCGGCATCTAAGAATGCGCGGCGTGCTTGCATCATTTCCGCGCTGTCGccc                                                                                                                         <  2:101719/233‑1 (MQ=255)
                                                                                                                         ggatagctgttccgcgccccTTTTAACGCTACCTTCGAACCATCCAGTCCAACCGTGGTGATTTCGGGCAACCCATCGGCAAATGCTGGTTGGTAATACCCTTCCCCACAGCCAATATCCAGCATCGCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTGGGCGACAATTGCATCACGCAGCGGCTGATAATGTCCGGCATCTAAGAATGCGCGGCGTGCTTGCATCATTTCCGCGCTGTCGCCCGGATc                                                                                                                    <  2:345307/237‑1 (MQ=255)
                                                                                                                                  ttccctgccgcTGTTATCGCTACCTTTGAACCCTCCAGACCAAACGTGGTGATTTCGGGGACCGCATCGGCAAATTCGTGTGTGTAATTCCCTTCACCACAGCCAATGCACAGCATCGCCGTGGCCTTATCAACAAGTCGTTCCCTCAGGTGGGCGACAATTGCATCACGCAGCGGCGGATAATGGCCGGCATCTAAGAATGCGCGGCGTGCTTGCATCGTTTCCGCGCGGTCGCCCGGATCACGATAcc                                                                                                            <  2:539179/244‑1 (MQ=255)
                                                                                                                                         ccgtctgtctcgcctcctgcGAAACATCCAGACCAAACGGGGGGATTTCGGGCAACGTAGCGGCAAATGCGGTTGTGTAATACCCTTCACCACAGCCAATATCCAGCATCGCCGGGGCCTTATCTTCAAGCCGTTCCCTCAGTTGGGCGACAATTGCATCACGCAGCGGCTGATAATGTCCGGCATCTAAGAATGCGCGGCGTGCTTTCATCATTTCCGCGCTGTCGTCCGGATCACGAGTCCTTTTATGc                                                                                                    <  2:507313/232‑1 (MQ=255)
                                                                                                                                               ttATCGCTACCTTCGAAACATCCAGACCAAACGTGGTGATTTCGGGCAACGCATCGGCAAATGCGTGTGTGTAATACCCTTCACCACAGCCAATATCCAGCATCGCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTGGGCGACAATTGCATCACGCAGCGGCTGATAATGTCCGGCATCTAAGAATGCGCGGCGTGCTTGCATCATTTCCGCGCTGTCGCCCGGATCACGAGACCGTTTATGCTGAACg                                                                                              >  2:423441/1‑251 (MQ=255)
                                                                                                                                               ttATCGCTACCTTCGAAACAACCAAACCACACGGGGGGATTTCGGGCAACGCATCGGCAAATGCGTGTGTTTAATCCCCTTCACCCCAGCCAATATCCAGCATCGCCGGGGCCTTATCATCAAGCCGTTCCCTCAGTTGGGCGACAATTGCATCACGCAGCGGCTGATAATGTCCGGCATCTAAGAATGCGCGGCGTGCTTGCATCATTTCCGCGCTGTCGCCCGGATCACGAGACCGTTTATGCTGAAc                                                                                               <  2:303776/250‑1 (MQ=255)
                                                                                                                                                tATCGCTACCTTCGAAACATCCAGACCAAACGTGGTGATTTCGGGCAACGCATCGGCAAATGCGTGTGTGTAATACCCTTCACCACAGCCAATATCCAGCATCGCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTGGGCGACAATTGCATCACGCAGCGGCTGATAATGTCCGGCATCTAAGAATGCGCGGCGTGCTTGCATCATTTCCGCGCTGTCGCCCGGATCACGAGACCGTTTATGCTGAACgg                                                                                             >  2:476307/1‑251 (MQ=255)
                                                                                                                                                 aTCGCTACCTTCGAAACATCCAGACCAAACGTGGTGATTTCGGGCAACGCATCGGCAAATGCGTGTGTGTAATACCCTTCACCACAGCCAATATCCAGCATCGCCGTGGCCTTATCATCAAGCCGTTCCCTCCGTTGGGCGACAATTGCATCACGCAGCGGCTGCTAATGCCCGGCATCTAAGAATGCGCGGCGTGCTTGCATCATTTCCGCGCTGTCGCCCGGATCCCGAGACCGTTTATGCTGAACggg                                                                                            >  1:36512/1‑251 (MQ=255)
                                                                                                                                                   cgctaccagcgaaaatcccAGACCAACCGTGGTGATTCCGGCCAACGCATCGGAAAATGGGTGGGTGTAATACCCTTCCCCACAGCCAATATCCAGCATCGCCGTGGCCTTATCATCAAGCCGTTCCCCCAGTTGGGCGACAATTGCATCACGCAGCGGCTGATAATGTCCGGCATCTAAGAATGCGCGGCGTGCTTGCATCATTTCCGCGCTGTCGCCCGGATCACGAGACCGTTTATGCTGAACGGgca                                                                                          <  2:246736/234‑1 (MQ=255)
                                                                                                                                                   cGCTACCTTCGAAACATCCAGACCAAACGTGGTGATTTCGGGCAACGCATCGGCAAATGCGTGTGTGTAATACCCTTCACCACAGCCAATATCCAGCATCGCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTGGGCGACAATTGCATCACGCAGCGGCTGATAATGTCCGGCATCTAAGAATGCGCGGCGTGCTTGCATCATTTCCGCGCTGTCGCCCGGATCACGAGACCGTTTATGCTGAACGGgc                                                                                           <  1:9426/250‑1 (MQ=255)
                                                                                                                                                       aCCTTCGAAACATCCAGACCAAACGTGGTGATTTCGGGCAACGCATCGGCAAATGCGTGTGTGTAATACCCTTCACCACAGCCAATATCCAGCATCGCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTGGGCGACAATTGCATCACGCAGCGGCTGATAATGTCCGGCATCTAAGAATGCGCGGCGTGCTTGCATCATTTCCGCGCTGTCGCCCGGATCACGAGACCGTTTATGCTGAACGGgcagca                                                                                       <  1:423441/250‑1 (MQ=255)
                                                                                                                                                       aCCTTCGAAACATCCAGACCAAACGTGGTGATGTCGGGCAACGCATCGTCAAAGGCGTGTGTGTAATTCCCTTCACCACAGCCAATATCCAGCATCGCCGTGGCCTTATCATCAAGCCGTGCCCTCAGTTGGGCGACAATTGCATCACGCAGCTGCTGATAATGTCCGGCATCTAAGAATGCGCGGCGTGCTTGCACCATTTCCGCGCTGTCGCCCGGATCACGAGACCGTTTATGCTGAACGGGcagcag                                                                                      <  1:559992/251‑1 (MQ=255)
                                                                                                                                                             gAAACATCCAGACCAAACGTGGTGATTTCGGGCAACGCATCGGCAAATGCGTGTGTGTAATACCCTTCACCACAGCCAATATCCAGCATCGCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTGGGCGACAATTGCATCACGCAGCGGCTGATAATGTCCGGCATCTAAGAATGCGCGGCGTGCTTGCATCATTTCCGCGCTGTCGCCCGGATCACGAGACCGTTTATGCTGAACGGGCAGCAGCTTgc                                                                                 >  2:252424/1‑249 (MQ=255)
                                                                                                                                                                    ccAGACCAAACGTGGTGATTTCGGGCAACGCATCGGCAAATGCGTGTGTGTAATACCCTTCACCACAGCCAATATCCAGCATCGCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTGGGCGACAATTGCATCACGCAGCGGCTGATAATGTCCGGCATCTAAGAATGCGCGGCGTGCTTGCATCATTTCCGCGCTGTCGCCCGGATCACGAGACCGTTTATGCTGAACGGGCAGCAGATTGACATACCct                                                                         <  1:395032/251‑1 (MQ=255)
                                                                                                                                                                     cataCCAAACGTGGATATTCGGGGCAACGCATCGGCAAATGCGTGTGTGTAATACCCTTCACCACAGCCAATATCCAGCATCGCCGTGGCCATATCAACAAGCCGTTCCCTCAGTTGGGCGACAATTGCATCACGCAGCGGCTGATAATGTCCGGCAGCTAAGAATGCGCGGCGTGCTTGCATCATTTCCGCGCTGTCGCCCGGATCACGAGACCGTTTATGCGGAACGGGCAGCAGATTGACATACCctt                                                                        <  2:67243/248‑1 (MQ=255)
                                                                                                                                                                       gACCAAACGTGGTGATTTCGGGCAACGCATCGGCAAATGCGTGTGTGTAATACCCTTCACCACAGCCAATATCCAGCATCGCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTGGGCGACAATTGCATCACGCAGCGGCTGATAATGTCCGGCATCTAAGAATGCGCGGCGTGCTTGCATCATTTCCGCGCTGTCGCCCGGATCACGAGACCGTTTATGCTGAACGGGCAGCAGATTGACATACCcttct                                                                      <  1:476307/251‑1 (MQ=255)
                                                                                                                                                                         ccaaactgggtgATTTCGGCCAACTCATCGGCAATTGCGTGTTTGTAATCCCCTTAACCACAGCCATTATCCAGCACCGCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTGGGCGACAATTGCATCACGCAGCGGCTGATAATGTCCGGCATCTAAGAATGCGCGGCGTGCTTGCATCATTTCCGCGCTGTCGCCTGGATCACGAGACCGTTTATGCTGAACGGGCAGCAGATTGACATACCCTTCttt                                                                    <  1:507544/243‑1 (MQ=255)
                                                                                                                                                                                     aTTTCGGGCAACGCATCGGCAAATGCGTGTGTGTAATACCCTTCACCACAGCCAATATCCAGCATCGCCGTGGCGTTATCATCAAGCCGTTCCCTCAGTTGGGCGACAATTGCATCACGCAGCGGCTGATAATGTCCGGCATCTAAGAATGCGCGGCGTGCTTGCATCATTTCCGCGCTGTCGCCCGGATCACGAGACCGTTTATGCTGAACGGGCAGCAGATTGACATACCCTTCTTTCGCCATATCaaa                                                        <  1:249167/251‑1 (MQ=255)
                                                                                                                                                                                       ttCGGGCAACGCATCGGCAAATGCGTGTGTGTAATACCCTTCACCACAGCCAATATCCAGCATCGCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTGGGCGACAATTGCATCACGCAGCGGCTGATAATGTCCGGCATCTAAGAATGCGCGGCGTGCTTGCATCATTTCCGCGCTGTCGCCCGGATCACGAGACCGTTTATGCTGAACGGGCAGCAGATTGACATACCCTTCTTTCGCCATATCAAAc                                                       <  1:53603/250‑1 (MQ=255)
                                                                                                                                                                                        ccGGGCAACGCACCGGAAAAGGCGTTGGTGTAATACCCTTCCCCCCAGCCAATATCCAGCATCGCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTGGGCGACAATTGCATCACGCAGCGGCTGATAATGTCCGGCATCTAAGAATGCGCGGCGTGCTTGCATCATTTCCGCGCTGTCGCCCGGATCACGAGACCGTTTATGCTGAACGGGCAGCAGATTGACATACCCTTCTTTCGCCATATCAAACTg                                                     <  2:539824/250‑1 (MQ=255)
                                                                                                                                                                                             caccGCATCGGCAAATGCGTGTGTGTAATACCCTTCACCACAGCCAATATCCAGCATCGCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTGGGCGACAATTGCATCACGCAGCGGCTGATAATGTCCGGCATCTAAGAATGCGCGGCGTGCTTGCATCATTTCCGCGCTGTCGCCCGGATCACGAGACCGTTTATGCTGAACGGGCAGCAGATTGACATACCCTTCTTTCGCCATATCAAACTGATGTc                                                <  1:354176/248‑1 (MQ=255)
                                                                                                                                                                                                 gCATCGGCAAATGCGTGTGTGTAATACCCTTCACCACAGCCAATATCCAGCATCGCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTGGGCGACAATTGCATCACGCAGCGGCTGATAATGTCCGGCATCTAAGAATGCGCGGCGTGCTTGCATCATTTCCGCGCTGTCGCCCGGATCACGAGACCGTTTATGCTGAACGGGCAGCAGATTGACATACCCTTCTTTCGCCATATCAAACTGATGTCGCTg                                            <  1:231291/251‑1 (MQ=255)
                                                                                                                                                                                                   aTCGGCAAATGCGTGTGTGTAATACCCTTCACCACAGCCAATATCCAGCATCGCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTGGGCGACAATTGCATCACGCAGCGGCTGATAATGTCCGGCATCTAAGAATGCGCGGCGTGCTTGCATCATTTCCGCGCTGTCGCCCGGATCACGAGACCGTTTATGCTGAACGGGCAGCAGATTGACATACACTTCTTTCGCCATATCAAACTGATGTCGCTggg                                          >  2:373436/1‑251 (MQ=255)
                                                                                                                                                                                                     cGGCAAATGCGTGTGTGTAATACCCTTCACCACAGCCAATATCCAGCATCGCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTGGGCGACAATTGCATCACGCAGCGGCTGATAATGTCCGGCATCTAAGAATGCGCGGCGTGCTTGCATCATTTCCGCGCTGTCGCCCGGATCACGAGACCGTTTATGCTGAACGGGCAGCAGATTGACATACCCTTCTTTCGCCATATCAAACTGATGTCGCTgggg                                         <  1:325827/250‑1 (MQ=255)
                                                                                                                                                                                                        cAAATGCGTGTGTGTAATACCCTTCACCACAGCCAATATCCAGCATCGCCGTGGCCTTATCATCAAGCCGTGCCCTCAGTTGGGCGACAATTGCATCACGCAGCGGCTGATAATGTCCGGCATCTAAGAATGCGCGGCGTGCTTGCATCATTTCCGCGCTGTCGCCCGGATCACGAGACCGTTTATGCTGAACGGGCAGCAGATTGACATCCCCTTCTTTCGCCATATCAAACTGATGTCGCTGGGGACAg                                     <  1:205730/251‑1 (MQ=255)
                                                                                                                                                                                                            tGCGTGAGTGTAATACCCTTCACCACAGCCAATATCCAGCATCGCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTGGGCGACAATTGCATCACGCAGCGGCTGATAATGTCCGACATCTAAGAATGCGCGGCGTGCTTGCATCATTTCCGCGCTGTCGCCCGGATCACGCGACCGTTTATGCTGAACGGGCAGCAGATTACCATACCCTTCTTTCGCCATATCAAACTGATGTCGCTGGGGACAGatat                                 >  2:352292/1‑251 (MQ=255)
                                                                                                                                                                                                              cGTGTGTGGAATACCCTTCACCACAGCCAATATCCAGCATCGCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTGGCCGACAATTGCATCACGCAGCGGCTGATAATGTCCGGCCTCTAAGCATGCGCGGCGTGCTTGCATCATTTCCGCGCTGACGCCCGGATCACGAGACCGTTTATGATGAACGGGCAGCAGCTTGACATACCATTCTTTCGCCATATCAAACTGATGTCGCTGGGGACCGatctct                               >  2:317850/1‑247 (MQ=255)
                                                                                                                                                                                                                   gtgtAATACCCTTCACCACAGCCAATATCCAGCATCGCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTGGGCGACAATTGCATCACGCAGCGGCTGATAATGTCCGGCATCTACGAATGCGCGGCGTGCTTGCATCATTTCCGCGCTGTCGCCCGGATCACGAGACCGTTTATGCTGAACGGGCAGCAGATTGACATACCCTTCTTTCGCCATATCAAACTGATGCCGCTGGGGACAGATATAGCTGtt                          >  2:219444/1‑251 (MQ=255)
                                                                                                                                                                                                                        atccccTTCACCACAGCCAATATCCAGCATCGCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTGGGCGACAATTGCATCACGCAGCGGCTGATAATGTCCGGCATCTAAGAATGCGCGGCGTGCTTGCATCATTTCCGCGCTGTCGCCCGGATCACGAGACCGTTTATGCTGAACGGGCAGCAGATTGACATACCCTTCTTTCGCCATATCAAACTGATGTCGCTGGGGACAGATATAGCTGttttttt                     <  1:306909/248‑1 (MQ=255)
                                                                                                                                                                                                                                caccacAGCCAATATCCAGCATCGCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTGGGCGACAATTGCATCACGCAGCGGCTGATAATGTCCGGCATCTAAGAATGCGCGGCGTGCTTGCATCATTTCCGCGCTGTCGCCCGGATCACGAGACCGTTTATGCTGAACGGGCAGCAGATTGACATACCCTTCTTTCGCCATATCAAACTGATGTCGCTGGGGACAGATATAGCTGTTTTTTTCACGCGa              <  1:252424/250‑1 (MQ=255)
                                                                                                                                                                                                                                           atatCCAGCATCGCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTGGGCGACAATTGCATCACGCAGCGGCTGATAATGTCCGGCATCTAAGAATGCGCGGCGTGCTTGCATCATTTCCGCGCTGTCGCCCGGATCACGAGACCGTTTATGCTGAACGGGCAGCAGATTGACATACCCTTCTTTCGCCATATCAAACTGATGTCGCTGGGGACAGATATAGCTGTTTTTTTCACGCGAAAGAGGCTGATg  <  1:13643/251‑1 (MQ=255)
                                                                                                                                                                                                                                                 aGCATCGCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTGGGCGACAATTGCATCACGCAGCGGCTGATAATGTCCGGCATCTAAGAATGCGCGGCGTGCTTGCATCAtt                                                                                                                                        >  1:529119/1‑111 (MQ=255)
                                                                                                                                                                                                                                                 aGCATCGCCGTGGCCGTATCATCAAGCCGTTCCCTCAGTTGGGCGACAATTGCATCACGCAGCGGCTGATAATGTCCGGCATCTAAGAATGCGCGGCGTGCTTGCATCAtt                                                                                                                                        <  2:529119/111‑1 (MQ=255)
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                
GGGCTTCACTACTCGTGCTAATTCTTCTGCTTTACACGGCGCGTAAATACGTATTATGGCGTCCATACTGGTATCGGAAAACGGCAAACGGTGGCTGGAAGCGACACAAAAAGTGACCTGCGGATAGCGTTTCGCCGCCGCTTTTATCGCTACCTTCGAAACATCCAGACCAAACGTGGTGATTTCGGGCAACGCATCGGCAAATGCGTGTGTGTAATACCCTTCACCACAGCCAATATCCAGCACCGCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTGGGCGACAATTGCATCACGCAGCGGCTGATAATGTCCGGCATCTAAGAATGCGCGGCGTGCTTGCATCATTTCCGCGCTGTCGCCCGGATCACGAGACCGTTTATGCTGAACGGGCAGCAGATTGACATACCCTTCTTTCGCCATATCAAACTGATGTCGCTGGGGACAGATATAGCTGTTTTTTTCACGCGAAAGAGGCTGATG  >  NC_000913/1906551‑1907036

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