Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,261,359 T→A 15.7% intergenic (+53/‑310) 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,3590TA15.7% 272.4 / 40.2 110intergenic (+53/‑310)yjbM/dusAuncharacterized protein/tRNA‑dihydrouridine synthase A
Reads supporting (aligned to +/- strand):  ref base T (51/41);  new base A (8/9);  total (60/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.06e-01

Only 100 of 123 total aligned reads displayed.

CTTCTAAATTAATTTATGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATAT‑AAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTTTATTTATCGTCGAACCTGGATTGTTTA  >  NC_000913/4261260‑4261458
                                                                                                    |                                                                                                   
cTTCTAAATTAATTTATGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATAT‑AAAAGGACCCAAtatt                                                                                                     >  2:952192/1‑100 (MQ=255)
  tCTAAATTAATTTATGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAATTTCTTCGATAGTATGGAAGGCATTATAT‑AAAAGGACCCAAtattta                                                                                                   <  2:1130800/100‑1 (MQ=255)
  tCTAAATTAATTTATGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATAT‑AAAAGGACCCAAtattta                                                                                                   >  2:1220250/1‑100 (MQ=255)
   cTAAATTAATTTATGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATAT‑AAAAGGACCCAAtatttat                                                                                                  >  2:220149/1‑100 (MQ=255)
    tAAATTAATTTATGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATAT‑AAAAGGACCCAAtatttatt                                                                                                 >  2:2436657/1‑100 (MQ=255)
     aaaTTAATTTATGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATAT‑AAAAGGACCCAATATTTATTg                                                                                                <  1:635644/100‑1 (MQ=255)
       attaattTATGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATAT‑AAAAGGACCCAATATTTATTggg                                                                                              >  1:518716/1‑100 (MQ=255)
       attaattTATGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATAT‑AAAAGGACCCAATATTTATTggg                                                                                              <  1:2062386/100‑1 (MQ=255)
        ttaattTATTATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATAT‑AAAAGGACCCAATATTTATTGGGt                                                                                             <  1:2442256/100‑1 (MQ=255)
          aattTATGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATAT‑AAAAGGACCCAATATTTATTGGGTTc                                                                                           <  2:784972/100‑1 (MQ=255)
           attTATGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATAT‑AAAAGAACCCAATAAATATTGGGTCCt                                                                                          <  1:179174/100‑1 (MQ=38)
            tttATGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATAT‑AAAAGGACCCAATATTTATTGGGTTCtt                                                                                         >  1:800013/1‑100 (MQ=255)
             ttATGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAAGAACCCAATAAATATTGGGTCCtt                                                                                         <  2:1003458/100‑1 (MQ=37)
             ttATGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATAT‑AAAAGAACCCAATAAATATTGGGTCCttt                                                                                        <  1:1814811/100‑1 (MQ=38)
             ttATGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATAT‑AAAAGAACCCAATAAATATTGGGTCCttt                                                                                        <  1:558050/100‑1 (MQ=38)
             ttATGATTTATTCGGGGTTAAATCGGAGAAATTATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATAT‑AAAAGGACCCATTATTTATTGGGTTCt                                                                                          <  2:1827984/98‑1 (MQ=255)
              tATGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATAT‑AAAAGGACCCAATAAATATTGGGTCCtttt                                                                                       <  1:2234942/100‑1 (MQ=39)
               aTGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATAT‑AAAAGGACTCAATATTTATTGGGTTCttttt                                                                                      >  2:1938558/1‑100 (MQ=255)
                 gATTTGTTCTTAATTAAATCGGAAAAATAATATATGAAGTTCTTCGATAGTATGGAAGTCATTATAT‑AAAAGGACCCAATATTTATTGGGTTCtt                                                                                         <  2:1474783/95‑1 (MQ=255)
                 gATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATAT‑AAAAGGACCCAATATTTATTGGGTTCTTTTTtc                                                                                    <  2:386298/100‑1 (MQ=255)
                   tttGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATAT‑AAAAGAACCCAATAAATATTGGGTCCttttatata                                                                                  <  1:2014307/100‑6 (MQ=25)
                   tttGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGCTCTTCGATAGTATGGAAGGCATTATAT‑AAAAGGACACAATATTTATTGGGTTCTTTTTtctc                                                                                  >  2:1369876/1‑100 (MQ=255)
                    ttGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATAT‑AAAAGAACCCAATAAATATTGGGTCCttttatataa                                                                                 <  1:2438707/100‑7 (MQ=25)
                       ttCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATAT‑AAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATc                                                                              <  1:286439/100‑1 (MQ=255)
                          ggggTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATAT‑AAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAAt                                                                           >  1:1018087/1‑100 (MQ=255)
                          ggggTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATAT‑AAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAAt                                                                           >  1:1511149/1‑100 (MQ=255)
                          ggggTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATAT‑AAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAAt                                                                           >  1:2549396/1‑100 (MQ=255)
                           gggTTAATTCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAGGTCATTATAT‑AAAATGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATa                                                                          <  2:1579641/100‑1 (MQ=255)
                            ggTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATAT‑AAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATAc                                                                         <  2:2505399/100‑1 (MQ=255)
                             gTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATAT‑AAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACt                                                                        >  2:1606565/1‑100 (MQ=255)
                              ttAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATAT‑AAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTa                                                                       >  1:627401/1‑100 (MQ=255)
                                aaaTCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATAT‑AAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTAtt                                                                     <  2:2473777/100‑1 (MQ=255)
                                aaaTCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATAT‑AAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTAtt                                                                     <  2:736530/100‑1 (MQ=255)
                                 aaTCTGAGAAATAATATGTGAAGTTCTTCGATATTGTGGTAGGCATTATAT‑AATAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTa                                                                    <  2:1261751/100‑1 (MQ=255)
                                  aTCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATAT‑AAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAg                                                                   >  1:863623/1‑100 (MQ=255)
                                   tCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATAT‑AAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGc                                                                  <  2:2241945/100‑1 (MQ=255)
                                      gagaAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATAT‑AAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAgg                                                               >  2:471022/1‑100 (MQ=255)
                                        gaAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATAT‑AAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGga                                                             >  1:52731/1‑100 (MQ=255)
                                        gaAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATAT‑AAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGga                                                             <  2:2016693/100‑1 (MQ=255)
                                         aaaTAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATAT‑AAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGgag                                                            <  2:463896/100‑1 (MQ=255)
                                          aataatATGTGAAGTTCTTCGATAGTATGGAAGGCATTATAT‑AAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGgaga                                                           >  2:1945331/1‑100 (MQ=255)
                                          aataatATGTGAAGTTCTTCGATAGTATGGAAGGCATTATAT‑AAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGgaga                                                           >  1:687842/1‑100 (MQ=255)
                                            taataTGTGAAGTTCTTCGATAGTATGGAAGGCATTATAT‑AAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGata                                                         <  1:2502019/100‑1 (MQ=255)
                                            taataTGTGAAGTTCTTCGATAGTATGGAAGGCATTATAT‑AAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGata                                                         <  1:929157/100‑1 (MQ=255)
                                             aataTGTGAAGGTTTTCGATAGTATGTAAGTCATTATAT‑AAAAGGACCTAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGatat                                                        <  2:1232646/100‑1 (MQ=255)
                                              atatGTGAAGTTCTTCGATAGTATGGAAGGCATTATAT‑AAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGatata                                                       <  1:279546/100‑1 (MQ=255)
                                              atatGTGAAGTTCTTCGATAGTATGGAAGGCATTATAT‑AAAAGGACCCAATATCTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGatata                                                       >  1:1226934/1‑100 (MQ=255)
                                               tatGTGAAGTTCTTCGATAGTATGGAAGGCATTATAT‑AAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGatatat                                                      >  1:684240/1‑100 (MQ=255)
                                                atGTGAAGTTCTTCGATAGTATGGAAGGCATTATAT‑AAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATc                                                     >  2:161562/1‑100 (MQ=255)
                                                  gtgAAGTTCTTCGATAGTATGGAAGGCATTATAT‑AAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCAc                                                   >  2:794701/1‑100 (MQ=255)
                                                    gAAGTTCTTCGATAGTATGGAAGGCATTATAT‑AAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCAcc                                                  >  2:1704745/1‑99 (MQ=255)
                                                      aGTTCTTCGATAGTATGGAAGGCATTATAT‑AAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCaga                                               <  1:1923392/100‑1 (MQ=255)
                                                       gTTCTTCGATAGTATGGAAGGCATTATAT‑AAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCagag                                              <  2:518709/100‑1 (MQ=255)
                                                           ttcGATAGTATCGAAGGCATTATAT‑AAAAGGACCCAATATTTACTGGGTTCTTTTTCCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGt                                             <  2:140276/97‑1 (MQ=255)
                                                            tcGATAGTATGGAAGGCATTATAT‑AAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTaa                                         >  1:1288551/1‑100 (MQ=255)
                                                             cGATAGTATGGAAGGCATTATAT‑AAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAAt                                        <  2:2235455/100‑1 (MQ=255)
                                                              gATAGTATGGAAGGCATTATAT‑AAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTATCAGGGAGATATATCACCAGAGTTTAAtg                                       >  2:623930/1‑100 (MQ=255)
                                                              gATAGTATGGAAGGCATTATAT‑AAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAAtg                                       >  2:1903543/1‑100 (MQ=255)
                                                              gATAGTATGGAAGGCATTATAT‑AAAAGGACCCAAGATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAAtg                                       <  1:2281729/100‑1 (MQ=255)
                                                                 aGTATGGAAGGCATTATAT‑AAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAAtgtgt                                    >  1:1176105/1‑100 (MQ=255)
                                                                   tATGGAAGGCATTATAT‑AAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGa                                  >  1:1460904/1‑100 (MQ=255)
                                                                   tATGGAAGGCATTATAT‑AAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGa                                  >  2:342662/1‑100 (MQ=255)
                                                                    aTGGAAGGCATTATAT‑AAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGAt                                 >  2:90381/1‑100 (MQ=255)
                                                                     tGGAAGGCATTATAT‑AAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGAtt                                >  1:833195/1‑100 (MQ=255)
                                                                       naaGGCATTATAT‑AAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGAttt                               >  2:141424/2‑99 (MQ=255)
                                                                           gCATTATAT‑AAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTtttat                          <  1:695629/100‑1 (MQ=255)
                                                                           gCATTATAT‑AAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTtttat                          <  2:1270632/100‑1 (MQ=255)
                                                                            cATTATAT‑AAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTtttatt                         <  1:401737/100‑1 (MQ=255)
                                                                             aTTATAT‑AAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTtttat                          <  2:2218083/98‑1 (MQ=255)
                                                                              ttATAT‑AAAAGGACCTAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGTTGATATATCACCAGAGTTTAATGTGTGATTTtttatt                         <  2:129001/98‑1 (MQ=255)
                                                                              ttATAT‑AAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTtttattta                       <  1:1363606/100‑1 (MQ=255)
                                                                              ttATAT‑AAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTtttattta                       >  1:1679779/1‑100 (MQ=255)
                                                                               agagaaaaaaGAACCCAATAAATATTGGGTCCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTtttattta                       >  1:1700268/7‑100 (MQ=25)
                                                                                   t‑aAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTTTATTTAtcgtc                  <  1:2436671/100‑1 (MQ=255)
                                                                                   t‑aAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATAACCAGAGTTTAATTTGTGATTTTTTATTTAtc                     >  2:1024630/1‑97 (MQ=255)
                                                                                   aaaaaaGAACCCAATAAATATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTTTATTTAtcgt                   >  2:307179/3‑100 (MQ=38)
                                                                                    aaaaaGAACCCAATAAATATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGATTAATGTGTGATTTtttattt                        >  2:568568/2‑94 (MQ=38)
                                                                                     aaaaGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTTTATTTAtcgtcg                 <  2:1728994/100‑1 (MQ=255)
                                                                                     aaaaGAACCCAATAAATATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTTTGATTTTTTATTTAtc                     >  2:2190317/1‑96 (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)
                                                                                            cccAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTTTATTTATCGTCGAACCTgg          >  2:1471021/1‑100 (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)
                                                                                              cAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTTTATTTATCGTCGAACCTgg          >  2:1328790/1‑98 (MQ=255)
                                                                                              cAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTTTATTTATCGTCGAACCTGGAt        >  1:1935459/1‑100 (MQ=255)
                                                                                               aaTATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTTTATTTATCGTCGAACCTGGAtt       <  2:420922/100‑1 (MQ=255)
                                                                                                 tatttattGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGTTATATATCACGAGAGTTTAATGTGTGATTTTTTACTTATCGTCGAACCTGGATTGt     <  2:1620869/100‑1 (MQ=255)
                                                                                                 tatttattGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTTTATTTATCGTCGAACCTGGATTGt     >  1:880733/1‑100 (MQ=255)
                                                                                                   tttattGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTTTATTTATCGTCGAACCTGGATTGttt   <  1:317480/100‑1 (MQ=255)
                                                                                                   tttattGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTTTATTTATCGTCGAACCTGGATTGttt   <  2:633779/100‑1 (MQ=255)
                                                                                                   tttattGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTTTATTTATCGTCGAACCTGGATTGttt   >  1:518368/1‑100 (MQ=255)
                                                                                                    |                                                                                                   
CTTCTAAATTAATTTATGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATAT‑AAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTTTATTTATCGTCGAACCTGGATTGTTTA  >  NC_000913/4261260‑4261458

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