Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,261,358 T→A 18.1% 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,3580TA18.1% 322.2 / 58.8 132intergenic (+52/‑311)yjbM/dusAuncharacterized protein/tRNA‑dihydrouridine synthase A
Reads supporting (aligned to +/- strand):  ref base T (53/55);  new base A (14/10);  total (67/65)
Fisher's exact test for biased strand distribution p-value = 5.00e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.86e-01

Only 100 of 134 total aligned reads displayed.

CTTCTAAATTAATTTATGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTTTATTTATCGTCGAACCTGGATTGTTT  >  NC_000913/4261260‑4261457
                                                                                                  |                                                                                                   
     aaaTTAATTTATGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTg                                                                                               <  2:909916/100‑1 (MQ=255)
      aattaattTATGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTgg                                                                                              >  1:1324892/1‑100 (MQ=255)
        ttaattTATGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGt                                                                                            <  2:198086/100‑1 (MQ=255)
        ttaattTATGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGt                                                                                            <  2:1490562/100‑1 (MQ=255)
        ttaattTATGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGt                                                                                            <  2:2304117/100‑1 (MQ=255)
           attTATGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCt                                                                                         <  1:1497497/100‑1 (MQ=255)
           attTATGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGAACCCAATAAATATTGGGTCCt                                                                                         <  2:584364/100‑1 (MQ=38)
             ttATGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCttt                                                                                       >  1:696853/1‑100 (MQ=255)
             ttATGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATAAATATTGGGTCCttt                                                                                       <  1:61594/100‑1 (MQ=39)
              tATGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATAAATATTGGGTCCtttt                                                                                      <  2:199818/100‑1 (MQ=39)
              tATGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGAACCCAATAAATATTGGGTCCtttt                                                                                      <  1:1901813/100‑1 (MQ=38)
              tATGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGAACCCAATAAATATTGGGTCCt                                                                                         <  1:2580190/97‑1 (MQ=39)
                tgtttttttGGGGGTTAAATCGGAGAAATATTATGTGAAGTTCTTCGATAGTATTGAAGGCATTATATAAAAGCACCCAATAAATATTGGGTCCttt                                                                                       <  1:2681724/94‑1 (MQ=25)
                tGATTTGTTCGGGGTTAAATCGGATAAATAATATTTAAAGTGCTTCGATAGTATGGAAGGCATGATATAAAAGGACCCAATATTTATTGGGTTCttttt                                                                                     <  1:527377/99‑1 (MQ=255)
                tGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCtttttt                                                                                    >  2:1779551/1‑100 (MQ=255)
                tGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCtttttt                                                                                    >  1:2142060/1‑100 (MQ=255)
                tGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCtttttn                                                                                    <  1:2298903/100‑2 (MQ=255)
                 gATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATAAATATTGGGTCCttttata                                                                                   <  2:913913/100‑4 (MQ=38)
                  aTTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTtct                                                                                  <  2:2619744/100‑1 (MQ=255)
                  aTTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTtct                                                                                  <  2:1300511/100‑1 (MQ=255)
                  aTTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTtct                                                                                  >  1:1834302/1‑100 (MQ=255)
                  aTTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGAACCCAATAAATATTGGGTCCttttatat                                                                                  <  1:1492341/100‑5 (MQ=37)
                    ttGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTctct                                                                                <  2:430930/100‑1 (MQ=255)
                     tGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTa                                                                               >  2:1199992/1‑100 (MQ=255)
                      gTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGAACCCAATAAATATTGGGTCCttttatataatg                                                                              <  2:2851225/100‑9 (MQ=25)
                       ttCGGGGTTAAATCGGAGAAATAATATGTGTAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGCCCCAATATTTATTGGGTTCTTTTTTCTCTATc                                                                             <  1:1614382/100‑1 (MQ=255)
                       ttCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATc                                                                             <  2:1317839/100‑1 (MQ=255)
                       ttCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATc                                                                             <  1:1689932/100‑1 (MQ=255)
                        tCGGGGGTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGAACCCAATAAATATCGGGTCCttttatataatgcc                                                                            <  1:1776379/100‑11 (MQ=21)
                         cGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCaa                                                                           <  2:1643876/100‑1 (MQ=255)
                         cGGGGTTAAATCGGAGAAATAATACGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCaa                                                                           >  1:2089707/1‑100 (MQ=255)
                          ggggTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAAt                                                                          >  2:702935/1‑100 (MQ=255)
                          gaggTTAAATCGGAGAGATAATATGTGAATTTCTTCGATAGTATTGAAGTCATTATATAAAATGACCCAATATTTATTGGGTGCTTTTTTCTCTATCAAt                                                                          <  1:1822686/98‑1 (MQ=255)
                             gTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACt                                                                       >  1:1188282/1‑100 (MQ=255)
                               tAAATCGGAGAAATAATATGTGGAGTTCTTCGATAGTATGGAAGGCATTATATAAAAAGACCCAATATTTATTGGTTTCTTTTTTCTCTATCACTACTAt                                                                     >  1:172401/1‑100 (MQ=255)
                                 aaTCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTa                                                                   >  2:333469/1‑100 (MQ=255)
                                 aaTCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTa                                                                   <  2:15106/100‑1 (MQ=255)
                                  aTCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAg                                                                  >  2:57128/1‑100 (MQ=255)
                                  aTCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAg                                                                  >  2:2301119/1‑100 (MQ=255)
                                   tCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGc                                                                 <  2:193003/100‑1 (MQ=255)
                                   tCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGc                                                                 <  2:1378421/100‑1 (MQ=255)
                                    cGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGc                                                                 <  2:130701/99‑1 (MQ=255)
                                     ggAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAg                                                               >  2:1556169/1‑100 (MQ=255)
                                      gagaAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAgg                                                              <  1:2550936/100‑1 (MQ=255)
                                        gaAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGga                                                            >  2:2879885/1‑100 (MQ=255)
                                        gaAATAATATGTGAAGTTCTTAGATAGTATGGAGGCCATTATATAAATGGATACAATATTAACTGGGTTCTTTTTTATCTATATATACTact                                                                    >  1:1004448/1‑90 (MQ=255)
                                         aaaTAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGgag                                                           <  1:2899202/100‑1 (MQ=255)
                                         aaaTAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGgag                                                           <  2:1938037/100‑1 (MQ=255)
                                            taataTGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGata                                                        <  1:2885011/100‑1 (MQ=255)
                                            taataTGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGata                                                        <  2:1165099/100‑1 (MQ=255)
                                             aataTGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGatat                                                       <  1:159689/100‑1 (MQ=255)
                                              atatGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGatata                                                      <  2:1834310/100‑1 (MQ=255)
                                                      aGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATAGTACCaga                                              <  1:2026612/100‑1 (MQ=255)
                                                           ttcGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTa                                         >  2:599192/1‑100 (MQ=255)
                                                            tcGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATTACCagag                                             >  1:1741732/1‑95 (MQ=255)
                                                               aTAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAAtgt                                     >  2:2811265/1‑100 (MQ=255)
                                                               aTAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAAtgt                                     >  1:1362023/1‑100 (MQ=255)
                                                               aTAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATGGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAAtgt                                     >  1:1365809/1‑100 (MQ=255)
                                                                tAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAAtgtg                                    <  2:850402/100‑1 (MQ=255)
                                                                tAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAAtgtg                                    <  2:1115473/100‑1 (MQ=255)
                                                                 aGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAAtgtgt                                   >  1:1904/1‑100 (MQ=255)
                                                                  gTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAAtgtgtg                                  >  2:1174851/1‑100 (MQ=255)
                                                                  gTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAAtgtgtg                                  >  2:675575/1‑100 (MQ=255)
                                                                    aTGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGAt                                >  2:938078/1‑100 (MQ=255)
                                                                      ggAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGAttt                              >  2:361047/1‑100 (MQ=255)
                                                                       gAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGAtttt                             >  1:2819220/1‑100 (MQ=255)
                                                                       gAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGAtttt                             <  1:507941/100‑1 (MQ=255)
                                                                         aGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGAtttttt                           <  1:2783278/100‑1 (MQ=255)
                                                                          ggCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTttta                          >  2:2086469/1‑100 (MQ=255)
                                                                          ggCATTATATAAAAGGACACAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGAtttttt                           <  1:2762141/99‑1 (MQ=255)
                                                                           gatagagaaaaaaGAACCCAATAAATATTGGGTCCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTtttat                         >  2:113517/10‑100 (MQ=21)
                                                                           gCATTATATACACGGCCCCAATATTTATTGGGTTCTTTTTTCTATATCAATAATATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTttta                          <  2:133445/99‑1 (MQ=255)
                                                                           gCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGAtttttttt                         <  2:1971223/100‑3 (MQ=255)
                                                                           gCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTtttat                         <  1:1058935/100‑1 (MQ=255)
                                                                             tagagaaaaaaGAACCCAATAAATATTGGGTCATTTTTTCTGTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTtttattt                       >  1:1586796/8‑100 (MQ=25)
                                                                             aTTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTtttattt                       >  2:2887032/1‑100 (MQ=255)
                                                                               tatataAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGATATATATCACCAGAGTTTAATGTGTGATTTtttattta                      <  1:773439/99‑1 (MQ=255)
                                                                               tatataAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTttatttat                     <  1:2336992/100‑1 (MQ=255)
                                                                               tatataAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTttatttat                     <  1:2337108/100‑1 (MQ=255)
                                                                               gagaaaaaaGAACCCAATAAATATTGGGTCCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTttatttat                     >  1:88239/6‑100 (MQ=25)
                                                                                atataAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTTTATTTAtc                    <  2:1362028/100‑1 (MQ=255)
                                                                                agaaaaaaGAACCCAATAAATATTGGGTCCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATAACCAGAGTTTCATGTGTGATTTTTTATTTAtc                    >  2:2470354/5‑100 (MQ=25)
                                                                                 tataAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCCGGGAGCTATATCCCCAGAATTTAATGTGTGATTTTTTATTTAtc                    >  1:1335654/1‑99 (MQ=255)
                                                                                 gaaaaaataaCCCAATAAATATTGGGTCCTTTTTTCTCTATCAATACCATTAGCAGGGAGATATATAACCAGAGTTTAATGTGTGATTTTTTATTTAtcg                   >  1:1011799/10‑100 (MQ=37)
                                                                                  ataAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTTTATTTAtcgt                  >  2:2564938/1‑100 (MQ=255)
                                                                                   taAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTTTATTTAtcgtc                 >  1:1184648/1‑100 (MQ=255)
                                                                                   aaaaaGAACCCAATAAATATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTTTATTTAtcgtc                 >  2:877793/2‑100 (MQ=39)
                                                                                   aaaaaGAACCCAATAAATATTGGGTCCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTTTATTTAtcgtc                 >  2:2281055/2‑100 (MQ=38)
                                                                                    aaaaGAACCCAATAAATATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTTTATTTAtcgtcg                >  2:370463/1‑100 (MQ=39)
                                                                                     aaaGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTTTATTTATCGTCGa               >  1:2618465/1‑100 (MQ=255)
                                                                                      aaGAACCCAATAAATATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTTTATTTATCGTCGaa              >  2:1489940/1‑100 (MQ=255)
                                                                                        ggACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTTTATTTATCGTCGAAcc            <  1:688974/100‑1 (MQ=255)
                                                                                          aCCCAATATTTATTGGGTTCTTGTTTCTCTATCAATACTATTAGCAGTGAGATATATCACCAGAGTTTAATGTGTGATTTTGTATTTATCGTCGAACCTg          <  1:939716/100‑1 (MQ=255)
                                                                                           cccAATATTTATTGGGTTCTTTTTTTTTTTTTTTTCCTATTAGCAGGGCGATATATCACCTTAGTTTAATGTGTGATTTTTTATTTAtcctccccc             >  1:67549/1‑92 (MQ=255)
                                                                                           cccAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTTTATTTATCGTCGAACCTgg         >  2:572571/1‑100 (MQ=255)
                                                                                           cccAATATTTATTGGGTTCTTTATTATCTATCAATACAATTAGCAGGGAAATATATCACCAGAGTTTACTGTGAGATTTTTTATTTATAGTCGAACCTg          >  1:149883/1‑99 (MQ=255)
                                                                                            ccAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTTTATTTATCGTCGAACCTGGa        <  1:2160940/100‑1 (MQ=255)
                                                                                             cAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTTTATTTATCGTCGAACCTGGAt       >  2:1866265/1‑100 (MQ=255)
                                                                                              aaTATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTTTATTTATCGTCGAACCTGGAtt      <  2:1414258/100‑1 (MQ=255)
                                                                                                  tttattGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTTTATTTATCGTCGAACCTGGATTGttt  <  2:2738789/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 < 22 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: