Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,197,365 G→A 28.5% pseudogene (892/2847 nt) yehH → DUF4132 family pseudogene; defective yehI paralog; molybdate metabolism regulator, first fragment; molybdate metabolism regulator, second fragment 2; yehI' molR mutations have phenotypes and the N‑terminal molR fragment of yehI' may be expressed and may have retained or evolved MolR function

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,197,3650GA28.5% 181.7 / 88.6 126pseudogene (892/2847 nt)yehHDUF4132 family pseudogene; defective yehI paralog; molybdate metabolism regulator, first fragment; molybdate metabolism regulator, second fragment 2; yehI' molR mutations have phenotypes and the N‑terminal molR fragment of yehI' may be expressed and may have retained or evolved MolR function
Reads supporting (aligned to +/- strand):  ref base G (48/42);  new base A (26/10);  total (74/52)
Fisher's exact test for biased strand distribution p-value = 7.10e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.64e-01

Only 100 of 130 total aligned reads displayed.

AGCTACTTGGAGCAATATTCCAGATTTGAATGACGCCTACGTAAACATTTATCACTGACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTGCCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGCTTTGTTATTACCCGAAAAACCAGAAATTGCACATGAAATTGCCCAGCGTCTACTGG  >  NC_000913/2197266‑2197464
                                                                                                   |                                                                                                   
aGCTACTTGGAGCAATATTCCAGATTTGAATGACGCCTACGTAAACATTTATCACTGACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTg                                                                                                     <  1:2105184/100‑1 (MQ=255)
   tACTTGGAGCAATATTCCAGATTTGAATGACGCCTACGTAAACATTTATCACTGACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTGCCg                                                                                                  <  2:2059765/100‑1 (MQ=255)
     cTTGGAGCAATATTCCAGATTTGAATGACGCCTACGTAAACATTTATCACTGACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTGCCGCa                                                                                                >  1:3899117/1‑100 (MQ=255)
      ttGGAGCAATATTCCAGATTTGAATGACGCCTACGTAAACATTTATCACTGACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTACCGCaa                                                                                               >  2:557940/1‑100 (MQ=255)
       tGGAGCAATATTCCAGATTTGAATGACGCCTACGTAAACATTTATCACTGACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTACCGCAAt                                                                                              >  2:757656/1‑100 (MQ=255)
         gAGCAATATTCCAGATTTGAATGACGCCTACGTAAACATTTATCACTGACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTGCCGCAATTc                                                                                            <  2:1591343/100‑1 (MQ=255)
         gAGCAATATTCCAGATTTGAATGACGCCTACGTAAACATTTATCACTGACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTACCGCAATTc                                                                                            >  2:2042022/1‑100 (MQ=255)
          aGCAATATTCCAGATTTGAATGACGCCTACGTAAACATGTATCACTGACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTGCCGCAATTcc                                                                                           <  2:2194701/100‑1 (MQ=255)
           gCAATATTCCAGATTTGAATGACGCCTACGTAAACATTTATCACTGACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTGCCGCAATTCca                                                                                          >  1:2941569/1‑100 (MQ=255)
           gCAATATTCCAGATTTGAATGACGCCTACGTAAACATTTATCACTGACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTACCTCAATTCca                                                                                          >  2:794858/1‑100 (MQ=255)
           gCAATATTCCAGATTTGAATGACGCCTACGGAAACATTTATCACTGACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTGCCGCAATTCca                                                                                          >  1:3532306/1‑100 (MQ=255)
            cAATATTCCAGATTTGAATGACGCCTACGTAAACATTTATCACTGACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTGCCGCAATTCcac                                                                                         <  2:380493/100‑1 (MQ=255)
            cAATATTCCAGATTTGAATGACGCCTACGTAAACATTTATCACTGACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTACCGCAATTCcac                                                                                         <  1:1650113/100‑1 (MQ=255)
             aaTATTCCAGATTTGAATGACGCCTACGTAAACATTTATCACTGACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTGCCGCAATTCcaca                                                                                        <  2:1654621/100‑1 (MQ=255)
               tatTCCAGATTTGAATGACGCCTACGTAAACATTTATCACTGACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTGCCGCAATTCCACata                                                                                      >  1:2592384/1‑100 (MQ=255)
                atTCCAGATTTGAATGACGCCTACGTAAACATTTATCACTGACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTGCCGCAATTCCACatat                                                                                     >  1:3566040/1‑100 (MQ=255)
                  tCCAGATTTGAATGACGCCTACGTAAACATTTATCACTGACAGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTGCCGCAATTCCACATATTc                                                                                   >  2:2254117/1‑100 (MQ=255)
                   ccAGATTTGAATGACGCCTACGTAAACATTTATCACTGACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTACCGCAATTCCACATATTcc                                                                                  >  1:801984/1‑100 (MQ=255)
                    cAGATTTGAATGACGCCTACCAAAACATTTATCACTAACGAATGAGTCGCTATGGCAAAGCTTTGTACAAAAACTTATCGCCGCCATGCCACATATTc                                                                                   >  2:1201663/1‑98 (MQ=255)
                     aGATTTGAATGACGCCTACGTAAACATTTATCACTGACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTGCCGCAATTCCACATATTCCAg                                                                                >  1:3807657/1‑100 (MQ=255)
                       aTTTGAATGACGCCTACGTAAACATTTATCACTGACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTGCCGCAATTCCACATATTCCAGaa                                                                              >  1:1798806/1‑100 (MQ=255)
                        tttGAATGACGCCTACGTAAACATTTATCACTGACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTGCCGCAATTCCACATATTCCAGAAt                                                                             >  1:2321498/1‑100 (MQ=255)
                        tttGAATGACGCCTACGTAAACATTTATCACTGACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTACCGCAATTCCACATATTCCAGAAt                                                                             >  1:3648484/1‑100 (MQ=255)
                         ttGAATGACGCCTACGTAAACATTTATCACTGACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTGCCGCAATTCCACATATTCCAGAATg                                                                            >  1:558944/1‑100 (MQ=255)
                          tGAATGACGCCTACGTAAACATTTATCACTGACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTACCGCAATTCCACATATTCCAGAATgg                                                                           >  2:2412710/1‑100 (MQ=255)
                           gAATGACGCCTACGTAAACATTTATCACTGACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTGCCGCAATTCCACATATTCCAGAATggc                                                                          >  1:2736286/1‑100 (MQ=255)
                             aTGACGCCTACGTAAACATTTATCACTGACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTGCCGCAATTCCACATATTCCAGAATggcgg                                                                        <  1:574088/100‑1 (MQ=255)
                             aTGACGCCTACGTAAACATTTATCACTGACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTGCCGCAATTCCACATATTCCAGAATggcg                                                                         >  2:4014281/1‑99 (MQ=255)
                             aTGACGCCTACGTAAACATTTATCACTGACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTACCGCAATTCAACATATTCCAGAATggcgg                                                                        <  2:3765170/100‑1 (MQ=255)
                              tGACGCCTACGTAAACATTTATCACTGACGGAAGAGTCTCTCTGGCAAAAATGTTCACAAAAGCTTATTGTCGCAATTCTACATATTCCAGTATggcggc                                                                       <  2:178383/100‑1 (MQ=255)
                              tGACGCCTACGTAAACATTTATCACTGACGGAAGAGTATCTCTGGCAAAAATGTGCACAAAAACTTATTGCCGCAATTCCACATATTCCAGAATggcggc                                                                       <  1:1847933/100‑1 (MQ=255)
                               gACGCCTACGTAAACATTTATCACTGACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTGCCGCAATTCCACATATTCCAGAATGGCGGCa                                                                      >  2:2595234/1‑100 (MQ=255)
                                aCGCCTACGTAAACATTTATCACTGACGGAAGAGCCTCTCTGGCAAAAATGTGCACAAAAACTTATTACCGCAATTCCACATATTCCAGAATGGCGGCaa                                                                     >  2:295617/1‑100 (MQ=255)
                                 cGCCTACGTAAACATTTATCACTGACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTGCCGCAATTCCACATATTCCAGAATGGCGGCAAc                                                                    <  2:2430709/100‑1 (MQ=255)
                                 cGCCTACGTAAACATTTATCACTGACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTGCCGCAATTCCACATATTCCAGAATGGCGGCAAc                                                                    <  2:104304/100‑1 (MQ=255)
                                   ccTACGTAAACATTTATCACTGACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTGCCGCAATTCCACATATTCCAGAATGGCGGCAACCa                                                                  <  1:4035597/100‑1 (MQ=255)
                                       cGTAAACATTTATCACTGACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTGCCGCAATTCCACATATTCCAGAATGGCGGCAACCATTaa                                                              >  1:1813941/1‑100 (MQ=255)
                                       cGTAAACATTTATCACTGACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTGCCGCAATTCCACATATTCCAGAATGGCGGCAACCATTaa                                                              >  1:3553659/1‑100 (MQ=255)
                                        gTAAACATTTATCACTGACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTGCCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAAt                                                             >  2:2128831/1‑100 (MQ=255)
                                         tAAACATTTATCACTGACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTACCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAAtt                                                            >  2:4088548/1‑100 (MQ=255)
                                          aaaCATTTATCACTGACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTGCCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTg                                                           >  1:1532411/1‑100 (MQ=255)
                                           aaCATTTATCACTGACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTGCCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGc                                                          <  1:1262078/100‑1 (MQ=255)
                                            aCATTTATCACTGACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTGCCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGCt                                                         >  2:541878/1‑100 (MQ=255)
                                             cATTTATCACTGACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTGCCGCAATTCCACATATTCCAGAATGGCGGCGACCATTAATTGCt                                                         >  1:482915/1‑99 (MQ=255)
                                                  aTCACTGACTGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTGCCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGCTTTGtta                                                   <  1:3654954/100‑1 (MQ=255)
                                                  aTCACTGACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTGCCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGCTTTGtta                                                   <  2:2668250/100‑1 (MQ=255)
                                                   tCACTGACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTGCCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGCTTTGttat                                                  <  2:134239/100‑1 (MQ=255)
                                                    cACTGACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTGCCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGCTTTGttatt                                                 >  1:3087065/1‑100 (MQ=255)
                                                    cACTGACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTACCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGCTTTGttatt                                                 >  1:3139055/1‑100 (MQ=255)
                                                     aCTGACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTACCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGCTTTGttatta                                                >  2:3860141/1‑100 (MQ=255)
                                                      cTGACGGAAGAGTGTCTCTGGCAAAAATGTGTACAAAAACTTATTACCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGCTTTGTTATTAc                                               >  1:2792966/1‑100 (MQ=255)
                                                      cTGACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTGCCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGCTTTGTTATTAc                                               >  1:1084183/1‑100 (MQ=255)
                                                      cTGACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTGCCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGCTTTGTTATTAc                                               <  1:2254130/100‑1 (MQ=255)
                                                       tGACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTGCCGCTATTCCACATATTCCAGAATGGCGGCAACCATTAATTGCTTTGTTATTAcc                                              <  2:2944456/100‑1 (MQ=255)
                                                       tGACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTGCCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGCTTTGTTATTAcc                                              <  2:3665181/100‑1 (MQ=255)
                                                       tGACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTGCCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGCTTTGTTATTAcc                                              >  2:2770515/1‑100 (MQ=255)
                                                        gACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTGCCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGCTTTGTTATTAccc                                             <  1:3609247/100‑1 (MQ=255)
                                                        gACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTGCCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGCTTTGTTATTAccc                                             <  2:1739539/100‑1 (MQ=255)
                                                        gACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTGCCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGCTTTGTTATTAccc                                             <  2:3532283/100‑1 (MQ=255)
                                                          cGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTGCCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGCTTTGTTATTACCCGa                                           >  1:1008304/1‑100 (MQ=255)
                                                              agagTCTCTCTGGCAAAAATGTGCACAAAAACTTATTGCCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGCTTTGTTATTACCCGaaaaa                                       >  2:1778481/1‑100 (MQ=255)
                                                                agTCTCTCTGGCAAAAATGTGCACAAAAACTTATTGCCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGCTTTGTTATTACCCGAAAAAcc                                     <  1:2323530/100‑1 (MQ=255)
                                                                agTCTCTCTGGCAAAAATGTGCACAAAAACTTATTACCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGCTTTGTTATTACCCGAAAAAcc                                     <  1:2042035/100‑1 (MQ=255)
                                                                agTCTCTCTGGCAAAAATGTGCACAAAAACTTATTACCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGCTTTGTTATTACCCGAAAAAcc                                     >  1:3731003/1‑100 (MQ=255)
                                                                  tctctcTGGCAAAAATGTGCACAAAAACTTATTGCCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGTTTTGTTATTACCCGAAAAACCAg                                   <  1:3213898/100‑1 (MQ=255)
                                                                   ctctctGGCAAAAATGTGCACAAAAGCTTATTGCCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGCTTTGTTATTACCCGAAAAACCAGa                                  <  2:2888593/100‑1 (MQ=255)
                                                                   ctctctGGCAAAAATGTGCACAAAAACTTATTGCCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGCTTTGTTATTACCCGACAAACCAGa                                  >  2:561660/1‑100 (MQ=255)
                                                                   ctctctGGCAAAAATGTGCACAAAAACTTATTGCCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGCTTTGTTATTACCCGAAAAACCAGa                                  >  1:1398827/1‑100 (MQ=255)
                                                                    tctctGGCAAAAATGTGCACAAAAACTTATTGCCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGCTTTGTTATTACCCGAAAAACCAGaa                                 >  2:3006593/1‑100 (MQ=255)
                                                                    tctctGGCAAAAATGTGCACAAAAACTTATTACCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGCTTTGTTATTACCCGAAAAACCAGaa                                 >  1:1985423/1‑100 (MQ=255)
                                                                     ctctGGCAAAAATGTGCACAAAAACTTATTGCCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGCTTTGTTATTACCCGAAAAACCAGaaa                                >  1:1020663/1‑100 (MQ=255)
                                                                       ctGGCAAAAATGTGCACAAAAACTTATTGCCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGCTTTGTTATTACCCGAAAAACCAGAAAtt                              <  1:2541318/100‑1 (MQ=255)
                                                                       ctGGCAAAAATGTGCACAAAAACTTATTACCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGCTTTGTTATTACCCGAAAAACCAGAAAtt                              <  1:1123412/100‑1 (MQ=255)
                                                                        tgtcAAAAATGTGCACAACAACTTATTGCCGTACTTCCACATATTCCAGACTGGCGGCAACCATTAATTGCTTTGTTATCACCCGACAAACCAGAAATTg                             <  2:821439/97‑1 (MQ=255)
                                                                        tGGCAAAAATGTGCACAAAAACTTATTGCCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGCTTTGTTATTACCCGAAAAACCAGAAATTg                             >  2:801079/1‑100 (MQ=255)
                                                                         ggCAAAAATGTGCACAAAAACTTATTGCCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGCTTTGTTATTACCCGAAAAACCAGAAATTGc                            >  1:2343622/1‑100 (MQ=255)
                                                                         ggCAAAAATGTGCACAAAAACTTATTGCCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGCTTTGTTATTACCCGAAAAACCAGAAATTGc                            >  1:336642/1‑100 (MQ=255)
                                                                         ggCAAAAATGTGCACAAAAACTTATTGCCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGCTATGTTATTACCCGAAAAACCAGAAATTGc                            >  1:2560029/1‑100 (MQ=255)
                                                                         ggCAAAAATGTGCACAAAAACTTATTACCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGCTTTGTTATTACCCGAAAAACCAGAAATTGc                            >  1:3851970/1‑100 (MQ=255)
                                                                          gCAAAAATGTGCACAAAAACTTATTGCCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGCTTTGTTATTACCCGAAAAACCAGAAATTGca                           <  1:3301640/100‑1 (MQ=255)
                                                                          gCAAAAATGTGCACAAAAACTTATTACCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGCTTTGTTATTACCCGAAAAACCAGAAATTGca                           >  1:1615712/1‑100 (MQ=255)
                                                                           cAAAAATGTGCACAAAAACTTATTACCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGCTTTGTTATTACCCGAAAAACCAGAAATTGca                           >  1:1787357/1‑99 (MQ=255)
                                                                             aaaaTGTGCACAAAAACTTATTACCGCAATTCCACATATTCCAGAATGGCGGCAACCACTAATTGCTTTGTTATAACCCGAAAGACCAGaaa                                >  2:1454830/1‑92 (MQ=255)
                                                                              aaaTGTGCACAAAAACTTATTGCCGCAATTCCACATATCCCAGAATGGCGGCAACCATTAATTGCTTTGTAATAAACCGAAAAACCAGAAATTGCACATg                       <  2:2987065/100‑1 (MQ=255)
                                                                                  gtgCACAAAAACTTATTACCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGCTTTGTTATTACCCGAAAAACCAGAAATTGCACATGAAAt                   <  2:396131/100‑1 (MQ=255)
                                                                                    gCACAAAAACTTATTACCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGCTTTGTTATTACCCGAAAAACCAGAAATTGCACATGAAATTg                 >  1:2775307/1‑100 (MQ=255)
                                                                                     cacaAAAACTTATTGCCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGCTTTGTTATTACCCGAAAAACCAGAAATTGCACATGAAATTGc                >  2:86956/1‑100 (MQ=255)
                                                                                     cacaAAAACTTATTACCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGCTTTGTTATTACCCGAAAAACCAGAAATTGCACATGAAATTGc                <  2:3543432/100‑1 (MQ=255)
                                                                                      acaAAACCTTATTACCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGCTTTGTTATTACCCGAAAAACCAGAAATTGCACATGAAATTGcc               >  1:3435653/1‑100 (MQ=255)
                                                                                      acaAAAACTTATTGCCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGCTTTGTTATTACCCGAAAAACCAGAAATTGCACATGAAATTGcc               <  1:913232/100‑1 (MQ=255)
                                                                                      acaAAAACTTATTGCCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGCTTTGTTATTACCCGAAAAACCAGAAATTGCACATGAAATTGc                >  1:4006618/1‑99 (MQ=255)
                                                                                      aaaaaaaCTTATTACCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGCTTTGTTATTACCCGAAAAACCAGAAATTGCACATGAAATTg                 <  2:1504613/96‑1 (MQ=255)
                                                                                       caAAAACTTATTGCCGCAATTCCACATATTCGAGAATGGCTGCAACCATTAATTGCTTTGTTATTACCCGAAAAACCAGAAATTGCACATGAAATTGcc               >  2:1219561/1‑99 (MQ=255)
                                                                                        aaaaaCTTATTACCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGCTTAAAAAAAAAAAAAAAAACCAGAAATTGCACATGAAATTGccc              <  2:3619628/99‑1 (MQ=255)
                                                                                         aaaaCTTATTGCCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGCTTTGTTATTACCCGAAAAACCAGAAATTGCACATGAAATTGCCCAg            <  1:4337/100‑1 (MQ=255)
                                                                                            aCTTATTGCCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGCTTTGTTATTACCCGAAAAACCAGAAATTGCACATGAAATTGCCCAGCGt         <  2:93559/100‑1 (MQ=255)
                                                                                             cTTATTGCCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGCTTTGTTATTACCCGAAAAACCAGAAATTGCACATGAAATTGCCCAGCGTc        <  2:3639226/100‑1 (MQ=255)
                                                                                              ttATTGCCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGCTTTGTTATTACCCGAAAAACCAGAAATTGCACATGAAATTACCCAGCGTCt       >  1:1372436/1‑100 (MQ=255)
                                                                                                   gCCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGCTTTGTTATTACCCGAAAAACCAGAAATTGCACATGAAATTGCCCAGCGTCTACTgg  >  1:788629/1‑100 (MQ=255)
                                                                                                   gCCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGCTTTGTTATTAACCGAAAAACCAGAAATTGCACATGAAATTGCCCAGCGTCTACTgg  >  2:1945077/1‑100 (MQ=255)
                                                                                                   |                                                                                                   
AGCTACTTGGAGCAATATTCCAGATTTGAATGACGCCTACGTAAACATTTATCACTGACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTGCCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGCTTTGTTATTACCCGAAAAACCAGAAATTGCACATGAAATTGCCCAGCGTCTACTGG  >  NC_000913/2197266‑2197464

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: