Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,197,365 G→A 72.0% 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,3650GA72.0% 165.9 / 75.9 110pseudogene (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 (11/20);  new base A (56/23);  total (67/43)
Fisher's exact test for biased strand distribution p-value = 9.88e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 111 total aligned reads displayed.

CTACTTGGAGCAATATTCCAGATTTGAATGACGCCTACGTAAACATTTATCACTGACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTGCCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGCTTTGTTATTACCCGAAAAACCAGAAATTGCACATGAAATTGCCCAGCGTCTACTGG  >  NC_000913/2197268‑2197464
                                                                                                 |                                                                                                   
cTACTTGGAGCAATATTCCAGATTTGAATGACGCCTACGTAAACATTTATCACTGACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTAcc                                                                                                   <  1:879232/100‑1 (MQ=255)
cTACTTGGAGCAATATTCCAGATTTGAATGACGCCTACGTAAACATTTATCACTGACGGAAGAGAAAAAAAAAAAAAAATGTGCACAAAAACTTATTGcc                                                                                                   <  1:1232420/100‑1 (MQ=255)
  aCTTGGAGCAATATTCCAGATTTGAATGACGCCTACGTAAACATTTATCACTGACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTACCGc                                                                                                 <  2:2908882/100‑1 (MQ=255)
  aCTTGGAGCAATATTCCAGATTTGAATGACGCCTACGTAAACATTTATCACTGACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTACCGc                                                                                                 <  1:2271459/100‑1 (MQ=255)
  aCTTGGAGCAATATTCCAGATTTGAATGACGCCTACGTAAACATTTATCACTGACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTACCGc                                                                                                 >  1:548075/1‑100 (MQ=255)
   cTTGGAGCAATATTCCAGATTTGAATGACGCCTACGTAAACATTTATCACTGACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTGCCGCa                                                                                                <  1:333647/100‑1 (MQ=255)
       gAGCAATATTCCAGATTTGAATGACGCCTACGTAAACATTTATCACTGACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTACCGCAATTc                                                                                            >  1:2381203/1‑100 (MQ=255)
       gAGCAATATTCCAGATTTGAATGACGCCTACGTAAACATTTATCACTGACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTACCGCAATTc                                                                                            >  2:2800364/1‑100 (MQ=255)
        aGCAATATTCCAGATTTGAATGACGCCTACGTAAACATTTATCACTGACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTACCGCAATTcc                                                                                           >  1:1291950/1‑100 (MQ=255)
         gCAATATTCCAGATTTGAATGACGCCTACGTAAACATTTATCACTGACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTACCGCAATTCca                                                                                          >  2:810766/1‑100 (MQ=255)
         gCAATATTCCAGATTTGAATGACGCCTACGTAAACATTTATCACTGACGGAAGAGGCTCTCTGGCAAAAATCTGAACAAAAACTTATTACCGCAATTc                                                                                            >  1:1257611/1‑98 (MQ=255)
          cAATATTCCAGATTTGAATGACGCCTACGTAAACATTTATCACTGACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTGCCGCAATTCcac                                                                                         >  2:1211329/1‑100 (MQ=255)
           aaTATTCCAGATTTGAATGACGCCTACGTAAACATTTATCACTGACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTGCCGCAATTCcaca                                                                                        <  2:2801982/100‑1 (MQ=255)
           aaTATTCCAGATTTGAATGACGCCTACGTAAACATTTATCACTGACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTACCGCAATTCcaca                                                                                        <  2:2633729/100‑1 (MQ=255)
            atatTCCAGATTTGAATGACGCCTACGTAAACATTTATCACTGACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTACCGCAATTCCACat                                                                                       <  2:69609/100‑1 (MQ=255)
            atatTCCAGATTTGAATGACGCCTACGTAAACATTTATCACTGACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTACCGCAATTCCACat                                                                                       >  2:55689/1‑100 (MQ=255)
               ttCCAGATTTGAATGACGCCTACGTAAACATTTATCACTGACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTGCCGCAATTCCACATAtt                                                                                    >  1:2946561/1‑100 (MQ=255)
                tCCAGATTTGAATGACGCCTACGTAAACATTTATCACTGACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTACCGCAATTCCACATATTc                                                                                   >  1:318779/1‑100 (MQ=255)
                 ccAGATTTGAATGACGCCTACGTAAACATTTATCACTGACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTGCCGCAATTCCACATAttac                                                                                  <  1:1967213/100‑3 (MQ=255)
                 ccAGATTTGAATGACGCCTACGTAAACATTTATCACTGACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTGCCGCAATTCCACATATTcc                                                                                  >  2:3261460/1‑100 (MQ=255)
                 ccAGATTTGAATGACGCCTACGTAAACATTTATCACTGACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTACCGCAATTCCACATATTcc                                                                                  >  2:2649703/1‑100 (MQ=255)
                 ccAGATTTGAATGACGCCTACGTAAACATTTATCACTGACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTACCGCAATTCCACATATTcc                                                                                  >  1:1158439/1‑100 (MQ=255)
                  cAGATTTGAATGACGCCTACGTAAACATTTATCACTGACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTACCGCAATTCCACATATTCCa                                                                                 <  2:1326940/100‑1 (MQ=255)
                    gATTTGAATGACGCCTACGTAAACATTTATCACTGACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTGCCGCAGTTCCACATATTCCAGa                                                                               <  1:365543/100‑1 (MQ=255)
                      tttGAATGACGCCTACGTAAACATTTATCACTGACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTACCGCAATTCCACATATTCCAGAAt                                                                             >  2:64667/1‑100 (MQ=255)
                         gAATGACGCCTACGTAAACATTTATCACTGACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTACCGCAATTCCACATATTCCAGAATggc                                                                          >  1:2698416/1‑100 (MQ=255)
                         gAATGACGCCTACGTAAACATTTATCACTGACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTACCGCAATTCCACATATTCCAGAATggc                                                                          >  2:1220259/1‑100 (MQ=255)
                           aTGACGCCTACGTAAACATTTATCACTGACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTGCCGCAATTCCACATATTCCAGAATggcgg                                                                        >  2:114522/1‑100 (MQ=255)
                           aTGACGCCTACGTAAACATTTATCACTGACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTGCCGCAATTCCACATATTCCAGAATggcgg                                                                        <  2:159121/100‑1 (MQ=255)
                              aCGCCTACGTAAACATTTATCACTGACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTACCGCAATTCCACATATTCCAGAATGGCGGCaa                                                                     >  2:2759079/1‑100 (MQ=255)
                               cGCCTACGTAAACATTTATCACTGACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTACCGCAATTCCACATATTCCAGAATGGCGGCAAc                                                                    >  1:3176834/1‑100 (MQ=255)
                               cGCCTACGTAAACATTTATCACTGACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTACCGCAATTCCACATATTCCAGAATGGCGGCAAc                                                                    >  1:1950107/1‑100 (MQ=255)
                                 ccTACGTAAACATTTATCACTGACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTACCGCAATTCCACATATTCCAGAATGGCGGCAACCa                                                                  >  2:2524335/1‑100 (MQ=255)
                                 ccTACGTAAACATTTATCACTGACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTACCGCAATTCCACATATTCCAGAATGGCGGCAACCa                                                                  >  1:1605670/1‑100 (MQ=255)
                                 ccTACGTAAACATTTATCACTGACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTACCGCAATTCCACATATTCCAGAATGGCGGCAACCa                                                                  <  1:39395/100‑1 (MQ=255)
                                     cGTAAACATTTATCACTGACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTACCGCAATTCCACATATTCCAGAATGGCGGCAACCATTaa                                                              >  2:1034217/1‑100 (MQ=255)
                                     cGTAAACATTTATCACTGACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTACCGCAATTCCACATATTCCAGAATGGCGGCAACCATTaa                                                              >  2:1429318/1‑100 (MQ=255)
                                      gaacaaaattACCACTGACAAAAGAGAAACTCTGGCAAAAATGTGCACACAAACTTATTACCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAAt                                                             >  1:1530142/9‑100 (MQ=255)
                                         aaCATTTATCACTGACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTACCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGc                                                          <  1:2383208/100‑1 (MQ=255)
                                         aaCATTTATCACTGACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTACCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGc                                                          >  2:1835533/1‑100 (MQ=255)
                                          aCATTTATCACTGACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTACCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGCt                                                         >  1:573307/1‑100 (MQ=255)
                                          aCATTTATCACTGACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTACCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGCt                                                         >  2:495858/1‑100 (MQ=255)
                                          aCATTTATCACTGACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTACCGCAATTCCACATATTCCAGAATGGCGGCAACCATGAATTGCt                                                         >  1:1818350/1‑100 (MQ=255)
                                           cATTTATCACTGACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTACCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGCtt                                                        >  2:223977/1‑100 (MQ=255)
                                            aTTTATCACTGACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTGCCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGCttt                                                       <  2:2123666/100‑1 (MQ=255)
                                             tttATCACTGACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTGCCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGCTTTg                                                      >  1:2628611/1‑100 (MQ=255)
                                             tttATCACTGACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTACCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGCTTTg                                                      <  2:497232/100‑1 (MQ=255)
                                             tttATCACTGACGGAAGAGGCTCTGTGGCAAATATGTGCATAAAGACTTATTGCCGCAATTCCAGATATTCCAGAATGGCGGCAACCATTAATTGAttt                                                       >  1:2122379/1‑99 (MQ=255)
                                               tATCACTGACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTACCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGCTTTGtt                                                    <  1:740641/100‑1 (MQ=255)
                                               tATCACTGACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTACCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGCTTTGtt                                                    <  2:1056934/100‑1 (MQ=255)
                                                 tCACTGACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTACCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGCTTTGttat                                                  <  2:864452/100‑1 (MQ=255)
                                                   aatGAAGGAAGAGTCTCTCTGGCAAAAATGTGCAAAAAAAATTATAGCCGCAATTCCACATATTCCAGAATGACGGCACCCATTAATTGCTTTGttatta                                                <  1:2195374/98‑1 (MQ=255)
                                                      gACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTACCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGCTTTGTTATTAc                                               >  2:3271757/1‑98 (MQ=255)
                                                      ctcGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTACCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGCTTTGTTATTAccc                                             <  2:2711036/98‑1 (MQ=255)
                                                         ggAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTGCCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGCTTTGTTATTACCCGaa                                          <  1:1956907/100‑1 (MQ=255)
                                                          gAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTACCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGCTTTGTTATTACCCGaaa                                         >  2:652772/1‑100 (MQ=255)
                                                              agTCTCTCTGGCAAAAATGTGCACAAAAACTTATTGCCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGCTTTGTTATTACCCGAAAAAcc                                     <  1:50862/100‑1 (MQ=255)
                                                               gTCTCTCTGGCAAAAATGTGCACAAAAACTTATTGCCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGCTTTGTTATTACCCGAAAAACCa                                    <  2:2243682/100‑1 (MQ=255)
                                                               gTCTCTCTGGCAAAAATGTGCACAAAAACTTATTACCGCAATTCCACATATTCCACAATGGCGGCAACCATTAATTGCTTTGTTATTACCCGAAAAAcc                                     >  1:1115065/1‑99 (MQ=255)
                                                                tctctcTGGCAAAAATGTGCACAAAAACTTATTACCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGCTTTGTTATTACTCGAAGAAc                                      >  1:179063/1‑97 (MQ=255)
                                                                tctctcTGGCAAAAATGTGCACAAAAACTTATTACCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGCTTTGTTATTACCCGAAAAACCAg                                   <  2:1603968/100‑1 (MQ=255)
                                                                 ctctctGGCAAAAATGTGCACAAAAACTTATTACCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGCTTTGTTATTACCCGAAAAACCAGa                                  >  1:2351369/1‑100 (MQ=255)
                                                                 ctctctGGCAAAAATGTGCACAAAAACTTATTACCGCAATTCCACATATTCCAGAATGGCGGACACAATTAATTGCTTTGTTCTTACACGAAAAACCAg                                   >  1:2015549/1‑99 (MQ=255)
                                                                   ctctGGCAAAAATGTGCACAAAAACTTATTACCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGCTTTGTTATTACCCGAAAAACCAGaaa                                >  2:40047/1‑100 (MQ=255)
                                                                      tGGCAAAAATGTGCACAAAAACTTATTACCGCAATTCCACATATTCCAGAATGGCGGCAACCGTTAATTGCTTTGTTATTACCCGAAAAACCAGAAATTg                             >  1:2278198/1‑100 (MQ=255)
                                                                      tGGCAAAAATGTGCACAAAAACTTATTACCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGCTTTGTTATTACCCGAAAAACCAGAAATTg                             >  2:1597693/1‑100 (MQ=255)
                                                                       ggCAAAAATGTGCACAAAAACTTATTGCCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGCTTTGTTATTACCCGAAAAACCAGAAATTGc                            <  1:428564/100‑1 (MQ=255)
                                                                       ggCAAAAATGTGCACAAAAACTTATTACCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGCTTTGTTATTACCCGAAAAACCAGAAATTGc                            >  2:2614693/1‑100 (MQ=255)
                                                                        gCAAAAATGTGCACAAAAACTTATTACCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGCTTTGTTATTACCCGAAAAACCAGAAATTGca                           >  1:135796/1‑100 (MQ=255)
                                                                        gCAAAAATGTGCACAAAAACTTATTACCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGCTTTGTTATTACCCGAAAAACCAGAAATTGca                           <  2:2853790/100‑1 (MQ=255)
                                                                        gCAAAAATGTGCACAAAAACTTATTACCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGCTTTGTTATTACCCGAAAAACCAGAAATTGca                           >  2:3288543/1‑100 (MQ=255)
                                                                         cAAAAATGTGCACAAAAACTTATTGCCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGCTTTGTTATTACCCGAAAAACCAGAAATTGcac                          >  1:2839021/1‑100 (MQ=255)
                                                                         cAAAAATGTGCACAAAAACTTATTACCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGCTTTGTTATTACCCGAAAAACCAGAAATTGcac                          <  2:548078/100‑1 (MQ=255)
                                                                         cAAAAATGTGCACAAAAACTTATTACCGCAATGCCACATATTCCAGAATGGCGGCAACCATTAATTCCTTTGTTATTACCCGAAAAACCAGAAATTGcac                          >  2:3212974/1‑100 (MQ=255)
                                                                             aaTGTGCACAAAAACTTATTACCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGCTTTGTTATTACCCGAAAAACCAGAAATTGCACATGa                      <  2:1878381/100‑1 (MQ=255)
                                                                              aTGTGCACAAAAACTTATTGCCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGCTTTGTTATTACCCGAAAAACCAGAAATTGCACATGaa                     <  2:2493423/100‑1 (MQ=255)
                                                                                  gCACAAAAACTTATTGCCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGCTTTGTTATTACCCGAAAAACCAGAAATTGCACATGAAATTg                 >  2:765954/1‑100 (MQ=255)
                                                                                  gCACAAAAACTTATTACCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGCTTTGTTATTACCCGAAAAACCAGAAATTGCACATGAATTTg                 >  1:2770462/1‑100 (MQ=255)
                                                                                  gCACAAAAACTTATTACCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGCTTTGTTATTACCCGAAAAACCAGAAATTGCACATGAAATTg                 >  2:1431534/1‑100 (MQ=255)
                                                                                   cacaAAAACTTATTGCCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGCTTTGTTATTACCCGAAAAACCAGAAATTGCACATGAAATTGc                <  2:2759248/100‑1 (MQ=255)
                                                                                   cacaAAAACTTATTACCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGCTTTGTTATTACCCGAAAAACCAGAAATTGCACATGAAATTg                 >  1:3173241/1‑99 (MQ=255)
                                                                                    acaAAAACTTATTGCCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGCTTTGTTATTACCCGAAAAACCAGAAATTGCACATGAAATTGcc               <  2:2814104/100‑1 (MQ=255)
                                                                                    acaAAAACTTATTACCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGCTTTGTTATTACCCGAAAAACCAGAAATTGCACATGAAATTGcc               <  1:2615757/100‑1 (MQ=255)
                                                                                    acaAAAACTTATTACCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGCTTTGTTATTACCCGAAAAACCAGAAATTGCACATGAAATTGcc               <  1:2877163/100‑1 (MQ=255)
                                                                                    acaAAAACTTATTACCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGCTTTGTTATTACCCGAAAAACCAGAAATTGCACATGAAATTGcc               >  1:981799/1‑100 (MQ=255)
                                                                                    acaAAAACTTATTACCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGCTTTGTTATTACCCGAAAAACCAGAAATTGCACATGAAATTGcc               >  1:2628096/1‑100 (MQ=255)
                                                                                    acaAAAACTTATTACCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGCTTTGTTATTACCCGAAAAACCAGAAATTGCACATGAAATTGcc               <  2:2153420/100‑1 (MQ=255)
                                                                                    acaAAAACTTATTACCGCAATCCCACATATTCCAGAATGGCGGCAACCATTAATTGCTTTGTTATTAACCGAAAAACCAGACATTGCACATGAAATTGcc               >  2:687702/1‑100 (MQ=255)
                                                                                     caAAAACTTATTACCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGCTTTGTTATTACCCGAAAAACCAGAAATTGCACATGAAATTGccc              >  2:849226/1‑100 (MQ=255)
                                                                                     caAAAACTTATTACCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGCTTTGTTATTACCCGAAAAACCAGAAATTGCACATGAAATTGccc              >  2:2584244/1‑100 (MQ=255)
                                                                                       aaaaCTTATTGACGCAATTCAACATATTCCAGAATGGCGGCAACCATTAATTGCTTTGTTATTAACCGAAAAACCAGAAATTGAACATGAAATTGACCAg            >  1:2015832/1‑100 (MQ=255)
                                                                                       aaaaCTTATTACCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGCTTTGTTATTACCCGAAAAACCAGAAATTGCACATGAAATTGCCCAg            <  2:920730/100‑1 (MQ=255)
                                                                                       aaaaCTTATTACCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGCTTTGTTATTACCCGAAAAACCAGAAATTGCACATGAAATTGCCCAg            <  1:2936413/100‑1 (MQ=255)
                                                                                       aaaaCTTATTACCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGCCTTGTTATTACCCGAAAAACCAGAAATTGCACATGAAATTGCCCa             >  1:979162/1‑99 (MQ=255)
                                                                                            ttATTGCCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGCTTTGTTATTACCCGAAAAACCAGAAATTGCACATGAAATTGCCCAGCGTCt       <  1:114522/100‑1 (MQ=255)
                                                                                              aTTGCCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGCTTTGTTATTACCCGAAAAACCAGAAATTGCACATGAAATTGCCCAGCGTCTAc     <  1:79219/100‑1 (MQ=255)
                                                                                              aTTACCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGCTTTGTTATTACCCGAAAAACCAGAAATTGCACATGAAATTGCCCAGCGTCTAc     <  1:1429312/100‑1 (MQ=255)
                                                                                               ttGCCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGCTTTGTTATTACCCGAAAAACCAGAAATTGCACATGAAATTGCCCAGCGTCTACt    <  1:847207/100‑1 (MQ=255)
                                                                                                tGCCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGCTTTGTTATTACCCGAAAAACCAGAAATTGCACATGAAATTGCCCAGCGTCTACTg   >  2:2873321/1‑100 (MQ=255)
                                                                                                 gCCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGCTTTGTTATTACCCGAAAAACCAGAAATTGCACATGAAATTGCCCAGCGTCTACTgg  <  2:2946587/100‑1 (MQ=255)
                                                                                                 |                                                                                                   
CTACTTGGAGCAATATTCCAGATTTGAATGACGCCTACGTAAACATTTATCACTGACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTGCCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGCTTTGTTATTACCCGAAAAACCAGAAATTGCACATGAAATTGCCCAGCGTCTACTGG  >  NC_000913/2197268‑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: