Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,197,365 G→A 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,3650GA98.5% 452.2 / ‑3.0 137pseudogene (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 (1/1);  new base A (74/61);  total (75/62)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.46e-01

Only 100 of 137 total aligned reads displayed.

TACTTGGAGCAATATTCCAGATTTGAATGACGCCTACGTAAACATTTATCACTGACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTGCCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGCTTTGTTATTACCCGAAAAACCAGAAATTGCACATGAAATTGCCCAGCGTC  >  NC_000913/2197269‑2197458
                                                                                                |                                                                                             
tACTTGGAGCAATATTCCAGATTTGAATGACGCCTACGTAAACATTTATCACTGACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTACCg                                                                                            >  2:1654890/1‑100 (MQ=255)
  cTTGGAGCAATATTCCAGATTTGAATGACGCCTACGTAAACATTTATCACTGACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTACCGc                                                                                           <  2:1699435/99‑1 (MQ=255)
   ttGGAGCAATATTCCAGATTTGAATGACGCCTACGTAAACATTTATCACTGACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTACCGCa                                                                                          <  1:453660/99‑1 (MQ=255)
     ggAGCAATATTCCAGATTTGAATGACGCCTACGTAAACATTTATCACTGACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTACCGCAAtt                                                                                       >  1:2776443/1‑100 (MQ=255)
      gAGCAATATTCCAGATTTGAATGACGCCTACGTAAACATTTATCACTGACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTACCGCAATTc                                                                                      >  2:2891150/1‑100 (MQ=255)
        gCAATATTCCAGATTTGAATGACGCCTACGTAAACATTTATCACTGACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTACCGCAATTCca                                                                                    >  2:926118/1‑100 (MQ=255)
         cAATATTCCAGATTTGAATGACGCCTACGTAAACATTTATCACTGACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTACCGCAATTCcac                                                                                   <  2:2402354/100‑1 (MQ=255)
          aaTATTCCAGATTTGAATGACGCCTACGTAAACATTTATCACTGACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTACCGCAATTCcaca                                                                                  >  2:1725228/1‑100 (MQ=255)
          aaTATTCCAGATTTGAATGACGCCTACGTAAACATTTATCACTGACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTACCGCAATTCcaca                                                                                  <  2:1416534/100‑1 (MQ=255)
           atatTCCAGATTTGAATGACGCCTACGTAAACATTTATCACTGACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTACCGCAATTCCACat                                                                                 <  2:2625456/100‑1 (MQ=255)
           atatTCCAGATTTGAATGACGCCTACGTAAACATTTATCACTGACGGAAGAGTCTCTCTGGCAAAAATATGCACAAAGACTTATTACCGCAATTCcaca                                                                                  >  1:2348175/1‑99 (MQ=255)
             atTCCAGATTTGAATGACGCCTACGTAAACATTTATCACTGAGGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTACCGCAATTCCACatat                                                                               >  2:1965743/1‑100 (MQ=255)
             atTCCAGATTTGAATGACGCCTACGTAAACATTTATCACTGACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTACCGCAATTCCACatat                                                                               >  2:724094/1‑100 (MQ=255)
              ttCCAGATTTGAATGACGCTTACGTAAACATTTATCACTGACGGAAGAGTCTCTCTGGCAAAAATGTGTAAAAAAACTTATTACCGCAATGCCACATAtt                                                                              >  2:1230863/1‑100 (MQ=255)
               tCCAGATTTGAATGACGCCTACGTAAACATTTATCACTGACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTACCGCAATTCCACATATTc                                                                             >  2:2406948/1‑100 (MQ=255)
                ccAGATTTGAATGACGCCTACGTAAACATTTATCACTGACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTACCGCAATTCCACATATTcc                                                                            >  2:2298677/1‑100 (MQ=255)
                 cAGATTTGAATGACGCCTACGTAAACATTTATCACTGACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTACCGCAATTCCACATATTCCa                                                                           <  2:2247385/100‑1 (MQ=255)
                     tttGAATGACGCCTACGTAAACATTTATCACTGACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTACCGCAATTCCACATATTCCAGAAt                                                                       >  2:1026347/1‑100 (MQ=255)
                       tGAATGACGCCTACGTAAACATTTATCACTGACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTACCGCAATTCCACATATTCCAGAATgg                                                                     >  1:1218373/1‑100 (MQ=255)
                        gAATGACGCCTACGTAAACATTTATCACTGACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTACCGCAATTCCACATATTCCAGAATggc                                                                    >  2:513736/1‑100 (MQ=255)
                         aaTGACGCCTACGTAAACATTTATCACTGACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTACCGCAATTCCACATATTCCAGAATggcg                                                                   >  2:1995530/1‑100 (MQ=255)
                           tGACGCCTACGTAAACATTTATCACTGACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTACCGCAATTCCACATATTCCAGAATggcggc                                                                 >  2:1864043/1‑100 (MQ=255)
                            gACGCCTACGTAAACATTTATCACTTACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTACCGCAAGACCACATATTCCAGAATggcgtc                                                                 >  2:2252759/1‑97 (MQ=255)
                            gACGCCTACGTAAACATTTATCACTGACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTACCGCAATTCCACATATTCCAGAATGGCGGCa                                                                <  1:1213367/100‑1 (MQ=255)
                              cGCCTACGTAAACATTTATCACTGACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTACCGCAATTCCACATATTCCAGAATGGCGGCAAc                                                              >  1:1686538/1‑100 (MQ=255)
                              cGCCTACGTAAACATTTATCACTGACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTACCGCAATTCCACATATTCCAGAATGGCGGCAAc                                                              >  2:1483474/1‑100 (MQ=255)
                               gCCTACGTAAACATTTATCACTGACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTACCGCAATTCCACATATTCCAGAATGGCGGCAAcc                                                             <  2:1787542/100‑1 (MQ=255)
                                ccTACGTAAACATTTATCACTGACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTACCGCAATTCCACATATTCCAGAATGGCGGCAACCa                                                            >  1:1546777/1‑100 (MQ=255)
                                ccTACGTAAACATTTATCACTGACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTACCGCAATTCCACATATTCCAGAATGGCGGCAACCa                                                            >  1:3047652/1‑100 (MQ=255)
                                 cTACGTAAACATTTATCACTGACGGAAGAGTCTCTCTGGCAAAAATGTGCACAACAACTTATTACCGCAATTCCACATATTCCAGAATGGCGGCAACCAt                                                           >  1:1177899/1‑100 (MQ=255)
                                  tACGTAAACATTTATCACTGACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTACCGCAATTCCACATATTCCAGAATGGCGGCAACCAtt                                                          <  2:341666/100‑1 (MQ=255)
                                      tAAACATTTATCACTGACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTACCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAAtt                                                      >  2:2364339/1‑100 (MQ=255)
                                      tAAACATTTATCACTGACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTACCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAAta                                                      >  2:2957540/1‑99 (MQ=255)
                                       aaaCATTTATCAGTGACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTACCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAAtt                                                      <  2:3272100/99‑1 (MQ=255)
                                       aaaCATTTATCACTGACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTACCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTg                                                     >  2:1966704/1‑100 (MQ=255)
                                        aaCATTTATCACTGACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTACCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGc                                                    <  2:698338/100‑1 (MQ=255)
                                          cATTTATCACTGACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTACCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGCtt                                                  >  1:1416600/1‑100 (MQ=255)
                                            tttATCACTGACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTACCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGCTTTg                                                <  2:479200/100‑1 (MQ=255)
                                            tttATCACTGACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTACCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGCTTTg                                                <  2:502864/100‑1 (MQ=255)
                                            tttATCACTGACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTACCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGCTTTg                                                >  1:2117749/1‑100 (MQ=255)
                                            tttATCACTGACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTACCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGCTTTg                                                <  2:70823/100‑1 (MQ=255)
                                            tttATCACTGACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTACCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGCTTTg                                                <  1:1518747/100‑1 (MQ=255)
                                            tttATCACTGACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTACCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGCTTTg                                                <  1:1120808/100‑1 (MQ=255)
                                             ttATCACTGACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTGCCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGCTTTGt                                               <  2:815201/100‑1 (MQ=255)
                                             ttATCACTGACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTACCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGCTTTGt                                               <  2:948110/100‑1 (MQ=255)
                                              tATCACTGACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTACCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGCTTTGtt                                              <  2:2658559/100‑1 (MQ=255)
                                              tATCACTGACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTACCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGCTTTGtt                                              >  2:188906/1‑100 (MQ=255)
                                               aTCACTGACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTACCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGCTTTGtta                                             <  2:1281446/100‑1 (MQ=255)
                                               aTCACTGACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTACCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGCTTTGtta                                             <  1:2922310/100‑1 (MQ=255)
                                               aTCACTGACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTACCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGCTTTGtta                                             >  2:1757894/1‑100 (MQ=255)
                                                tCACTGACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTACCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGCTTTGttat                                            <  1:2832263/100‑1 (MQ=255)
                                                   cTGACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTACCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGCTTTGTTATTAc                                         >  2:2633306/1‑100 (MQ=255)
                                                    tGACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTACCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGCTTTGTTATTAcc                                        <  1:1350016/100‑1 (MQ=255)
                                                     gACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTACCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGCTTTGTTATTAccc                                       >  2:2773345/1‑100 (MQ=255)
                                                     gACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTACCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGCTTTGTTATTAccc                                       <  2:2117780/100‑1 (MQ=255)
                                                     gACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTACCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGCTTTGTTATTAccc                                       <  1:1068095/100‑1 (MQ=255)
                                                     gACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTACCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGCTTTGTTATTAccc                                       >  2:1669629/1‑100 (MQ=255)
                                                     gACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTACCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGCTTTGTTATTAccc                                       >  2:1571409/1‑100 (MQ=255)
                                                        ggAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTACCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGCTTTGTTATTACCCGaa                                    <  2:2380146/100‑1 (MQ=255)
                                                          aaGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTACCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGCTTTGTTATTACCCGaaaa                                  <  1:439203/100‑1 (MQ=255)
                                                             agTCTCTCTGGCAAAAATGTGCACAAAAACTTATTACCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGCTTTGTTATTACCCGAAAAAcc                               >  1:530532/1‑100 (MQ=255)
                                                             agTCTCTCTGGCAAAAATGTGCACAAAAACTTATTACCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGCTTTGTTATTACCCGAAAAAcc                               <  2:1577854/100‑1 (MQ=255)
                                                             agTCTCTCTGGCAAAAATGTGCACAAAAACTTATTACCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGCTTTGTTATTACCCGAAAAAcc                               <  2:1815903/100‑1 (MQ=255)
                                                               tctctcTGGCAAAAATGTGCACAAAAACTTATTACCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGCTTTGTTATTACCCGAAAAACCAg                             >  2:1696645/1‑100 (MQ=255)
                                                                ctctctGGCAAAAATGTGCACAAAAACTTATTACCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGCTTTGTTATTACCCGAAAAACCAGa                            >  2:1746493/1‑100 (MQ=255)
                                                                ctctctGGCAAAAATGTGCACAAAAACTTATTACCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGCTTTGTTATTACCCGAAAAACCAGa                            >  2:2889032/1‑100 (MQ=255)
                                                                 tctctGGCAAAAATGTGCACAAAAACTTATTACCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGCTTTGTTATTACCCGAAAAACCAGaa                           >  2:836711/1‑100 (MQ=255)
                                                                  ctctGGCAAAAATGTGCACAAAAACTTATTACCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGCTTTGTTATTACCCGAAAAACCAGaaa                          >  2:1456033/1‑100 (MQ=255)
                                                                  ctctGGCAAAAATGTGCACAAAAACTTATTACCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGCTTTGTTATTACCCGAAAAACCAGaaa                          >  1:2492237/1‑100 (MQ=255)
                                                                    ctGGCAAAAATGTGCACAAAAACTTATTACCGCAATTGCACATATTCCAGAATGGCGGCAACCATTAATTGCTTTGTTATTACCCGAAAAACCAGAAAtt                        >  1:1851247/1‑100 (MQ=255)
                                                                      ggCAAAAATGTGCACAAAAACTTATTACCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGCTTTGTTATTACCCGAAAAACCAGAAATTGc                      >  2:1916724/1‑100 (MQ=255)
                                                                       gCAAAAATGTGCACAAAAACTTATTACCGCAATTCCACATATTCCATAATGGCGGCAACCATTAATTGCTTTGTTATTACCCGAAAAACCAGAAATTGca                     <  1:990677/100‑1 (MQ=255)
                                                                       gCAAAAATGTGCACAAAAACTTATTACCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGCTTTGTTATTACCCGAAAAACCAGAAATTgaa                     >  1:2108508/1‑98 (MQ=255)
                                                                       gCAAAAATGTGCACAAAAACTTATTACCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGCTTTGTTATTACCCGAAAAACCAGAAATTGca                     <  1:2399062/100‑1 (MQ=255)
                                                                        cgaaaaTGTGCACAAAAACTTATTGCCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGCTTTGTTATTACCCGAAAAACCAGAAATTGcac                    >  1:2969659/3‑100 (MQ=255)
                                                                        cAAAAATGTGCACAAAAACTTATTACCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGCTTTGTTATTACCCGAAAAACCAGAAATTGcac                    <  2:1218388/100‑1 (MQ=255)
                                                                        cAAAAATGTGCACAAAAACTTATTACCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGCTTTGTTATTACCCGAAAAACCAGAAATTGcac                    >  2:2868278/1‑100 (MQ=255)
                                                                        cAAAAATGTGCACAAAAACTTATTACCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGCTTTGTTATTACCCGAAAAACCAGAAATTGca                     <  1:1031291/99‑1 (MQ=255)
                                                                          aaaaTGTGCACAAAAACTTATTACCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGCTTTGTTATTACCCGAAAAACCAGAAATTGcaca                   >  2:1416500/1‑99 (MQ=255)
                                                                           aaaTGTGCACAAAAACTTATTACCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGCTTTGTTATTACCCGAAAAACCAGAAATTGCACATg                 >  2:2150352/1‑100 (MQ=255)
                                                                             aTGTGCACAAAAACTTATTACCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGCTTTGTTATTACCCGAAAAACCAGAAATTGCACATGaa               <  2:3306067/100‑1 (MQ=255)
                                                                              tgtgCACAAAAACTTATTACCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGCTTTGTTATTACCCGAAAAACCAGAAATTGCACATGaaa              >  1:300603/1‑100 (MQ=255)
                                                                                tgCACAAAAACTTATTACCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGCTTTGTTATTACCCGAAAAACCAGAAATTGCACATGAAAtt            >  2:2758955/1‑100 (MQ=255)
                                                                                tgCACAAAAACTTATTACCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGCTTTGTTATTACCCGAAAAACCAGAAATTGCACATGAAAtt            >  2:2465510/1‑100 (MQ=255)
                                                                                 gCACAAAAACTTATTACCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGCTTTGTTATTACCCGAAAAACCAGAAATTGCACATGAAATTg           >  2:1493789/1‑100 (MQ=255)
                                                                                 gCACAAAAACTTATTACCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGCTTTGTTATTACCCGAAAAACCAGAAATTGCACATGAAATTg           <  2:774307/100‑1 (MQ=255)
                                                                                 gCACAAAAACTTATTACCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGCTTTGTTATAACCCGAAAAACCAGAAATTGCACATGAAATTg           >  2:436913/1‑100 (MQ=255)
                                                                                  cacaAAAACTTATTACCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGCTTTGTTATTACCCGAAAAACCAGAAATTGCACATGAAATTg           >  1:2182659/1‑99 (MQ=255)
                                                                                  cacaAAAACTTATTACCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGCTTTGTTATTACCCGAAAAACCAGAAATTGCACATGAAATTGc          >  1:498030/1‑100 (MQ=255)
                                                                                   acaAAAACTTATTACCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGCTTTGTTATTACCCGAAAACCCAGAAATTGCACATGAAATTccc         >  1:1605478/1‑100 (MQ=255)
                                                                                    caAAAACTTATTACCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGCTTTGTTATTACCCGAAAAACCAGAAATTTCACATGAAATTGccc        >  1:1715018/1‑100 (MQ=255)
                                                                                    caAAAACTTATTACCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGCTTTGTTATTACCCGAAAAACCAGAAATTGCACATGAAATTGccc        >  1:2064059/1‑100 (MQ=255)
                                                                                    caAAAACTTATTACCGCAATTCCACATATTCCAGAATAGCGGCAACCATTAATTGCTTTGTTATTACCCGAAAAACCAGAAATTGCACATGAAATTGccc        >  1:202719/1‑100 (MQ=255)
                                                                                      aaaaCTTATTACCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGCTTTGTTATTACCCGAAAAACCAGAAATTGCACATGAAATTGCCCAg      >  2:2959448/1‑100 (MQ=255)
                                                                                      aaaaCTTATTACCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGCTTTGTTATTACCCGAAAAACCAGAAATTGCACATGAAATTGCCCAg      <  2:3313649/100‑1 (MQ=255)
                                                                                       aaaCTTATTACCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGCTTTGTTATTACCCGAAAAACCAGAAATTGCACATGAAATTGccc        >  1:3232908/1‑97 (MQ=255)
                                                                                        aaCTTATTACCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGCTTTGTTATTACCCGAAAAACCAGAAATTGCACATGAAATTGCCCAGCg    <  2:1043637/100‑1 (MQ=255)
                                                                                        aaCTTATTACCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGCTTTGTTATTACCCGAAAAACCAGAAATTGCACATGAAATTGCCCAGCg    <  1:3400147/100‑1 (MQ=255)
                                                                                         aCTTATTACCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGCTTTGTTATTACCCGAAAAACCAGAAATTGCACATGAAATTGCCCAGCGt   <  2:2040531/100‑1 (MQ=255)
                                                                                          cTTATTACCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGCTTTGTTATTACCCGAAAAACCAGAAATTGCACATGAAATTGCCCAGCGTc  <  1:3366719/100‑1 (MQ=255)
                                                                                                |                                                                                             
TACTTGGAGCAATATTCCAGATTTGAATGACGCCTACGTAAACATTTATCACTGACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTGCCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGCTTTGTTATTACCCGAAAAACCAGAAATTGCACATGAAATTGCCCAGCGTC  >  NC_000913/2197269‑2197458

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: