Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,641,488 A→G T183A (ACA→GCA)  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,4880AG98.2% 167.9 / ‑4.5 55T183A (ACA→GCA) yjtDputative tRNA/rRNA methyltransferase YjtD
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base G (26/28);  total (26/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 = 3.22e-01

ATGGATGAGCCATGCCGCGCTGGTGTTTGGTCGCGAAGATTCCGGGTTGACTAACGAAGAGTTAGCGTTGGCTGACGTTCTTACTGGTGTGCCGATGGTGGCGGATTATCCTTCGCTCAATCTGGGGCAGGCGGTGATGGTCTATTGCTATCAATTAGCAACATTAATACAACAACCGGCGAAAAGTGATGCAACGGCAGACCAACATCAACTGCAAGCTTTACGCGAACGAGCCATGACATTGCTGACGACTCTGGCAGTGGCAGATGACATAAAA‑CTGGTCGACTGGTTACAACAACGCCTGGGGCTTTTAGAGCAACGAGACACGGCAATGTTGCACCGTTTGCTGCATGATATTGAAAAAAATATCACCAAATAAAAAACGCCTTAGTAAGTATTTTTC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_000913/4641250‑4641652
                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                               
aTGGATGAGCCATGCCGCGCTGGTGTTTGGTCGCGAAGATTCCGGGCTGACTAACGAAGAGTTAGCGTTGGCTGACGTTCTTACTGGTGTGCCGATGGTGGCAGATTATCCTTCGCTCAATCTGGGGCAGGCAGTGATGGTCTATTGCTATCAATTAGCAACATTAATACAACAACCGACGAAAAGTGACACAACGGCAGACCAACATCAACTGCAAGCTTTACGAGAACGAGTCATGGCATTGCTgtcg                                                                                                                                                                                                                                      >  2:347005/1‑247 (MQ=255)
        gCCATGCCGCGCTGGTGTTTGGTCGCGAAGATTCCGGGCTGACTAACGAAGAGTTAGCGTTGGCTGACGTTCTTACTGGTGTGCCGATGGTGGCAGATTATCCTTCGCTCAATCTGGGGCAGGCAGTGATGGTCTATTGCTATCAATTAGCAACATTAATACAACAACCGACGAAAAGTGACACAACGGCAGACCAACATCAACTGCAAGCTTTACGCGAACGAGTCATGGCATTGCTGACGActc                                                                                                                                                                                                                                  >  1:497378/1‑246 (MQ=255)
        gCCATGCCGCGCTGGTGTTTGGTCGCGAAGATTCCGGGCTGACTAACGAAGAGTTAGCGTTGGCTGACGTTCTTACTGGTGTGCCGATGGTGGCAGATTATCCTTCGCTCAATCTGGGGCAGGCAGTGATGGTCTATTGCTATCAATTAGCAACATTAATACAACAACCGACGAAAAGTGACACAACGGCAGACCAACATCAACTGCAAGCTTTACGCGAACGAGTCATGGCATTGCTGACGActc                                                                                                                                                                                                                                  <  2:497378/246‑1 (MQ=255)
                        gtTTGGTCGCGAAGATTCCGGGCTGACTAACGAAGAGTTAGCGTTGGCTGACGTTCTTACTGGTGTGCCGATGGTGGCAGATTATCCTTCGCTCAATCTGGGGCAGGCAGTGATGGTCTATTGCTATCAATTAGCAACATTAATACAACAAACGACGAAAAGTGACACAACGGCAGACCAACATCAACTGCAAGCTTTACGCGAACGAGTCATGGCATTGCTGACGACTCTGGCAGTGGCAGATGACATa                                                                                                                                                                                                              >  1:20907/1‑250 (MQ=255)
                                         cctggctGACGAACGAAGAGTTAGCGATGGCTGTCGTTCTTGCTGGGGTGCCGATGGTGGCAGATTATCCTTCGCTCAATCTGGGGCAGGCAGTGATGGTCTATTGCTATCAATTAGCAACATTAATACAACAACCGACGAAAAGTGACACAACGGCAGACCAACATCAACTGCAAGCTTTACGCGAACGAGTCATGGCATTGCTGACGACTCTGGCAGTGGCAGATGACATAAAA‑CTGGTCGACTGGtt                                                                                                                                                                                            <  2:334224/244‑1 (MQ=255)
                                          cGGGCTGACTAACGAAGAGTTAGCGTTGGCTGACGTTCTTACTGGTGTGCCGATGGTGGCAGATTATCCTTCGCTCAATCTGGGGCAGGCAGTGATGGTCTATTGCTATCAATTAGCAACATTAATACAACAACCGACGAAAAGTGACACAACGGCAGACCAACATCAAATGCAAGCTTTACGAGAACGAGGCATGGCATTGCTGACGACTCTGGCAGTGGCAGATGA‑ATAAAAACTGGTCGACTGGTTa                                                                                                                                                                                           >  2:308384/1‑250 (MQ=255)
                                             gctGACTAACGAAGAGTTAGCGTTGGCTGACGTTCTTACTGGTGTGCCGATGGAGGCAGATTATCCTTCGCTCAATCTGGGGCAGGCAGTGATGGTCTATTGCTATCAATTAGCAACATTAATACAACAACCGACGAAAAGTGACACAACGGCAGACCAACATCAACTGCAAGCTTTACGCGAACGAGTCATGGCATTGCTGACGACTCTGGCAGTGGCAGATGACATAAAA‑CTGGTCGACTGGTTac                                                                                                                                                                                          >  2:276830/3‑248 (MQ=255)
                                              ctGACTAACGAAGAGTTAGCGTTGGCTGACGTTCTTACTGGTGTGCCGATGGTGGCAGATTATCCTTCGCTCAATCTGGGGCAGGCAGTGATGGTCTATTGCTATCAATTAGCAACATTAATACAACAACCGACGAAAAGTGACACAACGGCAGACCAACATCAACTGCAAGCTTTACGCGAACGAGTCATGGCATTGCTGACGACTCTGGCAGTGGCAGATGACATAAAA‑CTGGTCGACTGGTTacagc                                                                                                                                                                                       >  2:259067/2‑248 (MQ=255)
                                                                     ggCTGACGTTCTTACTGGTGTGCCGATGGTGGCAGATTATCCTTCGCTCAATCTGGGGCAGGCAGTGATGGTCTATTGCTATCAATTAGCAACATTAATACAACAACCGACGAAAAGTGACACAACGGCAGACCAACATCAACTGCAAGCTTTACGCGAACGAGTCATGGCATTGCTGACGACTCTGGCAGTGGCAGATGACATAAAA‑CTGGTCGACTGGTTACAGCAACGCCTGGGGCTTTTAGAGCAAc                                                                                                                                                               <  1:123695/251‑1 (MQ=255)
                                                                            gTTCTTACTGGTGTGCCGATGGTGGCAGATTATCCTTCGCTCAATCTGGGGCAGGCAGTGATGGTCTATTGCTATCAATTAGCAACATTAATACAACAACCGACGAAAAGTGACACAACGGCAGACCAACATCAACTGCAAGCTTTACGCGAACGAGTCATGGCATTGCTGACGACTCTGGCAGTGGCAGATGACATAAAA‑CTGGTCGACTGGTTACAGCAACGCCTGGGGCTTTTAGAGCAACGAGacac                                                                                                                                                        <  1:214402/251‑1 (MQ=255)
                                                                                 tACTGGTGTGCCGATGGTGGCAGATTATCCTTCGCTCAATCTGGGGCAGGCAGTGATGGTCTATTGCTATCAATTAGCAACATTAATACAACAACCGACGAAAAGTGACACAACGGCAGACCAACATCAACTGCAAGCTTTACGCGAACGAGTCATGGCATCGCTGACGACTCTGGCAGTGGCAGATGACATAAAA‑CTGGTCGACTGGTTACAGCAACGCCTGGGGCTTTTAGAGCAACGAGACACgg                                                                                                                                                      <  1:251473/248‑1 (MQ=255)
                                                                                      gtgtGCCGATGGTGGCAGATTATCCTTCTCCCAATCTGGGGAAGGCAGTGATGGTCTATTGTTATCAATTAGCAACATTAATACAACAACCGACGAAAAGTGACACAACGGCAGACCAACATCAACTGCAAGCTTTTCGCGAACGAGTTATGGCATTGCTGACGACCCTGGCAGTGGCAGATGACATAAAA‑CTGGTCGACTGGTTACAGCAACGCCAGGGGCTTTTAGAGCAACGAGACACGGCAATGTTg                                                                                                                                              <  2:54510/251‑1 (MQ=255)
                                                                                          gCCGATGGTGGCAGATTATCCTTCGCTCAATCTGGGGCAGGCAGTGATGGTCTATTGCTATCAATTAGCAACATTAATACAACAACCGACGAAAAGTGACACAACGGCAGACCAACATCAACTGCAAGCTTTACGCGAACGAGTCATGGCATTGCTGACGACTCTGGCAGTGGCAGATGACATAAAA‑CTGGTCGACTGGTTACAGCAACGCCTGGGGCTTTTAGAGCAACGAGACACGGCAATGTTGCAcc                                                                                                                                          <  1:168249/251‑1 (MQ=255)
                                                                                            cGATGGTGGCAGATTATCCTTCGCTCAATCTGGGGCAGGCAGTGATTGTCTATTGCTATCAATTAGCAACATTAATACAACAACCGACGAAAAGTGACACAACGGCAGACCAACATCAACTGCAAGCTTTACGCGAACGAGTCATGGCATTGCTGACGACTCTGGCAGTGGCAGATGACATAAAA‑CTGGTCGACTGGTTACAGCAACGCCTGGGGCTTTTAGAGCAACGAGACACGGCAATGTTGCACCGt                                                                                                                                        <  1:456719/251‑1 (MQ=255)
                                                                                               ggggggcagATTATCCTTCGCTCAATCTGGGGCAGGCAGTGATGGTCTATTGCTATCAATTAGCCACATTAATACAACAACCGACGAAAAGTGACACAACGGCAGACCAACATCAACTGCAAGCTTTACGCGAACGTGACAAGGCATTGCTGTAGACTCGAGAAGTGGCAGATGGCCAAAAG‑CGGGACGACTGATAACATCAACGCCTGGGGCTTTTCAAGCAACGAGACACGGCAATATGGCACCGttt                                                                                                                                      >  2:120036/9‑250 (MQ=255)
                                                                                                ggtggCAGATTATCCTTCGCTCAATCTAGGGCAGGCAGTGATGGTCTATTGCTATCAATTAGCAGCATTAATACAACAACCGACGAAAAGTGACACAACGGCAGACCAACATCAACTGCAAGCTTTACGCGAACGAGTCATGGCATTGCTGACGACTCTGGCAGTGGCAGATGACATAAAA‑CTGGTCGACTGGTTACAGCAACGCCTGGGGCTTTTAGAGCAACGAGACACGGCAATGTTGCACCGTTtgc                                                                                                                                    <  2:338181/251‑1 (MQ=255)
                                                                                                   ggcagattatactgcgacaaaTCTGGGGCAGGCAGTGATGGTCTATGTCTATCAATTAGCAACATTAATACAACAACCGACGAAAAGTGACACAACGGCAGACCAACATCAACTGCAAGCTTTACGCGAACGAGTCATGGCATTGCTGACGACTCTGGCAGTGGCAGATGACATAAAA‑CTGGTCGACTGGTTACAGCAACGCCTGGGGCTTTTAGAGCAACGAGACACGGCAATGTTGCACCGTTtgcgg                                                                                                                                  <  2:224/231‑3 (MQ=255)
                                                                                                   ggcagATTATCCTTCGCTCAATCTGGGGCAGGCAGTGATGGTCTATTGCTATCAATTAGCAACATTAATACAACAACCGACGAAAAGTGACACAACGGCAGACCAACATCAACTGCAAGCTTTACGCGAACGAGTCATGGCATTGCTGACGACTCTGGCAGTGGCAGATGACATAAAA‑CTGGTCGACTGGTTACAGCAACGCCTGGGGCTTTTAGAGCAACGAGACACGGCAATGTTGCACCGTTtgctg                                                                                                                                  <  1:1428/246‑1 (MQ=255)
                                                                                                     cagATTATCCTTCGCTCAATCTGGGGCAGGCAGTGATGGTCTATTGCTATCAATTAGCAACATTAATACAACAACCGACGAAAAGTGACACAACGGCAGACCAACAGCAACTGCAAGCTTTACGCGAACGAGTCATGGCATTGCTGACGACTCTGGCAGTGGCAGATGACATAAAA‑CTGGTCGACTGGTTACAGCAACGCCTGGGGCTTTTAGAGCAACGAGACACGGCAATGTTGCACCGTTTGCTGCAt                                                                                                                               <  2:232419/249‑1 (MQ=255)
                                                                                                        aTTATCCTTCGCTCAATCGGGGACAGGCAGTGATGGTCTATTGCTAACAATTAGCAACATTAATACAACAACCGACGAAAAGTGACACAACGGCAGACCAACATCAACTGCAAGCTTTACGCGAACGAGTCATGGCATTGCTGACGACTCTGGCAGTGGCAGATGACATAAAA‑CTGGTCGACTGGTTACAGCAACGCCTGGGGCTTTTAGAGCAACGAGACACGGCAATGTTGCACCGTTTGCTGCATGat                                                                                                                            <  2:20907/251‑1 (MQ=255)
                                                                                                              cTTCGCTCAATCTGGGGCAGGCAGTGATGGTCTATTGCTATCAATTAGCAACATTAATACAACAACCGACGAAAAGTGACACAACGGCAGACCAACATCAACTGCAAGCTTTACGCGAACGAGTCATGGCATTGCTGACGACTCTGGCAGTGGCAGATGACATAAAA‑CTGGTCGACTGGTTACAGCAACGCCTGGGGCTTTTAGAGCAACGAGACACGGCAATGTTGCACCGTTTGCTGCATGATATTGa                                                                                                                       >  1:413138/1‑250 (MQ=255)
                                                                                                                     cAATCTGGGGCAGGCAGTGATGGTCTATTGCTATCAATTAGCAACATTAATACAACAACCGACGAAAAGTGACACAACGGCAGACCAACATCAACTGCAAGCTTTACGCGAACGAGTCATGGCATTGCTGACGACTCTGGCAGTGGCAGATGACATAAAA‑CTGGTCGACTGGTTACAGCAACGCCTGGGGCTTTTAGAGCAACGAGACACGGCAATGTGGCACCGTTTGCTGCATGATAATGaaaaaaa                                                                                                                 <  2:235414/249‑1 (MQ=255)
                                                                                                                          tGGGGCAGGCAGTGATGGTCTATTGCTATCAATTAGCAACATGAATACAACAACCGACGAAAAGTGACACAACGGCAGACCAACATCAACTGCAAGCTTTACGCGAACGAGTCATGGCATTGCTGACGACTCTGGCAGTGGCAGATGACATAAAA‑CTGGTCGACTGGTTACAGCAACGCCTGGGGCTTTTAGAGCAACGAGACACGGCAATGTTGCACCGTTTGCTGCATGATATTGAAAAAAATATCAc                                                                                                           <  2:94365/250‑1 (MQ=255)
                                                                                                                                aggcaatGATGGTCTATTGCTCTCTAGTCGCCACATGAATACAACAACCGTCGAAAAGCGACACAACGGCAGACCAACATCAACTGCAAGCATTACGCGAACGAGTCATGGCATTGCCGACGACTCTGGCAGTGGCATATGACATAAAA‑CGGGTCGACTGGTTACAGCAACGCTTGGGGCTTTTAGAGCAACGAGACACGGCAATGTTGCACCGTTTGCTGCATGATCGTGAAAAAAATACCACCCAATaa                                                                                                    <  2:179066/245‑1 (MQ=255)
                                                                                                                                    agTGATGGTCTATTGCTATCCATTAGCAACATTAATACAACAACCGACGAAAAGTGACACAACGGCAGACCAACATCAACTGCAAGCTTTAAGCGAACGAGTCATGGCATTGCTGACGACTCTGGCAGTGGCAGATGACATAAAA‑CTGGTCGACTGGTTACAGCAACGCCTGGGGCTTTAAGAGCAAAGAGACACGGCACTGTTGCACCGTTTGCTGCATGATATTGAAAAAAATATCACCAAATaaaaa                                                                                                 >  2:42340/2‑250 (MQ=255)
                                                                                                                                       gATGGTCTATTGCTATCAATTAGCAACATTAATACAACAACCGACGAAAAGTGACACAACGGCAGACCAACATCAACTGCAAGCTTTACGCGAACGAGTCATGGCATTGCTGACGACTCTGGCAGTGGCAGATGACATAAAA‑CTGGTCGACTGGTTACAGCAACGCCTGGGGCTTTTAGAGCAACGAGACACGGCAATGTTGCACCGTTTGCTGCATGatat                                                                                                                          >  1:254058/1‑222 (MQ=255)
                                                                                                                                       gATGGTCTATTGCTATCAATTAGCAACATTAATACAACAACCGACGAAAAGTGACACAACGGCAGACCAACATCAACTGCAAGCTTTACGCGAACGAGTCATGGCATTGCTGACGACTCTGGCAGTGGCAGATGACATAAAA‑CTGGTCGACTGGTTACAGCAACGCCTGGGGCTTTTAGAGCAACGAGACACGGCAATGTTGCACCGTTTGCTGCATGatat                                                                                                                          <  2:254058/222‑1 (MQ=255)
                                                                                                                                       gATGGTCTATTGCTATCAATTAGCAACATTAATACAACAACCGACGAAAAGTGACACAACGGCAGACCAACATCAACTGCAAGCTTTACGCGAACGAGTCATGGCATTGCTGACGACTCTGGCAGTGGCAGATGACATAAAA‑CTGGTCGACTGGTTACAGCAACGCCTGGGGCTTTTAGAGCAACGAGACACGGCAATGTTGCACCGTTTGCTGCATGATATTGAAAAAAATATCACCAAATAAAAAACGc                                                                                             >  1:146385/1‑251 (MQ=255)
                                                                                                                                           gTCTATTGCTATCAATTAGCAACATTAATACAACAACCGACGAAAAGTGACACAACGGCAGACCAACATCAACTGCAAGCTTTACGCGAACGAGTCATGGCATTGCTGACGACTCTGGCAGTGGCAGATGACATAAAA‑CTGGTCGACTGGTTACAGCAACGCCTGGGGCTTTTAGAGCAACGAGACACGGCAATGTTGCACCGTTTGCTGCATGATATTGAAAAAAATATCACCAAATAAAAAACGCCTTa                                                                                         <  1:259067/251‑1 (MQ=255)
                                                                                                                                            tCTATTGCTATCAATTAGCAACATTAATACAACAACCGACGAAAAGTGACACAACGGCAGACCAACATCAACTGCAAGCTTTACGCGAACGAGTCATGGCATTGCTGACGACTCTGGCAGTGGCAGATGACATAAAA‑CTGGTCGACTGGTTACAGCAACGCCTGGGGCTTTTAGAGCAACGAGACACGGCAATGTTGCACCGTTTGCTGCATGATATTGAAAAAAATATCACCAAATAAAAAACGCCtt                                                                                          >  1:149121/1‑249 (MQ=255)
                                                                                                                                               atttcTATCAATTAGCAACATGAATCAAACAACCGGCGAAAAGTGATGCAACGGCAGACCAACATCAACTGCAAGCTTGACGCGAACGAGCCATGACATTGCTGACGACTCTGGCAGTGGCAGATGACATAAAA‑CTGGTCGACTGGTTACAACAACGCCTGGGGCTTTTAGAGCAACGAGACACGGCAATGTTGCACCGTTTGCTGCATGATATTGAAAAAAATATCACCAAATAAAAAACGCCTTagt                                                                                       <  1:186998/245‑1 (MQ=255)
                                                                                                                                                ttGCTATCAATTAGCAACATTAATCCAACACCCGACGACAAGTGACACAACGGCAGACCAACATCAACTGCAAGATATACGCGCAAGAGTCATCGCATTGCTGACGACTCTGGCCGTGGCAGATGACATAAAA‑CTGGTCGACTGGTTACAGCAACGCCGGGGGCTTTTAGAGCAACGCGACAAGGCACTGTTGCACCGATTGCTGCCTGATATTGAAAAAAATATAACCAAATAAAAAACGCATTagtaa                                                                                     >  2:468293/1‑250 (MQ=255)
                                                                                                                                                 tGCTATCAATTGTCAACATTAATCCAACAACCGAAGAAAAGTGACACAACGGCAGACCAACATCAACTGCAAGCTTTACGCGAACGAGTCATGGCATTGCTGACGACTCTGGCAGTGGCAGATGACATAAAA‑CGGGTCGACTGGTTACAGCAACGCCTGGGGCTTTTAGAGCAACGAGACACGGCAATGTTGCACCGTTTGCTGCATGATATTGAAAAAAATATCACCAAATAAAAAACGACTTagtaag                                                                                    <  2:371669/250‑1 (MQ=255)
                                                                                                                                                 tGCTATCAATTAGCAACATTAATACAAAAACCGTCGAAAAGTGACACAACGGCAGACCAACATCACCTGCAAGCTTTACGCGAACGAGTCATGGCATTGCTGACGACTCTGGCAGTGGCAGATGACATAAAA‑CTGGTCGAATGGTTACAGCAACGCCAGGGGCTTTTAGAGAAACGAGACACGGCAATGTGGCACCGTTTGCTGCATGGTATTGAAAAAAATATACCCAAATAAAAAAAGCATTagta                                                                                      >  2:391447/1‑248 (MQ=255)
                                                                                                                                                   cTATCAATTAGAAACATTAATACAACAACCGACGAAAAGTGACACAACGGCAGACCAACATCAACTGCAAGCTTTACGCGAACGAGTCATGGCATTGCTGACGACTCTGGAAGTGGCAGATGACATAAAA‑CTGGTCGACTGGTTACAGCAACGCCTGGGGCTTTTAGAGCAACGAGAAACGGCAATGTTGCACCGTTTGCTGCCTGATATTGAAAAAACTATAACCAAATAAAAAAAGCATTAGTAAGTAt                                                                                 >  2:462414/1‑251 (MQ=255)
                                                                                                                                                    tATCAATTAGCAACATTAATCCAACAACCGACGCAAAGTGACACAACGGCAGACCAACATCAACTGCAAGCTTTACGCGAACGAGTCATGGCATTGCTGACTACTCTGGCAGTGGCAGCGGACATAACA‑CTGGTCGACTGGTTACAGCAACGCCTGGGGCTTTTAGAGCAACGAGACACGGCAATGTTGCACCGTTTGCTGCATGCTATTGACAACAATATCACCAAATAAAAACTGCCTTAGTAAGTAtt                                                                                >  1:263225/1‑251 (MQ=255)
                                                                                                                                                       aaaTTAGCAACATTAATACAACAACCGACGAAAAGTGACACAACGGCAGACCAACATCAACTGCAAGCTTTACGCGAACGAGTCATGGCATTGCTGACGACTCTGGCAGTGGCAGATGACATAAAA‑CTGGTCGACTGGTTACAGCAACGCCTGGGGCTTTTAGAGCAACGAGACACGGCAATGTTGCACCGTTTGCTGCATGATATTGAAAAAAATATCACCAAATAAAAAACGCCTTAGTAAGTAttttt                                                                             <  2:149233/250‑1 (MQ=255)
                                                                                                                                                        aaTTAGCAACATTCATACAACAACCGGCGAAAAGTGAAACAACGGCAGACCAACATCAACTGCAATCTTTACGCGAACGAGTCATGGCAGTGCTGACGACTCTGGCAGTGGCAGATGACATAAAA‑CTGGTCGACGGTTTACAGCAACGTCTGTGGCTTTTAGAGCAACGAGACACGGCAATGTTGCACCGTTTGCTGCATGATATTGAAAAAAATATCCCCAAATAAAAAACGCCTTAGTAAGCATTTTTc                                                                            <  2:433125/251‑1 (MQ=255)
                                                                                                                                                         aTTAGCAACATTAATACAACAACCGACGAAAAGTGACACAACGGCAGACCAACATCAACTGCAAGCTTTACGCGAACGAGTCATGGCATTGCTGACGACTCTGGCAGTGGCAGATGACATAAAA‑CTGGTCGACTGGTTACAGCAACGCCTGGGGCTTTTAGAGCAACGAGACACGGCAATGTTGCACCGTTTGCTGCATGATATTGAAAAAAATATCACCAAATAAAAAACGCCTTAGTAAGTATTTTTCa                                                                           >  1:240347‑M2/1‑250 (MQ=255)
                                                                                                                                                            aGCTACAGTAATACAACAACCGAAGAAAAGTGACACAACGGCAGACCAACATCAACTGCAAGCTTTACGCGAACGAGTCATGGCATTGCTTACGACTCTGGCAGTGGCAGATGACATAAAA‑CTGGTCGACTGGTTACAGCAACGCCTGGGGCGTTTAGAGCAACGAGACACGGCAATGTTGCACCGTTTGCTGCATGATATTGAAAAAAATATCACCAAATAAAAAACGCCTTAGTAAGTATTTTTCagc                                                                         >  1:468481‑M2/1‑247 (MQ=255)
                                                                                                                                                              cAACATTAATACAACAACCGACGAAAAGTGACACAACGGCAGACCAACATCCACTGCAAGATTTACGCGAACGAGTCATGGCATTGCTGACGACTATGGCAGTCGCAGATGACATAAAA‑CTGGTCGACTGGTTACAGCAACGCCTGGGGCTTTTAGAGCAACGCGACCCGGCAATGTTGCACCGTTTGCTGCATGATATTGAAAAAAATATCACCAAATAAAAAACGCCTTAGTAAGTATTTTTCagcttt                                                                      <  1:498161‑M2/251‑7 (MQ=255)
                                                                                                                                                                 aatTAATACAACAACCGACGAAAAGTGACACAACGGCAGACCAACATCAACTGCAAGCTGTACGCGAACGAGTCATGGCATTGCTGACGACTCTGGCAGTGGCAGATGACATAAAA‑CTGGTCGACTGGTTACAGCAACGCCTGGGGCTTTTAGAGCAACGAGACACGGCAATGTTGCACCGTTTGCTGCATGATATTGAAAAAAATATCACCAAATAAAAAACGCCTTAGTAAGTATTTTTCagcttttca                                                                   >  2:146385‑M2/2‑242 (MQ=255)
                                                                                                                                                                     aaTACAACAACCGACGAAAAGTGACACAACGGCAGACCAACATCAACTGCAAGCTTTACGCGAACGAGTCATGGCATTGCTGACGACTCTGGCAGTGGCAGATGACATAAAA‑CTGGTCGACTGGTTACAGCAACGCCTGGGGCTTTTAGAGCAACGAGACCCGGCAATGTTGCACCGtt                                                                                                                                       >  2:366734/1‑179 (MQ=255)
                                                                                                                                                                     aaTACAACAACCGACGAAAAGTGACACAACGGCAGACCAACATCAACTGCAAGCTTTACGCGAACGAGTCATGGCATTGCTGACGACTCTGGCAGTGGCAGATGACATAAAA‑CTGGTCGACTGGTTACAGCAACGCCTGGGGCTTTTAGAGCAACGAGACACGGCAATGTTGCACCGtt                                                                                                                                       <  1:366734/179‑1 (MQ=255)
                                                                                                                                                                          aacaacCGACGAAAAGGGAGCCAACGGCAGACCAACATCAACTGCACGCTTTACGCGAACGATTCATGGCATTGCTGACGACGCTGGCAGTGGCAGATGACATAAAT‑CTGGTCGACTGGTTATAGCAACGCCTGGGGCTTTTAGAGCAACGAGACACGGCAATGTTGCACCGTTTGCTGCATGATATTGAAAAATTTATCACCAAATAAAAAACGCCTGAGTAATTATTTTTCagcttttcattctgactg                                                          >  2:502048‑M2/1‑233 (MQ=255)
                                                                                                                                                                                         gtgacacaacggcAGACCAACATCAACTGCAAGCTTTACGCGAACGAGTCATGGCATTGCTGACGACTCTGGCAGTGGCAGATGACATAAAA‑CTGGTCGACTGGTTACAGCAACGCCTGGGGCTTTTAGAGCAACGAGACACGGCAATGTTGCACCGTTTGCTGCATGATATTGAAAAAAATATCACCAAATAAAAAACGCCTTAGTAAGTATTTTTCagcttttcattctgactgcaacgggcaatatgt                                           <  1:51908‑M2/245‑34 (MQ=255)
                                                                                                                                                                                         gtgacacaacggcAGACCAACATCAACTGCAAGCTTTACGCGAACGAGTCATGGCATTGCTGACGACTCTGGCAGTGGCAGATGACATAAAA‑CTGGTCGACTGGTTACAGCAACGCCTGGGGCTTTTAGAGCAACGAGACACGGCAATGTTGCACCGTTTGCTGCATGATATTGAAAAAAATATCACCAAATAAAAAACGCCTTAGTAAGTATTTTTCagcttttcattctg                                                              <  1:472914‑M2/226‑15 (MQ=255)
                                                                                                                                                                                         gtgacacaacggcAGACCAACAGCAACTGCAAGCTTTACGCGAACGAGTCATGGCATTGCTGACGACTCTGGCAGTGGCAGATGACATAAAA‑CTGGTCGACTGGTTACAGCAACGCCTGGGGCTTTTAGAGCAACGAGACACGGCAATGTTGCACCGTTTGCTGCATGATATTGAAAAAAATAACACCAAATAAAAAACGCCTTAGTAAGTATTTTTCagcttttcattctg                                                              >  2:472914‑M2/7‑218 (MQ=255)
                                                                                                                                                                                                    gCAGACCAACATCAACTGCAAGCTTTACGCGAACGAGTCATGGCATTGCTGACGACTCTGGCAGTGGCAGATGACATAAAA‑CTGGTCGACTGGTTACAGCAACGCCTGGGGCTTTTAGAGCAACGAGACACGGCAATGTTGCACCGTTTGCTGCATGATATTGAAAAAAATATCACCAAATAAAAAACGCCTTAGAAAGTATTTTTCagcttttcattctgactgcaacgggcaatatgtctctgtgtgg                                 <  2:330122‑M2/250‑44 (MQ=255)
                                                                                                                                                                                                     cAGACCAACATCAACTGCAAGCTTTACGCGAACGAGTCATGGCATTGCTGACGACTCTGGCAGTGGCAGATGACATAAAA‑CTGGTCGACTGGTTACAGCAACGCCTGGGGCTTTTAGAGCAACGAGACACGGCAATGTTGCACCGTTTGCTGCATGATATTGAAAAAAATATCACCAAATAAAAAACGCCTTAGTAAGTATTTTTCagcttttcattctgactgcaacgggcaatatgtttctgtgtggat                               >  1:276830‑M2/1‑206 (MQ=255)
                                                                                                                                                                                                      aGACCAACATCAACTGCAAGCTTTACGCGAACGAGTCATGGCATTGCTGACGACTCTGGCAGTGGCAGATGACATAAAA‑CTGGTCGACTGGTTACAGCAACGCCTGGGGCTTTTAGAGCAACGAGACACGGCAATGTTGCACCGTTTGCTGCATGATATTGAAAAAAATATCACCAAATAAAAAACGCCTTAGTAAGTATTTTTCagcttttcattctgactgcaacgggcaatatgtctctgtgtggatt                              >  1:330122‑M2/1‑205 (MQ=255)
                                                                                                                                                                                                      aGACCAACATCAACTGCAAGCTTTACGCGAACGAGTCATGGCATTGCTGACGACTCTGGCAGTGGCAGATGACATAAAA‑CTGGTCGACTGGTTACAGCAACGCCTGGGGCTTTTAGAGCAACGAGACACGGCAATGTTGCACCGTTTGCTGCATGATATTGAAAAAAATATCACCAAATAAAAAACGCCTTAGTAAGTATTTTTCagcttttcattctgactgcaacgggcaatatgtctctgtgtggatt                              >  1:113151‑M2/1‑205 (MQ=255)
                                                                                                                                                                                                                cAACTGCAAGCTTTACGCGAACGAGTCATGGCATTGCTGACGACTCTGGCAGTGGCAGATGACATAAAA‑CTGGTCGACTGGTTACAGCAACGCCTGGGGCTTTTAGAGCAACGAGACACGGCAATGTTGCACCGTTTGCTGCATGATATTGAAAAAAATATCACCAAATAAAAAACGCCTTAGTAAGTATTTTTCagcttttcattctgactgcaacgggcaatatgtctctgtgtggattaaaaaaaga                     >  2:337024‑M2/1‑195 (MQ=255)
                                                                                                                                                                                                                    tGCAAGCTTTACGCGAACGAGTCATGGCATTGCTGACGACTCTGGCAGTGGCAGATGACATAAAA‑CTGGTCGACTGGTTACAGCAACGCCTGGGGCTTTTAGAGCAACGAGACACGGCAATGTTGCACCGTTTGCTGCATGATATTGAAAAAAATATCACCAAATAAAAAACGCCTTAGTAAGTATTTTTCagcttttcattctgactgcaacgggcaatatgtctctgtgtggattaaaaaaagagtgtc                <  2:266953‑M2/251‑61 (MQ=255)
                                                                                                                                                                                                                                  gAACGAGTCATGGCATTGCTGACGACTCTGGCAGTGGCAGATGACATAAAA‑CTGGTCGACTGGTTACAGCAACGCCTGGGGCTTTTAGAGCAACGAGACACGGCAATGTTGCACCGTTTGCTGCATGATATTGAAAAAAATATCACCAAATAAAAAACGCCTTAGTAAGTATTTTTCagcttttcattctgactgcaacgggcaatatgtctctgtgtggattaaaaaaagagtgtctgatagcagcttct  <  1:492670‑M2/251‑75 (MQ=255)
                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                               
ATGGATGAGCCATGCCGCGCTGGTGTTTGGTCGCGAAGATTCCGGGTTGACTAACGAAGAGTTAGCGTTGGCTGACGTTCTTACTGGTGTGCCGATGGTGGCGGATTATCCTTCGCTCAATCTGGGGCAGGCGGTGATGGTCTATTGCTATCAATTAGCAACATTAATACAACAACCGGCGAAAAGTGATGCAACGGCAGACCAACATCAACTGCAAGCTTTACGCGAACGAGCCATGACATTGCTGACGACTCTGGCAGTGGCAGATGACATAAAA‑CTGGTCGACTGGTTACAACAACGCCTGGGGCTTTTAGAGCAACGAGACACGGCAATGTTGCACCGTTTGCTGCATGATATTGAAAAAAATATCACCAAATAAAAAACGCCTTAGTAAGTATTTTTC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_000913/4641250‑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: ‑