Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,410,053 (T)9→8 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.96.2% 654.8 / 18.1 157intergenic (+39/‑88)rlmB/yjfI23S rRNA mG2251 2'‑O‑ribose methyltransferase, SAM‑dependent/DUF2170 family protein
Reads supporting (aligned to +/- strand):  ref base T (6/0);  new base . (73/77);  total (80/77)
Fisher's exact test for biased strand distribution p-value = 2.84e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.80e-01

Only 100 of 172 total aligned reads displayed.

CCTGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGAATAATATATGACATGGA  >  NC_000913/4409946‑4410142
                                                                                                   |                                                                                                 
ccTGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGt                                                                                                   >  1:820781/1‑100 (MQ=255)
  tGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGttt                                                                                                 <  2:2517066/100‑1 (MQ=255)
        tttCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGttttttttg                                                                                           <  2:1404643/100‑2 (MQ=255)
         ttCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGttttttttgt                                                                                          >  1:605462/1‑98 (MQ=255)
         ttCGGTTGCGACCAGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATATGAAAAGCCATCCAGATTTGGATGGttttttttgtc                                                                                         <  2:1758860/100‑4 (MQ=255)
          tCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGttttttttgtc                                                                                         <  2:557301/100‑4 (MQ=255)
          tCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGttttttttgtc                                                                                         <  2:1524687/100‑4 (MQ=255)
            ggTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCta                                                                                      <  1:2275765/100‑1 (MQ=255)
                gCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGc                                                                                  <  2:3070040/100‑1 (MQ=255)
                   aCCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTgg                                                                               <  1:1385190/100‑1 (MQ=255)
                     cGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTa                                                                             <  1:2668428/100‑1 (MQ=255)
                      ggAATTTGCTTATTTGAAGCGGTGCGACAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTaa                                                                            >  1:1878729/1‑100 (MQ=255)
                       gAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAg                                                                           >  2:1636287/1‑100 (MQ=255)
                       gAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAg                                                                           <  1:1520293/100‑1 (MQ=255)
                       gAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAg                                                                           >  1:511304/1‑100 (MQ=255)
                          tttGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATa                                                                        <  2:2285103/100‑1 (MQ=255)
                             gCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTa                                                                    >  2:2685956/1‑100 (MQ=255)
                              cttaatcaaaGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGCTTTGGATGGTTTTTTTTGTCTATACCTGGTAACATAAATa                                                                    >  2:1367312/9‑100 (MQ=255)
                              cTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTa                                                                    >  1:1593785/1‑100 (MQ=255)
                              cTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTa                                                                    <  2:131949/100‑1 (MQ=255)
                              cTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTa                                                                    >  2:1270788/1‑100 (MQ=255)
                                tATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACg                                                                  >  1:3056411/1‑100 (MQ=255)
                                 aTTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGt                                                                 <  2:1985398/100‑1 (MQ=255)
                                 aTTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGt                                                                 >  2:947252/1‑100 (MQ=255)
                                  tttGAAGCGGTGCTCCAGCGCAGCTAATTTCTCACAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTATATATCTGGTAAGATAATTACGTa                                                                >  2:2886638/1‑100 (MQ=255)
                                   ttGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTAt                                                               <  2:1514675/100‑1 (MQ=255)
                                     gAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTg                                                             >  1:2599618/1‑100 (MQ=255)
                                      aaGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGc                                                            <  2:1665527/100‑1 (MQ=255)
                                      aaGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGc                                                            <  1:1840125/100‑1 (MQ=255)
                                        gCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCaa                                                          <  2:2580211/100‑1 (MQ=255)
                                             gcgcCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATAtt                                                     >  1:2712761/1‑100 (MQ=255)
                                                 cAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCAtt                                                 <  1:2730783/100‑1 (MQ=255)
                                                  aGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCAttt                                                >  1:441480/1‑100 (MQ=255)
                                                      cAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTttatt                                            <  1:774847/100‑1 (MQ=255)
                                                      cAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTttatt                                            <  1:2548016/100‑1 (MQ=255)
                                                      cAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTAGTAAGATAATTACGTATTGCAAATATTCATTttatt                                            >  2:2686505/1‑100 (MQ=255)
                                                      cAGCTAATTTCTCAGAATTATGAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTttatta                                           <  1:2685981/100‑1 (MQ=255)
                                                        gCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTaa                                          >  1:1302798/1‑100 (MQ=255)
                                                          tAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAata                                        <  1:1562831/100‑1 (MQ=255)
                                                          tAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAat                                         <  2:2057095/99‑1 (MQ=255)
                                                           aaTTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTTTTTTTTTCTGGTAAGATAATTACGTATTGCAAATATTTATTTTATTAatat                                       >  2:210599/1‑100 (MQ=255)
                                                           aaTTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAatat                                       <  2:2068361/100‑1 (MQ=255)
                                                            aTTTCTCAGATTTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAatat                                       <  2:2616229/99‑1 (MQ=255)
                                                            aTTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATAtt                                      >  1:45083/1‑100 (MQ=255)
                                                               tctcAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTa                                    <  2:2621179/99‑1 (MQ=255)
                                                                ctcAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTaaa                                  <  1:1041071/100‑1 (MQ=255)
                                                                ctcAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTaaa                                  >  1:1991802/1‑100 (MQ=255)
                                                                ctcAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTaaa                                  >  2:2343721/1‑100 (MQ=255)
                                                                  cAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACt                                <  1:1661301/100‑1 (MQ=255)
                                                                  cAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACt                                >  2:2597367/1‑100 (MQ=255)
                                                                   aGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACtt                               >  1:38620/1‑100 (MQ=255)
                                                                      gacaaaaaaaACCATCCAAATCTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAAt                            >  2:1248133/7‑100 (MQ=37)
                                                                      aTTATGAAAAGTCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTa                              <  2:1578730/98‑1 (MQ=255)
                                                                      aTTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAAt                            >  1:559650/1‑100 (MQ=255)
                                                                      aTTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAACTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAAt                            <  2:1423969/100‑1 (MQ=255)
                                                                       ttATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATc                           <  2:699089/100‑1 (MQ=255)
                                                                       ttATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATc                           >  2:1447791/1‑100 (MQ=255)
                                                                        caaaaaaaACCATCCAAATCTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCt                          >  1:331165/5‑100 (MQ=37)
                                                                         aagAAAAGCCATCAAGATTTGGATGGTTTTTATTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCta                         >  1:2082053/3‑100 (MQ=255)
                                                                         aaaaaaaACCATCCAAATCTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTTCAAATATTTATTTTATTAATATTTAAACTTAATc                           >  2:653140/4‑98 (MQ=37)
                                                                         aTGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCta                         >  1:86719/1‑100 (MQ=255)
                                                                          tGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCtat                        >  1:2050646/1‑100 (MQ=255)
                                                                          tGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCtat                        >  2:1357434/1‑100 (MQ=255)
                                                                          aaaaaaACCATCCAAATCTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCtat                        >  2:1711823/3‑100 (MQ=38)
                                                                           gAAAAGCCATCGAGATTTGGATGGTTTTGTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCtata                       <  2:1065542/100‑1 (MQ=255)
                                                                           gAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCtata                       >  2:436570/1‑100 (MQ=255)
                                                                           gAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCtata                       >  2:103967/1‑100 (MQ=255)
                                                                           aaaaaACCATCCAAATCTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCtata                       >  1:501028/2‑100 (MQ=255)
                                                                           aaaaaACCATCCAAATCTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACATAATCtata                       >  2:1188147/2‑100 (MQ=38)
                                                                            aaaaACCATCCAAATCTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTatat                      >  2:1405128/1‑100 (MQ=255)
                                                                             aaaaCCATCCAAATCTGGATGGTTTTTTTTTTTTTTATCTGGTAAGATAATTACGTATTGCAAATATTTTTTTTTTTAATATTTAAATTTCATCtattt                      >  2:908623/1‑97 (MQ=255)
                                                                              aaaccATCCAAATCTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATgg                    >  1:2053431/4‑100 (MQ=255)
                                                                              aaaccATCCAAATCTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATgg                    >  1:2747198/4‑100 (MQ=255)
                                                                              aaaccATCCAAATCTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATgg                    >  1:1777146/4‑100 (MQ=255)
                                                                              aaGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATgg                    >  2:1890998/1‑100 (MQ=255)
                                                                              aaGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATgg                    <  1:2642660/100‑1 (MQ=255)
                                                                              aaGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATgg                    <  1:322830/100‑1 (MQ=255)
                                                                               aGCCATTCAGATTTGGATGGTTTTTTTTGTCTATAGCTTGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATg                     <  2:1533368/98‑1 (MQ=255)
                                                                                gtcATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGaa                  >  1:1191645/3‑100 (MQ=255)
                                                                                accATCCAAATCTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGaa                  >  2:271582/2‑100 (MQ=255)
                                                                                 ccATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGaat                 <  2:1593768/100‑1 (MQ=255)
                                                                                      cAGATTTGGATGGTTTTTTTTGTTTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGAATAatat            <  2:1144005/100‑1 (MQ=255)
                                                                                      cAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGAATAatat            <  2:824433/100‑1 (MQ=255)
                                                                                      cAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGAATAatat            <  2:42421/100‑1 (MQ=255)
                                                                                       aGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGAATAatata           >  1:1450044/1‑100 (MQ=255)
                                                                                        gATTTGGATGGTTTTTTTTGTCTATAGCTTGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGAATAatatat          <  2:484931/100‑1 (MQ=255)
                                                                                         aTTTGGATGGTTTATTTTGTCTATAGCTTTTAAAATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGAataata             <  2:605460/96‑1 (MQ=255)
                                                                                          ttttGATGGTTTTTTTTGTCTATAGCTTGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGAATAATATATGa        <  2:2860470/100‑1 (MQ=255)
                                                                                          tttGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGAATAATATATGa        >  2:324970/1‑100 (MQ=255)
                                                                                           ttGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGAATAATATATGAc       <  1:1268483/100‑1 (MQ=255)
                                                                                           ttGGATGGTTTTTTGTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGAATAATATATGAc       <  1:1357446/100‑1 (MQ=255)
                                                                                           ttGAATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTTCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGAATAATATATGAc       <  2:721248/100‑1 (MQ=255)
                                                                                            tGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGAATAATATATGAc       <  2:1926038/99‑1 (MQ=255)
                                                                                             ggATGGCTTTTTTTTTTTTTTTTCTGGTAAGATAATTACGTATTGTAAATTTTTTTTTTATTAATATTTAAATTTAATCTATATGGAATAatatattttt      >  2:1745137/1‑96 (MQ=255)
                                                                                             ggATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGAATAATATATGACAt     >  1:2848453/1‑100 (MQ=255)
                                                                                             ggATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGAATAATATATGACAt     <  1:2587192/100‑1 (MQ=255)
                                                                                             ggATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGAATAATATATGACAt     <  1:1479822/100‑1 (MQ=255)
                                                                                              gATGGTTTTTTTTTTTTTTTTCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGAATAATATATTACAt     >  2:2385515/1‑99 (MQ=255)
                                                                                              gATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGAATAATATATGACATg    <  1:148082/100‑1 (MQ=255)
                                                                                                tGGTTTTTTTTTTTTTTTGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGAATAATATATGACATTGa  >  2:224439/1‑100 (MQ=255)
                                                                                                   |                                                                                                 
CCTGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGAATAATATATGACATGGA  >  NC_000913/4409946‑4410142

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