Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,261,349 G→A 15.5% intergenic (+43/‑320) 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,3490GA15.5% 294.6 / 39.9 124intergenic (+43/‑320)yjbM/dusAuncharacterized protein/tRNA‑dihydrouridine synthase A
Reads supporting (aligned to +/- strand):  ref base G (46/58);  new base A (13/6);  total (59/65)
Fisher's exact test for biased strand distribution p-value = 7.90e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.05e-01

Only 100 of 125 total aligned reads displayed.

TACAGACACTTCTAAATTAATTTATGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTTTATTTATCGTCGAACC  >  NC_000913/4261252‑4261447
                                                                                                 |                                                                                                  
 aCAGACACTTCTAAATTAATTTATGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGAcc                                                                                                 <  1:2217215/100‑1 (MQ=255)
  cAGACACTTCTAAATTAATTTATGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGAccc                                                                                                <  2:600733/100‑1 (MQ=255)
   aGACACTTCTAAATTAATTTATGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCa                                                                                               >  1:1660012/1‑100 (MQ=255)
    gACACTTCTAAATTAATTTATGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCaa                                                                                              <  2:824926/100‑1 (MQ=255)
        cTTCTAAATTAATTTATGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAAtatt                                                                                          >  2:2248963/1‑100 (MQ=255)
             aaaTTAATTTATGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTg                                                                                     <  2:909916/100‑1 (MQ=255)
              aattaattTATGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTgg                                                                                    >  1:1324892/1‑100 (MQ=255)
               attaattTATGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTggg                                                                                   <  2:1506981/100‑1 (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)
                     ttATGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCttt                                                                             >  2:1329302/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)
                       aTGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCttttt                                                                           <  1:1390100/100‑1 (MQ=255)
                        tgtttttttGGGGGTTAAATCGGAGAAATATTATGTGAAGTTCTTCGATAGTATTGAAGGCATTATATAAAAGCACCCAATAAATATTGGGTCCttt                                                                             <  1:2681724/94‑1 (MQ=25)
                        tGATTTGTTCGGGGTTAAATCGGATAAATAATATTTAAAGTGCTTCGATAGTATGGAAGGCATGATATAAAAGGACCCAATATTTATTGGGTTCttttt                                                                           <  1:527377/99‑1 (MQ=255)
                        tGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCtttttt                                                                          >  1:2142060/1‑100 (MQ=255)
                        tGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCtttttt                                                                          >  2:1779551/1‑100 (MQ=255)
                        tGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCtttttn                                                                          <  1:2298903/100‑2 (MQ=255)
                         gATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTtc                                                                         <  1:1329297/100‑1 (MQ=255)
                         gATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATAAATATTGGGTCCttttata                                                                         <  2:913913/100‑4 (MQ=38)
                         gATTTGTTCGGGGTTAAATCGGAGAAATAATATATGAAGGTCTTCGGTAGTATGGAATGCATTATATAAACGGACGCAGCATTTATTGGGTTCtttct                                                                           >  1:1862378/1‑96 (MQ=255)
                          aTTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTtct                                                                        <  2:2619744/100‑1 (MQ=255)
                          aTTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTtct                                                                        <  2:1300511/100‑1 (MQ=255)
                          aTTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGAACCCAATAAATATTGGGTCCttttatat                                                                        <  1:1492341/100‑5 (MQ=37)
                            ttGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTctct                                                                      <  2:430930/100‑1 (MQ=255)
                             tGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTa                                                                     >  2:1199992/1‑100 (MQ=255)
                               ttCGGGGTTAAATCGGAGAAATAATATGTGTAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGCCCCAATATTTATTGGGTTCTTTTTTCTCTATc                                                                   <  1:1614382/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)
                                   gggTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATa                                                               >  2:592318/1‑100 (MQ=255)
                                     gTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACt                                                             >  1:730390/1‑100 (MQ=255)
                                       tAAATCGGAGAAATAATATGTGGAGTTCTTCGATAGTATGGAAGGCATTATATAAAAAGACCCAATATTTATTGGTTTCTTTTTTCTCTATCACTACTAt                                                           >  1:172401/1‑100 (MQ=255)
                                        aaaTCGGAGAAATAATATGTGAAGTTGTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTAtt                                                          >  1:1668309/1‑100 (MQ=255)
                                         aaTCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTa                                                         >  2:333469/1‑100 (MQ=255)
                                         aaTCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTa                                                         <  2:15106/100‑1 (MQ=255)
                                          aTCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAg                                                        >  2:2301119/1‑100 (MQ=255)
                                           tCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGc                                                       <  2:1378421/100‑1 (MQ=255)
                                            cGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGc                                                       <  2:130701/99‑1 (MQ=255)
                                               agaAATAATATGTGACGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAggg                                                   <  2:2872754/100‑1 (MQ=255)
                                                gaAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGga                                                  >  2:2879885/1‑100 (MQ=255)
                                                gaAATAATATGTGAAGTTCTTAGATAGTATGGAGGCCATTATATAAATGGATACAATATTAACTGGGTTCTTTTTTATCTATATATACTact                                                          >  1:1004448/1‑90 (MQ=255)
                                                 aaaTAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGgag                                                 <  2:1938037/100‑1 (MQ=255)
                                                    taataTGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGata                                              <  2:1165099/100‑1 (MQ=255)
                                                    taataTGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGata                                              <  1:2885011/100‑1 (MQ=255)
                                                    taataTGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGata                                              >  1:2666202/1‑100 (MQ=255)
                                                     aataTGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGatat                                             <  1:159689/100‑1 (MQ=255)
                                                      atatGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGatata                                            <  2:1834310/100‑1 (MQ=255)
                                                              aGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATAGTACCaga                                    <  1:2026612/100‑1 (MQ=255)
                                                                  cttcGATAGTATGGAAGGCATTATATAAAAGGACCTAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGttt                                <  1:2789181/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)
                                                                        tAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAAtgtg                          <  2:1115473/100‑1 (MQ=255)
                                                                        tAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAAtgtg                          <  2:850402/100‑1 (MQ=255)
                                                                         aGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAAtgtgt                         >  1:1904/1‑100 (MQ=255)
                                                                          gTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAAtgtgtg                        >  2:1163922/1‑100 (MQ=255)
                                                                            aTGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGAt                      >  2:938078/1‑100 (MQ=255)
                                                                              ggAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGAttt                    >  2:361047/1‑100 (MQ=255)
                                                                               gAAGGCATTATATACAAGGACCCAATATTTATTGGGTTCTTGTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTTTGAtttt                   <  1:401651/100‑1 (MQ=255)
                                                                               gAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGAtttt                   >  1:2819220/1‑100 (MQ=255)
                                                                               gAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGAtttt                   <  1:507941/100‑1 (MQ=255)
                                                                                 aGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGAtttttt                 <  1:2783278/100‑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)
                                                                                    aatgatataAAATGAGCCAATATTTATTGGGGTCTTTTTTCTCTGTCGATACTTTTAGCAGGGAGATATAGTACCAGAGTTTAATGTGTGATTTtttatt              >  1:243273/5‑100 (MQ=255)
                                                                                     aTTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGCGATATATCACCAGAGTTTAATGTGTGATTTtttattt             <  1:1554163/100‑1 (MQ=255)
                                                                                     aTTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTtttattt             >  2:1249037/1‑100 (MQ=255)
                                                                                     aTTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTtttattt             >  2:2887032/1‑100 (MQ=255)
                                                                                      ttATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTtttattta            >  2:838449/1‑100 (MQ=255)
                                                                                      agagaaaaaaGAACCCAATAAATATTGGGTCCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTtttattta            >  2:1894478/7‑100 (MQ=25)
                                                                                       tatataAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGATATATATCACCAGAGTTTAATGTGTGATTTtttattta            <  1:773439/99‑1 (MQ=255)
                                                                                       tatataAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTttatttat           <  2:913206/100‑1 (MQ=255)
                                                                                       tatataAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTttatttat           <  1:2336992/100‑1 (MQ=255)
                                                                                       tatataAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTttatttat           <  1:2337108/100‑1 (MQ=255)
                                                                                       gagaaaaaaGAACCCAATAAATATTGGGTCCTTTTTTCTCTATCAATACTATTAGCAGGGCGATATATCACCAGAGTTTAATGTGTGATTTTttatttat           >  2:537659/6‑100 (MQ=25)
                                                                                       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)
                                                                                         gaaaaaaGAACCCAATAAATATTGGGTCCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTTTATTTAtcg         >  1:2197555/4‑100 (MQ=37)
                                                                                          ataAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTTTATTTAtcgt        >  2:2564938/1‑100 (MQ=255)
                                                                                           taAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTTTATTTAtcgtc       >  1:1184648/1‑100 (MQ=255)
                                                                                           aaaaaGAACCCAATAAATATTGGGTCCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTTTATTTAtcgtc       >  2:2281055/2‑100 (MQ=38)
                                                                                            aaaaGAACCCAATAAATATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTTTATTTAtcgtcg      >  2:370463/1‑100 (MQ=39)
                                                                                             aaaGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTTTATTTATCGTCGa     >  1:2618465/1‑100 (MQ=255)
                                                                                             aaaGAACCCAATAAATATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCATCCGTGTTTAATGTGTTATTTTTTATTTAtcgtgca     >  1:1168322/1‑97 (MQ=37)
                                                                                              aaGAACCCAATAAATATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTTTATTTATCGTCGaa    >  2:1489940/1‑100 (MQ=255)
                                                                                                 |                                                                                                  
TACAGACACTTCTAAATTAATTTATGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTTTATTTATCGTCGAACC  >  NC_000913/4261252‑4261447

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: