Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,261,358 T→A 15.6% 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,3580TA15.6% 267.9 / 42.1 109intergenic (+52/‑311)yjbM/dusAuncharacterized protein/tRNA‑dihydrouridine synthase A
Reads supporting (aligned to +/- strand):  ref base T (51/41);  new base A (8/9);  total (59/50)
Fisher's exact test for biased strand distribution p-value = 6.01e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.87e-01

Only 100 of 124 total aligned reads displayed.

ACTTCTAAATTAATTTATGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTTTATTTATCGTCGAACCTGGATTGTTT  >  NC_000913/4261259‑4261457
                                                                                                   |                                                                                                   
acTTCTAAATTAATTTATGATTTGTTCGGGGTTAAATCGGAGAAATAATCTGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAatat                                                                                                     >  1:2071354/1‑100 (MQ=255)
acTTCTAAATTAATTTATGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAatat                                                                                                     >  2:1209146/1‑100 (MQ=255)
 cTTCTAAATTAATTTATGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAAtatt                                                                                                    >  2:952192/1‑100 (MQ=255)
   tCTAAATTAATTTATGATTTGTTCGGGGTTAAATCGGAGTAATACTATGAGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAAtattta                                                                                                  >  1:1901932/1‑100 (MQ=255)
   tCTAAATTAATTTATGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAATTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAAtattta                                                                                                  <  2:1130800/100‑1 (MQ=255)
    cTAAATTAATTTATGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAAtatttat                                                                                                 >  2:220149/1‑100 (MQ=255)
      aaaTTAATTTATGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTg                                                                                               <  1:635644/100‑1 (MQ=255)
        attaattTATGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTggg                                                                                             <  1:2062386/100‑1 (MQ=255)
        attaattTATGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTggg                                                                                             >  1:518716/1‑100 (MQ=255)
         ttaattTATGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGt                                                                                            >  1:463294/1‑100 (MQ=255)
           aattTATGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTc                                                                                          >  2:1716887/1‑100 (MQ=255)
           aattTATGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTc                                                                                          <  2:784972/100‑1 (MQ=255)
            attTATGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCt                                                                                         >  2:307233/1‑100 (MQ=255)
            attTATGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGAACCCAATAAATATTGGGTCCt                                                                                         <  1:179174/100‑1 (MQ=38)
             tttATGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCtt                                                                                        >  1:800013/1‑100 (MQ=255)
              ttATGATTTATTCGGGGTTAAATCGGAGAAATTATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCATTATTTATTGGGTTCt                                                                                         <  2:1827984/98‑1 (MQ=255)
                 tGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGAACCCAATAAATATTGGGTCCttttat                                                                                    <  2:2305030/100‑3 (MQ=37)
                  gATTTGTTCTTAATTAAATCGGAAAAATAATATATGAAGTTCTTCGATAGTATGGAAGTCATTATATAAAAGGACCCAATATTTATTGGGTTCtt                                                                                        <  2:1474783/95‑1 (MQ=255)
                  gATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTtc                                                                                   >  1:607413/1‑100 (MQ=255)
                  gATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTtc                                                                                   <  2:386298/100‑1 (MQ=255)
                    tttGTTCGGGGTTAAGTCGGAGAAATAATATGTGAAGTTCTTCCATAGTATGGAAGGCATTATATAAAAGGACCCAGTATTTATTGGGTTCTTTTTtctc                                                                                 <  2:1567833/100‑1 (MQ=255)
                    tttGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGCTCTTCGATAGTATGGAAGGCATTATATAAAAGGACACAATATTTATTGGGTTCTTTTTtctc                                                                                 >  2:1369876/1‑100 (MQ=255)
                     ttGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGAACCCAATAAATATTGGGTCCttttatataa                                                                                <  1:2438707/100‑7 (MQ=25)
                       gTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGAACCCAATAAATATTGGGTCCttttatataatg                                                                              <  1:723952/100‑9 (MQ=25)
                        ttCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATc                                                                             <  1:286439/100‑1 (MQ=255)
                           ggggTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAAt                                                                          >  1:2549396/1‑100 (MQ=255)
                           ggggTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAAt                                                                          >  1:1018087/1‑100 (MQ=255)
                           ggggTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAAt                                                                          >  1:1511149/1‑100 (MQ=255)
                            gggTTAATTCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAGGTCATTATATAAAATGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATa                                                                         <  2:1579641/100‑1 (MQ=255)
                             ggTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATAc                                                                        <  2:2505399/100‑1 (MQ=255)
                               ttAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTa                                                                      >  1:627401/1‑100 (MQ=255)
                                 aaaTCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTAtt                                                                    <  2:736530/100‑1 (MQ=255)
                                 aaaTCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTAtt                                                                    <  2:2473777/100‑1 (MQ=255)
                                  aaTCTGAGAAATAATATGTGAAGTTCTTCGATATTGTGGTAGGCATTATATAATAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTa                                                                   <  2:1261751/100‑1 (MQ=255)
                                   aTCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAg                                                                  >  1:863623/1‑100 (MQ=255)
                                    tCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGc                                                                 <  2:2241945/100‑1 (MQ=255)
                                       gagaAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAgg                                                              >  2:471022/1‑100 (MQ=255)
                                         gaAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGga                                                            >  1:52731/1‑100 (MQ=255)
                                         gaAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGga                                                            <  2:2016693/100‑1 (MQ=255)
                                           aataatATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGgaga                                                          <  2:1407996/100‑1 (MQ=255)
                                           aataatATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGgaga                                                          >  2:1945331/1‑100 (MQ=255)
                                             taataTGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGata                                                        <  1:929157/100‑1 (MQ=255)
                                             taataTGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGata                                                        <  1:2502019/100‑1 (MQ=255)
                                              aataTGTGAAGGTTTTCGATAGTATGTAAGTCATTATATAAAAGGACCTAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGatat                                                       <  2:1232646/100‑1 (MQ=255)
                                               atatGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGatata                                                      <  1:279546/100‑1 (MQ=255)
                                               atatGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATCTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGatata                                                      >  1:1226934/1‑100 (MQ=255)
                                                tatGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGatatat                                                     >  1:684240/1‑100 (MQ=255)
                                                 atGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATc                                                    >  2:161562/1‑100 (MQ=255)
                                                   gtgAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCAc                                                  >  2:794701/1‑100 (MQ=255)
                                                     gAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCAcc                                                 >  2:1704745/1‑99 (MQ=255)
                                                     gAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCa                                                <  2:191508/100‑1 (MQ=255)
                                                       aGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCaga                                              <  1:1923392/100‑1 (MQ=255)
                                                        gTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCagag                                             >  2:460652/1‑100 (MQ=255)
                                                        gTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCagag                                             <  2:518709/100‑1 (MQ=255)
                                                            ttcGATAGTATCGAAGGCATTATATAAAAGGACCCAATATTTACTGGGTTCTTTTTCCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGt                                            <  2:140276/97‑1 (MQ=255)
                                                             tcGATAGTATGGAAGGCATTATATAAATGTACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTaa                                        <  2:1073115/100‑1 (MQ=255)
                                                             tcGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTaa                                        >  1:1288551/1‑100 (MQ=255)
                                                              cGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAAt                                       <  2:2235455/100‑1 (MQ=255)
                                                               gATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTATCAGGGAGATATATCACCAGAGTTTAAtg                                      >  2:623930/1‑100 (MQ=255)
                                                               gATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAAtg                                      >  2:1903543/1‑100 (MQ=255)
                                                               gATAGTATGGAAGGCATTATATAAAAGGACCCAAGATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAAtg                                      <  1:2281729/100‑1 (MQ=255)
                                                                 tAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAAtatg                                    >  2:2160371/1‑100 (MQ=255)
                                                                  aGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAAtgtgt                                   >  1:1176105/1‑100 (MQ=255)
                                                                    tATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGa                                 >  1:1460904/1‑100 (MQ=255)
                                                                     aTGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGAt                                >  2:90381/1‑100 (MQ=255)
                                                                        naaGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGAttt                              >  2:141424/2‑99 (MQ=255)
                                                                          aGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGAtttttt                           >  1:2314116/1‑100 (MQ=255)
                                                                            gCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTtttat                         <  2:1270632/100‑1 (MQ=255)
                                                                            gCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTtttat                         <  1:695629/100‑1 (MQ=255)
                                                                              aTTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTtttat                         <  2:2218083/98‑1 (MQ=255)
                                                                               ttATATAAAAGGACCTAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGTTGATATATCACCAGAGTTTAATGTGTGATTTtttatt                        <  2:129001/98‑1 (MQ=255)
                                                                               ttATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTtttattta                      <  1:1363606/100‑1 (MQ=255)
                                                                               ttATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTtttattta                      >  1:1679779/1‑100 (MQ=255)
                                                                               agagaaaaaaGAACCCAATAAATATTGGGTCCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTtttattta                      >  1:1700268/7‑100 (MQ=25)
                                                                                 agaaaaaaGAACCCAATAAATATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTTTATTTAtc                    >  1:879443/5‑100 (MQ=38)
                                                                                   aaaaaaGAACCCAATAAATATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTTTATTTAtcgt                  >  2:307179/3‑100 (MQ=38)
                                                                                    taAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTTTATTTAtcgtc                 <  1:2436671/100‑1 (MQ=255)
                                                                                    taAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATAACCAGAGTTTAATTTGTGATTTTTTATTTAtc                    >  2:1024630/1‑97 (MQ=255)
                                                                                    aaaaaGAACCCAATAAATATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGATTAATGTGTGATTTtttattt                       >  2:568568/2‑94 (MQ=38)
                                                                                      aaaGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTTTATTTATCGTCGa               >  2:2147898/1‑100 (MQ=255)
                                                                                      aaaGAACCCAATAAATATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTTTATTTATCGTCGa               >  1:1495111/1‑100 (MQ=255)
                                                                                       aaGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTTTATTTATCGTCGaa              >  2:1821010/1‑100 (MQ=255)
                                                                                        agaaCCCAATAAATATTGGGTCCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTTTATTTATCGTCGAAc             >  2:592434/4‑100 (MQ=38)
                                                                                         gaaCCCAATAAATATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTTTATTTAtcgt                  >  2:407801/3‑94 (MQ=255)
                                                                                           aCCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTTTATTTATCGTCGAACCTg          >  1:1084206/1‑100 (MQ=255)
                                                                                           aCCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTTTATTTATCGTCGAACCTg          <  1:808946/100‑1 (MQ=255)
                                                                                            cccAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTTTATTTATCGTCGAACCt           >  2:1285472/1‑98 (MQ=255)
                                                                                             ccAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATGACCAGAGTTTAATGTGTGATTTTTTATTTATTGTCGAACCTgg         >  1:1981013/1‑99 (MQ=255)
                                                                                             ccAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTTTATTTATCGTCGAACCTGGa        >  1:2110506/1‑100 (MQ=255)
                                                                                             ccAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTTTATGTATCGTCGAACCTGGa        >  2:1496854/1‑100 (MQ=255)
                                                                                              cAATATTTATTGTGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTTTGATTTTTTATTTAtcgtcg                <  2:64043/91‑1 (MQ=255)
                                                                                              cAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTTTATTTATCGTCGAACCTggct       >  1:1323515/1‑98 (MQ=255)
                                                                                              cAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTTTATTTATCGTCGAACCTGGAt       >  1:1935459/1‑100 (MQ=255)
                                                                                               aaTATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTTTATTTATCGTCGAACCTGGAtt      <  2:420922/100‑1 (MQ=255)
                                                                                                atatTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTTTATTTATCGTCGAACCTGGATTg     >  1:202363/1‑100 (MQ=255)
                                                                                                 tatttattGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGTTATATATCACGAGAGTTTAATGTGTGATTTTTTACTTATCGTCGAACCTGGATTGt    <  2:1620869/100‑1 (MQ=255)
                                                                                                 tatttattGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTTTATTTATCGTCGAACCTGGATTGt    >  1:880733/1‑100 (MQ=255)
                                                                                                   tttattGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTTTATTTATCGTCGAACCTGGATTGttt  <  2:633779/100‑1 (MQ=255)
                                                                                                   tttattGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTTTATTTATCGTCGAACCTGGATTGttt  <  1:317480/100‑1 (MQ=255)
                                                                                                   tttattGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTTTATTTATCGTCGAACCTGGATTGttt  >  1:518368/1‑100 (MQ=255)
                                                                                                   |                                                                                                   
ACTTCTAAATTAATTTATGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTTTATTTATCGTCGAACCTGGATTGTTT  >  NC_000913/4261259‑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: