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.94.7% 606.3 / 26.0 151intergenic (+39/‑88)rlmB/yjfI23S rRNA mG2251 2'‑O‑ribose methyltransferase, SAM‑dependent/DUF2170 family protein
Reads supporting (aligned to +/- strand):  ref base T (6/2);  new base . (59/84);  total (65/86)
Fisher's exact test for biased strand distribution p-value = 7.56e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 179 total aligned reads displayed.

CCTGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGAATAATATATGACATGG  >  NC_000913/4409946‑4410141
                                                                                                    |                                                                                                 
ccTGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGt                                                                                                   >  1:3011621/1‑100 (MQ=255)
ccTGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGt                                                                                                   >  1:1770527/1‑100 (MQ=255)
   gAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAACCATCCAAATCTGGATGGcttt                                                                                                <  1:231185/100‑1 (MQ=38)
      cgatctggTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGttttttt                                                                                             <  2:278523/94‑1 (MQ=255)
       gTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTATTTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGtttttttt                                                                                            <  2:1817913/100‑1 (MQ=255)
       gTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAACCATCCAAATCTGGATGGcttttcat                                                                                            <  2:2638610/100‑4 (MQ=37)
        tttCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCTTCCAGATTTGGATGGttttttttg                                                                                           <  1:3006889/100‑2 (MQ=255)
        tttCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGttttttttg                                                                                           <  2:1218794/100‑2 (MQ=255)
         ttCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGttttttttgt                                                                                          >  1:1184634/1‑98 (MQ=255)
          tCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGttttttttgtc                                                                                         <  2:855816/100‑4 (MQ=255)
          tCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAACCATCCAAATCTGGATGGcttttcataat                                                                                         <  1:1844988/100‑7 (MQ=25)
          tCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAACCATCCAAATCTGGATGGcttttcataat                                                                                         <  1:2633491/100‑7 (MQ=37)
          tCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAACCATCCAAATCTGGATGGcttttcataat                                                                                         <  2:374708/100‑7 (MQ=37)
          tCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAACCATCCAAATCTGGATGGcttttcataat                                                                                         <  1:1848354/100‑7 (MQ=25)
          tCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAACCATCCAAATCTGGATGGcttttcataat                                                                                         <  2:1495260/100‑7 (MQ=37)
           cGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGttttttttgtct                                                                                        >  1:81683/1‑96 (MQ=255)
              ttGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCtata                                                                                    <  2:79561/100‑1 (MQ=255)
              ttGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCtata                                                                                    <  2:2675000/100‑1 (MQ=255)
                gCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCtat                                                                                     >  2:346965/1‑97 (MQ=255)
                gCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGc                                                                                  >  1:1892787/1‑100 (MQ=255)
                 cGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCt                                                                                 <  2:2349714/100‑1 (MQ=255)
                  gACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTg                                                                                <  2:2466211/100‑1 (MQ=255)
                  gACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTg                                                                                <  2:59957/100‑1 (MQ=255)
                  gACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTg                                                                                >  1:730314/1‑100 (MQ=255)
                    ccGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGt                                                                              >  2:3094683/1‑100 (MQ=255)
                    ccGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGt                                                                              <  1:2360478/100‑1 (MQ=255)
                       gAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAg                                                                           >  1:944440/1‑100 (MQ=255)
                       gAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAg                                                                           <  2:3172124/100‑1 (MQ=255)
                       gAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAg                                                                           <  2:2496909/100‑1 (MQ=255)
                        aaTTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGa                                                                          <  2:1535222/100‑1 (MQ=255)
                        aaTTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGa                                                                          >  1:294784/1‑100 (MQ=255)
                         aTTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGAt                                                                         >  2:2656548/1‑100 (MQ=255)
                          tttGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCTGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATa                                                                        <  2:569954/100‑1 (MQ=255)
                           ttGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATaa                                                                       >  1:2468022/1‑100 (MQ=255)
                              cTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTa                                                                    <  2:935506/100‑1 (MQ=255)
                                 aTTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGt                                                                 >  1:1606414/1‑100 (MQ=255)
                                 aTTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGt                                                                 <  2:1504184/100‑1 (MQ=255)
                                 aTTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTCTTTTGTCTATAGCTGGTAAGATAATTACGt                                                                 >  1:1097432/1‑100 (MQ=255)
                                  tttGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTa                                                                >  1:1937061/1‑100 (MQ=255)
                                   ttGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTAt                                                               <  1:2724373/100‑1 (MQ=255)
                                   ttGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAAATATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTAt                                                               <  1:3111369/100‑1 (MQ=255)
                                    tGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTAtt                                                              >  2:1873857/1‑100 (MQ=255)
                                    tGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTAtt                                                              <  1:128142/100‑1 (MQ=255)
                                      aaGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGc                                                            <  1:2682217/100‑1 (MQ=255)
                                      aaGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGc                                                            <  1:2265517/100‑1 (MQ=255)
                                        gCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCaa                                                          <  2:2612630/100‑1 (MQ=255)
                                         cGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCaaa                                                         >  1:73765/1‑100 (MQ=255)
                                         cGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCaaa                                                         <  1:2501361/100‑1 (MQ=255)
                                          ggTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAat                                                        <  2:2717587/100‑1 (MQ=255)
                                            tGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAatat                                                     >  1:1971465/1‑100 (MQ=255)
                                             gcgcCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTCCAAATAtt                                                    <  2:1389726/100‑1 (MQ=255)
                                                tccgatctgCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCAt                                                 <  2:2304876/92‑1 (MQ=255)
                                                    cgcAGCTAATTTCTCAGAATTATGAAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCAtttt                                              <  2:1695740/100‑1 (MQ=255)
                                                      cAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTTTTCAATAGCTGGAAGGAAATTAACGTATTCAAATAATTCAtttttt                                            >  1:1520517/1‑98 (MQ=255)
                                                      cAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTttatt                                           <  2:2048743/100‑1 (MQ=255)
                                                       aGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTttatt                                           <  1:272957/100‑1 (MQ=255)
                                                       aGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTCCGTATTGCAAATATTCATTttatta                                          >  2:257057/1‑100 (MQ=255)
                                                       aGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTttatta                                          <  1:82598/100‑1 (MQ=255)
                                                       aGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTttatta                                          >  2:2861450/1‑100 (MQ=255)
                                                       aGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTttatta                                          <  2:528063/100‑1 (MQ=255)
                                                       aGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTttatta                                          <  1:100037/100‑1 (MQ=255)
                                                       aGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTttatt                                           >  2:953861/1‑99 (MQ=255)
                                                        gCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTaa                                         >  1:2885478/1‑100 (MQ=255)
                                                          tAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTttatt                                           <  2:2794362/96‑1 (MQ=255)
                                                            aTTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATAtt                                     <  2:2951344/100‑1 (MQ=255)
                                                                ctcAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTaaa                                 >  2:1149989/1‑100 (MQ=255)
                                                                ctcAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTaaa                                 >  1:2471816/1‑100 (MQ=255)
                                                                 tcAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTaaa                                 <  2:1436463/99‑1 (MQ=255)
                                                                 tcAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAAc                                >  2:2827564/1‑100 (MQ=255)
                                                                 tcAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAAc                                <  2:1770505/100‑1 (MQ=255)
                                                                  cAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAAc                                <  1:94491/99‑1 (MQ=255)
                                                                  cAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACt                               >  2:1471706/1‑100 (MQ=255)
                                                                     aaTTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTaa                            <  2:1892764/100‑1 (MQ=255)
                                                                      agacaaaaaaaACCATCCAAATCTGGATGGCTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTaa                            >  2:46959/8‑100 (MQ=25)
                                                                      aTTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAAt                           >  2:1933644/1‑100 (MQ=255)
                                                                       ttATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATc                          >  2:1264596/1‑100 (MQ=255)
                                                                       ttATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATc                          >  1:1423947/1‑100 (MQ=255)
                                                                        tATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCt                         <  2:117513/100‑1 (MQ=255)
                                                                         caaaaaaaACCATCCAAATCTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCt                         >  2:1979962/5‑100 (MQ=38)
                                                                         aTGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCta                        <  2:1754192/100‑1 (MQ=255)
                                                                           gAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCtatn                      <  1:686232/100‑2 (MQ=255)
                                                                           gAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCtata                      >  1:1522780/1‑100 (MQ=255)
                                                                           aaaaaaACCATCCAAATCTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCtat                       >  2:479716/3‑100 (MQ=38)
                                                                             aaaaGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTatat                     <  2:1096003/100‑1 (MQ=255)
                                                                             aaaaGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTatat                     <  1:1471078/100‑1 (MQ=255)
                                                                               aaaccATCCAAATCTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATgg                   >  1:416467/4‑100 (MQ=255)
                                                                                aGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATgg                   <  2:2207320/99‑1 (MQ=255)
                                                                                 gCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAATCTTAATCTATATGGaa                 >  2:335709/1‑100 (MQ=255)
                                                                                  ccATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGaat                <  1:1675/100‑1 (MQ=255)
                                                                                     tCCAAATTTGGATGGTTTTTTTTGTCTATAGTTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGaataat             <  2:738224/100‑1 (MQ=255)
                                                                                      ccAGATTTGGATGGTTTTTTTTTTTTTTTTTTTGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGAataata            >  2:679359/1‑100 (MQ=255)
                                                                                      ccAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGAataata            <  1:1167331/100‑1 (MQ=255)
                                                                                       caaatctGGATGGTTTTTTTTTTTTTTAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAATCTTAATCTATATGGaataa              >  2:1290504/7‑97 (MQ=255)
                                                                                       cAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGAataata            <  2:3011585/99‑1 (MQ=255)
                                                                                       cAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGAATAatat           >  1:1038152/1‑100 (MQ=255)
                                                                                        aGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGAATAatata          <  2:789487/100‑1 (MQ=255)
                                                                                         gATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGAATAatatat         <  1:2820478/100‑1 (MQ=255)
                                                                                          aTTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGAATAATATATg        <  1:1022419/100‑1 (MQ=255)
                                                                                               gATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGAATAATATATGACATg   <  2:485424/100‑1 (MQ=255)
                                                                                               gATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGAATAATATATGACATg   <  1:2053455/100‑1 (MQ=255)
                                                                                                    |                                                                                                 
CCTGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGAATAATATATGACATGG  >  NC_000913/4409946‑4410141

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