Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,410,053 (T)9→8 100% intergenic (+47/‑80) rlmB → / → yjfI 23S rRNA mG2251 2'‑O‑ribose methyltransferase, SAM‑dependent/DUF2170 family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,410,0450T.89.2% 522.6 / 58.6 148intergenic (+39/‑88)rlmB/yjfI23S rRNA mG2251 2'‑O‑ribose methyltransferase, SAM‑dependent/DUF2170 family protein
Reads supporting (aligned to +/- strand):  ref base T (7/8);  new base . (73/59);  total (81/67)
Fisher's exact test for biased strand distribution p-value = 5.91e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 183 total aligned reads displayed.

CCTGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGAATAATATATGACATG  >  NC_000913/4409946‑4410140
                                                                                                   |                                                                                               
ccTGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGt                                                                                                 >  1:2722110/1‑100 (MQ=255)
  tGAACGTTTCTGTTGTGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAACCATCCAAATCTGGATGGctt                                                                                               <  1:947580/100‑1 (MQ=38)
  tGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGttt                                                                                               <  1:2595306/100‑1 (MQ=255)
   gAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGtttt                                                                                              >  2:2292458/1‑100 (MQ=255)
    aaCGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGTCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAAATCTGGATGGctttt                                                                                             <  1:1082301/100‑1 (MQ=255)
     aCGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGtttttt                                                                                            >  1:239912/1‑100 (MQ=255)
     aCGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTACGAAAAGCCATCCAGATTTGGATGGtttttt                                                                                            >  1:2922308/1‑100 (MQ=255)
      cGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGttttttt                                                                                           <  1:528562/100‑1 (MQ=255)
       gTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAACCATCCAAATCTGGATGGcttttcat                                                                                          <  1:3697208/100‑4 (MQ=37)
        tttCGGTTGCGACCGGAATTTGCTTATTTGATGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGttttttttg                                                                                         <  2:862849/100‑2 (MQ=255)
        tttCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGttttttttg                                                                                         <  2:1566699/100‑2 (MQ=255)
         ttCGGTTGCGACCGGAATTTGCTTATTTTAAGCGGTGCGCCAGCTCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGtttttttttn                                                                                        <  1:3279562/100‑2 (MQ=255)
          tCGTTTGCGACTGTAATTTGCTTATTTGATGCGGTGCGTCAGCGTAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGttttttttgtc                                                                                       <  1:960604/100‑4 (MQ=255)
          tCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGttttttttgtc                                                                                       <  1:716875/100‑4 (MQ=255)
          tCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGttttttttgtc                                                                                       <  1:2664347/100‑4 (MQ=255)
          tCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCCCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGttttttttgtc                                                                                       <  2:2509549/100‑4 (MQ=255)
            ggTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAACCATCCAAATCTGGATGGcttttcataattc                                                                                     <  1:1392079/100‑9 (MQ=25)
              ttGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCtata                                                                                  <  2:119036/100‑1 (MQ=255)
                gtgatcGGAATTTGCTTATTTGAAGCGGTGCGTCAGCTCAGCTAATTTTTCATAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAgn                                                                                <  2:2703606/95‑2 (MQ=255)
                 cGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCt                                                                               >  1:876127/1‑100 (MQ=255)
                  gACCGGAATTTGCTTATTTGAAGCGGTGCGTCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTg                                                                              <  1:1105644/100‑1 (MQ=255)
                  gACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTATCAGAATTATGAAAAGCCATCAAGATTTGGATGGTTTTTTTTGTTTATAgtt                                                                               >  2:3307969/1‑97 (MQ=255)
                   aCCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGtttttttttttttttgtt                                                                                >  2:1265120/1‑89 (MQ=255)
                     cgaaaTTTGATTATTTGAAGCGGTGCGTCAGCGCAGTTAATTGTTCAGAATTATGAAAAGCAATCCAGATTTGGATGGTTTTTTTTTGTCTATAGCTGGt                                                                            <  1:2292479/97‑1 (MQ=255)
                      ggAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTaa                                                                          <  1:2891778/100‑1 (MQ=255)
                       gAATTTGCTTATTTGAAGCTGTGCGCCAGCGCAGTTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTg                                                                              <  1:255990/95‑1 (MQ=255)
                         aTTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGan                                                                       <  2:387733/100‑2 (MQ=255)
                          tttGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATa                                                                      <  2:1781187/100‑1 (MQ=255)
                          tttGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATa                                                                      <  1:33341/100‑1 (MQ=255)
                             gCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAAtt                                                                   >  1:1260461/1‑100 (MQ=255)
                              cTTATTTGAAGCGGTGCGCTAACGCATCTAATATCTCAGAATTATGAAAAGCCATACAGATTTGTATGGTTGTTATTTTTTATAGTTGTTAAGATaagt                                                                   >  2:3789338/1‑97 (MQ=255)
                              cTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATTGTTTTTTTTGTATATAGCTGGTAAGATAATTa                                                                  >  1:1891112/1‑100 (MQ=255)
                              cTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTa                                                                  >  1:2009710/1‑100 (MQ=255)
                              cTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGACATCCAGAGTTGGATGGTTTTTTTTTTCTATAGTTGGTAAGATAttt                                                                   >  2:260761/1‑99 (MQ=255)
                              cTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAACCATCCAGATATGGATGGATTTTTTTGTCTATAGCTGGTAAGATAAtt                                                                   >  1:36351/1‑99 (MQ=255)
                               ttATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTTGTCTATAGCTGGTAAGATAATTa                                                                  >  1:3416967/1‑100 (MQ=255)
                               ttATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTAc                                                                 <  2:754617/100‑1 (MQ=255)
                                    tGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTAtt                                                            >  2:2603528/1‑100 (MQ=255)
                                    tGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTAtt                                                            >  1:1202205/1‑100 (MQ=255)
                                      aaGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGc                                                          >  2:186843/1‑100 (MQ=255)
                                         agctgcgctggcgcAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCaaa                                                       >  2:2553312/11‑100 (MQ=255)
                                          ggTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAat                                                      <  2:3625795/100‑1 (MQ=255)
                                             gcgcCAGCGCAGCTAATTTCTCAGTATTATGAATAGCCATCCATATTTTGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATAtt                                                   <  1:68207/100‑1 (MQ=255)
                                             gcgcCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAata                                                     >  1:2580599/1‑98 (MQ=255)
                                              cgcCAGCGCAGCTAATTTCTCATAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTc                                                  <  1:699506/100‑1 (MQ=255)
                                              cgcCAGCGCAGCTAATTTACAAAAAAAATAAAAAGCCAACAACATTTGGATCGTTTTTTTTGTCTATACCTTGTAACCTAATTACATATTGCAAat                                                      >  2:3157091/1‑96 (MQ=255)
                                               gcCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCa                                                 >  1:1414663/1‑100 (MQ=255)
                                               gcCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCa                                                 <  2:1994718/100‑1 (MQ=255)
                                                ccAGCGCAGCTAATTTCTTAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCAt                                                <  1:406361/100‑1 (MQ=255)
                                                 cAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCAtt                                               >  1:2518813/1‑100 (MQ=255)
                                                  aGCGCAGCTAATTTTTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCAtt                                               <  2:600977/99‑1 (MQ=255)
                                                  aGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCAttt                                              >  1:2310409/1‑100 (MQ=255)
                                                   gcgcAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGATGGTAAGATAATTACGTATTGCAAATATTCAtttt                                             >  2:3848903/1‑100 (MQ=255)
                                                     gcAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTTGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTttat                                           <  2:2774286/100‑1 (MQ=255)
                                                      cAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTttat                                           >  2:1172196/1‑99 (MQ=255)
                                                      cAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTttagt                                          <  2:1462053/100‑3 (MQ=255)
                                                       aGCTAATTACTCAGAATTATGAAAAGACATCCAGATTTGGATGGTTTTTTTTGTCTATAGATGGTAAGATAATTACGTATTGCAAATATTCATTttatta                                         >  2:257365/1‑100 (MQ=255)
                                                        gttAATTTCTTAGAATTATGAAAAGTTATCTTGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTaa                                        <  1:2349212/98‑1 (MQ=255)
                                                        gCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTaa                                        >  1:3816006/1‑100 (MQ=255)
                                                          tAATTTTTCAGAATTATGTAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCAttt                                              <  1:825322/92‑1 (MQ=255)
                                                          tAATTTCTCAGTATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAata                                      <  1:1483337/100‑1 (MQ=255)
                                                          tAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAata                                      <  1:3314283/100‑1 (MQ=255)
                                                           aaTTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAatat                                     >  2:897505/1‑100 (MQ=255)
                                                            aTTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAatact                                    >  2:2372614/1‑98 (MQ=255)
                                                              ttCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTa                                  >  1:3768926/1‑100 (MQ=255)
                                                                atctgAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTaaa                                <  2:2142706/96‑1 (MQ=255)
                                                                  cAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACt                              >  1:2553185/1‑100 (MQ=255)
                                                                  cAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACt                              >  2:1611964/1‑100 (MQ=255)
                                                                   aGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACtt                             >  2:347241/1‑100 (MQ=255)
                                                                   aGAATTATGAAAAGCCATCCAAATTTGAATGGATTTTTTTGTCTATAGCTGATAAGATAATTAAGTATTGCAAATATTCATTTTATTAATATTTAAAcat                             >  2:347619/1‑98 (MQ=255)
                                                                    tagacaaaaaaaACCATCCAAATCTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTa                            >  2:1018256/9‑100 (MQ=37)
                                                                    gAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTaaaagat                            >  1:2142725/1‑96 (MQ=255)
                                                                    gAATTATGAAAAGCCATCCAGATTTAGATGGTTTTATTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACt                              >  2:471162/1‑98 (MQ=255)
                                                                     agacaaaaaaaACCATCCAAATCTGGATGGCTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTaa                           >  2:2922617/8‑100 (MQ=25)
                                                                      gacaaaaaaaACCATCCAAATCTGGATGGCTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTaa                           >  1:421108/7‑99 (MQ=25)
                                                                      gacaaaaaaaACCATCCAAATCTGGATGGCTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAAt                          >  1:1754783/7‑100 (MQ=25)
                                                                      gacaaaaaaaACCATCCAAATCTGCATGGTTTTTTTTGTATATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAAt                          >  2:3020934/7‑100 (MQ=255)
                                                                       acaaaaaaaACCATCCAAATCTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATc                         >  2:3079807/6‑100 (MQ=37)
                                                                        tataAAAGGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAAATAATTACGTATTGCAAATATTCTTTTTATTAATATTTAAACTTAATCt                        <  2:2456495/100‑1 (MQ=255)
                                                                        caaaaaaaACTATTCAGTTTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATc                         <  1:1778510/95‑1 (MQ=255)
                                                                        caaaaaaaACCATCCAAATCTGAAAGGCTTTTTTTTGTCTATAGCTGGTAAAATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAAt                          >  2:1110475/5‑99 (MQ=37)
                                                                           gAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCtata                     >  1:3899425/1‑100 (MQ=255)
                                                                           gAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCtata                     >  2:3397321/1‑100 (MQ=255)
                                                                           aaaaaACCATCCAAATCTGTATGGTTTTTTTTGTCTATAGTTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCtat                      >  2:3365937/2‑99 (MQ=255)
                                                                             aaaGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATg                   <  2:1260447/100‑1 (MQ=255)
                                                                             aaaGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATg                   <  2:2036148/100‑1 (MQ=255)
                                                                                gtcATTCAGATTTGGATTGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATg                   <  1:3553941/95‑1 (MQ=255)
                                                                                  cATCCAGTTTTGGATGGTTTTTTTTGTCTATAGCTGGTATGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGaat               <  1:1044730/99‑1 (MQ=255)
                                                                                   aTCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGaataa             >  2:281016/1‑100 (MQ=255)
                                                                                   aTCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGAAAATTACTCATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGaataa             >  2:1338000/1‑100 (MQ=255)
                                                                                     ccAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGAataata           >  1:283922/1‑100 (MQ=255)
                                                                                     ccAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGAataata           <  1:462327/100‑1 (MQ=255)
                                                                                      cAGATTTGGATGGTTTTTTTTTGTCTATAGCTTGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGAataata           <  1:503880/100‑1 (MQ=255)
                                                                                       aGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGAGTAatata         <  1:2054309/100‑1 (MQ=255)
                                                                                        gATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGAATAatatat        >  1:594680/1‑100 (MQ=255)
                                                                                          tttGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGAATAATATATGa      <  1:1894921/100‑1 (MQ=255)
                                                                                             ggATGGTTTTTTTTTTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGAATAATATATGACAt   >  2:2383582/1‑100 (MQ=255)
                                                                                             ggATGGTTTTTTTTTTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGAATAATATATGACAt   >  2:424427/1‑100 (MQ=255)
                                                                                             ggATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGAATAATATATGACAt   <  1:3222995/100‑1 (MQ=255)
                                                                                             ggATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGAATAATATATGACAt   >  1:1585356/1‑100 (MQ=255)
                                                                                                   |                                                                                               
CCTGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGAATAATATATGACATG  >  NC_000913/4409946‑4410140

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