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.

CTTCTAAATTAATTTATGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTTTATTTATCGTCGAACCTGGATTGTTT  >  NC_000913/4261260‑4261457
                                                                                                   |                                                                                                   
cTTCTAAATTAATTTATGATTTGTTCGGGGTAAAATCGGAGAAATAATATGTGAAGTTCTACGATAGTATGGAAGGCATTATATAAAAGGACCCAAtatt                                                                                                    <  1:1387365/100‑1 (MQ=255)
     aaaTTAATTTATGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTg                                                                                               <  1:1038170/100‑1 (MQ=255)
      aattaattTATGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTgg                                                                                              <  2:1983942/100‑1 (MQ=255)
       attaattTATGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTggg                                                                                             >  2:1758350/1‑100 (MQ=255)
          aattTATGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGAACCCAATAAATATTGGGtcc                                                                                          <  2:687256/100‑3 (MQ=39)
           attTATGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCt                                                                                         <  1:607752/100‑1 (MQ=255)
           attTATGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATAAATATTGGGTCCt                                                                                         <  2:3603825/100‑1 (MQ=255)
           attTATGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGAACCCAATAAATATTGGGTCCt                                                                                         <  2:195297/100‑1 (MQ=38)
              tATGGTTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCtttt                                                                                      <  1:240231/100‑1 (MQ=255)
                tGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCtttttt                                                                                    <  2:3799504/100‑1 (MQ=255)
                 gATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCTATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTtc                                                                                   <  1:1144095/100‑1 (MQ=255)
                 gATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTtc                                                                                   <  1:660810/100‑1 (MQ=255)
                 gATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTtc                                                                                   <  2:2138617/100‑1 (MQ=255)
                 gATTTGGTCGGAGTTAAATCGGAGAAATAATATGTGAAGTTTGGCGATAGTATGTAAGGCATTATATAAGAGGAGGCTATAGTAATTGAATTCttttt                                                                                     <  1:3268733/98‑1 (MQ=255)
                   tttGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTtctc                                                                                 <  2:840665/100‑1 (MQ=255)
                   tttGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTtctc                                                                                 <  2:1018907/100‑1 (MQ=255)
                   tttGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTtctc                                                                                 <  2:1687325/100‑1 (MQ=255)
                    ttGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGAACCCAATAAATATTGGGTCCttttatataa                                                                                <  2:3533432/100‑7 (MQ=25)
                     tGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTa                                                                               <  1:1223349/100‑1 (MQ=255)
                     tGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTa                                                                               <  2:3844237/100‑1 (MQ=255)
                      gTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTAt                                                                              <  2:2327928/100‑1 (MQ=255)
                      gTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTAt                                                                              >  2:459722/1‑100 (MQ=255)
                      gTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTAt                                                                              >  2:2829015/1‑100 (MQ=255)
                       ttCGGTGTTAAATCGGAGAAATAATATGTAAAGGTATTCGATAGTATGGAAGTCATTATATAAAAGAACCCAATAAATATTGGGTCCttttatataatg                                                                              <  1:1861033/99‑9 (MQ=18)
                       ttCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATc                                                                             <  2:2434754/100‑1 (MQ=255)
                        tCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAAGAACCCAATAAATATTGGGTCCttttatataatgc                                                                             <  1:3452187/100‑10 (MQ=18)
                        tCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGAACCCAATAAATATTGGGTCCttttatataatgcc                                                                            <  1:3809902/100‑11 (MQ=21)
                           gggTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATa                                                                         >  1:3057428/1‑100 (MQ=255)
                           gggTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATa                                                                         <  2:3210406/100‑1 (MQ=255)
                            ggTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATAc                                                                        <  1:80365/100‑1 (MQ=255)
                            ggTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATAc                                                                        >  2:3932481/1‑100 (MQ=255)
                              ttAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTa                                                                      <  1:3997480/100‑1 (MQ=255)
                               tAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTAt                                                                     <  2:1625821/100‑1 (MQ=255)
                                 aaTCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTa                                                                   >  1:3740934/1‑100 (MQ=255)
                                  aTCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAg                                                                  >  2:3067505/1‑100 (MQ=255)
                                   tCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAg                                                                  >  2:881247/1‑99 (MQ=255)
                                    cGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCa                                                                >  1:517991/1‑100 (MQ=255)
                                    cGGAGAAATAATATGTGAAGTTATTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCa                                                                >  1:2449603/1‑100 (MQ=255)
                                     ggAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAg                                                               <  2:4070836/100‑1 (MQ=255)
                                      gagaAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAgg                                                              <  2:1778312/100‑1 (MQ=255)
                                      gagaAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAgg                                                              >  2:168044/1‑100 (MQ=255)
                                          aataatATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGgaga                                                          <  1:4043245/100‑1 (MQ=255)
                                          aataatATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGgaga                                                          <  2:3199360/100‑1 (MQ=255)
                                           ataataTGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGat                                                         >  1:284280/1‑100 (MQ=255)
                                            taataTGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGata                                                        <  2:2024437/100‑1 (MQ=255)
                                            taataTGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGata                                                        <  2:1963289/100‑1 (MQ=255)
                                            taataTGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGata                                                        >  2:2511996/1‑100 (MQ=255)
                                            taataTGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGata                                                        <  2:793704/100‑1 (MQ=255)
                                             aataTGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGatat                                                       <  2:3965149/100‑1 (MQ=255)
                                             aataTGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGatat                                                       >  1:4219992/1‑100 (MQ=255)
                                              atatGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGatata                                                      <  1:168043/100‑1 (MQ=255)
                                              atatGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGatata                                                      <  2:2550986/100‑1 (MQ=255)
                                                 tttgAAGTTCTTCAATAGTATGGAAGGCATTATATAAAAGGAGCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCa                                                   >  2:2983304/3‑100 (MQ=255)
                                                   tgAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCAcc                                                 <  2:3379784/100‑1 (MQ=255)
                                                    gAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCa                                                >  1:4237871/1‑100 (MQ=255)
                                                       gTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCagag                                             >  2:4221390/1‑100 (MQ=255)
                                                          cttcGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGttt                                          <  1:3819383/100‑1 (MQ=255)
                                                          cttcGATAGTATGAAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGttt                                          <  2:3483225/100‑1 (MQ=255)
                                                              gATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGtata                                         >  1:4168743/1‑97 (MQ=255)
                                                               aTAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATCCTATTAGCCGTGAGATATATCACCAGAGTTTAAtgt                                     >  2:56475/1‑100 (MQ=255)
                                                                 aGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAAtgtgt                                   <  1:1804309/100‑1 (MQ=255)
                                                                 aGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAAtgtgt                                   <  1:4058218/100‑1 (MQ=255)
                                                                  gTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAAtgtgtg                                  >  2:2187681/1‑100 (MQ=255)
                                                                  gTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAAtgtgtg                                  >  2:2792284/1‑100 (MQ=255)
                                                                    aTGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGAt                                <  1:2751282/100‑1 (MQ=255)
                                                                    aTGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGAt                                <  2:1501025/100‑1 (MQ=255)
                                                                    aTGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGAGATATCACCAGAGTTTAATGTGTGAt                                <  1:2329762/100‑1 (MQ=255)
                                                                     tGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGAtt                               >  2:1002495/1‑100 (MQ=255)
                                                                       gAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGAtttt                             >  2:1599971/1‑100 (MQ=255)
                                                                        aaGGCATTATTTATCAGGACCCAATATGTATTGGGTTCTTTTTTCTCTATCACTACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGAttttt                            <  1:320650/100‑1 (MQ=255)
                                                                         aGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGAtttttt                           >  2:768181/1‑100 (MQ=255)
                                                                          ggCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGTAGGGAGATATATCACCAGAGTTTAATGTGTGAtttttt                           <  1:2702476/99‑1 (MQ=255)
                                                                          ggCATTATATAAAAGGACCCAAGATTTAGTGGGTTCTTTTTCCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTttta                          <  2:1435896/100‑1 (MQ=255)
                                                                           tgatagagaaaaaaGAACCCAATAAATATTGGGTCCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTttta                          >  1:2499827/11‑100 (MQ=21)
                                                                           gCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTtttat                         <  2:128844/100‑1 (MQ=255)
                                                                             atagagaaaaaaGAACCCAATAAATATTGGGTCCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTtttatt                        >  2:225254/9‑100 (MQ=25)
                                                                             aTTATATAAAAGGACCCAATATATATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCATCAGAGTTTAATGTGTGATGTttgattt                       >  1:1533882/1‑100 (MQ=255)
                                                                              tagagaaaaaaGAACCCAATAAATATTGGGTCCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTtttattt                       >  1:192596/8‑100 (MQ=25)
                                                                                gagaaaaaaGAACCCAATAAATATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGATTTTAATGTGTGATTTTttatttat                     >  1:1094815/6‑100 (MQ=37)
                                                                                atataAAAGGACCCAATATTTATTGGGCTCTTTCTTCTCTATAAAAACTATTAACAGGGCGATATATCACCCGCCTTCAACGTGTGATTTTTTCTTtctc                    >  1:1495080/1‑100 (MQ=255)
                                                                                 tataAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTTTATTTAtc                    >  1:458443/1‑99 (MQ=255)
                                                                                   taAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTTTATTTAtcgtc                 <  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)
                                                                                                   |                                                                                                   
CTTCTAAATTAATTTATGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTTTATTTATCGTCGAACCTGGATTGTTT  >  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: