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)
              ttGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTGTCtata                                                                                  <  2:119036/100‑1 (MQ=255)
                gtgatcGGAATTTGCTTATTTGAAGCGGTGCGTCAGCTCAGCTAATTTTTCATAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTGTCTATAgn                                                                                <  2:2703606/95‑2 (MQ=255)
                 cGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTGTCTATAGCt                                                                               >  1:876127/1‑100 (MQ=255)
                  gACCGGAATTTGCTTATTTGAAGCGGTGCGTCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTGTCTATAGCTg                                                                              <  1:1105644/100‑1 (MQ=255)
                  gACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTATCAGAATTATGAAAAGCCATCAAGATTTGGATGG‑TTTTTTTTGTTTATAgtt                                                                               >  2:3307969/1‑97 (MQ=255)
                   aCCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGtttttttttttttttgtt                                                                                >  2:1265120/1‑89 (MQ=255)
                     cgaaaTTTGATTATTTGAAGCGGTGCGTCAGCGCAGTTAATTGTTCAGAATTATGAAAAGCAATCCAGATTTGGATGGTTTTTTTTTGTCTATAGCTGGt                                                                            <  1:2292479/97‑1 (MQ=255)
                      ggAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTGTCTATAGCTGGTaa                                                                          <  1:2891778/100‑1 (MQ=255)
                       gAATTTGCTTATTTGAAGCTGTGCGCCAGCGCAGTTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTGTCTATAGCTg                                                                              <  1:255990/95‑1 (MQ=255)
                         aTTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTGTCTATAGCTGGTAAGan                                                                       <  2:387733/100‑2 (MQ=255)
                          tttGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTGTCTATAGCTGGTAAGATa                                                                      <  2:1781187/100‑1 (MQ=255)
                          tttGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTGTCTATAGCTGGTAAGATa                                                                      <  1:33341/100‑1 (MQ=255)
                             gCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTGTCTATAGCTGGTAAGATAAtt                                                                   >  1:1260461/1‑100 (MQ=255)
                              cTTATTTGAAGCGGTGCGCTAACGCATCTAATATCTCAGAATTATGAAAAGCCATACAGATTTGTATGGTTGTTATTT‑TTTATAGTTGTTAAGATaagt                                                                   >  2:3789338/1‑97 (MQ=255)
                              cTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATTG‑TTTTTTTTGTATATAGCTGGTAAGATAATTa                                                                  >  1:1891112/1‑100 (MQ=255)
                              cTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTGTCTATAGCTGGTAAGATAATTa                                                                  >  1:2009710/1‑100 (MQ=255)
                              cTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGACATCCAGAGTTGGATGGTTTTTTTTT‑TCTATAGTTGGTAAGATAttt                                                                   >  2:260761/1‑99 (MQ=255)
                              cTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAACCATCCAGATATGGATGG‑ATTTTTTTGTCTATAGCTGGTAAGATAAtt                                                                   >  1:36351/1‑99 (MQ=255)
                               ttATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTTGTCTATAGCTGGTAAGATAATTa                                                                  >  1:3416967/1‑100 (MQ=255)
                               ttATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTGTCTATAGCTGGTAAGATAATTAc                                                                 <  2:754617/100‑1 (MQ=255)
                                    tGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTGTCTATAGCTGGTAAGATAATTACGTAtt                                                            >  2:2603528/1‑100 (MQ=255)
                                    tGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTGTCTATAGCTGGTAAGATAATTACGTAtt                                                            >  1:1202205/1‑100 (MQ=255)
                                      aaGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGc                                                          >  2:186843/1‑100 (MQ=255)
                                         agctgcgctggcgcAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCaaa                                                       >  2:2553312/11‑100 (MQ=255)
                                          ggTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAat                                                      <  2:3625795/100‑1 (MQ=255)
                                             gcgcCAGCGCAGCTAATTTCTCAGTATTATGAATAGCCATCCATATTTTGATGG‑TTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATAtt                                                   <  1:68207/100‑1 (MQ=255)
                                             gcgcCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAata                                                     >  1:2580599/1‑98 (MQ=255)
                                              cgcCAGCGCAGCTAATTTCTCATAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTc                                                  <  1:699506/100‑1 (MQ=255)
                                              cgcCAGCGCAGCTAATTTACAAAAAAAATAAAAAGCCAACAACATTTGGATCG‑TTTTTTTTGTCTATACCTTGTAACCTAATTACATATTGCAAat                                                      >  2:3157091/1‑96 (MQ=255)
                                               gcCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCa                                                 >  1:1414663/1‑100 (MQ=255)
                                               gcCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCa                                                 <  2:1994718/100‑1 (MQ=255)
                                                ccAGCGCAGCTAATTTCTTAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCAt                                                <  1:406361/100‑1 (MQ=255)
                                                 cAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCAtt                                               >  1:2518813/1‑100 (MQ=255)
                                                  aGCGCAGCTAATTTTTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCAtt                                               <  2:600977/99‑1 (MQ=255)
                                                  aGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCAttt                                              >  1:2310409/1‑100 (MQ=255)
                                                   gcgcAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTGTCTATAGATGGTAAGATAATTACGTATTGCAAATATTCAtttt                                             >  2:3848903/1‑100 (MQ=255)
                                                     gcAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTTGATGG‑TTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTttat                                           <  2:2774286/100‑1 (MQ=255)
                                                      cAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTttat                                           >  2:1172196/1‑99 (MQ=255)
                                                      cAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTttagt                                          <  2:1462053/100‑3 (MQ=255)
                                                       aGCTAATTACTCAGAATTATGAAAAGACATCCAGATTTGGATGG‑TTTTTTTTGTCTATAGATGGTAAGATAATTACGTATTGCAAATATTCATTttatta                                         >  2:257365/1‑100 (MQ=255)
                                                        gttAATTTCTTAGAATTATGAAAAGTTATCTTGATTTGGATGG‑TTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTaa                                        <  1:2349212/98‑1 (MQ=255)
                                                        gCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTaa                                        >  1:3816006/1‑100 (MQ=255)
                                                          tAATTTTTCAGAATTATGTAAAGCCATCCAGATTTGGATGG‑TTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCAttt                                              <  1:825322/92‑1 (MQ=255)
                                                          tAATTTCTCAGTATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAata                                      <  1:1483337/100‑1 (MQ=255)
                                                          tAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAata                                      <  1:3314283/100‑1 (MQ=255)
                                                           aaTTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAatat                                     >  2:897505/1‑100 (MQ=255)
                                                            aTTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAatact                                    >  2:2372614/1‑98 (MQ=255)
                                                              ttCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTa                                  >  1:3768926/1‑100 (MQ=255)
                                                                atctgAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTaaa                                <  2:2142706/96‑1 (MQ=255)
                                                                  cAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACt                              >  1:2553185/1‑100 (MQ=255)
                                                                  cAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACt                              >  2:1611964/1‑100 (MQ=255)
                                                                   aGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACtt                             >  2:347241/1‑100 (MQ=255)
                                                                   aGAATTATGAAAAGCCATCCAAATTTGAATGG‑ATTTTTTTGTCTATAGCTGATAAGATAATTAAGTATTGCAAATATTCATTTTATTAATATTTAAAcat                             >  2:347619/1‑98 (MQ=255)
                                                                    tagacaaaaaaaACCATCCAAATCTGGATGG‑TTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTa                            >  2:1018256/9‑100 (MQ=37)
                                                                    gAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTaaaagat                            >  1:2142725/1‑96 (MQ=255)
                                                                    gAATTATGAAAAGCCATCCAGATTTAGATGG‑TTTTATTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACt                              >  2:471162/1‑98 (MQ=255)
                                                                     agacaaaaaaaACCATCCAAATCTGGATGG‑CTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTaa                           >  2:2922617/8‑100 (MQ=25)
                                                                      gacaaaaaaaACCATCCAAATCTGGATGG‑CTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTaa                           >  1:421108/7‑99 (MQ=25)
                                                                      gacaaaaaaaACCATCCAAATCTGGATGG‑CTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAAt                          >  1:1754783/7‑100 (MQ=25)
                                                                      gacaaaaaaaACCATCCAAATCTGCATGG‑TTTTTTTTGTATATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAAt                          >  2:3020934/7‑100 (MQ=255)
                                                                       acaaaaaaaACCATCCAAATCTGGATGG‑TTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATc                         >  2:3079807/6‑100 (MQ=37)
                                                                        tataAAAGGCCATCCAGATTTGGATGG‑TTTTTTTTGTCTATAGCTGGTAAAATAATTACGTATTGCAAATATTCTTTTTATTAATATTTAAACTTAATCt                        <  2:2456495/100‑1 (MQ=255)
                                                                        caaaaaaaACTATTCAGTTTTGGATGG‑TTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATc                         <  1:1778510/95‑1 (MQ=255)
                                                                        caaaaaaaACCATCCAAATCTGAAAGGCTTTTTTTTGTCTATAGCTGGTAAAATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAAt                          >  2:1110475/5‑99 (MQ=37)
                                                                           gAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCtata                     >  1:3899425/1‑100 (MQ=255)
                                                                           gAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCtata                     >  2:3397321/1‑100 (MQ=255)
                                                                           aaaaaACCATCCAAATCTGTATGG‑TTTTTTTTGTCTATAGTTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCtat                      >  2:3365937/2‑99 (MQ=255)
                                                                             aaaGCCATCCAGATTTGGATGG‑TTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATg                   <  2:1260447/100‑1 (MQ=255)
                                                                             aaaGCCATCCAGATTTGGATGG‑TTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATg                   <  2:2036148/100‑1 (MQ=255)
                                                                                gtcATTCAGATTTGGATTG‑TTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATg                   <  1:3553941/95‑1 (MQ=255)
                                                                                  cATCCAGTTTTGGATGG‑TTTTTTTTGTCTATAGCTGGTATGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGaat               <  1:1044730/99‑1 (MQ=255)
                                                                                   aTCCAGATTTGGATGG‑TTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGaataa             >  2:281016/1‑100 (MQ=255)
                                                                                   aTCCAGATTTGGATGG‑TTTTTTTTGTCTATAGCTGGTAAGAAAATTACTCATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGaataa             >  2:1338000/1‑100 (MQ=255)
                                                                                     ccAGATTTGGATGG‑TTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGAataata           >  1:283922/1‑100 (MQ=255)
                                                                                     ccAGATTTGGATGG‑TTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGAataata           <  1:462327/100‑1 (MQ=255)
                                                                                      cAGATTTGGATGGTTTTTTTTTGTCTATAGCTTGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGAataata           <  1:503880/100‑1 (MQ=255)
                                                                                       aGATTTGGATGG‑TTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGAGTAatata         <  1:2054309/100‑1 (MQ=255)
                                                                                        gATTTGGATGG‑TTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGAATAatatat        >  1:594680/1‑100 (MQ=255)
                                                                                          tttGGATGG‑TTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGAATAATATATGa      <  1:1894921/100‑1 (MQ=255)
                                                                                             ggATGGTTTTTTTTT‑TCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGAATAATATATGACAt   >  2:2383582/1‑100 (MQ=255)
                                                                                             ggATGGTTTTTTTTT‑TCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGAATAATATATGACAt   >  2:424427/1‑100 (MQ=255)
                                                                                             ggATGG‑TTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGAATAATATATGACAt   <  1:3222995/100‑1 (MQ=255)
                                                                                             ggATGG‑TTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGAATAATATATGACAt   >  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: ‑