Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,641,428 G→A A163T (GCG→ACG)  yjtD → putative tRNA/rRNA methyltransferase YjtD

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,641,4280GA97.8% 142.1 / NA 45A163T (GCG→ACG) yjtDputative tRNA/rRNA methyltransferase YjtD
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base A (20/24);  total (20/25)
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.99e-01

CGGGCGAAATATCATTACTACGCCACGCCAGTTGAACTGGTGCCGCTGTTAGAGGAAAAATCTTCATGGATGAGCCATGCCGCGCTGGTGTTTGGTCGCGAAGATTCCGGGTTGACTAACGAAGAGTTAGCGTTGGCTGACGTTCTTACTGGTGTGCCGATGGTGGCGGATTATCCTTCGCTCAATCTGGGGCAGGCGGTGATGGTCTATTGCTATCAATTAGCAACATTAATACAACAACCGGCGAAAAGTGATGCAACGGCAGACCAACATCAACTGCAAGCTTTACGCGAACGAGCCATGACATTGCTGACGACTCTGGCAGTGGCAGATGACATAAAACTGGTCGACTGGTTACAACAACGCCTGGGGCTTTTAGAGCAACGAGACACGGCAATGTTGCACCGTTTGCTGCATGATATTGAAAAAAATATCACCAAATAAAAAACGCCTTAGTAAGTATTTTTC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_000913/4641185‑4641652
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                    
cGGGCGAAATATCATTACTACGCTACGCCAGTAGAACTGGTACCGCTGTTAGAGGAAAAATCTTCATGGATGAGCCATGCCGCGCTGGTGTTTGGTCGCGAAGATTCCGGGCTGACTAACGAAGAGTTAGCGTTGGCTGACGTTCTTACTGGTGTGCCGATGGTGGCAGATTATCCTTCGCTCAATCTGGGGCAGGCAGTGATGGTCTATTGCTATCAATTAGCAACATTAATAAAACAACCGACGaaa                                                                                                                                                                                                                                                 >  2:341912/1‑249 (MQ=255)
  ggCGAAATATCATTACTACGCTACGCCAGTAGAACTGGTACCGCTGTTAGAGGAAAAATCTTCATGGATGAGCCATGCCGCGCTGGTGTTTGGTCGCGAAGATTCCGGGCTGACTAACGAAGAGTTAGCGTTGGCTGACGTTCTTACTGGTGTGCCGATGGTGGCAGATTATCCTTCGCTCAATCTGGGTCAGGCAGTGATGGTCTATTGCTATCAATTAGCAACATTAATACAACAACCGACGAAAAGTg                                                                                                                                                                                                                                             >  2:105406/1‑251 (MQ=255)
                          gccaggagAACTGGTACCGCTGTTAGAGGAAAAATCTTCATGGATGAGCCATGCCGCGCTGGTGTTTGGTCGCGAAGATTCCGGGCTGACTAACGAAGAGTTAGCGTTGGCTGACGTTCTTACTGGTGTGCCGATGGTGGCAGATTATCCTTCGCTCAATCTGGGGCAGGCAGTGATGGTCTATTGCTATCAATTAGCAACATTAATACAACAACCGACGAAAAGTGACACAACGGCAGACCAACATCAAc                                                                                                                                                                                                                     <  2:45333/244‑1 (MQ=255)
                            caGTAGAACTGGTACCGCTGTTAGAGGAAAAATCTTCATGGATGAGCCATGCCGCGCTGGTGTTTGGTCGCGAAGATTCCGGGCTGACTAACGAAGAGTTAGCGTTGGCTGACGTTCTTACTGGTGTGCCGATGGTGGCAGATTATCCTTCGCTCAATCTGGGGCAGGCAGTGATGGTCTATTGCTATCAATTAGCAACATTAATACAACAACCGACGAAAAGTGACACAACGGCAGACCAACATCAACTg                                                                                                                                                                                                                   <  1:303379/251‑1 (MQ=255)
                                agAACTGGTACCGCTGTTAGAGGAAAAAGCTTCATGGATGAGCCATGCCGCGCTGGTGTTTGGTCGCGAAGATTCCGGGCTGACTAACGAAGAGTTAGCGTTGGCTGACGTTCTTACTGGTGTGCCGATGGTGGCAGATTATCCTTCGCTCAATCTGGGGCAGGTAGTGATGGTCTATTGCTATCAATTAGCAACATGAATACAACAACCGACGAAAAGTGACACAACGGCAGACCAACATCAACTGCAAg                                                                                                                                                                                                               <  1:267184/250‑1 (MQ=255)
                                   acaggtaccGCTGTTAGAGGAAAAATCTTCATGGATGAGCCATGCCGCGCTGGTGTTTGGTCGCGAAGATTCCGGGCTGACTAACGAAGAGTTAGCGTTGGCTGACGTTCTTACTGGTGTGCCGAGGGTGGCAGATTATCCTTCGCTCAATCTGGGGCAGGCAGTGATGGTCTATTGCTATCAATTAGCAACATTAATACAACAACCGACGAAAAGTGACACAACGGCAGACCAACAACAACTGCAAGCtt                                                                                                                                                                                                            <  2:19796/244‑1 (MQ=255)
                                    cTGGTACCGCTGTTAGAGGAAAAATCTTCATGGATGAGCCATGCCGCGCTGGTGTTTGGTCGCGAAGATTCCGGGCTGACTAACGAAGAGTTAGCGTTGGCTGACGTTCTTACTGGTGTGCCGATGGTGGCAGATTATCCTTCGCTCAATCTGGGGCAGGCAGTGATGGTCTATTGCTATCAATTAGCAACATTAATACAACAACCGACGAAAAGTGACACAACGGCAGACCAACATCAACTGCAAGCttt                                                                                                                                                                                                           >  1:452131/1‑251 (MQ=255)
                                    cTGGTACCGCTGTTAGAGGAAAAATCTTCATGGATGAGCCATGCCGCGATGGTGTTTGGTCGCGAAGATACCGGGCTGACTAACGAAGAGTTAGCGTTGGCTGACGTTCTTACTGGTGTGACGATGGTGGCAGATTATCCTTCGCTCAATCTGGGGCAGGCAGTGATGGTCTATTGCTATCAATTAGCAACATTAATACAACAACAGAAGAAAAGTGACACAAAGGCAGGCCAACATCAACAGAAAGCttt                                                                                                                                                                                                           >  1:228804/1‑251 (MQ=255)
                                           cGCTGTTAGAGGAAAAATCTTCATGGATGAGCCATGCCGCGCTGGTGTTTGGTCGCGAAGATTCCGGGCTGACTAACGAAGAGTTAGCGTTGGCTGACGTTCTTACTGGTGTGCCGATGGTGGCAGATTATCCTTCGCTCAATCTGGGGCAGGCAGTGATGGTCTATTGCTATCAATTAGCAACATTAATACAACAACCGACGAAAAGTGACACAACGGCAGACCAACATCAACTGCAAGCTTTACGCGa                                                                                                                                                                                                     >  1:321078/1‑250 (MQ=255)
                                           cGCTGTTAGAGGAAAAATCTTCATGGATGAGCCATGCCGCGCTGGTGTTTGGTCGCGAAGATTCCGGGCTGACTAACGAAGAGTTAGCGTTGGCTGACGTTCTTACTGGTGTGCCGATGGTGGCAGATTATCCTTCGCTCAATCTGGGGCAGGCAGTGATGGTCTATTGCTATCAATTAGCAACATTAATACAACAACCGACGAAAAGTGACACAACGGCAGACCAACATCAACTGCAAGCTTTACGCGa                                                                                                                                                                                                     >  2:184913/1‑250 (MQ=255)
                                                 tAGAGGAAAAATCTTCATGGATGAGCCATGCCGCGCTGGTGTTTGGTCGCGAAGATTCCGGGCTGACTAACGAAGAGTTAGCGTTGGCTGCCGTTCTTACTGGTGTGCCGATGGTGGCAGATTATCCTTCGCTCAATCTGGGGCAGGCAGTGATGGTCTATTGCTATCAATTAGCAACATTAATACAACAACCGACGAAAAGTGACACAACGGCAGACCAACATCAACTGCAAGCTTTACGCGAACGAgtc                                                                                                                                                                                              >  1:298562/1‑249 (MQ=255)
                                                   aagGAATAAGCTTTATGGATGAGCCATGCAGCGATGTTTTGTGGGCGCTTAGATTCCGGGATGACGAACGAAGAGTTAGCGTGGGCTGACGTTCTTACTGGTGTGCCGATGGTGGCAGATTATCCTTCGCTCAATCTGGGGAAGGCATTGATGGTCTATTGCTAGCAATTAGCAACATTAATACAACAACCGGCGAAAAGTGACACAACGGCAGACCAACATCAACTGCAAGATTTGCTCGAACGAGTCAt                                                                                                                                                                                            <  2:301738/250‑1 (MQ=255)
                                                    gggAAAAATCTTCATGGATGAGCCATGCCGCGCTGGTGTTTGGTCGCGAAGATGCCGGGCTGACTAACGAAGAGTTAGCGTTGGCTGACGTTCTTACTGGTGTGCCGATGGTGGCAGATTATCCTTCGCTCAATCTGGGGCAGGCAGTGATGGTCTATTGCTATCAATTAGCAACATTAATACAACAACCGACGAAAAGTGACACAACGGCAGACCAACATCAACTGCAAGCTTTACGCGAACGAGTCATg                                                                                                                                                                                           <  1:438371/250‑1 (MQ=255)
                                                               tCATGGATGAGCCATGCCGCGCTGGTGTTTGGTCGCGAAGATTCCGGGCTGACTAACGAAGAGTTAGCGTTGGCTGACGTTCTTACTGGTGTGCCGATGGTGGCAGATTATCCTTCGCTCAATCTGGGGCAGGCAGTGATGGTCTATTGCTATCAATTAGCAACATTAATACAACAACCGACGAAAAGTGACACAACGGCAGACCAACGTCAACTGCAAGCTTTACGCGAACGAGTCATGGCATTGc                                                                                                                                                                                    <  2:321078/247‑1 (MQ=255)
                                                               tCATGGATGAGCCATGCCGCGCTGGTGTTTGGTCGCGAAGATTCCGGGCTGACTAACGAAGAGTTAGCGTTGGCTGACGTTCTTACTGGTGTGCCGATGGTGGCAGATTATCCTTCGCTCAATCTGGGGCAGGCAGTGATGGTCTATTGCTATCAATTAGCAACATTAATACAACAACCGACGAAAAGTGACACAACGGCAGACCAACATCAACTGCAAGCTTTACGCGAACGAGTCATGGCATTGCt                                                                                                                                                                                   <  1:403231/248‑1 (MQ=255)
                                                                cATGGATGAGCCATGCCGCGCTGGTGTTTGGTCGCGAAGATTCCGGGCTGACTAACGAAGAGTTAGCGTTGGCTGACGTTCTTACTGGTGTGCCGATGGTGGCAGATTATCCTTCGCTCAATCTGGGGCAGGCAGTGATGGTCTATTGCTATCAATTAGCAACATTAATACAACAACCGACGAAAAGTGACACAACGGCAGACCAACATCAACTGCAAGCTTTACGCGAACGAGTCATGGCATTGCTgacg                                                                                                                                                                               <  1:418016/251‑1 (MQ=255)
                                                                                    gtGGTGTTTGGTCGCGAAGATTCCGGGCTGACTAACGAAGAGTTAGCGTTGGCTGACGTTCTTACTGGTGTGCCGATGGTGGCAGATTATCCTTCGCTCAATCTGGGGCAGGCAGTGATGGTCTATTGCTATCAATTAGCAACATTAATACAACAACCGACGAAAAGTGACACAACGGCAGACCAACATCAACTGCAAGCTTTACGCGAACGAGTCATGGCATTGCTGACGACTCTGGCAGTGGCAGATGa                                                                                                                                                           <  2:142156/250‑1 (MQ=255)
                                                                                          tttGGTCGCGAAGATTCCGGGCTGACTAACGAAGAGTTAGCGTTGGCTGACGTTCTTACTGGTGTGCCGATGGTGGCAGATTATCCTTCGCTCAATCTGGGGCAGGCAGTGATGGTCTATTGCTATCAATTAGCAACATTAATACAACAACCGACGAAAAGTGACACAACGGCAGACCAACATCAACTGCAAGCTTTACGCGAACGAGTCATGGCATTGCTGACGACTCTGGCAGTGGCAGATGACATaa                                                                                                                                                      >  1:136533/1‑250 (MQ=255)
                                                                                           ttGGTCGCGAAGATTCCGGGATGACTAACGAAGAGTGAGCGTTGGCTGACGTTCTTACGGGTGGGCCGATGGTGGCAGATTATCCTTCGCTCAATCTGGGGCAGGCAGTGATGGTCTATTGCTATCAATTAGCAACATTAATACAACAACCGACGAAAAGTGACACAACGGCAGACCAACATCAACTGCAAGCTTTACGCGAACGAGTCATGGCATTGCTGACGACTCTGGCAGTGGCAGATGACATaaa                                                                                                                                                     <  2:59834/250‑1 (MQ=255)
                                                                                              gTCGCGAAGATTCCGGGCTGACTAACGAAGAGTTAGCGTTGGCTGACGTTCTTACTGGTGTGACGATGGGGGCAGATTATCCTTCGCTCAATCTGGGGCAGGCAGTGATGGTCTATTGCTATCAATTAGCAACATTAATACAACAACCGACGAAAAGTGACACAACGGCAGACCAACATCAACTGCAAGCTTTACGCGAACGAGTCATGGCATTGCTGACGACTCTGGCAGTGGCAGATGACATAAAACTg                                                                                                                                                 >  1:275515/1‑251 (MQ=255)
                                                                                                  cgAAGATTCCGGGCTGACTAACGAAGAGTTAGCGTTGGCTGACGTTCTTACTGGTGTGCCGATGGTGGCAGATTATCCTTCGCTCAATCTGGGGCAGGCAGTGATGGTCTATTGCTATCAATTAGCAACATTAATACAACAACCGACGAAAAGTGACACAACGGCAGACCAACATCAACTGCAAGCTTTACGCGAACGAGTCATGGCATTGCTGACGACTCTGGCAGTGGCAGATGACATAAAACTGGTCg                                                                                                                                             >  1:461196/1‑251 (MQ=255)
                                                                                                                 gACTAACGAAGAGTTAGCGTTGGCTGACGTTCTTACTGGTGTGCCGATGGTGGCAGATTATCCTTCGCTCAATCTGGGGCAGGCAGTGATGGTCTATTGCTATCAATTAGCAACATTAATACAACAACCGACGAAAAGTGACACAACGGCAGACCAACATCAACTGCAAGCTTTACGCGAACGAGTCATGGCATTGCTGACGACTCTGGCAGTGGCAGATGACATAAAACTGGTCGACTGGTTacagca                                                                                                                                <  1:105406/249‑4 (MQ=255)
                                                                                                                                             gTTCTTACGGGTGTGCCGATGGTGGCAGATTATCCTTCGCTCAATCTGGGGCAGGCAGTGATGGTCTATTGCTATCAATTAGCAACATTAATACAACAACCGACGAAAAGTGACACAACGGCAGACCAACATCAACTGCAAGCTTTACGCGAACGAGTCATGGCATTGCTGACGACTCTGGCAGTGGCAGATGACATAAAACTGGTCGACTGGTTACAGCAACGCCTGGGGCTTTTAGAGCAACGAGacac                                                                                                  <  2:136533/251‑1 (MQ=255)
                                                                                                                                                  tACTGGGGTGCCGATGGTGGCAGATTATCCTTCGCTCAATCTGGGGCAGGCAGTGATGGTCTATTGCTATCAATTAGCAACATTAATACAACAACCGACGAAAAGTGACACAACGGCAGACCAACATCAACTGCAAGCTTTACGCGAACGAGTCATGGCATTGCTGACGACTCTGGCAGTGGCAGATGACATAAAACTGGTCGACTGGTTACAGCAACGCCTGGGGCTTTTAGAGCAACGAGACACGGCa                                                                                              <  1:310796/250‑1 (MQ=255)
                                                                                                                                                      ggTGTGCCGATGGTGGCAGATTATCCTTCGCTCAATCTGGGGCAGGCAGTGATGGTCTATTGCTATCAATTAGCAACATTAATACAACAACCGACGAAAAGTGACACAACGGCAGACCAACATCAACTGCAAGCTTTACGCGAACGAGTCATGGCATTGCTGACGACTCTGGCAGTGGCAGATGACATAAAACTGGTCGACTGGTTACAGCAACGCCTGGGGCTTTTAGAGCAACGAGACACGGCAATg                                                                                           <  1:82576/249‑1 (MQ=255)
                                                                                                                                                         gtgCCGATGGTGGCAGATTATCCTTCGCTCAATCTGGGGCAGGCAGTGATGGTCTATTGCTATCAATTAGCAACATTAATACAACAACCGACGAAAAGTGACACAACGGCAGACCAACATCAACTGCAAGCTTTACGCGAACGAGTCATGGCATTGCTGACGACTCTGGCAGTGGCAGATGACATAAAACTGGTCGACTGGTTACAGCAACGCCTGGGGCTTTTAGAGCAACGAGACACGGCAATGTTGc                                                                                       >  1:477668/1‑250 (MQ=255)
                                                                                                                                                          ggCCGATGGTGGCAGATTATCCTTCGCTCAATCTGGGGCAGGCAGTGATGGTCTATTGCTATCAATTAGCAACATTAATACAACAACCGACGAAAAGTGACACAACGGCAGACCAACATCAACTGCAAGCTTTACGCGAACGAGTCATGGCATTGCTGACGACTCTGGCAGTGGCAGATGACATAAAACTGGTCGACTGGTTACAGCAACGCCTGGGGCTTTTAGAGCAACGAGACACGGCAATGTTGCAc                                                                                     <  1:453603/250‑1 (MQ=255)
                                                                                                                                                            ccGATGGTGGCAGATTATCCTTCGCTCAATCTGGGGCAGGCAGTGATGGTCTATTGCTATCAATTAGCAACATTAATACAACAACCGACGAAAAGTGACACAACGGCAGACCAACATCAACTGCAAGCTTTACGCGAACGAGTCATGGCATTGCTGACGACTCTGGCAGTGGCAGATGACATAAAACTGGTCGACTGGTTACAGCAACGCCTGGGGCTTTTAGAGCAACGAGACACGGCAATGTTGCACCg                                                                                   <  1:27132/251‑1 (MQ=255)
                                                                                                                                                                      aagATTATCCGTCGCTCAATCTGGGGCAGGTAGTGATGGTCTATTGCTAGCAATTAGCAACATTAATACAACAACCGACGAAAAGTGACACAACGGCAGACCAACATCAACTGCAAGCTTTACGCGAACGAGTCATGGCATTGCTGACGACTCTGGCAGTGGCAGATGACATAAAACTGGTCGACTGGTTACAGCAACGCCTGGGGCTTTTAGAGCAACGAGACACGGCAGTGTTGCACCGTTTGCTGCAt                                                                         <  2:298562/249‑1 (MQ=255)
                                                                                                                                                                        gATTATCCTTCGCTCAATCTGGGGCAGGCAGTGATGGTCTATTGCTATCAATTAGCAACATTAATACAACAACCGACGAAAAGTGACACAACGGCAGACCAACATCAACTGCAAGCTTTACGCGAACGAGTCATGGCATTGCTGACGACTCTGGCAGTGGCAGATGACATAAAACTGGTCGACTGGTTACAGCAACGCCTGGGGCTTTTAGAGCAACGAGACACGGCAATGTTGCACCGTTTGCTGCATg                                                                        >  1:56955/1‑250 (MQ=255)
                                                                                                                                                                        gATTATCCTTCGCTCAATCTGGGGCAGGCAGTGATGGTCTATTGCTATCAATTAGCAACATTAATACAACAACCGACGAAAAGTGACACAACGGCAGACCAACATCAACTGCAAGCTTTACGCGAACGAGTCATGGCATTGCTGACGACTCTGGCAGTGGCAGATGACATAAAACTGGTCGACTGGTTACAGAAACGCCTGGGGCTTTTAGAGCAACGGGACACGGCAATGTTGCACCGTTTGCTGCATg                                                                        >  1:326432/1‑250 (MQ=255)
                                                                                                                                                                                          cTGGGGCAGGCAGTGATGGTCTATTGCTCTCAATTAGCAACATTAATACAACAACCGACGAAAAGTGACACAACGGCAGACCAACATCAACTGCAAGCTTTACGCGAACGAGTAATGGCATTGCTGACGACTCTGGCAGTGGCAGATGACATAAAACTGGTCGACTGGTTACAGCAACGCCTGGGGCTTTTAGAGAAACGAGACACGGCAATGTTGCACCGTTTGCTGCATGATATTGAAAAAAATATCAc                                                     >  1:85341/1‑251 (MQ=255)
                                                                                                                                                                                             gggCAGGCAGTGATGGTCTATTGCTATCAATTTGCAACATTAATACAACAACCGACGAAAAGTGACACAACGGCAGACCAACATCAACTGCAAGCTTTACGCGAACGAGTCATGGCATTGCTGACGACTCTGGCAGTGGCAGATGACATAAAACTGGTCGACTGGTTACAGCAACGCCTGGGGCTTTTAGAGCAACGAGACACGGCAATGTTGCACCGTTtgct                                                                             >  2:314221/1‑224 (MQ=255)
                                                                                                                                                                                             gggCAGGCAGTGATGGTCTATTGCTATCAATTTGCAACATTAATACAACAACCGACGAAAAGTGACACAACGGCAGACCAACATCAACTGCAAGCTTTACGCGAACGAGTCATGGCATTGCTGACGACTCTGGCAGTGGCAGATGACATAAAACTGGTCGACTGGTTACAGCAACGCCTGGGGCTTTTAGAGCAACGAGACACGGCAATGTTGCACCGTTtgct                                                                             <  1:314221/224‑1 (MQ=255)
                                                                                                                                                                                                  ggcagTGATGGTCTATTGCTATCAATTAGCAACATTAATACAACAACCGACGAAAAGTGACACAACGGCAGACCAACATCAACTGCAAGCTTTACGCGAACGAGTCATGGCATTGCTGACGACTCTGGCAGTGGCAGATGACATAAAACTGGTCGACTGGTTACAGCAACGCCTGGGGCTTTTAGAGCAACGAGACACGGCAATGTTGCACCGTTTGCTGCATGATATTGAAAAAAATATCACCAAATaa                                              >  1:126039/5‑250 (MQ=255)
                                                                                                                                                                                                             tCTATTGCTATCAATTAGCAACATTAATACAACAACCGACGAAAAGTGACACAACGGCAGACCAACATCAACTGCAAGCTTTACGCGAACGAGTCATGGCATTGCTGACGACTCTGGCAGTGGCAGATGACATAAAACTGGTCGACTGGTTACAGCAACGCCTGGGGCTTTTAGAGCAACGAGACACGGCAATGTTGCACCGTTTGCTGCATGATATTGAAAAAAATATCACCAAATAAAAAACGCCtt                                    >  1:435184/1‑249 (MQ=255)
                                                                                                                                                                                                                aTTGCTATCAATTAGCAACATTAATACAACAACCGACGAAAAGTGACACAACGGCAGACCAACATCAACTGCAAGCTTTACGCGAACGAGTCATGGCATTGCTGACGACTCtgg                                                                                                                                                                        >  1:419353/1‑114 (MQ=255)
                                                                                                                                                                                                                aTTGCTATCAATTAGCAACATTAATACAACAACCGACGAAAAGTGACACAACGGCAGACCAACATCAACTGCAAGCTTTACGCGAACGAGTCATGGCATTGCTGACGACTCtgg                                                                                                                                                                        <  2:419353/114‑1 (MQ=255)
                                                                                                                                                                                                                       acaagtAGCAACATTAATACAACAACCGACGAAAAGTGACACAACGGCAGACCAACATCAACTGCAAGCTTTACGCGAACGAGTCATGGCATTGCTGACGACTCTGGCAGTGGCAGATGACATAAAACTGGTCGACTGGTTACAGCAACGCCTGGGGCTTTTAGAGCAACGAGACACGGCAATGTTGCACCGTTTGCTGCATGATATTGAAAAAAATATCACCAAATAAAAAACGCCTTAGTAAGTAtttt                        <  2:360633/246‑1 (MQ=255)
                                                                                                                                                                                                                              gCAACATTAATACCACAACTGACGAAAAGTGACACAACGGCAGACCAACATTAACTGCAAGCTTTACTCGAACCAGTCATGGCATGTCTTACGACTCTGGCAGTGGCAGATGACATAAAACTGGTCGACTTGTTACAGCAGCGCCTGGTGCTTTTAGAGCAACGAGACACGGCAATGTTGCACCTTTTGCTGTATGATATTTAAAAAAATATCACCAAAAAAAAAACGCCTTAGTAAGTATTTTTcagctt                 <  2:283627/251‑6 (MQ=255)
                                                                                                                                                                                                                              gCAACATTAATACAACAACCGACGAAAAGTGACACAACGGCAGACCAACATCAACTGCAAGCTTTACGCGAACGAGTCATGGCATTGCTGACGACTCTGGCAGTGGCAGATGACATAAAACTGGTCGACTGGTTACAGCAACGCCTGGGGCTTTTAGAGCAACGAGACACGGCAATGTTGCACCGTTTGCTGCATGATATTGAAAAAAATATCACCAAATAAAAAACGCCTTAGTAAGTATTTTTCagctt                 >  1:289724‑M2/1‑246 (MQ=255)
                                                                                                                                                                                                                              gCAACATTAATACAACAACCGACGAAAAGTGACACAACGGCAGACCAACATCAACTGCAAGCTTTACGCGAACGAGTCATGGCATTGCTGACGACTCTGGCAGTGGCAGATGACATAAAACTGGTCGACTGGTTACAGCAACGCCTGGGGCTTTTAGAGCAACGAGACACGGCAATGTTGCACCGTTTGCTGCATGATATTGAAAAAAATATCACCAAATAAAAAACGCCTTAGTAAGTATTTTTCagctt                 >  2:435184‑M2/1‑246 (MQ=255)
                                                                                                                                                                                                                                   aTTAATACAACAACCGACGAAAAGTGACACAACGGCAGACCAACATCAACAGAAAGCTTTACGCGAACGAGTCATGGCATTGCTGACGACTCTGGCAGTGGAAGATGACATAAAACTGGTCGACTGGTTACAGCAACGCCTGGGGCTTTTAGAGCAACGAGACACGGCAATGTTGCACCGTTTGCTGCATTATATTGAAAAAAATATCAACAACTAAAAAAAGCCTTAGCAAGTATTTTTCagcttttgat            <  2:217485‑M2/251‑11 (MQ=255)
                                                                                                                                                                                                                                              aacCGACGAAAAGTGACACAACGGCAGACCAACATCAACTGCAAGCTTTACGCGAACGCGTCATGGCATTGCTGACGACTCTGGCAGTGGCAGATGACATAAAACTGGTCGACTGGTTACAGCAACGCCTGGGGCTTTTAGAGCAACGAGACACGGCAATGTTGCACCGTTTGCTGCATGATATTGAAAAAAATATCACCAAATAAAAAACGCCTTAGTAAGTATTTTTCagcttttcattctgactgca  <  1:64239‑M2/250‑21 (MQ=255)
                                                                                                                                                                                                                                              aacCGACGAAAAGTGACACAACGGCAGACCAACATCAACTGCAAGCTTTACGCGAACGAGTCATGGCATTGCTGACGACTCTGGCAGTGGCAGATGACATAAAACTGGTCGACTGGTTACAGCAACGCCTGGGGCTTTTAGAGCAACGAGACACGGCAATGTTGCACCGTTTGCTGCATGATATTGAAAAAAATATCACCAAATAAAAAACGCCTTAGTAAGTATTTTTCagcttttcattctgactgca  <  1:277703‑M2/250‑21 (MQ=255)
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                    
CGGGCGAAATATCATTACTACGCCACGCCAGTTGAACTGGTGCCGCTGTTAGAGGAAAAATCTTCATGGATGAGCCATGCCGCGCTGGTGTTTGGTCGCGAAGATTCCGGGTTGACTAACGAAGAGTTAGCGTTGGCTGACGTTCTTACTGGTGTGCCGATGGTGGCGGATTATCCTTCGCTCAATCTGGGGCAGGCGGTGATGGTCTATTGCTATCAATTAGCAACATTAATACAACAACCGGCGAAAAGTGATGCAACGGCAGACCAACATCAACTGCAAGCTTTACGCGAACGAGCCATGACATTGCTGACGACTCTGGCAGTGGCAGATGACATAAAACTGGTCGACTGGTTACAACAACGCCTGGGGCTTTTAGAGCAACGAGACACGGCAATGTTGCACCGTTTGCTGCATGATATTGAAAAAAATATCACCAAATAAAAAACGCCTTAGTAAGTATTTTTC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_000913/4641185‑4641652

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