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.

CCTGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATG‑AAAAGCCATCCAGATTTGGATGGTTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAAT‑ATTCATTTTATTAATATTTAAACTTAATCTATATGGAATAATATATGACATGG  >  NC_000913/4409946‑4410141
                                                                                                    |                                                                                                 
ccTGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATG‑AAAAGCCATCCAGATTTGGATGGt                                                                                                   >  1:3011621/1‑100 (MQ=255)
ccTGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATG‑AAAAGCCATCCAGATTTGGATGGt                                                                                                   >  1:1770527/1‑100 (MQ=255)
   gAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATG‑AAAAACCATCCAAATCTGGATGGcttt                                                                                                <  1:231185/100‑1 (MQ=38)
      cgatctggTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATG‑AAAAGCCATCCAGATTTGGATGGttttttt                                                                                             <  2:278523/94‑1 (MQ=255)
       gTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTATTTTCTCAGAATTATG‑AAAAGCCATCCAGATTTGGATGGtttttttt                                                                                            <  2:1817913/100‑1 (MQ=255)
       gTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATG‑AAAAACCATCCAAATCTGGATGGcttttcat                                                                                            <  2:2638610/100‑4 (MQ=37)
        tttCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATG‑AAAAGCCTTCCAGATTTGGATGGttttttttg                                                                                           <  1:3006889/100‑2 (MQ=255)
        tttCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATG‑AAAAGCCATCCAGATTTGGATGGttttttttg                                                                                           <  2:1218794/100‑2 (MQ=255)
         ttCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATG‑AAAAGCCATCCAGATTTGGATGGttttttttgt                                                                                          >  1:1184634/1‑98 (MQ=255)
          tCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATG‑AAAAGCCATCCAGATTTGGATGGttttttttgtc                                                                                         <  2:855816/100‑4 (MQ=255)
          tCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATG‑AAAAACCATCCAAATCTGGATGGcttttcataat                                                                                         <  1:1844988/100‑7 (MQ=25)
          tCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATG‑AAAAACCATCCAAATCTGGATGGcttttcataat                                                                                         <  1:2633491/100‑7 (MQ=37)
          tCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATG‑AAAAACCATCCAAATCTGGATGGcttttcataat                                                                                         <  2:374708/100‑7 (MQ=37)
          tCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATG‑AAAAACCATCCAAATCTGGATGGcttttcataat                                                                                         <  1:1848354/100‑7 (MQ=25)
          tCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATG‑AAAAACCATCCAAATCTGGATGGcttttcataat                                                                                         <  2:1495260/100‑7 (MQ=37)
           cGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATG‑AAAAGCCATCCAGATTTGGATGGttttttttgtct                                                                                        >  1:81683/1‑96 (MQ=255)
              ttGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATG‑AAAAGCCATCCAGATTTGGATGG‑TTTTTTTTGTCtata                                                                                    <  2:79561/100‑1 (MQ=255)
              ttGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATG‑AAAAGCCATCCAGATTTGGATGG‑TTTTTTTTGTCtata                                                                                    <  2:2675000/100‑1 (MQ=255)
                gCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATG‑AAAAGCCATCCAGATTTGGATGG‑TTTTTTTTGTCtat                                                                                     >  2:346965/1‑97 (MQ=255)
                gCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATG‑AAAAGCCATCCAGATTTGGATGG‑TTTTTTTTGTCTATAGc                                                                                  >  1:1892787/1‑100 (MQ=255)
                 cGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATG‑AAAAGCCATCCAGATTTGGATGG‑TTTTTTTTGTCTATAGCt                                                                                 <  2:2349714/100‑1 (MQ=255)
                  gACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATG‑AAAAGCCATCCAGATTTGGATGG‑TTTTTTTTGTCTATAGCTg                                                                                <  2:2466211/100‑1 (MQ=255)
                  gACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATG‑AAAAGCCATCCAGATTTGGATGG‑TTTTTTTTGTCTATAGCTg                                                                                <  2:59957/100‑1 (MQ=255)
                  gACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATG‑AAAAGCCATCCAGATTTGGATGG‑TTTTTTTTGTCTATAGCTg                                                                                >  1:730314/1‑100 (MQ=255)
                    ccGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATG‑AAAAGCCATCCAGATTTGGATGG‑TTTTTTTTGTCTATAGCTGGt                                                                              >  2:3094683/1‑100 (MQ=255)
                    ccGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATG‑AAAAGCCATCCAGATTTGGATGG‑TTTTTTTTGTCTATAGCTGGt                                                                              <  1:2360478/100‑1 (MQ=255)
                       gAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATG‑AAAAGCCATCCAGATTTGGATGG‑TTTTTTTTGTCTATAGCTGGTAAg                                                                           >  1:944440/1‑100 (MQ=255)
                       gAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATG‑AAAAGCCATCCAGATTTGGATGG‑TTTTTTTTGTCTATAGCTGGTAAg                                                                           <  2:3172124/100‑1 (MQ=255)
                       gAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATG‑AAAAGCCATCCAGATTTGGATGG‑TTTTTTTTGTCTATAGCTGGTAAg                                                                           <  2:2496909/100‑1 (MQ=255)
                        aaTTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATG‑AAAAGCCATCCAGATTTGGATGG‑TTTTTTTTGTCTATAGCTGGTAAGa                                                                          <  2:1535222/100‑1 (MQ=255)
                        aaTTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATG‑AAAAGCCATCCAGATTTGGATGG‑TTTTTTTTGTCTATAGCTGGTAAGa                                                                          >  1:294784/1‑100 (MQ=255)
                         aTTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATG‑AAAAGCCATCCAGATTTGGATGG‑TTTTTTTTGTCTATAGCTGGTAAGAt                                                                         >  2:2656548/1‑100 (MQ=255)
                          tttGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCTGAATTATG‑AAAAGCCATCCAGATTTGGATGG‑TTTTTTTTGTCTATAGCTGGTAAGATa                                                                        <  2:569954/100‑1 (MQ=255)
                           ttGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATG‑AAAAGCCATCCAGATTTGGATGG‑TTTTTTTTGTCTATAGCTGGTAAGATaa                                                                       >  1:2468022/1‑100 (MQ=255)
                              cTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATG‑AAAAGCCATCCAGATTTGGATGG‑TTTTTTTTGTCTATAGCTGGTAAGATAATTa                                                                    <  2:935506/100‑1 (MQ=255)
                                 aTTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATG‑AAAAGCCATCCAGATTTGGATGG‑TTTTTTTTGTCTATAGCTGGTAAGATAATTACGt                                                                 >  1:1606414/1‑100 (MQ=255)
                                 aTTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATG‑AAAAGCCATCCAGATTTGGATGG‑TTTTTTTTGTCTATAGCTGGTAAGATAATTACGt                                                                 <  2:1504184/100‑1 (MQ=255)
                                 aTTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATG‑AAAAGCCATCCAGATTTGGATGG‑TTTCTTTTGTCTATAGCTGGTAAGATAATTACGt                                                                 >  1:1097432/1‑100 (MQ=255)
                                  tttGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATG‑AAAAGCCATCCAGATTTGGATGG‑TTTTTTTTGTCTATAGCTGGTAAGATAATTACGTa                                                                >  1:1937061/1‑100 (MQ=255)
                                   ttGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATG‑AAAAGCCATCCAGATTTGGATGG‑TTTTTTTTGTCTATAGCTGGTAAGATAATTACGTAt                                                               <  1:2724373/100‑1 (MQ=255)
                                   ttGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAAATATG‑AAAAGCCATCCAGATTTGGATGG‑TTTTTTTTGTCTATAGCTGGTAAGATAATTACGTAt                                                               <  1:3111369/100‑1 (MQ=255)
                                    tGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATG‑AAAAGCCATCCAGATTTGGATGG‑TTTTTTTTGTCTATAGCTGGTAAGATAATTACGTAtt                                                              >  2:1873857/1‑100 (MQ=255)
                                    tGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATG‑AAAAGCCATCCAGATTTGGATGG‑TTTTTTTTGTCTATAGCTGGTAAGATAATTACGTAtt                                                              <  1:128142/100‑1 (MQ=255)
                                      aaGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATG‑AAAAGCCATCCAGATTTGGATGG‑TTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGc                                                            <  1:2682217/100‑1 (MQ=255)
                                      aaGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATG‑AAAAGCCATCCAGATTTGGATGG‑TTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGc                                                            <  1:2265517/100‑1 (MQ=255)
                                        gCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATG‑AAAAGCCATCCAGATTTGGATGG‑TTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCaa                                                          <  2:2612630/100‑1 (MQ=255)
                                         cGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATG‑AAAAGCCATCCAGATTTGGATGG‑TTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCaaa                                                         >  1:73765/1‑100 (MQ=255)
                                         cGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATG‑AAAAGCCATCCAGATTTGGATGG‑TTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCaaa                                                         <  1:2501361/100‑1 (MQ=255)
                                          ggTGCGCCAGCGCAGCTAATTTCTCAGAATTATG‑AAAAGCCATCCAGATTTGGATGG‑TTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAat                                                        <  2:2717587/100‑1 (MQ=255)
                                            tGCGCCAGCGCAGCTAATTTCTCAGAATTATG‑AAAAGCCATCCAGATTTGGATGG‑TTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAat‑at                                                     >  1:1971465/1‑100 (MQ=255)
                                             gcgcCAGCGCAGCTAATTTCTCAGAATTATG‑AAAAGCCATCCAGATTTGGATGG‑TTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTCCAAAT‑Att                                                    <  2:1389726/100‑1 (MQ=255)
                                                tccgatctgCTAATTTCTCAGAATTATG‑AAAAGCCATCCAGATTTGGATGG‑TTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAAT‑ATTCAt                                                 <  2:2304876/92‑1 (MQ=255)
                                                    cgcAGCTAATTTCTCAGAATTATGAAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAAT‑ATTCAtttt                                              <  2:1695740/100‑1 (MQ=255)
                                                      cAGCTAATTTCTCAGAATTATG‑AAAAGCCATCCAGATTTGGATGGTTTTTTTTTTTCAATAGCTGGAAGGAAATTAACGTATT‑CAAATAATTCAtttttt                                            >  1:1520517/1‑98 (MQ=255)
                                                      cAGCTAATTTCTCAGAATTATG‑AAAAGCCATCCAGATTTGGATGG‑TTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAAT‑ATTCATTttatt                                           <  2:2048743/100‑1 (MQ=255)
                                                       aGCTAATTTCTCAGAATTATG‑AAAAGCCATCCAGATTTGGATGGTTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAAT‑ATTCATTttatt                                           <  1:272957/100‑1 (MQ=255)
                                                       aGCTAATTTCTCAGAATTATG‑AAAAGCCATCCAGATTTGGATGG‑TTTTTTTTGTCTATAGCTGGTAAGATAATTCCGTATTGCAAAT‑ATTCATTttatta                                          >  2:257057/1‑100 (MQ=255)
                                                       aGCTAATTTCTCAGAATTATG‑AAAAGCCATCCAGATTTGGATGG‑TTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAAT‑ATTCATTttatta                                          <  1:82598/100‑1 (MQ=255)
                                                       aGCTAATTTCTCAGAATTATG‑AAAAGCCATCCAGATTTGGATGG‑TTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAAT‑ATTCATTttatta                                          >  2:2861450/1‑100 (MQ=255)
                                                       aGCTAATTTCTCAGAATTATG‑AAAAGCCATCCAGATTTGGATGG‑TTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAAT‑ATTCATTttatta                                          <  2:528063/100‑1 (MQ=255)
                                                       aGCTAATTTCTCAGAATTATG‑AAAAGCCATCCAGATTTGGATGG‑TTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAAT‑ATTCATTttatta                                          <  1:100037/100‑1 (MQ=255)
                                                       aGCTAATTTCTCAGAATTATG‑AAAAGCCATCCAGATTTGGATGG‑TTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAAT‑ATTCATTttatt                                           >  2:953861/1‑99 (MQ=255)
                                                        gCTAATTTCTCAGAATTATG‑AAAAGCCATCCAGATTTGGATGG‑TTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAAT‑ATTCATTTTATTaa                                         >  1:2885478/1‑100 (MQ=255)
                                                          tAATTTCTCAGAATTATG‑AAAAGCCATCCAGATTTGGATGG‑TTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAAT‑ATTCATTttatt                                           <  2:2794362/96‑1 (MQ=255)
                                                            aTTTCTCAGAATTATG‑AAAAGCCATCCAGATTTGGATGG‑TTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAAT‑ATTCATTTTATTAATAtt                                     <  2:2951344/100‑1 (MQ=255)
                                                                ctcAGAATTATG‑AAAAGCCATCCAGATTTGGATGG‑TTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAAT‑ATTCATTTTATTAATATTTaaa                                 >  2:1149989/1‑100 (MQ=255)
                                                                ctcAGAATTATG‑AAAAGCCATCCAGATTTGGATGG‑TTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAAT‑ATTCATTTTATTAATATTTaaa                                 >  1:2471816/1‑100 (MQ=255)
                                                                 tcAGAATTATG‑AAAAGCCATCCAGATTTGGATGG‑TTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAAT‑ATTCATTTTATTAATATTTaaa                                 <  2:1436463/99‑1 (MQ=255)
                                                                 tcAGAATTATG‑AAAAGCCATCCAGATTTGGATGG‑TTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAAT‑ATTCATTTTATTAATATTTAAAc                                >  2:2827564/1‑100 (MQ=255)
                                                                 tcAGAATTATG‑AAAAGCCATCCAGATTTGGATGG‑TTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAAT‑ATTCATTTTATTAATATTTAAAc                                <  2:1770505/100‑1 (MQ=255)
                                                                  cAGAATTATG‑AAAAGCCATCCAGATTTGGATGG‑TTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAAT‑ATTCATTTTATTAATATTTAAAc                                <  1:94491/99‑1 (MQ=255)
                                                                  cAGAATTATG‑AAAAGCCATCCAGATTTGGATGG‑TTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAAT‑ATTCATTTTATTAATATTTAAACt                               >  2:1471706/1‑100 (MQ=255)
                                                                     aaTTATG‑AAAAGCCATCCAGATTTGGATGG‑TTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAAT‑ATTCATTTTATTAATATTTAAACTTaa                            <  2:1892764/100‑1 (MQ=255)
                                                                      agacaaaaaaaACCATCCAAATCTGGATGG‑CTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAAT‑ATTCATTTTATTAATATTTAAACTTaa                            >  2:46959/8‑100 (MQ=25)
                                                                      aTTATG‑AAAAGCCATCCAGATTTGGATGG‑TTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAAT‑ATTCATTTTATTAATATTTAAACTTAAt                           >  2:1933644/1‑100 (MQ=255)
                                                                       ttATG‑AAAAGCCATCCAGATTTGGATGG‑TTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAAT‑ATTCATTTTATTAATATTTAAACTTAATc                          >  2:1264596/1‑100 (MQ=255)
                                                                       ttATG‑AAAAGCCATCCAGATTTGGATGG‑TTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAAT‑ATTCATTTTATTAATATTTAAACTTAATc                          >  1:1423947/1‑100 (MQ=255)
                                                                        tATG‑AAAAGCCATCCAGATTTGGATGG‑TTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAAT‑ATTCATTTTATTAATATTTAAACTTAATCt                         <  2:117513/100‑1 (MQ=255)
                                                                         caaaaaaaACCATCCAAATCTGGATGG‑TTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAAT‑ATTCATTTTATTAATATTTAAACTTAATCt                         >  2:1979962/5‑100 (MQ=38)
                                                                         aTG‑AAAAGCCATCCAGATTTGGATGG‑TTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAAT‑ATTCATTTTATTAATATTTAAACTTAATCta                        <  2:1754192/100‑1 (MQ=255)
                                                                           g‑AAAAGCCATCCAGATTTGGATGG‑TTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAAT‑ATTCATTTTATTAATATTTAAACTTAATCtatn                      <  1:686232/100‑2 (MQ=255)
                                                                           g‑AAAAGCCATCCAGATTTGGATGG‑TTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAAT‑ATTCATTTTATTAATATTTAAACTTAATCtata                      >  1:1522780/1‑100 (MQ=255)
                                                                           aaaaaaACCATCCAAATCTGGATGG‑TTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAAT‑ATTCATTTTATTAATATTTAAACTTAATCtat                       >  2:479716/3‑100 (MQ=38)
                                                                             aaaaGCCATCCAGATTTGGATGG‑TTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAAT‑ATTCATTTTATTAATATTTAAACTTAATCTatat                     <  2:1096003/100‑1 (MQ=255)
                                                                             aaaaGCCATCCAGATTTGGATGG‑TTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAAT‑ATTCATTTTATTAATATTTAAACTTAATCTatat                     <  1:1471078/100‑1 (MQ=255)
                                                                               aaaccATCCAAATCTGGATGG‑TTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAAT‑ATTCATTTTATTAATATTTAAACTTAATCTATATgg                   >  1:416467/4‑100 (MQ=255)
                                                                                aGCCATCCAGATTTGGATGG‑TTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAAT‑ATTCATTTTATTAATATTTAAACTTAATCTATATgg                   <  2:2207320/99‑1 (MQ=255)
                                                                                 gCCATCCAGATTTGGATGG‑TTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAAT‑ATTCATTTTATTAATATTTAATCTTAATCTATATGGaa                 >  2:335709/1‑100 (MQ=255)
                                                                                  ccATCCAGATTTGGATGG‑TTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAAT‑ATTCATTTTATTAATATTTAAACTTAATCTATATGGaat                <  1:1675/100‑1 (MQ=255)
                                                                                     tCCAAATTTGGATGG‑TTTTTTTTGTCTATAGTTGGTAAGATAATTACGTATTGCAAAT‑ATTCATTTTATTAATATTTAAACTTAATCTATATGGaataat             <  2:738224/100‑1 (MQ=255)
                                                                                      ccAGATTTGGATGGTTTTTTTTT‑TTTTTTTTTTGTAAGATAATTACGTATTGCAAAT‑ATTCATTTTATTAATATTTAAACTTAATCTATATGGAataata            >  2:679359/1‑100 (MQ=255)
                                                                                      ccAGATTTGGATGG‑TTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAAT‑ATTCATTTTATTAATATTTAAACTTAATCTATATGGAataata            <  1:1167331/100‑1 (MQ=255)
                                                                                       caaatctGGATGGTTTTTTTTT‑TTTTTAGCTGGTAAGATAATTACGTATTGCAAAT‑ATTCATTTTATTAATATTTAATCTTAATCTATATGGaataa              >  2:1290504/7‑97 (MQ=255)
                                                                                       cAGATTTGGATGG‑TTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAAT‑ATTCATTTTATTAATATTTAAACTTAATCTATATGGAataata            <  2:3011585/99‑1 (MQ=255)
                                                                                       cAGATTTGGATGG‑TTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAAT‑ATTCATTTTATTAATATTTAAACTTAATCTATATGGAATAatat           >  1:1038152/1‑100 (MQ=255)
                                                                                        aGATTTGGATGG‑TTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAAT‑ATTCATTTTATTAATATTTAAACTTAATCTATATGGAATAatata          <  2:789487/100‑1 (MQ=255)
                                                                                         gATTTGGATGG‑TTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAAT‑ATTCATTTTATTAATATTTAAACTTAATCTATATGGAATAatatat         <  1:2820478/100‑1 (MQ=255)
                                                                                          aTTTGGATGG‑TTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAAT‑ATTCATTTTATTAATATTTAAACTTAATCTATATGGAATAATATATg        <  1:1022419/100‑1 (MQ=255)
                                                                                               gATGG‑TTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAAT‑ATTCATTTTATTAATATTTAAACTTAATCTATATGGAATAATATATGACATg   <  2:485424/100‑1 (MQ=255)
                                                                                               gATGG‑TTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAAT‑ATTCATTTTATTAATATTTAAACTTAATCTATATGGAATAATATATGACATg   <  1:2053455/100‑1 (MQ=255)
                                                                                                    |                                                                                                 
CCTGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATG‑AAAAGCCATCCAGATTTGGATGGTTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAAT‑ATTCATTTTATTAATATTTAAACTTAATCTATATGGAATAATATATGACATGG  >  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: ‑