Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,410,037 T→C 15.4% intergenic (+31/‑96) 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,0370TC15.4% 403.3 / 63.5 181intergenic (+31/‑96)rlmB/yjfI23S rRNA mG2251 2'‑O‑ribose methyltransferase, SAM‑dependent/DUF2170 family protein
Reads supporting (aligned to +/- strand):  ref base T (76/76);  new base C (10/18);  total (86/95)
Fisher's exact test for biased strand distribution p-value = 2.17e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.11e-02

Only 100 of 188 total aligned reads displayed.

TTTCTTCCCTGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TT‑TTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGAATAATATATGA  >  NC_000913/4409939‑4410136
                                                                                                  |                                                                                                     
tttCTTCCCTGAACGTTTCGGTTGCGACCGGAATTTTATTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGAtt                                                                                                       >  1:2357038/1‑99 (MQ=255)
tttCTTCCCTGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGAttt                                                                                                      <  2:1088173/100‑1 (MQ=255)
 ttcttcCCTGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTg                                                                                                     <  1:2870306/100‑1 (MQ=255)
 ttcttcCCTGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTg                                                                                                     <  2:2910756/100‑1 (MQ=255)
 ttcttcCCTGAACGTCTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTg                                                                                                     <  2:4335069/100‑1 (MQ=255)
  tcttcCCTGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTgg                                                                                                    <  1:903226/100‑1 (MQ=255)
   cttcCCTGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGa                                                                                                   <  1:108780/100‑1 (MQ=255)
    ttcCCTGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGAt                                                                                                  <  2:1617400/100‑1 (MQ=255)
       cttGAACGTTTCGATTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑t                                                                                              >  2:2990829/3‑100 (MQ=255)
       ccTGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑t                                                                                              >  1:2003809/1‑100 (MQ=255)
       ccTGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑t                                                                                              >  2:3208043/1‑100 (MQ=255)
        cTGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAACAGCCATCCAGATTTGGATGG‑tt                                                                                             >  2:1773254/1‑100 (MQ=255)
           aaCGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAAATCTGGATGG‑ct‑ttt                                                                                         <  2:2468866/100‑1 (MQ=255)
            acttttCGGTTGCGACCGGAATTTTCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAAATCTGGATGG‑ct‑tttc                                                                                        <  2:3188900/97‑2 (MQ=255)
              gTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑tt‑tttttt                                                                                      <  1:210235/100‑1 (MQ=255)
                ttCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TT‑TTTTTTTg                                                                                    <  1:76920/100‑1 (MQ=255)
                ttCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAACCATCCAAATCTGGATGG‑ct‑tttcat                                                                                      <  1:1175051/98‑4 (MQ=37)
                tccGGTGGCGACCTGAATTTGCTTACTTAAAGCGGTGCGCCAGCGCAGCTAATGTCTCAGAATTATGAAAAGCCATCCAAATCTGGATGG‑ct‑tttc                                                                                        <  1:3883929/94‑2 (MQ=37)
                 tCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAACCATCCAAATCTGGATGG‑ct‑tttcataat                                                                                   <  2:4066807/100‑7 (MQ=37)
                  cGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TT‑TTTTTTTGTc                                                                                  <  2:2829659/100‑1 (MQ=255)
                   ggTTGCTGCCGGAATTTGCTTATTGGAAGCGGTGCGCCAGCGCAGCTAATTTCTCTGAATTATGAAAAGCCATCCAGATTTGGATGG‑TT‑TTTTTTTGTCt                                                                                 <  1:3740468/100‑1 (MQ=255)
                   ggTTGCCACAGGAATTTGCTTATTGAAAGCGTTGCGCCAGCCCAGCCAATTTCTCAGAATTATGAAAAACCATCCAACTCTGGAGGG‑ct‑tttca                                                                                       <  1:327368/94‑3 (MQ=255)
                      tGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TT‑TTTTTTTGTCtata                                                                              >  2:1197168/1‑100 (MQ=255)
                        cGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TT‑TTTTTTTGTCTATAGc                                                                            >  2:1837039/1‑100 (MQ=255)
                          aCCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TT‑TTTTTTTGTCTATAGCTg                                                                          <  2:4250961/100‑1 (MQ=255)
                           ccGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TT‑TTTTTTTGTCTATAGCTgg                                                                         >  2:19576/1‑100 (MQ=255)
                           ccGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TT‑TTTTTTTGTCTATAGCTgg                                                                         >  2:2236388/1‑100 (MQ=255)
                             ggAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TT‑TTTTTTTGTCTATAGCTGGTa                                                                       <  1:3262269/100‑1 (MQ=255)
                                 tttGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TT‑TTTTTTTGTCTATAGCTGGTAAGAt                                                                   <  1:3379288/100‑1 (MQ=255)
                                   tGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCAAGATTTGGATGG‑AT‑TTTTTTTGTTTATTGCTTTTAAGTTaa                                                                 >  1:2607343/1‑100 (MQ=255)
                                    gCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TT‑TTTTTTTGTCTATAGCTGGTAAGATAAt                                                                >  2:2874697/1‑100 (MQ=255)
                                     cTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TT‑TTTTTTTGTCTATAGCTGGTAAGATAAtt                                                               >  2:1567241/1‑100 (MQ=255)
                                       tATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TT‑TTTTTTTGTCTATAGCTGGTAAGATAATTAc                                                             >  2:521193/1‑100 (MQ=255)
                                       tATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TT‑TTTTTTTGTCTATAGCTGGTAAGATAATTAc                                                             >  2:943707/1‑100 (MQ=255)
                                        aTTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TT‑TTTTTTTGTCTATAGCTGGTAAGATAATTACg                                                            >  2:1463177/1‑100 (MQ=255)
                                         tttGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TT‑TTTTTTTGTCTATAGCTGGTAAGATAATTACg                                                            <  2:619526/99‑1 (MQ=255)
                                         tttGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TT‑TTTTTTTGTCTATAGCTGGTAAGATAATTACGt                                                           <  2:3272653/100‑1 (MQ=255)
                                         tttGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TT‑TTTTTTTGTCTATAGCTGGTAAGATAATTACGt                                                           >  1:3225088/1‑100 (MQ=255)
                                          ttGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGCTTTGGATGG‑TT‑TTTTTTTGTCTATAGATGGTAAGATAATTACGTa                                                          >  1:3520256/1‑100 (MQ=255)
                                           tGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TT‑TTTTTTTGTCTATAGCTGGTAAGATAATTACGTAt                                                         >  2:3594132/1‑100 (MQ=255)
                                            ggaGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TT‑TTTTTTTGTCTATAGCTGGTAAGATAATTACGTAt                                                         <  2:581739/97‑1 (MQ=255)
                                            gAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TT‑TTTTTTTGTCTATAGCTGGTAAGATAATTACGTAtt                                                        >  2:1037555/1‑100 (MQ=255)
                                              aGCGGTGCTCCAGCGCAGCTAATGTCTCAGAATTATGAAAAGCCCTCCAGATGTGCATGG‑TT‑TTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGc                                                      <  1:342922/100‑1 (MQ=255)
                                              aGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TT‑TTTTTTGTTCAAAGGCTGGTAAGTTAATTActtttt                                                        >  2:3377874/1‑93 (MQ=255)
                                               gCGGGGCCCCGGCGCACATATTTTTTCATATTTATAAAAAGCCACCCCATTTTGGTTGT‑TT‑TTTTTTTGTCTATAGCTGGTAAGATAATTACGTAtt                                                        <  2:312763/97‑1 (MQ=255)
                                                  gTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TT‑TTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAat                                                  >  2:2583957/1‑100 (MQ=255)
                                                   tGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TT‑TTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAata                                                 <  2:2122878/100‑1 (MQ=255)
                                                   tGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TT‑TTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAata                                                 <  1:3054162/100‑1 (MQ=255)
                                                    gcgcCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TT‑TTTTTTTTTCTTTTTTTGATAAGATAATTACGTATTGCAAatat                                                >  1:1572820/1‑100 (MQ=255)
                                                     cgcCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TT‑TTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATAtt                                               >  2:3485655/1‑100 (MQ=255)
                                                      gcCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TT‑TTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTc                                              >  2:3026518/1‑100 (MQ=255)
                                                      gcCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑‑T‑TTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCa                                             <  2:1921275/100‑1 (MQ=255)
                                                             cAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TT‑TTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTttat                                       >  2:1586504/1‑100 (MQ=255)
                                                             cAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TT‑TTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTttat                                       <  1:1042749/100‑1 (MQ=255)
                                                              aGCTAATTTCGCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TT‑TTTTTTTGTTCATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTttat                                       <  1:2055782/99‑1 (MQ=255)
                                                               gCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TT‑TTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTttatta                                     <  1:3026501/100‑1 (MQ=255)
                                                               gCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TT‑TTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTttatta                                     <  2:48684/100‑1 (MQ=255)
                                                                    tttCTCATAATTATGGAAAGACAGCAATATTTGGATTG‑TT‑TTTTTGTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAata                                  <  1:84203/98‑1 (MQ=255)
                                                                    tttCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TT‑TTTTTTTGTTTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATAtt                                >  1:677698/1‑100 (MQ=255)
                                                                    tttCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TT‑TTTTTTTGTATATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAata                                  >  2:4274488/1‑98 (MQ=255)
                                                                     tgttttgAATTATGTAAAGCCATTTTTTTTTTTATGG‑TTGTTTTTTTGTCTATACCTGGTAATATAATTACGTATTGCAAATATTCATTTTATTAata                                  <  1:1337252/92‑1 (MQ=255)
                                                                      tctcAGATTTTGGAACAGCCATCGAGATTTTGATTG‑TT‑TTTTTTTGTCTATAACTTCTACGATAATTACGTATTGCAAATATTCATTTTATTAATATTTa                              <  1:4157355/100‑1 (MQ=255)
                                                                       ctcAGAATTATGAAAAGCCATCTAGATTTGGATGG‑TT‑TTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTaa                             <  1:2519844/100‑1 (MQ=255)
                                                                       ctcAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TT‑TTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTaa                             >  2:1126133/1‑100 (MQ=255)
                                                                        tcttaaTTATGATAAGCCATCTATATTTTGATGG‑TT‑TTTTTTTGTATATAGCTTTTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTaaa                            <  1:2175525/96‑1 (MQ=255)
                                                                        tcAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TT‑TTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTACAAATATTCATTTTATTAATATTTaaa                            <  2:3506557/100‑1 (MQ=255)
                                                                         cAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TT‑TTTTTTTGTCTTTAGCTGGTATGATAATTACGTATTGCAAATATTCATTTTATTGATATTTaaa                            >  1:2813267/1‑99 (MQ=255)
                                                                         cAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TT‑TTTTTTTGTCTATAGCTGTTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAAc                           <  1:2449926/100‑1 (MQ=255)
                                                                          aGAATTATGAAAAGCCATCCAGATTTGGATGG‑TT‑TTTTTTTTTTTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACt                          >  1:2192982/1‑100 (MQ=255)
                                                                            gacaaaaaaaaACCATCCAAATCTGGATGG‑CT‑TTTTTTTGTCTATAGCGGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTa                        >  2:400535/8‑100 (MQ=25)
                                                                            aaTTATGAAAAGCCATCCAGATTTGGATGG‑TT‑TTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTa                        >  2:2321419/1‑100 (MQ=255)
                                                                             acaaaaaaaaACCATCCAAATCTGGATGG‑CT‑TTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTaa                       >  2:2178370/7‑100 (MQ=25)
                                                                              ttATGAAAAGCCATCCAGATTTGGATGG‑TT‑TTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAAt                      >  2:1921543/1‑100 (MQ=255)
                                                                              caaaaaaaaACCATCCAAATCTGGATGGCTT‑TTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTaa                       >  2:1707264/6‑100 (MQ=37)
                                                                               tATGAAAAGCCATCCAGATTCGGATGG‑TT‑TTTTTTTTTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAAt                      >  1:2691511/1‑99 (MQ=255)
                                                                                ataaaaaACCATCAAAATCTGGATGG‑TT‑TTTGTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCt                    >  2:2963093/4‑100 (MQ=255)
                                                                                aaaaaaaACCATCCAAATCTGGATGG‑TT‑TTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCt                    >  2:1422079/4‑100 (MQ=38)
                                                                                aTGAAAAGCCATCCATATTTGAATGG‑TT‑TTTTTTTGTTTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCt                    <  1:1595290/100‑1 (MQ=255)
                                                                                aTGAAAAGCCATCCAGATTTGGATGG‑TT‑TTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCt                    >  1:2963225/1‑100 (MQ=255)
                                                                                aTGAAAAGCCATCCAGATTTGGATGG‑TT‑TTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCt                    >  1:3123232/1‑100 (MQ=255)
                                                                                 tGAAAAGCCATCCAGATTTGGATGG‑TT‑TTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCta                   >  2:1518622/1‑100 (MQ=255)
                                                                                 tGAAAAGCCATCCAGATTTGGATGG‑TT‑TTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCta                   >  2:1597970/1‑100 (MQ=255)
                                                                                 aaaaaaACCATCCAAATCTGGATGG‑TT‑TTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCta                   >  1:446610/3‑100 (MQ=38)
                                                                                  gAAAAGCCATCCAGATTTGGATGG‑TT‑TTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCtat                  >  2:3245511/1‑100 (MQ=255)
                                                                                  gAAAAGCCATCCAGATTTGGATGG‑TT‑TTTTTTTGTCTATAGATGGTATGATAATTAACTATTGCAAATATTCATTTTATTAATATTTAAACTTAATc                     >  1:3615956/1‑97 (MQ=255)
                                                                                   aaaaGCCATTCAGATTTGGATTG‑TT‑TTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCtata                 <  1:3447558/100‑1 (MQ=255)
                                                                                   aaaaGCCATCCAGATTTGGATGG‑TT‑TTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCtata                 <  2:3777852/100‑1 (MQ=255)
                                                                                   aaaaGCCATCCAGATTTGGATGG‑TT‑TTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCtata                 >  2:4356476/1‑100 (MQ=255)
                                                                                    aaaGCCATCCAGATTTGGATGG‑TT‑TTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTatat                <  1:2321408/100‑1 (MQ=255)
                                                                                     taGCTATCGAGATTTTGATGG‑TT‑TTTTTTTGTCTATATCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCtata                 <  1:4062276/97‑1 (MQ=255)
                                                                                     cgatcttCCAGATTTGGATTG‑TT‑TTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATg               <  1:2335078/94‑1 (MQ=255)
                                                                                     aaGCCATCCAGATTTGGATGG‑TT‑TTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATg               <  2:3535724/100‑1 (MQ=255)
                                                                                      aGCCATCCAGATTTGGATTG‑TT‑TTTTTTTGTCTATAGCTTGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATgg              <  1:1810030/100‑1 (MQ=255)
                                                                                       gCCATCCAGATTTGGATGG‑TT‑TTTTTTTTTTTTTATCTGGTAAGATAATTACGTATTTCAAATATTTATTTTATTAATATTTAAACTTTATCTatat                >  1:148692/1‑97 (MQ=255)
                                                                                        ccATCCAGATTTGGATGG‑TT‑TTTTTTTGTCTATCGCTATTAAGATAATTACGTTTTGCAAATATTCATTTTATTAATATTTAAACTTAATCtata                 <  1:2394570/95‑1 (MQ=255)
                                                                                        ccATCCAGATTTGGATGG‑TT‑TTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGaa            >  2:3894214/1‑100 (MQ=255)
                                                                                           tCCAGATTTGGATGG‑TT‑TTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATgagatcg         >  2:2335507/1‑94 (MQ=255)
                                                                                           tCCAGATTTGGATGG‑TT‑TTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATgagatcg         >  2:2335089/1‑94 (MQ=255)
                                                                                             cAGATTTGGATGG‑TT‑TTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGAataata       <  2:891797/100‑1 (MQ=255)
                                                                                                  ttGGATGG‑TT‑TTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGAATAATATATGa  <  1:1921537/100‑1 (MQ=255)
                                                                                                  |                                                                                                     
TTTCTTCCCTGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TT‑TTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGAATAATATATGA  >  NC_000913/4409939‑4410136

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