Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,261,358 T→A 19.2% intergenic (+52/‑311) yjbM → / → dusA uncharacterized protein/tRNA‑dihydrouridine synthase A

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,261,3580TA19.2% 501.0 / 95.7 207intergenic (+52/‑311)yjbM/dusAuncharacterized protein/tRNA‑dihydrouridine synthase A
Reads supporting (aligned to +/- strand):  ref base T (74/92);  new base A (20/20);  total (94/113)
Fisher's exact test for biased strand distribution p-value = 5.97e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 217 total aligned reads displayed.

CTTCTAAATTAATTTATGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATAT‑AAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTTTATTTATCGTCGAACCTGGATTGTTT  >  NC_000913/4261260‑4261457
                                                                                                   |                                                                                                   
cTTCTAAATTAATTTATGATTTGTTCGGGGTAAAATCGGAGAAATAATATGTGAAGTTCTACGATAGTATGGAAGGCATTATAT‑AAAAGGACCCAAtatt                                                                                                    <  1:1387365/100‑1 (MQ=255)
     aaaTTAATTTATGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATAT‑AAAAGGACCCAATATTTATTg                                                                                               <  1:1038170/100‑1 (MQ=255)
      aattaattTATGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATAT‑AAAAGGACCCAATATTTATTgg                                                                                              <  2:1983942/100‑1 (MQ=255)
       attaattTATGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATAT‑AAAAGGACCCAATATTTATTggg                                                                                             >  2:1758350/1‑100 (MQ=255)
          aattTATGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATAT‑AAAAGAACCCAATAAATATTGGGtcc                                                                                          <  2:687256/100‑3 (MQ=39)
           attTATGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATAT‑AAAAGGACCCAATATTTATTGGGTTCt                                                                                         <  1:607752/100‑1 (MQ=255)
           attTATGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATAT‑AAAAGGACCCAATAAATATTGGGTCCt                                                                                         <  2:3603825/100‑1 (MQ=255)
           attTATGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATAT‑AAAAGAACCCAATAAATATTGGGTCCt                                                                                         <  2:195297/100‑1 (MQ=38)
              tATGGTTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATAT‑AAAAGGACCCAATATTTATTGGGTTCtttt                                                                                      <  1:240231/100‑1 (MQ=255)
                tGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATAT‑AAAAGGACCCAATATTTATTGGGTTCtttttt                                                                                    <  2:3799504/100‑1 (MQ=255)
                 gATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCTATAGTATGGAAGGCATTATAT‑AAAAGGACCCAATATTTATTGGGTTCTTTTTtc                                                                                   <  1:1144095/100‑1 (MQ=255)
                 gATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATAT‑AAAAGGACCCAATATTTATTGGGTTCTTTTTtc                                                                                   <  1:660810/100‑1 (MQ=255)
                 gATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATAT‑AAAAGGACCCAATATTTATTGGGTTCTTTTTtc                                                                                   <  2:2138617/100‑1 (MQ=255)
                 gATTTGGTCGGAGTTAAATCGGAGAAATAATATGTGAAGTTTGGCGATAGTATGTAAGGCATTATAT‑AAGAGGAGGCTATAGTAATTGAATTCttttt                                                                                     <  1:3268733/98‑1 (MQ=255)
                   tttGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATAT‑AAAAGGACCCAATATTTATTGGGTTCTTTTTtctc                                                                                 <  2:840665/100‑1 (MQ=255)
                   tttGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATAT‑AAAAGGACCCAATATTTATTGGGTTCTTTTTtctc                                                                                 <  2:1018907/100‑1 (MQ=255)
                   tttGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATAT‑AAAAGGACCCAATATTTATTGGGTTCTTTTTtctc                                                                                 <  2:1687325/100‑1 (MQ=255)
                    ttGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATAT‑AAAAGAACCCAATAAATATTGGGTCCttttatataa                                                                                <  2:3533432/100‑7 (MQ=25)
                     tGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATAT‑AAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTa                                                                               <  1:1223349/100‑1 (MQ=255)
                     tGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATAT‑AAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTa                                                                               <  2:3844237/100‑1 (MQ=255)
                      gTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATAT‑AAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTAt                                                                              <  2:2327928/100‑1 (MQ=255)
                      gTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATAT‑AAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTAt                                                                              >  2:459722/1‑100 (MQ=255)
                      gTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATAT‑AAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTAt                                                                              >  2:2829015/1‑100 (MQ=255)
                       ttCGGTGTTAAATCGGAGAAATAATATGTAAAGGTATTCGATAGTATGGAAGTCATTATAT‑AAAAGAACCCAATAAATATTGGGTCCttttatataatg                                                                              <  1:1861033/99‑9 (MQ=18)
                       ttCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATAT‑AAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATc                                                                             <  2:2434754/100‑1 (MQ=255)
                        tCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAAGAACCCAATAAATATTGGGTCCttttatataatgc                                                                             <  1:3452187/100‑10 (MQ=18)
                        tCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATAT‑AAAAGAACCCAATAAATATTGGGTCCttttatataatgcc                                                                            <  1:3809902/100‑11 (MQ=21)
                           gggTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATAT‑AAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATa                                                                         >  1:3057428/1‑100 (MQ=255)
                           gggTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATAT‑AAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATa                                                                         <  2:3210406/100‑1 (MQ=255)
                            ggTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATAT‑AAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATAc                                                                        <  1:80365/100‑1 (MQ=255)
                            ggTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATAT‑AAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATAc                                                                        >  2:3932481/1‑100 (MQ=255)
                              ttAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATAT‑AAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTa                                                                      <  1:3997480/100‑1 (MQ=255)
                               tAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATAT‑AAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTAt                                                                     <  2:1625821/100‑1 (MQ=255)
                                 aaTCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATAT‑AAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTa                                                                   >  1:3740934/1‑100 (MQ=255)
                                  aTCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATAT‑AAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAg                                                                  >  2:3067505/1‑100 (MQ=255)
                                   tCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATAT‑AAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAg                                                                  >  2:881247/1‑99 (MQ=255)
                                    cGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATAT‑AAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCa                                                                >  1:517991/1‑100 (MQ=255)
                                    cGGAGAAATAATATGTGAAGTTATTCGATAGTATGGAAGGCATTATAT‑AAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCa                                                                >  1:2449603/1‑100 (MQ=255)
                                     ggAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATAT‑AAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAg                                                               <  2:4070836/100‑1 (MQ=255)
                                      gagaAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATAT‑AAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAgg                                                              <  2:1778312/100‑1 (MQ=255)
                                      gagaAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATAT‑AAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAgg                                                              >  2:168044/1‑100 (MQ=255)
                                          aataatATGTGAAGTTCTTCGATAGTATGGAAGGCATTATAT‑AAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGgaga                                                          <  1:4043245/100‑1 (MQ=255)
                                          aataatATGTGAAGTTCTTCGATAGTATGGAAGGCATTATAT‑AAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGgaga                                                          <  2:3199360/100‑1 (MQ=255)
                                           ataataTGTGAAGTTCTTCGATAGTATGGAAGGCATTATAT‑AAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGat                                                         >  1:284280/1‑100 (MQ=255)
                                            taataTGTGAAGTTCTTCGATAGTATGGAAGGCATTATAT‑AAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGata                                                        <  2:2024437/100‑1 (MQ=255)
                                            taataTGTGAAGTTCTTCGATAGTATGGAAGGCATTATAT‑AAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGata                                                        <  2:1963289/100‑1 (MQ=255)
                                            taataTGTGAAGTTCTTCGATAGTATGGAAGGCATTATAT‑AAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGata                                                        >  2:2511996/1‑100 (MQ=255)
                                            taataTGTGAAGTTCTTCGATAGTATGGAAGGCATTATAT‑AAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGata                                                        <  2:793704/100‑1 (MQ=255)
                                             aataTGTGAAGTTCTTCGATAGTATGGAAGGCATTATAT‑AAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGatat                                                       <  2:3965149/100‑1 (MQ=255)
                                             aataTGTGAAGTTCTTCGATAGTATGGAAGGCATTATAT‑AAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGatat                                                       >  1:4219992/1‑100 (MQ=255)
                                              atatGTGAAGTTCTTCGATAGTATGGAAGGCATTATAT‑AAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGatata                                                      <  1:168043/100‑1 (MQ=255)
                                              atatGTGAAGTTCTTCGATAGTATGGAAGGCATTATAT‑AAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGatata                                                      <  2:2550986/100‑1 (MQ=255)
                                                 tttgAAGTTCTTCAATAGTATGGAAGGCATTATAT‑AAAAGGAGCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCa                                                   >  2:2983304/3‑100 (MQ=255)
                                                   tgAAGTTCTTCGATAGTATGGAAGGCATTATAT‑AAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCAcc                                                 <  2:3379784/100‑1 (MQ=255)
                                                    gAAGTTCTTCGATAGTATGGAAGGCATTATAT‑AAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCa                                                >  1:4237871/1‑100 (MQ=255)
                                                       gTTCTTCGATAGTATGGAAGGCATTATAT‑AAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCagag                                             >  2:4221390/1‑100 (MQ=255)
                                                          cttcGATAGTATGGAAGGCATTATAT‑AAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGttt                                          <  1:3819383/100‑1 (MQ=255)
                                                          cttcGATAGTATGAAAGGCATTATAT‑AAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGttt                                          <  2:3483225/100‑1 (MQ=255)
                                                              gATAGTATGGAAGGCATTATAT‑AAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGtata                                         >  1:4168743/1‑97 (MQ=255)
                                                               aTAGTATGGAAGGCATTATAT‑AAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATCCTATTAGCCGTGAGATATATCACCAGAGTTTAAtgt                                     >  2:56475/1‑100 (MQ=255)
                                                                 aGTATGGAAGGCATTATAT‑AAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAAtgtgt                                   <  1:1804309/100‑1 (MQ=255)
                                                                 aGTATGGAAGGCATTATAT‑AAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAAtgtgt                                   <  1:4058218/100‑1 (MQ=255)
                                                                  gTATGGAAGGCATTATAT‑AAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAAtgtgtg                                  >  2:2187681/1‑100 (MQ=255)
                                                                  gTATGGAAGGCATTATAT‑AAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAAtgtgtg                                  >  2:2792284/1‑100 (MQ=255)
                                                                    aTGGAAGGCATTATAT‑AAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGAt                                <  1:2751282/100‑1 (MQ=255)
                                                                    aTGGAAGGCATTATAT‑AAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGAt                                <  2:1501025/100‑1 (MQ=255)
                                                                    aTGGAAGGCATTATAT‑AAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGAGATATCACCAGAGTTTAATGTGTGAt                                <  1:2329762/100‑1 (MQ=255)
                                                                     tGGAAGGCATTATAT‑AAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGAtt                               >  2:1002495/1‑100 (MQ=255)
                                                                       gAAGGCATTATAT‑AAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGAtttt                             >  2:1599971/1‑100 (MQ=255)
                                                                        aaGGCATTATTT‑ATCAGGACCCAATATGTATTGGGTTCTTTTTTCTCTATCACTACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGAttttt                            <  1:320650/100‑1 (MQ=255)
                                                                         aGGCATTATAT‑AAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGAtttttt                           >  2:768181/1‑100 (MQ=255)
                                                                          ggCATTATAT‑AAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGTAGGGAGATATATCACCAGAGTTTAATGTGTGAtttttt                           <  1:2702476/99‑1 (MQ=255)
                                                                          ggCATTATAT‑AAAAGGACCCAAGATTTAGTGGGTTCTTTTTCCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTttta                          <  2:1435896/100‑1 (MQ=255)
                                                                           tgatagagaaaaaaGAACCCAATAAATATTGGGTCCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTttta                          >  1:2499827/11‑100 (MQ=21)
                                                                           gCATTATAT‑AAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTtttat                         <  2:128844/100‑1 (MQ=255)
                                                                             atagagaaaaaaGAACCCAATAAATATTGGGTCCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTtttatt                        >  2:225254/9‑100 (MQ=25)
                                                                             aTTATAT‑AAAAGGACCCAATATATATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCATCAGAGTTTAATGTGTGATGTttgattt                       >  1:1533882/1‑100 (MQ=255)
                                                                              tagagaaaaaaGAACCCAATAAATATTGGGTCCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTtttattt                       >  1:192596/8‑100 (MQ=25)
                                                                                gagaaaaaaGAACCCAATAAATATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGATTTTAATGTGTGATTTTttatttat                     >  1:1094815/6‑100 (MQ=37)
                                                                                atat‑aAAAGGACCCAATATTTATTGGGCTCTTTCTTCTCTATAAAAACTATTAACAGGGCGATATATCACCCGCCTTCAACGTGTGATTTTTTCTTtctc                    >  1:1495080/1‑100 (MQ=255)
                                                                                 tat‑aAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTTTATTTAtc                    >  1:458443/1‑99 (MQ=255)
                                                                                   t‑aAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTTTATTTAtcgtc                 <  2:3934822/100‑1 (MQ=255)
                                                                                    aaaaaGAACCCAATAAATATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCAGCAGAGTTTAATGTGTGATTTTTTGTTTAtcgtc                 >  2:2594037/2‑100 (MQ=38)
                                                                                    aaaaaGAACCCAATAAAAATTGGGTTGTTTTTTCTCTATCAATGCTTTTAACTGGGAGATATATTAGTAGAGTTGGATGTGAGATTTgttaattagttta                 >  1:4118727/2‑95 (MQ=255)
                                                                                       aaGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTTTATTTATCGTCGaa              <  1:3255552/100‑1 (MQ=255)
                                                                                         ggACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTTTATTTATCGTCGAAcc            >  1:221961/1‑100 (MQ=255)
                                                                                         gaaCCCAATAAATATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTTTATTTATCGTCGAAcc            >  2:3354621/3‑100 (MQ=255)
                                                                                           aCCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTTTATTTATCGTCGAACCTg          <  1:2113552/100‑1 (MQ=255)
                                                                                           aCCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTTTATTTATCGTCGAACCTg          >  2:4153577/1‑100 (MQ=255)
                                                                                            cccAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTTTATTTATCGTCGAACCTgg         <  2:2885465/100‑1 (MQ=255)
                                                                                             ccAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTTTATTTATCGTCGAACCTGGa        <  1:2478593/100‑1 (MQ=255)
                                                                                             ccAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTTTATTTATCGTCGAACCTGGa        >  2:541320/1‑100 (MQ=255)
                                                                                             ccAATAATTATTGAGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGCTATATCACCAGAGATTAATGTGTGATTTTTTATTTATCGTCGAATCt           >  1:1728418/1‑97 (MQ=255)
                                                                                              cAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTTTATTTATCGTCGAACCTg          >  1:1249463/1‑97 (MQ=255)
                                                                                              cAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTTTATTTATCGTCGAACCTGGAt       <  1:1556131/100‑1 (MQ=255)
                                                                                              cAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTTTATTTATCGTCGAACCTGGAt       >  2:2778628/1‑100 (MQ=255)
                                                                                                  atttattGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTTTATTTATCGTCGAACCTGGATTGtt   <  1:3029050/100‑1 (MQ=255)
                                                                                                  atttattGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTTTATTTATCGTCGAACCTGGATTGtt   <  1:1002491/100‑1 (MQ=255)
                                                                                                   tttattGGGTTCTTTTTTCTCTATCAATACTATTAGCTGTGATATATATCACGAGAGTTTAATGTTTGATTTTTTATTTATCGTCTAACCTGGATTGttt  <  1:970174/100‑1 (MQ=255)
                                                                                                   tttattGGGTTCTTTTTTCTCTATCAATACTATGAGCAGGGAGATATATCTCCAGAGTTTAATGTGTTATTTGTTATTTATCGTCGAACCTTTATTGttt  <  1:361572/100‑1 (MQ=255)
                                                                                                   |                                                                                                   
CTTCTAAATTAATTTATGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATAT‑AAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTTTATTTATCGTCGAACCTGGATTGTTT  >  NC_000913/4261260‑4261457

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