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.

TTTCTTCCCTGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGAATAATATATGA  >  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)
       cttGAACGTTTCGATTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGt                                                                                              >  2:2990829/3‑100 (MQ=255)
       ccTGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGt                                                                                              >  1:2003809/1‑100 (MQ=255)
       ccTGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGt                                                                                              >  2:3208043/1‑100 (MQ=255)
        cTGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAACAGCCATCCAGATTTGGATGGtt                                                                                             >  2:1773254/1‑100 (MQ=255)
           aaCGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAAATCTGGATGGctttt                                                                                         <  2:2468866/100‑1 (MQ=255)
            acttttCGGTTGCGACCGGAATTTTCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAAATCTGGATGGcttttc                                                                                        <  2:3188900/97‑2 (MQ=255)
              gTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGtttttttt                                                                                      <  1:210235/100‑1 (MQ=255)
                ttCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTTg                                                                                    <  1:76920/100‑1 (MQ=255)
                ttCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAACCATCCAAATCTGGATGGcttttcat                                                                                      <  1:1175051/98‑4 (MQ=37)
                tccGGTGGCGACCTGAATTTGCTTACTTAAAGCGGTGCGCCAGCGCAGCTAATGTCTCAGAATTATGAAAAGCCATCCAAATCTGGATGGcttttc                                                                                        <  1:3883929/94‑2 (MQ=37)
                 tCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAACCATCCAAATCTGGATGGcttttcataat                                                                                   <  2:4066807/100‑7 (MQ=37)
                  cGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTTGTc                                                                                  <  2:2829659/100‑1 (MQ=255)
                   ggTTGCTGCCGGAATTTGCTTATTGGAAGCGGTGCGCCAGCGCAGCTAATTTCTCTGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTTGTCt                                                                                 <  1:3740468/100‑1 (MQ=255)
                   ggTTGCCACAGGAATTTGCTTATTGAAAGCGTTGCGCCAGCCCAGCCAATTTCTCAGAATTATGAAAAACCATCCAACTCTGGAGGGcttttca                                                                                       <  1:327368/94‑3 (MQ=255)
                      tGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTTGTCtata                                                                              >  2:1197168/1‑100 (MQ=255)
                        cGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTTGTCTATAGc                                                                            >  2:1837039/1‑100 (MQ=255)
                          aCCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTTGTCTATAGCTg                                                                          <  2:4250961/100‑1 (MQ=255)
                           ccGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTTGTCTATAGCTgg                                                                         >  2:19576/1‑100 (MQ=255)
                           ccGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTTGTCTATAGCTgg                                                                         >  2:2236388/1‑100 (MQ=255)
                             ggAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTTGTCTATAGCTGGTa                                                                       <  1:3262269/100‑1 (MQ=255)
                                 tttGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTTGTCTATAGCTGGTAAGAt                                                                   <  1:3379288/100‑1 (MQ=255)
                                   tGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCAAGATTTGGATGGATTTTTTTTGTTTATTGCTTTTAAGTTaa                                                                 >  1:2607343/1‑100 (MQ=255)
                                    gCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTTGTCTATAGCTGGTAAGATAAt                                                                >  2:2874697/1‑100 (MQ=255)
                                     cTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTTGTCTATAGCTGGTAAGATAAtt                                                               >  2:1567241/1‑100 (MQ=255)
                                       tATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTTGTCTATAGCTGGTAAGATAATTAc                                                             >  2:521193/1‑100 (MQ=255)
                                       tATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTTGTCTATAGCTGGTAAGATAATTAc                                                             >  2:943707/1‑100 (MQ=255)
                                        aTTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTTGTCTATAGCTGGTAAGATAATTACg                                                            >  2:1463177/1‑100 (MQ=255)
                                         tttGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTTGTCTATAGCTGGTAAGATAATTACg                                                            <  2:619526/99‑1 (MQ=255)
                                         tttGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTTGTCTATAGCTGGTAAGATAATTACGt                                                           <  2:3272653/100‑1 (MQ=255)
                                         tttGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTTGTCTATAGCTGGTAAGATAATTACGt                                                           >  1:3225088/1‑100 (MQ=255)
                                          ttGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGCTTTGGATGGTTTTTTTTTGTCTATAGATGGTAAGATAATTACGTa                                                          >  1:3520256/1‑100 (MQ=255)
                                           tGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTAt                                                         >  2:3594132/1‑100 (MQ=255)
                                            ggaGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTAt                                                         <  2:581739/97‑1 (MQ=255)
                                            gAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTAtt                                                        >  2:1037555/1‑100 (MQ=255)
                                              aGCGGTGCTCCAGCGCAGCTAATGTCTCAGAATTATGAAAAGCCCTCCAGATGTGCATGGTTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGc                                                      <  1:342922/100‑1 (MQ=255)
                                              aGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTTCAAAGGCTGGTAAGTTAATTActtttt                                                        >  2:3377874/1‑93 (MQ=255)
                                               gCGGGGCCCCGGCGCACATATTTTTTCATATTTATAAAAAGCCACCCCATTTTGGTTGTTTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTAtt                                                        <  2:312763/97‑1 (MQ=255)
                                                  gTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAat                                                  >  2:2583957/1‑100 (MQ=255)
                                                   tGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAata                                                 <  2:2122878/100‑1 (MQ=255)
                                                   tGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAata                                                 <  1:3054162/100‑1 (MQ=255)
                                                    gcgcCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTTTTCTTTTTTTGATAAGATAATTACGTATTGCAAatat                                                >  1:1572820/1‑100 (MQ=255)
                                                     cgcCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATAtt                                               >  2:3485655/1‑100 (MQ=255)
                                                      gcCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTc                                              >  2:3026518/1‑100 (MQ=255)
                                                      gcCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCa                                             <  2:1921275/100‑1 (MQ=255)
                                                             cAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTttat                                       >  2:1586504/1‑100 (MQ=255)
                                                             cAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTttat                                       <  1:1042749/100‑1 (MQ=255)
                                                              aGCTAATTTCGCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTTGTTCATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTttat                                       <  1:2055782/99‑1 (MQ=255)
                                                               gCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTttatta                                     <  1:3026501/100‑1 (MQ=255)
                                                               gCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTttatta                                     <  2:48684/100‑1 (MQ=255)
                                                                    tttCTCATAATTATGGAAAGACAGCAATATTTGGATTGTTTTTTTGTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAata                                  <  1:84203/98‑1 (MQ=255)
                                                                    tttCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTTGTTTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATAtt                                >  1:677698/1‑100 (MQ=255)
                                                                    tttCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTTGTATATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAata                                  >  2:4274488/1‑98 (MQ=255)
                                                                     tgttttgAATTATGTAAAGCCATTTTTTTTTTTATGGTTGTTTTTTTGTCTATACCTGGTAATATAATTACGTATTGCAAATATTCATTTTATTAata                                  <  1:1337252/92‑1 (MQ=255)
                                                                      tctcAGATTTTGGAACAGCCATCGAGATTTTGATTGTTTTTTTTTGTCTATAACTTCTACGATAATTACGTATTGCAAATATTCATTTTATTAATATTTa                              <  1:4157355/100‑1 (MQ=255)
                                                                       ctcAGAATTATGAAAAGCCATCTAGATTTGGATGGTTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTaa                             <  1:2519844/100‑1 (MQ=255)
                                                                       ctcAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTaa                             >  2:1126133/1‑100 (MQ=255)
                                                                        tcttaaTTATGATAAGCCATCTATATTTTGATGGTTTTTTTTTGTATATAGCTTTTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTaaa                            <  1:2175525/96‑1 (MQ=255)
                                                                        tcAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTACAAATATTCATTTTATTAATATTTaaa                            <  2:3506557/100‑1 (MQ=255)
                                                                         cAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTTGTCTTTAGCTGGTATGATAATTACGTATTGCAAATATTCATTTTATTGATATTTaaa                            >  1:2813267/1‑99 (MQ=255)
                                                                         cAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTTGTCTATAGCTGTTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAAc                           <  1:2449926/100‑1 (MQ=255)
                                                                          aGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTTTTTTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACt                          >  1:2192982/1‑100 (MQ=255)
                                                                            gacaaaaaaaaACCATCCAAATCTGGATGGCTTTTTTTTGTCTATAGCGGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTa                        >  2:400535/8‑100 (MQ=25)
                                                                            aaTTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTa                        >  2:2321419/1‑100 (MQ=255)
                                                                             acaaaaaaaaACCATCCAAATCTGGATGGCTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTaa                       >  2:2178370/7‑100 (MQ=25)
                                                                              ttATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAAt                      >  2:1921543/1‑100 (MQ=255)
                                                                              caaaaaaaaACCATCCAAATCTGGATGGCTTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTaa                       >  2:1707264/6‑100 (MQ=37)
                                                                               tATGAAAAGCCATCCAGATTCGGATGGTTTTTTTTTTTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAAt                      >  1:2691511/1‑99 (MQ=255)
                                                                                ataaaaaACCATCAAAATCTGGATGGTTTTTGTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCt                    >  2:2963093/4‑100 (MQ=255)
                                                                                aaaaaaaACCATCCAAATCTGGATGGTTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCt                    >  2:1422079/4‑100 (MQ=38)
                                                                                aTGAAAAGCCATCCATATTTGAATGGTTTTTTTTTGTTTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCt                    <  1:1595290/100‑1 (MQ=255)
                                                                                aTGAAAAGCCATCCAGATTTGGATGGTTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCt                    >  1:2963225/1‑100 (MQ=255)
                                                                                aTGAAAAGCCATCCAGATTTGGATGGTTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCt                    >  1:3123232/1‑100 (MQ=255)
                                                                                 tGAAAAGCCATCCAGATTTGGATGGTTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCta                   >  2:1518622/1‑100 (MQ=255)
                                                                                 tGAAAAGCCATCCAGATTTGGATGGTTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCta                   >  2:1597970/1‑100 (MQ=255)
                                                                                 aaaaaaACCATCCAAATCTGGATGGTTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCta                   >  1:446610/3‑100 (MQ=38)
                                                                                  gAAAAGCCATCCAGATTTGGATGGTTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCtat                  >  2:3245511/1‑100 (MQ=255)
                                                                                  gAAAAGCCATCCAGATTTGGATGGTTTTTTTTTGTCTATAGATGGTATGATAATTAACTATTGCAAATATTCATTTTATTAATATTTAAACTTAATc                     >  1:3615956/1‑97 (MQ=255)
                                                                                   aaaaGCCATTCAGATTTGGATTGTTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCtata                 <  1:3447558/100‑1 (MQ=255)
                                                                                   aaaaGCCATCCAGATTTGGATGGTTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCtata                 <  2:3777852/100‑1 (MQ=255)
                                                                                   aaaaGCCATCCAGATTTGGATGGTTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCtata                 >  2:4356476/1‑100 (MQ=255)
                                                                                    aaaGCCATCCAGATTTGGATGGTTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTatat                <  1:2321408/100‑1 (MQ=255)
                                                                                     taGCTATCGAGATTTTGATGGTTTTTTTTTGTCTATATCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCtata                 <  1:4062276/97‑1 (MQ=255)
                                                                                     cgatcttCCAGATTTGGATTGTTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATg               <  1:2335078/94‑1 (MQ=255)
                                                                                     aaGCCATCCAGATTTGGATGGTTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATg               <  2:3535724/100‑1 (MQ=255)
                                                                                      aGCCATCCAGATTTGGATTGTTTTTTTTTGTCTATAGCTTGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATgg              <  1:1810030/100‑1 (MQ=255)
                                                                                       gCCATCCAGATTTGGATGGTTTTTTTTTTTTTTTATCTGGTAAGATAATTACGTATTTCAAATATTTATTTTATTAATATTTAAACTTTATCTatat                >  1:148692/1‑97 (MQ=255)
                                                                                        ccATCCAGATTTGGATGGTTTTTTTTTGTCTATCGCTATTAAGATAATTACGTTTTGCAAATATTCATTTTATTAATATTTAAACTTAATCtata                 <  1:2394570/95‑1 (MQ=255)
                                                                                        ccATCCAGATTTGGATGGTTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGaa            >  2:3894214/1‑100 (MQ=255)
                                                                                           tCCAGATTTGGATGGTTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATgagatcg         >  2:2335507/1‑94 (MQ=255)
                                                                                           tCCAGATTTGGATGGTTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATgagatcg         >  2:2335089/1‑94 (MQ=255)
                                                                                             cAGATTTGGATGGTTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGAataata       <  2:891797/100‑1 (MQ=255)
                                                                                                  ttGGATGGTTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGAATAATATATGa  <  1:1921537/100‑1 (MQ=255)
                                                                                                  |                                                                                                     
TTTCTTCCCTGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGAATAATATATGA  >  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: