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.95.6% 564.9 / 14.6 131intergenic (+39/‑88)rlmB/yjfI23S rRNA mG2251 2'‑O‑ribose methyltransferase, SAM‑dependent/DUF2170 family protein
Reads supporting (aligned to +/- strand):  ref base T (3/2);  new base . (56/69);  total (60/71)
Fisher's exact test for biased strand distribution p-value = 6.58e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 157 total aligned reads displayed.

CCTGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGAATAATATATGACATGGA  >  NC_000913/4409946‑4410142
                                                                                                   |                                                                                                 
ccTGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGt                                                                                                   >  1:1811370/1‑100 (MQ=255)
    aaCGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGttttt                                                                                               >  1:2152989/1‑100 (MQ=255)
      cGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGttttttt                                                                                             >  2:2663989/1‑100 (MQ=255)
      cGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGttttttt                                                                                             <  2:2384616/100‑1 (MQ=255)
      cGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGttttttt                                                                                             >  2:124829/1‑100 (MQ=255)
      cGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGttttttt                                                                                             >  2:2608934/1‑100 (MQ=255)
      cGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAACCATCCAAATCTGGATGGcttttca                                                                                             <  2:2876549/100‑3 (MQ=38)
      cGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAACCATCCAAATCTGGATGGcttttca                                                                                             <  2:901183/100‑3 (MQ=38)
       gTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGtttttttt                                                                                            <  2:1039786/100‑1 (MQ=255)
       gTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAACCATCCAAATCTGGATGGcttttcat                                                                                            <  1:949225/100‑4 (MQ=37)
        tttCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGttttttttg                                                                                           <  2:498319/100‑2 (MQ=255)
        tttCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGttttttttg                                                                                           <  1:1769351/100‑2 (MQ=255)
        tttCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGttttttttg                                                                                           >  1:2311282/1‑99 (MQ=255)
        tttCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGttttttttg                                                                                           <  2:394953/100‑2 (MQ=255)
        tttCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGttttttttg                                                                                           <  2:1575299/100‑2 (MQ=255)
         ttCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTTg                                                                                          <  2:2410794/100‑1 (MQ=255)
          tCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGttttttttgtc                                                                                         <  2:1452967/100‑4 (MQ=255)
          tCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAACCATCCAAATCTGGATGGcttttcataat                                                                                         <  2:1762140/100‑7 (MQ=37)
           cGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGttttttttgtct                                                                                        <  2:2990646/100‑5 (MQ=255)
             gTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCtat                                                                                     <  2:3343676/100‑1 (MQ=255)
                gCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGc                                                                                  <  2:1343594/100‑1 (MQ=255)
                 cGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCt                                                                                 >  1:214815/1‑100 (MQ=255)
                 cGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCt                                                                                 >  2:2074634/1‑100 (MQ=255)
                 cGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCt                                                                                 >  2:260279/1‑100 (MQ=255)
                 cGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCt                                                                                 >  2:765715/1‑100 (MQ=255)
                   aCCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTgg                                                                               <  2:1648420/100‑1 (MQ=255)
                    ccGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGt                                                                              >  2:408662/1‑100 (MQ=255)
                     cGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTa                                                                             <  1:2496471/100‑1 (MQ=255)
                      ggAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTaa                                                                            <  1:1706723/100‑1 (MQ=255)
                      ggAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTaa                                                                            <  1:1706939/100‑1 (MQ=255)
                      ggAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTaa                                                                            <  2:2501060/100‑1 (MQ=255)
                      ggAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTaa                                                                            >  1:1819792/1‑100 (MQ=255)
                       gAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAg                                                                           <  1:2354013/100‑1 (MQ=255)
                         aTTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGAt                                                                         >  2:1758774/1‑100 (MQ=255)
                            tGCTTATGTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCATAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATa                                                                        <  1:2873662/98‑1 (MQ=255)
                             gCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAAtt                                                                     >  1:2988259/1‑100 (MQ=255)
                              cTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCCTCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATa                                                                        >  1:760536/1‑96 (MQ=255)
                                tATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACg                                                                  <  1:2754910/100‑1 (MQ=255)
                                tATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACg                                                                  >  2:1993595/1‑100 (MQ=255)
                                 aTTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGt                                                                 >  1:1570700/1‑100 (MQ=255)
                                  tttGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTa                                                                >  1:562998/1‑100 (MQ=255)
                                    tGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTTTTTTTTGTTGGTAAGATAATTACGTAtt                                                              >  1:1179738/1‑100 (MQ=255)
                                    tGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTAtt                                                              <  2:2241195/100‑1 (MQ=255)
                                    tGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTAtt                                                              >  1:1592604/1‑100 (MQ=255)
                                      aaGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGc                                                            >  2:1224663/1‑100 (MQ=255)
                                        gCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCaa                                                          >  1:1201277/1‑100 (MQ=255)
                                           gTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAata                                                       >  1:514711/1‑100 (MQ=255)
                                            tGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAatat                                                      >  2:385291/1‑100 (MQ=255)
                                             gcgcCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATAtt                                                     <  2:1646460/100‑1 (MQ=255)
                                              cgcCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTc                                                    >  2:1530763/1‑100 (MQ=255)
                                                ccAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTTTCTTTTTCTTGTAAGATAATTACGTATTTCAAATATTCAt                                                  >  1:2051542/1‑100 (MQ=255)
                                                ccAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCAt                                                  <  1:1344342/100‑1 (MQ=255)
                                                  aGCGCAGCTAATTTCTCAGAATTATGTAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCAtt                                                 <  1:722997/99‑1 (MQ=255)
                                                   gcgcAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCAtttt                                               <  2:2043649/100‑1 (MQ=255)
                                                    cgcAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTtta                                              <  1:813359/100‑1 (MQ=255)
                                                     gcAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTttat                                             >  1:1027956/1‑100 (MQ=255)
                                                     gcAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTttat                                             <  2:1078581/100‑1 (MQ=255)
                                                       aGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGTTGGTAAGATAATTACGTATTGCAAATATTCATTttatt                                            >  1:930584/1‑99 (MQ=255)
                                                          tAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAata                                        <  2:1967911/100‑1 (MQ=255)
                                                           aaTTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAatat                                       <  2:214817/100‑1 (MQ=255)
                                                            aTTTTTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATAtt                                      <  1:1004042/100‑1 (MQ=255)
                                                              ttCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATAtt                                      <  1:3157533/98‑1 (MQ=255)
                                                                ctcAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTaaa                                  >  2:1444497/1‑100 (MQ=255)
                                                                ctcAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTaaa                                  >  2:76731/1‑100 (MQ=255)
                                                                ctcAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTaaa                                  <  2:802386/100‑1 (MQ=255)
                                                                ctcAGAATTATGAAAAACCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTaaa                                  >  1:463360/1‑100 (MQ=255)
                                                                 tcttatttttGAAAAGCAATCCAGATTTATATGGTTTTTTTTGTCTATATCTGGTAAAATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAAc                                 <  1:3427065/91‑1 (MQ=255)
                                                                 tcAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTTGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTaaa                                  <  1:544658/99‑1 (MQ=255)
                                                                 tcAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAAc                                 >  1:1928/1‑100 (MQ=255)
                                                                 tcAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACATATTGCAAATATTCATTTTATTAATATTTAAAc                                 >  2:2749656/1‑100 (MQ=255)
                                                                  cAGAATTATGAAAAGCCATCAATATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACt                                <  2:1177542/100‑1 (MQ=255)
                                                                   aGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACt                                <  2:42294/99‑1 (MQ=255)
                                                                      aTTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAAt                            >  2:3230073/1‑100 (MQ=255)
                                                                      aTTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAAt                            <  2:1305367/100‑1 (MQ=255)
                                                                       ttATGAAAAGCGATCCAGATTCGGATGGTTCTTTTTCTTTATAGTTTGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAAt                            <  1:2327206/99‑1 (MQ=255)
                                                                       ttATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAAc                                 <  1:3208028/94‑1 (MQ=255)
                                                                       ttATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACt                                <  1:3340801/95‑1 (MQ=255)
                                                                       ttATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATc                           <  2:186460/100‑1 (MQ=255)
                                                                        caaaaaaaACCATCCAAATCTGGATGGCTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCt                          >  1:1783896/5‑100 (MQ=37)
                                                                         aaaaaaaACCATCCAAATCTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCta                         >  2:2875275/4‑100 (MQ=38)
                                                                         aTGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCta                         >  2:1005090/1‑100 (MQ=255)
                                                                           gAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCtata                       <  2:2311300/100‑1 (MQ=255)
                                                                           gAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCtat                        >  1:2692715/1‑99 (MQ=255)
                                                                            aaaaGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTatat                      >  1:1961020/1‑100 (MQ=255)
                                                                             aaaGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATg                     <  1:1606124/100‑1 (MQ=255)
                                                                             aaaGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATg                     <  2:3114192/100‑1 (MQ=255)
                                                                             aaaGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATg                     <  2:464913/100‑1 (MQ=255)
                                                                                 ccATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGaat                 <  2:971102/100‑1 (MQ=255)
                                                                                   aTCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGaataa               >  2:3345769/1‑100 (MQ=255)
                                                                                    tCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGaataat              <  1:385288/100‑1 (MQ=255)
                                                                                     ccAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGAataata             <  2:1570678/100‑1 (MQ=255)
                                                                                     ccAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGAataata             >  1:1499158/1‑100 (MQ=255)
                                                                                      cAGATTTGGATGGTTTTTTTTGTCTATATCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGAATAatat            >  1:2598227/1‑100 (MQ=255)
                                                                                       aaatctGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGAATAatata           >  2:2493160/6‑100 (MQ=255)
                                                                                       aGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGAATAatata           <  1:2608915/100‑1 (MQ=255)
                                                                                         aTTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGAATAATATATg         <  2:528997/100‑1 (MQ=255)
                                                                                         aTTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGAATAATATATg         <  2:1027965/100‑1 (MQ=255)
                                                                                         aTTTGGATGGGTTTTTTTGTCTATAGCTTGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGAATAATATATg         <  1:2942645/100‑1 (MQ=255)
                                                                                          tttGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGAATAATATATGa        >  2:2774531/1‑100 (MQ=255)
                                                                                                tGGTTTTTTTTTTCTATAACTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGATTAATATATGAccttaa  >  1:2437857/1‑95 (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 < 32 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: