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.93.0% 617.6 / 37.4 156intergenic (+39/‑88)rlmB/yjfI23S rRNA mG2251 2'‑O‑ribose methyltransferase, SAM‑dependent/DUF2170 family protein
Reads supporting (aligned to +/- strand):  ref base T (8/2);  new base . (67/78);  total (75/81)
Fisher's exact test for biased strand distribution p-value = 5.08e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.23e-01

Only 100 of 187 total aligned reads displayed.

CCTGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGAATAATATATGACAT  >  NC_000913/4409946‑4410139
                                                                                                   |                                                                                              
ccTGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGt                                                                                                >  1:2908719/1‑100 (MQ=255)
ccTGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGt                                                                                                >  1:543395/1‑100 (MQ=255)
  tGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGGATTATGAAAAGCCATCCAGATTTGGATGGtct                                                                                              >  2:1323723/1‑98 (MQ=255)
  tGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGttt                                                                                              >  2:157228/1‑100 (MQ=255)
  tGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGttt                                                                                              <  1:2127082/100‑1 (MQ=255)
   gAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGtttt                                                                                             <  2:2494246/100‑1 (MQ=255)
   gAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAACTTCTCAGAATTATGAAAAACCATCCAAATCTGGATGGcttt                                                                                             <  2:3292167/100‑1 (MQ=38)
     aCGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGtttttt                                                                                           <  1:3706957/100‑1 (MQ=255)
     aCGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGtttttt                                                                                           <  2:3560614/100‑1 (MQ=255)
      cGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAGAGCCATCCAGATTTGGATGGttttttt                                                                                          >  1:786108/1‑100 (MQ=255)
      cGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGttttttt                                                                                          <  1:35071/100‑1 (MQ=255)
      cGTATCGGTTGCCACCGTAATTTGCTTATTTGAAGCCGTGCGCCAGCGCAGCTAATTTCTCCGTATTATGAAAAGCCATCCAGATTTGGATCGatttttt                                                                                          <  1:154557/100‑1 (MQ=255)
       gTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTGATTTCTCAGAATTATGAAAAACCATCCAAATCTGGATGGcttttcat                                                                                         <  1:2248047/100‑4 (MQ=37)
       gTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGtttttttt                                                                                         >  2:2316039/1‑100 (MQ=255)
        tttCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGttttttttg                                                                                        <  2:2414390/100‑2 (MQ=255)
        tttCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGttttttttg                                                                                        <  2:3485456/100‑2 (MQ=255)
        tttCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGttttttttg                                                                                        <  2:66537/100‑2 (MQ=255)
        tttCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGttttttt                                                                                          >  2:1840823/1‑98 (MQ=255)
          tCGGGTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGttttttttgtc                                                                                      <  1:2152697/100‑4 (MQ=255)
            ggTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTGTCta                                                                                   <  2:3537070/100‑1 (MQ=255)
                 cGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTGTCTATAGCt                                                                              >  2:795010/1‑100 (MQ=255)
                  gACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTGTCTATAGc                                                                               <  1:1492545/98‑1 (MQ=255)
                    ccGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTGTCTATAGCTGGt                                                                           >  2:1978786/1‑100 (MQ=255)
                     cGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTGTCTATAGCTGGTa                                                                          <  1:1956643/100‑1 (MQ=255)
                     cGGAATTTGCTTATTTGAAGCGGTACGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTGTCTATAGCTGGt                                                                           <  1:1087746/99‑1 (MQ=255)
                      ggAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTGTCTATATCTGGTaa                                                                         >  1:1127104/1‑100 (MQ=255)
                      ggAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTGTCTATAGCTGGTaa                                                                         >  2:1354572/1‑100 (MQ=255)
                      ggAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTGTCTATAGCTGGTaa                                                                         >  2:3431290/1‑100 (MQ=255)
                      ggAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTGTCTATAGCTGGTaa                                                                         >  2:1211323/1‑100 (MQ=255)
                       taaTTTGTTTATTTGTTGTGGTGCGCCAGCGCTGCTAATTTCTCTGTATTATGATAAGCCATTCAGATTTGGATGG‑TTTTTTTTGTTTATAGTTGGTa                                                                          <  1:1713789/97‑1 (MQ=255)
                           ttGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTGTCTATAGCTGGTAAGATaa                                                                    <  2:2910710/100‑1 (MQ=255)
                            tGCTTATTTGAAGCGGTGCGCCAGAGCAGCTAATTTCTCAGAATTATGAAAAGGCATCCTGAGGTGGATGG‑TTCTTTTTGTCTATAGc                                                                               >  1:3419299/1‑88 (MQ=255)
                             gCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTGTCTATAGCTGGTAAGATAAtt                                                                  >  2:658300/1‑100 (MQ=255)
                              cTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTGTCTATAGCTGGTAAGATAATTa                                                                 >  1:3192273/1‑100 (MQ=255)
                              cTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTGTCTATAGCTGGTAAGATAATTa                                                                 >  2:602269/1‑100 (MQ=255)
                              cTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTGTCTATAGCTGGTAAGATAATTa                                                                 >  2:3152387/1‑100 (MQ=255)
                                tATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTGTCTATAGCTGTTAAGATAATTAc                                                                >  1:2417093/1‑99 (MQ=255)
                                tATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTGTCTATAGCTGGTAAGATAATTACg                                                               <  2:525056/100‑1 (MQ=255)
                                tATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTGTCTATAGCTGGTAAGATAATTACg                                                               <  1:711666/100‑1 (MQ=255)
                                tATTAGAGGCGGTGAGCTTGTGCAGATAATTTCTCAGAATTATGAAGGCCTATCCAGATTTGG‑TTGTTTTTTTTTTTCTTTAGTTGGTAAGATAAGTAca                                                               <  1:515970/100‑2 (MQ=255)
                                  tttGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTGTCTATAGCTGGTAAGATAATTACGTa                                                             >  2:2807540/1‑100 (MQ=255)
                                  tttGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTGTCTATAGCTGGTAAGATAATTACGTa                                                             >  2:430464/1‑100 (MQ=255)
                                   ttGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTGTCTATAGCTGGTAAGATAATTACGTAt                                                            >  2:477798/1‑100 (MQ=255)
                                   ttGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTGTCTATAGCTGGTAAGATAATTACGTAt                                                            >  2:103176/1‑100 (MQ=255)
                                   ttGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTGTCTATAGCTGGTAAGATAATTACGTAt                                                            <  2:840756/100‑1 (MQ=255)
                                     gAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTg                                                          >  2:3231285/1‑100 (MQ=255)
                                      aaGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGc                                                         <  2:3406520/100‑1 (MQ=255)
                                      aaGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGc                                                         <  2:3453115/100‑1 (MQ=255)
                                       aGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTGTCTATATCTGGTAAGATAATTACGTATTg                                                          >  1:2128138/1‑98 (MQ=255)
                                         cGGTGCGCCAGCGCTGCTAATTTGTGAGAATTATGAAAAGACATCAAGATTTGGATGG‑TTTTTTTTGTCTATAGCTGGTAAGATAATTACGGATTGCaaa                                                      <  1:2590418/100‑1 (MQ=255)
                                         cGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCaaa                                                      <  2:3688623/100‑1 (MQ=255)
                                           gTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAata                                                    <  2:1722860/100‑1 (MQ=255)
                                            tGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAatat                                                   <  2:3209044/100‑1 (MQ=255)
                                              cgcCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTc                                                 <  2:233862/100‑1 (MQ=255)
                                                  aGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCAttt                                             >  2:570391/1‑100 (MQ=255)
                                                  aGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCAttt                                             >  2:2825451/1‑100 (MQ=255)
                                                   gcgcAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCAtttt                                            <  2:3335713/100‑1 (MQ=255)
                                                     gcAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTttat                                          <  2:3501338/100‑1 (MQ=255)
                                                     gcAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTttat                                          <  1:1786372/100‑1 (MQ=255)
                                                      cAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTGTCTATAGCTGGTTAGATAATTACGTATTGCAAATATTCATTttatt                                         <  1:2473862/100‑1 (MQ=255)
                                                        gCTATTTTCTCAGTATTATGATCAGCTATCCACATTTGGATCG‑TTGTTTTTGTCTTTAGCTCGTAAGATAATTACGTATTGCAAATATTCATTttat                                          <  1:2753632/97‑1 (MQ=255)
                                                        gCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTaa                                       <  2:2514849/100‑1 (MQ=255)
                                                          tAATTTCTCAGAATTATGAAAAGCTATCCAGATTTGGATAG‑TTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTttatt                                         <  1:2845083/96‑1 (MQ=255)
                                                            aTTTGTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATAtt                                   <  2:292029/100‑1 (MQ=255)
                                                             tttCTCAGAATTATGAAAAAAAATGAAGATTTGGATGG‑ATTTTTTTATCTATAGCTAGTAAGATAAATACGTATTGTAAATATTTATTTTATTAatat                                    >  1:253055/1‑98 (MQ=255)
                                                              ttCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTa                                 <  1:2314324/100‑1 (MQ=255)
                                                               tctcAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTGTCTATAGCTGGTAAGATAATTCCGTATTGCAAATATTCATTTTATTAATAttt                                  >  1:1452933/1‑98 (MQ=255)
                                                                ctcAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTT‑TTTTTTTTTGTTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTaa                                >  2:1085892/1‑99 (MQ=255)
                                                                 tcAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAAc                              <  1:660047/100‑1 (MQ=255)
                                                                 tcAGAATTATGAAAAGCCATCCAGAGTTGGATGG‑TTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTaaa                               <  1:222151/99‑1 (MQ=255)
                                                                  cAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTGTCTATAGTTTGTAAGATAATTACGTATTGAAAATATTCATTTTATTAATAGTTa                                 >  1:2054136/1‑96 (MQ=255)
                                                                  cAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACt                             >  2:92822/1‑100 (MQ=255)
                                                                  cAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACt                             >  2:2587201/1‑100 (MQ=255)
                                                                    gAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAATCTTa                           >  2:1382332/1‑100 (MQ=255)
                                                                     aaTTATGAAAAGACATCGAGATTTGGATAG‑TTTTTTTTATCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACt                             <  1:2099440/97‑1 (MQ=255)
                                                                      aTTATGAAAAGACATCCAGATTTAGATGG‑TTTTTTTTGTCTATAACTGGTAAGATAATTAAGTATTGCAAATATTCATTTTATTAATATTTAAACTTAAt                         <  1:497203/100‑1 (MQ=255)
                                                                       ttATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATc                        >  2:1660195/1‑100 (MQ=255)
                                                                       acaaaaaaaACCATCCAAATCTGGATGG‑TTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATc                        >  1:61252/6‑100 (MQ=37)
                                                                         aTGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCta                      <  2:2852564/100‑1 (MQ=255)
                                                                         aTGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCta                      >  2:3276994/1‑100 (MQ=255)
                                                                          aaaaaaACCATCCAAATCTGGATGG‑TTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCtat                     >  2:1939593/3‑100 (MQ=38)
                                                                           gAAAAGCCATCCAGATTTGGATGGTTTTTTTTT‑TTTTTTTTTTTTAAGATAATTACTTATTTCATTTATTCTTTTTTTTAATATTTAAACTTAATCtata                    >  1:3565144/1‑100 (MQ=255)
                                                                            aaaaGCCATCCAGATTTGGATGG‑TTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTatat                   >  1:1602419/1‑100 (MQ=255)
                                                                             aaaaCCATCCAAATCTGGATGG‑TTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTAATTTTATTAATATTTAAACTTAATCTatat                   >  1:2650186/1‑99 (MQ=255)
                                                                             aaaGCCATCCAGATTTGGATGG‑TTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATg                  <  1:3069102/100‑1 (MQ=255)
                                                                               aGCCATCCAGATTTGGATGG‑TTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGa                <  2:3419324/100‑1 (MQ=255)
                                                                                gCCATCCAGATTTTGATGG‑TTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTatat                   <  1:2568028/96‑1 (MQ=255)
                                                                                 ccATCCAGATTTGGATGG‑TTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAACATTTAAACTTAATCTATATGGaat              >  1:1873003/1‑100 (MQ=255)
                                                                                  cATCCAGATTTGGATGGTTTTTTTTT‑TTTTTTTTTTTTAAGATAATTACGTATTGCAAATTTTCATTTTATTGATATTTAAACTTAATTTATATggttt              >  1:2988882/1‑96 (MQ=255)
                                                                                   aTCCAGATTTGGATGG‑TTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGaataa            >  2:2979824/1‑100 (MQ=255)
                                                                                    tCCAGATTTGGATGGTTTTTTTTT‑TTTTTTTTTGGTAAGATAATTACGTATTGCAATTATTTATTTTATTAATATTTAAACTTAATCTATATTGaataat           >  1:826456/1‑100 (MQ=255)
                                                                                     ccAGATTTGGATGGTTTTTTTTT‑TCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGAataata          >  1:1333046/1‑100 (MQ=255)
                                                                                     ccAGATTTGGATGG‑TTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGAataata          >  2:1088135/1‑100 (MQ=255)
                                                                                      caaatctGGATGG‑TTTTTTTTGTCTTTAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGAATAatat         >  2:2191635/7‑100 (MQ=255)
                                                                                       aGATTTGGATGG‑TTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGAATAatata        <  1:1194883/100‑1 (MQ=255)
                                                                                        gATTTGGATGG‑TTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGAATAatatat       <  2:255844/100‑1 (MQ=255)
                                                                                         aTTTGGATGG‑TTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGAATAATATATg      <  2:2327285/100‑1 (MQ=255)
                                                                                          tttGGATGG‑TTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGAATAATATATGa     <  2:1602426/100‑1 (MQ=255)
                                                                                          tttGGATGG‑TTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGAATAATATATGa     >  2:2659181/1‑100 (MQ=255)
                                                                                            tGGATGG‑TTTTTTTTGTCTATAGCTGGGAAGATAATTACATATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGAataata          <  1:1054853/93‑1 (MQ=255)
                                                                                                   |                                                                                              
CCTGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGAATAATATATGACAT  >  NC_000913/4409946‑4410139

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: ‑