Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,641,544 A→G Q201Q (CAA→CAG 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,5440AG98.0% 157.6 / ‑4.4 51Q201Q (CAA→CAGyjtDputative tRNA/rRNA methyltransferase YjtD
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base G (21/28);  total (22/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.49e-01

GGCTGACGTTCTTACTGGTGTGCCGATGGTGGCGGATTATCCTTCGCTCAATCTGGGGCAGGCGGTGATGGTCTATTGCTATCAATTAGCAACATTAATACAACAACCGGCGAAAAGTGATGCAACGGCAGACCAACATCAACTGCAAGCTTTACGCGAACGAGCCATGACATTGCTGACGACTCTGGCAGTGGCAGATGACATAAAACTGGTCGACTGGTTACAACAACGCCTGGGGCTTTTAGAGCAACGAGACACGGCAATGTTGCACCGTTTGCTGCATGATATTGAAAAAAATATCACCAAATAAAAAACGCCTTAGTAAGTATTTTTC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_000913/4641319‑4641652
                                                                                                                                                                                                                                 |                                                                                                                                                                                                      
ggCTGACGTTCTTACTGGTGTGCCGATGGTGGCAGATTATCCTTCGCTCAATCTGGGGCAGGCAGTGATGGTCTATTGCTATCAATTAGCAACATTAATACAACAACCGACGAAAAGTGACACAACGGCAGACCAACATCAACTGCAAGCTTTACGCGAACGAGTCATGGCATTGCTGACGACTCTGGCAGTGGCAGATGACATAAAACTGGTCGACTGGTTACAGCAACGCCTGGGGCTTTTAGAGCAAc                                                                                                                                                                               <  1:123695/251‑1 (MQ=255)
       gTTCTTACTGGTGTGCCGATGGTGGCAGATTATCCTTCGCTCAATCTGGGGCAGGCAGTGATGGTCTATTGCTATCAATTAGCAACATTAATACAACAACCGACGAAAAGTGACACAACGGCAGACCAACATCAACTGCAAGCTTTACGCGAACGAGTCATGGCATTGCTGACGACTCTGGCAGTGGCAGATGACATAAAACTGGTCGACTGGTTACAGCAACGCCTGGGGCTTTTAGAGCAACGAGacac                                                                                                                                                                        <  1:214402/251‑1 (MQ=255)
            tACTGGTGTGCCGATGGTGGCAGATTATCCTTCGCTCAATCTGGGGCAGGCAGTGATGGTCTATTGCTATCAATTAGCAACATTAATACAACAACCGACGAAAAGTGACACAACGGCAGACCAACATCAACTGCAAGCTTTACGCGAACGAGTCATGGCATCGCTGACGACTCTGGCAGTGGCAGATGACATAAAACTGGTCGACTGGTTACAGCAACGCCTGGGGCTTTTAGAGCAACGAGACACgg                                                                                                                                                                      <  1:251473/248‑1 (MQ=255)
                 gtgtGCCGATGGTGGCAGATTATCCTTCTCCCAATCTGGGGAAGGCAGTGATGGTCTATTGTTATCAATTAGCAACATTAATACAACAACCGACGAAAAGTGACACAACGGCAGACCAACATCAACTGCAAGCTTTTCGCGAACGAGTTATGGCATTGCTGACGACCCTGGCAGTGGCAGATGACATAAAACTGGTCGACTGGTTACAGCAACGCCAGGGGCTTTTAGAGCAACGAGACACGGCAATGTTg                                                                                                                                                              <  2:54510/251‑1 (MQ=255)
                     gCCGATGGTGGCAGATTATCCTTCGCTCAATCTGGGGCAGGCAGTGATGGTCTATTGCTATCAATTAGCAACATTAATACAACAACCGACGAAAAGTGACACAACGGCAGACCAACATCAACTGCAAGCTTTACGCGAACGAGTCATGGCATTGCTGACGACTCTGGCAGTGGCAGATGACATAAAACTGGTCGACTGGTTACAGCAACGCCTGGGGCTTTTAGAGCAACGAGACACGGCAATGTTGCAcc                                                                                                                                                          <  1:168249/251‑1 (MQ=255)
                       cGATGGTGGCAGATTATCCTTCGCTCAATCTGGGGCAGGCAGTGATTGTCTATTGCTATCAATTAGCAACATTAATACAACAACCGACGAAAAGTGACACAACGGCAGACCAACATCAACTGCAAGCTTTACGCGAACGAGTCATGGCATTGCTGACGACTCTGGCAGTGGCAGATGACATAAAACTGGTCGACTGGTTACAGCAACGCCTGGGGCTTTTAGAGCAACGAGACACGGCAATGTTGCACCGt                                                                                                                                                        <  1:456719/251‑1 (MQ=255)
                          ggggggcagATTATCCTTCGCTCAATCTGGGGCAGGCAGTGATGGTCTATTGCTATCAATTAGCCACATTAATACAACAACCGACGAAAAGTGACACAACGGCAGACCAACATCAACTGCAAGCTTTACGCGAACGTGACAAGGCATTGCTGTAGACTCGAGAAGTGGCAGATGGCCAAAAGCGGGACGACTGATAACATCAACGCCTGGGGCTTTTCAAGCAACGAGACACGGCAATATGGCACCGttt                                                                                                                                                      >  2:120036/9‑250 (MQ=255)
                           ggtggCAGATTATCCTTCGCTCAATCTAGGGCAGGCAGTGATGGTCTATTGCTATCAATTAGCAGCATTAATACAACAACCGACGAAAAGTGACACAACGGCAGACCAACATCAACTGCAAGCTTTACGCGAACGAGTCATGGCATTGCTGACGACTCTGGCAGTGGCAGATGACATAAAACTGGTCGACTGGTTACAGCAACGCCTGGGGCTTTTAGAGCAACGAGACACGGCAATGTTGCACCGTTtgc                                                                                                                                                    <  2:338181/251‑1 (MQ=255)
                              ggcagattatactgcgacaaaTCTGGGGCAGGCAGTGATGGTCTATGTCTATCAATTAGCAACATTAATACAACAACCGACGAAAAGTGACACAACGGCAGACCAACATCAACTGCAAGCTTTACGCGAACGAGTCATGGCATTGCTGACGACTCTGGCAGTGGCAGATGACATAAAACTGGTCGACTGGTTACAGCAACGCCTGGGGCTTTTAGAGCAACGAGACACGGCAATGTTGCACCGTTtgcgg                                                                                                                                                  <  2:224/231‑3 (MQ=255)
                              ggcagATTATCCTTCGCTCAATCTGGGGCAGGCAGTGATGGTCTATTGCTATCAATTAGCAACATTAATACAACAACCGACGAAAAGTGACACAACGGCAGACCAACATCAACTGCAAGCTTTACGCGAACGAGTCATGGCATTGCTGACGACTCTGGCAGTGGCAGATGACATAAAACTGGTCGACTGGTTACAGCAACGCCTGGGGCTTTTAGAGCAACGAGACACGGCAATGTTGCACCGTTtgctg                                                                                                                                                  <  1:1428/246‑1 (MQ=255)
                                cagATTATCCTTCGCTCAATCTGGGGCAGGCAGTGATGGTCTATTGCTATCAATTAGCAACATTAATACAACAACCGACGAAAAGTGACACAACGGCAGACCAACAGCAACTGCAAGCTTTACGCGAACGAGTCATGGCATTGCTGACGACTCTGGCAGTGGCAGATGACATAAAACTGGTCGACTGGTTACAGCAACGCCTGGGGCTTTTAGAGCAACGAGACACGGCAATGTTGCACCGTTTGCTGCAt                                                                                                                                               <  2:232419/249‑1 (MQ=255)
                                   aTTATCCTTCGCTCAATCGGGGACAGGCAGTGATGGTCTATTGCTAACAATTAGCAACATTAATACAACAACCGACGAAAAGTGACACAACGGCAGACCAACATCAACTGCAAGCTTTACGCGAACGAGTCATGGCATTGCTGACGACTCTGGCAGTGGCAGATGACATAAAACTGGTCGACTGGTTACAGCAACGCCTGGGGCTTTTAGAGCAACGAGACACGGCAATGTTGCACCGTTTGCTGCATGat                                                                                                                                            <  2:20907/251‑1 (MQ=255)
                                         cTTCGCTCAATCTGGGGCAGGCAGTGATGGTCTATTGCTATCAATTAGCAACATTAATACAACAACCGACGAAAAGTGACACAACGGCAGACCAACATCAACTGCAAGCTTTACGCGAACGAGTCATGGCATTGCTGACGACTCTGGCAGTGGCAGATGACATAAAACTGGTCGACTGGTTACAGCAACGCCTGGGGCTTTTAGAGCAACGAGACACGGCAATGTTGCACCGTTTGCTGCATGATATTGa                                                                                                                                       >  1:413138/1‑250 (MQ=255)
                                                cAATCTGGGGCAGGCAGTGATGGTCTATTGCTATCAATTAGCAACATTAATACAACAACCGACGAAAAGTGACACAACGGCAGACCAACATCAACTGCAAGCTTTACGCGAACGAGTCATGGCATTGCTGACGACTCTGGCAGTGGCAGATGACATAAAACTGGTCGACTGGTTACAGCAACGCCTGGGGCTTTTAGAGCAACGAGACACGGCAATGTGGCACCGTTTGCTGCATGATAATGaaaaaaa                                                                                                                                 <  2:235414/249‑1 (MQ=255)
                                                     tGGGGCAGGCAGTGATGGTCTATTGCTATCAATTAGCAACATGAATACAACAACCGACGAAAAGTGACACAACGGCAGACCAACATCAACTGCAAGCTTTACGCGAACGAGTCATGGCATTGCTGACGACTCTGGCAGTGGCAGATGACATAAAACTGGTCGACTGGTTACAGCAACGCCTGGGGCTTTTAGAGCAACGAGACACGGCAATGTTGCACCGTTTGCTGCATGATATTGAAAAAAATATCAc                                                                                                                           <  2:94365/250‑1 (MQ=255)
                                                           aggcaatGATGGTCTATTGCTCTCTAGTCGCCACATGAATACAACAACCGTCGAAAAGCGACACAACGGCAGACCAACATCAACTGCAAGCATTACGCGAACGAGTCATGGCATTGCCGACGACTCTGGCAGTGGCATATGACATAAAACGGGTCGACTGGTTACAGCAACGCTTGGGGCTTTTAGAGCAACGAGACACGGCAATGTTGCACCGTTTGCTGCATGATCGTGAAAAAAATACCACCCAATaa                                                                                                                    <  2:179066/245‑1 (MQ=255)
                                                               agTGATGGTCTATTGCTATCCATTAGCAACATTAATACAACAACCGACGAAAAGTGACACAACGGCAGACCAACATCAACTGCAAGCTTTAAGCGAACGAGTCATGGCATTGCTGACGACTCTGGCAGTGGCAGATGACATAAAACTGGTCGACTGGTTACAGCAACGCCTGGGGCTTTAAGAGCAAAGAGACACGGCACTGTTGCACCGTTTGCTGCATGATATTGAAAAAAATATCACCAAATaaaaa                                                                                                                 >  2:42340/2‑250 (MQ=255)
                                                                  gATGGTCTATTGCTATCAATTAGCAACATTAATACAACAACCGACGAAAAGTGACACAACGGCAGACCAACATCAACTGCAAGCTTTACGCGAACGAGTCATGGCATTGCTGACGACTCTGGCAGTGGCAGATGACATAAAACTGGTCGACTGGTTACAGCAACGCCTGGGGCTTTTAGAGCAACGAGACACGGCAATGTTGCACCGTTTGCTGCATGatat                                                                                                                                          <  2:254058/222‑1 (MQ=255)
                                                                  gATGGTCTATTGCTATCAATTAGCAACATTAATACAACAACCGACGAAAAGTGACACAACGGCAGACCAACATCAACTGCAAGCTTTACGCGAACGAGTCATGGCATTGCTGACGACTCTGGCAGTGGCAGATGACATAAAACTGGTCGACTGGTTACAGCAACGCCTGGGGCTTTTAGAGCAACGAGACACGGCAATGTTGCACCGTTTGCTGCATGatat                                                                                                                                          >  1:254058/1‑222 (MQ=255)
                                                                  gATGGTCTATTGCTATCAATTAGCAACATTAATACAACAACCGACGAAAAGTGACACAACGGCAGACCAACATCAACTGCAAGCTTTACGCGAACGAGTCATGGCATTGCTGACGACTCTGGCAGTGGCAGATGACATAAAACTGGTCGACTGGTTACAGCAACGCCTGGGGCTTTTAGAGCAACGAGACACGGCAATGTTGCACCGTTTGCTGCATGATATTGAAAAAAATATCACCAAATAAAAAACGc                                                                                                             >  1:146385/1‑251 (MQ=255)
                                                                      gTCTATTGCTATCAATTAGCAACATTAATACAACAACCGACGAAAAGTGACACAACGGCAGACCAACATCAACTGCAAGCTTTACGCGAACGAGTCATGGCATTGCTGACGACTCTGGCAGTGGCAGATGACATAAAACTGGTCGACTGGTTACAGCAACGCCTGGGGCTTTTAGAGCAACGAGACACGGCAATGTTGCACCGTTTGCTGCATGATATTGAAAAAAATATCACCAAATAAAAAACGCCTTa                                                                                                         <  1:259067/251‑1 (MQ=255)
                                                                       tCTATTGCTATCAATTAGCAACATTAATACAACAACCGACGAAAAGTGACACAACGGCAGACCAACATCAACTGCAAGCTTTACGCGAACGAGTCATGGCATTGCTGACGACTCTGGCAGTGGCAGATGACATAAAACTGGTCGACTGGTTACAGCAACGCCTGGGGCTTTTAGAGCAACGAGACACGGCAATGTTGCACCGTTTGCTGCATGATATTGAAAAAAATATCACCAAATAAAAAACGCCtt                                                                                                          >  1:149121/1‑249 (MQ=255)
                                                                          atttcTATCAATTAGCAACATGAATCAAACAACCGGCGAAAAGTGATGCAACGGCAGACCAACATCAACTGCAAGCTTGACGCGAACGAGCCATGACATTGCTGACGACTCTGGCAGTGGCAGATGACATAAAACTGGTCGACTGGTTACAACAACGCCTGGGGCTTTTAGAGCAACGAGACACGGCAATGTTGCACCGTTTGCTGCATGATATTGAAAAAAATATCACCAAATAAAAAACGCCTTagt                                                                                                       <  1:186998/245‑1 (MQ=255)
                                                                           ttGCTATCAATTAGCAACATTAATCCAACACCCGACGACAAGTGACACAACGGCAGACCAACATCAACTGCAAGATATACGCGCAAGAGTCATCGCATTGCTGACGACTCTGGCCGTGGCAGATGACATAAAACTGGTCGACTGGTTACAGCAACGCCGGGGGCTTTTAGAGCAACGCGACAAGGCACTGTTGCACCGATTGCTGCCTGATATTGAAAAAAATATAACCAAATAAAAAACGCATTagtaa                                                                                                     >  2:468293/1‑250 (MQ=255)
                                                                            tGCTATCAATTGTCAACATTAATCCAACAACCGAAGAAAAGTGACACAACGGCAGACCAACATCAACTGCAAGCTTTACGCGAACGAGTCATGGCATTGCTGACGACTCTGGCAGTGGCAGATGACATAAAACGGGTCGACTGGTTACAGCAACGCCTGGGGCTTTTAGAGCAACGAGACACGGCAATGTTGCACCGTTTGCTGCATGATATTGAAAAAAATATCACCAAATAAAAAACGACTTagtaag                                                                                                    <  2:371669/250‑1 (MQ=255)
                                                                            tGCTATCAATTAGCAACATTAATACAAAAACCGTCGAAAAGTGACACAACGGCAGACCAACATCACCTGCAAGCTTTACGCGAACGAGTCATGGCATTGCTGACGACTCTGGCAGTGGCAGATGACATAAAACTGGTCGAATGGTTACAGCAACGCCAGGGGCTTTTAGAGAAACGAGACACGGCAATGTGGCACCGTTTGCTGCATGGTATTGAAAAAAATATACCCAAATAAAAAAAGCATTagta                                                                                                      >  2:391447/1‑248 (MQ=255)
                                                                              cTATCAATTAGAAACATTAATACAACAACCGACGAAAAGTGACACAACGGCAGACCAACATCAACTGCAAGCTTTACGCGAACGAGTCATGGCATTGCTGACGACTCTGGAAGTGGCAGATGACATAAAACTGGTCGACTGGTTACAGCAACGCCTGGGGCTTTTAGAGCAACGAGAAACGGCAATGTTGCACCGTTTGCTGCCTGATATTGAAAAAACTATAACCAAATAAAAAAAGCATTAGTAAGTAt                                                                                                 >  2:462414/1‑251 (MQ=255)
                                                                               tATCAATTAGCAACATTAATCCAACAACCGACGCAAAGTGACACAACGGCAGACCAACATCAACTGCAAGCTTTACGCGAACGAGTCATGGCATTGCTGACTACTCTGGCAGTGGCAGCGGACATAACACTGGTCGACTGGTTACAGCAACGCCTGGGGCTTTTAGAGCAACGAGACACGGCAATGTTGCACCGTTTGCTGCATGCTATTGACAACAATATCACCAAATAAAAACTGCCTTAGTAAGTAtt                                                                                                >  1:263225/1‑251 (MQ=255)
                                                                                  aaaTTAGCAACATTAATACAACAACCGACGAAAAGTGACACAACGGCAGACCAACATCAACTGCAAGCTTTACGCGAACGAGTCATGGCATTGCTGACGACTCTGGCAGTGGCAGATGACATAAAACTGGTCGACTGGTTACAGCAACGCCTGGGGCTTTTAGAGCAACGAGACACGGCAATGTTGCACCGTTTGCTGCATGATATTGAAAAAAATATCACCAAATAAAAAACGCCTTAGTAAGTAttttt                                                                                             <  2:149233/250‑1 (MQ=255)
                                                                                   aaTTAGCAACATTCATACAACAACCGGCGAAAAGTGAAACAACGGCAGACCAACATCAACTGCAATCTTTACGCGAACGAGTCATGGCAGTGCTGACGACTCTGGCAGTGGCAGATGACATAAAACTGGTCGACGGTTTACAGCAACGTCTGTGGCTTTTAGAGCAACGAGACACGGCAATGTTGCACCGTTTGCTGCATGATATTGAAAAAAATATCCCCAAATAAAAAACGCCTTAGTAAGCATTTTTc                                                                                            <  2:433125/251‑1 (MQ=255)
                                                                                    aTTAGCAACATTAATACAACAACCGACGAAAAGTGACACAACGGCAGACCAACATCAACTGCAAGCTTTACGCGAACGAGTCATGGCATTGCTGACGACTCTGGCAGTGGCAGATGACATAAAACTGGTCGACTGGTTACAGCAACGCCTGGGGCTTTTAGAGCAACGAGACACGGCAATGTTGCACCGTTTGCTGCATGATATTGAAAAAAATATCACCAAATAAAAAACGCCTTAGTAAGTATTTTTCa                                                                                           >  1:240347‑M2/1‑250 (MQ=255)
                                                                                       aGCTACAGTAATACAACAACCGAAGAAAAGTGACACAACGGCAGACCAACATCAACTGCAAGCTTTACGCGAACGAGTCATGGCATTGCTTACGACTCTGGCAGTGGCAGATGACATAAAACTGGTCGACTGGTTACAGCAACGCCTGGGGCGTTTAGAGCAACGAGACACGGCAATGTTGCACCGTTTGCTGCATGATATTGAAAAAAATATCACCAAATAAAAAACGCCTTAGTAAGTATTTTTCagc                                                                                         >  1:468481‑M2/1‑247 (MQ=255)
                                                                                         cAACATTAATACAACAACCGACGAAAAGTGACACAACGGCAGACCAACATCCACTGCAAGATTTACGCGAACGAGTCATGGCATTGCTGACGACTATGGCAGTCGCAGATGACATAAAACTGGTCGACTGGTTACAGCAACGCCTGGGGCTTTTAGAGCAACGCGACCCGGCAATGTTGCACCGTTTGCTGCATGATATTGAAAAAAATATCACCAAATAAAAAACGCCTTAGTAAGTATTTTTCagcttt                                                                                      <  1:498161‑M2/251‑7 (MQ=255)
                                                                                            aatTAATACAACAACCGACGAAAAGTGACACAACGGCAGACCAACATCAACTGCAAGCTGTACGCGAACGAGTCATGGCATTGCTGACGACTCTGGCAGTGGCAGATGACATAAAACTGGTCGACTGGTTACAGCAACGCCTGGGGCTTTTAGAGCAACGAGACACGGCAATGTTGCACCGTTTGCTGCATGATATTGAAAAAAATATCACCAAATAAAAAACGCCTTAGTAAGTATTTTTCagcttttca                                                                                   >  2:146385‑M2/2‑242 (MQ=255)
                                                                                                aaTACAACAACCGACGAAAAGTGACACAACGGCAGACCAACATCAACTGCAAGCTTTACGCGAACGAGTCATGGCATTGCTGACGACTCTGGCAGTGGCAGATGACATAAAACTGGTCGACTGGTTACAGCAACGCCTGGGGCTTTTAGAGCAACGAGACCCGGCAATGTTGCACCGtt                                                                                                                                                       >  2:366734/1‑179 (MQ=255)
                                                                                                aaTACAACAACCGACGAAAAGTGACACAACGGCAGACCAACATCAACTGCAAGCTTTACGCGAACGAGTCATGGCATTGCTGACGACTCTGGCAGTGGCAGATGACATAAAACTGGTCGACTGGTTACAGCAACGCCTGGGGCTTTTAGAGCAACGAGACACGGCAATGTTGCACCGtt                                                                                                                                                       <  1:366734/179‑1 (MQ=255)
                                                                                                     aacaacCGACGAAAAGGGAGCCAACGGCAGACCAACATCAACTGCACGCTTTACGCGAACGATTCATGGCATTGCTGACGACGCTGGCAGTGGCAGATGACATAAATCTGGTCGACTGGTTATAGCAACGCCTGGGGCTTTTAGAGCAACGAGACACGGCAATGTTGCACCGTTTGCTGCATGATATTGAAAAATTTATCACCAAATAAAAAACGCCTGAGTAATTATTTTTCagcttttcattctgactg                                                                          >  2:502048‑M2/1‑233 (MQ=255)
                                                                                                                    gtgacacaacggcAGACCAACATCAACTGCAAGCTTTACGCGAACGAGTCATGGCATTGCTGACGACTCTGGCAGTGGCAGATGACATAAAACTGGTCGACTGGTTACAGCAACGCCTGGGGCTTTTAGAGCAACGAGACACGGCAATGTTGCACCGTTTGCTGCATGATATTGAAAAAAATATCACCAAATAAAAAACGCCTTAGTAAGTATTTTTCagcttttcattctgactgcaacgggcaatatgt                                                           <  1:51908‑M2/245‑34 (MQ=255)
                                                                                                                    gtgacacaacggcAGACCAACATCAACTGCAAGCTTTACGCGAACGAGTCATGGCATTGCTGACGACTCTGGCAGTGGCAGATGACATAAAACTGGTCGACTGGTTACAGCAACGCCTGGGGCTTTTAGAGCAACGAGACACGGCAATGTTGCACCGTTTGCTGCATGATATTGAAAAAAATATCACCAAATAAAAAACGCCTTAGTAAGTATTTTTCagcttttcattctg                                                                              <  1:472914‑M2/226‑15 (MQ=255)
                                                                                                                    gtgacacaacggcAGACCAACAGCAACTGCAAGCTTTACGCGAACGAGTCATGGCATTGCTGACGACTCTGGCAGTGGCAGATGACATAAAACTGGTCGACTGGTTACAGCAACGCCTGGGGCTTTTAGAGCAACGAGACACGGCAATGTTGCACCGTTTGCTGCATGATATTGAAAAAAATAACACCAAATAAAAAACGCCTTAGTAAGTATTTTTCagcttttcattctg                                                                              >  2:472914‑M2/7‑218 (MQ=255)
                                                                                                                               gCAGACCAACATCAACTGCAAGCTTTACGCGAACGAGTCATGGCATTGCTGACGACTCTGGCAGTGGCAGATGACATAAAACTGGTCGACTGGTTACAGCAACGCCTGGGGCTTTTAGAGCAACGAGACACGGCAATGTTGCACCGTTTGCTGCATGATATTGAAAAAAATATCACCAAATAAAAAACGCCTTAGAAAGTATTTTTCagcttttcattctgactgcaacgggcaatatgtctctgtgtgg                                                 <  2:330122‑M2/250‑44 (MQ=255)
                                                                                                                                cAGACCAACATCAACTGCAAGCTTTACGCGAACGAGTCATGGCATTGCTGACGACTCTGGCAGTGGCAGATGACATAAAACTGGTCGACTGGTTACAGCAACGCCTGGGGCTTTTAGAGCAACGAGACACGGCAATGTTGCACCGTTTGCTGCATGATATTGAAAAAAATATCACCAAATAAAAAACGCCTTAGTAAGTATTTTTCagcttttcattctgactgcaacgggcaatatgtttctgtgtggat                                               >  1:276830‑M2/1‑206 (MQ=255)
                                                                                                                                 aGACCAACATCAACTGCAAGCTTTACGCGAACGAGTCATGGCATTGCTGACGACTCTGGCAGTGGCAGATGACATAAAACTGGTCGACTGGTTACAGCAACGCCTGGGGCTTTTAGAGCAACGAGACACGGCAATGTTGCACCGTTTGCTGCATGATATTGAAAAAAATATCACCAAATAAAAAACGCCTTAGTAAGTATTTTTCagcttttcattctgactgcaacgggcaatatgtctctgtgtggatt                                              >  1:330122‑M2/1‑205 (MQ=255)
                                                                                                                                 aGACCAACATCAACTGCAAGCTTTACGCGAACGAGTCATGGCATTGCTGACGACTCTGGCAGTGGCAGATGACATAAAACTGGTCGACTGGTTACAGCAACGCCTGGGGCTTTTAGAGCAACGAGACACGGCAATGTTGCACCGTTTGCTGCATGATATTGAAAAAAATATCACCAAATAAAAAACGCCTTAGTAAGTATTTTTCagcttttcattctgactgcaacgggcaatatgtctctgtgtggatt                                              >  1:113151‑M2/1‑205 (MQ=255)
                                                                                                                                           cAACTGCAAGCTTTACGCGAACGAGTCATGGCATTGCTGACGACTCTGGCAGTGGCAGATGACATAAAACTGGTCGACTGGTTACAGCAACGCCTGGGGCTTTTAGAGCAACGAGACACGGCAATGTTGCACCGTTTGCTGCATGATATTGAAAAAAATATCACCAAATAAAAAACGCCTTAGTAAGTATTTTTCagcttttcattctgactgcaacgggcaatatgtctctgtgtggattaaaaaaaga                                     >  2:337024‑M2/1‑195 (MQ=255)
                                                                                                                                               tGCAAGCTTTACGCGAACGAGTCATGGCATTGCTGACGACTCTGGCAGTGGCAGATGACATAAAACTGGTCGACTGGTTACAGCAACGCCTGGGGCTTTTAGAGCAACGAGACACGGCAATGTTGCACCGTTTGCTGCATGATATTGAAAAAAATATCACCAAATAAAAAACGCCTTAGTAAGTATTTTTCagcttttcattctgactgcaacgggcaatatgtctctgtgtggattaaaaaaagagtgtc                                <  2:266953‑M2/251‑61 (MQ=255)
                                                                                                                                                             gAACGAGTCATGGCATTGCTGACGACTCTGGCAGTGGCAGATGACATAAAACTGGTCGACTGGTTACAGCAACGCCTGGGGCTTTTAGAGCAACGAGACACGGCAATGTTGCACCGTTTGCTGCATGATATTGAAAAAAATATCACCAAATAAAAAACGCCTTAGTAAGTATTTTTCagcttttcattctgactgcaacgggcaatatgtctctgtgtggattaaaaaaagagtgtctgatagcagcttct                  <  1:492670‑M2/251‑75 (MQ=255)
                                                                                                                                                                            ttGCTGACGACTCTGGCAGTGGCAGATGACATAAAACGGGTCGACTGGTTACAGCAACGCCTGGGGCTTTAAGAGCAACGAGACACGGCAATGTTGCACCGTTAGCTGCATGATATTGAAAAAAATATCACCAAATAAAAAACGCCTTAGTAAGTATTTTTCagcttttcattctgactgcaacgggcaatatgtctctgtgtggagtaaaaaaagagtgtctgatagcagcttcggaactggttccctg    >  2:209184‑M2/1‑162 (MQ=255)
                                                                                                                                                                             tGCTGACGACTCTGGCAGTGGCAGATGACATAAAACTGGTCGACTGGTTACAGCAACGCCTGGGGCTTTTAGAGCAACGAGACACGGCAATGTTGCACCGTTTGCTGCATGATATTGAAAAAAATATCACCAAATAAAAAACGCCTTAGTAAGTATTTTTCagcttttcattctgactgcaacgggcaatatgtctctgtgtggattaaaaaaagagtgtctgatagcagcttctgaactggttacctgcc  <  1:446003‑M2/251‑91 (MQ=255)
                                                                                                                                                                                                             aaaCTGGTCGACTGGTTACAGCAACGCCTGGGGCTTTTAGAGCAACGAGACACGGCAATGTTGCACCGTTTGCTGCATGATATTGAAAAAAATATCACCAAATAAAAAACGCCTTAGTAAGTATTTTTCagcttttcattctgactgcaactggcaatatgtctctgtgtggattaaaaaaagagtgtctg                              >  1:132969‑M2/1‑129 (MQ=255)
                                                                                                                                                                                                             aaaCTGGACGACTGGTTACAGCAACGCCTGGGGCTTTTAGAGCAACGAGACACGGCAATGTTGCACCGTTTGCTGCATGATATTGAAAAAAATATCACCAAATAAAAAACGCCTTAGTAAGTATTTTTCagcttttcattctgactgcaacgggcaatatgtctctgtgtggattaaaaaaagagtgtctg                              <  2:132969‑M2/191‑63 (MQ=255)
                                                                                                                                                                                                                                 |                                                                                                                                                                                                      
GGCTGACGTTCTTACTGGTGTGCCGATGGTGGCGGATTATCCTTCGCTCAATCTGGGGCAGGCGGTGATGGTCTATTGCTATCAATTAGCAACATTAATACAACAACCGGCGAAAAGTGATGCAACGGCAGACCAACATCAACTGCAAGCTTTACGCGAACGAGCCATGACATTGCTGACGACTCTGGCAGTGGCAGATGACATAAAACTGGTCGACTGGTTACAACAACGCCTGGGGCTTTTAGAGCAACGAGACACGGCAATGTTGCACCGTTTGCTGCATGATATTGAAAAAAATATCACCAAATAAAAAACGCCTTAGTAAGTATTTTTC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_000913/4641319‑4641652

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