Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,199,587 G→A pseudogene (266/948 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,199,5870GA98.8% 280.6 / ‑4.4 82pseudogene (266/948 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 (0/1);  new base A (36/45);  total (36/46)
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.07e-01

CGAACGATGCGGTTAACCGCTACAAATTGCTGAAAAAAGATGCGCGTACCATCGCCGCCCAGCAGGTAGCAAGGCTGGAATCCGCCATGTGCCTGCGCCGCCGCTGGTCGCTGGAAAACTTCCAGCTCTTCCTGGTTGAGCATCCGCTGGTTCGTCACTTAACCCGCCGTCTGATTTGGGGCGTTTATAGC  >  NC_000913/2199491‑2199681
                                                                                                |                                                                                              
cGAACGATGCGGTTAACCGCTACAAATTGCTGAAAAAAGATGCGCGTACCATCGCCGCCCAGCAGGTAGCAAGGCTGGAATCCGCCATGTGCCTGcaccg                                                                                             >  1:596390/1‑100 (MQ=255)
 gAACGATGCGGTTAACCGCTACAAATTGCTGAAAAAAGATGCGCGTACCATCGCCGCCCAGCAGGTAGCAAGGCTGGAATCCGCCATGTGCCTGcaccgc                                                                                            >  2:2335802/1‑100 (MQ=255)
  aaCGATGCGGTTAACCGCTACAAATTGCTGAAAAAAGATGCGCGTACCATCGCCGCCCAGCAGGTAGCAAGGCTGGAATCCGCCATGTGCCTGcaccgcc                                                                                           <  2:2726693/100‑1 (MQ=255)
  aaCGATGCGGTTAACCGCTACAAATTGCTGAAAAAAGATGCGCGTACCATCGCCGCCCAGCAGGTAGCAAGGCTGGAATCCGCCATGTGCCTGcaccgcc                                                                                           >  1:2426062/1‑100 (MQ=255)
  aaCGATGCGGTTAACCGCTACAAATTGCTGAAAAAAGATGCGCGTACCATCGCCGCCCAGCAGGTAGCAAGGCTGGAATCCGCCATGTGCCTGcaccgcc                                                                                           <  1:1320112/100‑1 (MQ=255)
    cGATGCGGTTAACCGCTACAAATTGCTGAAAAAAGATGCGCGTACCATCGCCGCCCAGCAGGTAGCAAGGCTGGAATCCGCCATGTGCCTGcaccgccgc                                                                                         >  2:510037/1‑100 (MQ=255)
       tGCGGTTAACCGCTACAAATTGCTGAAAAAAGATGCGCGTACCATCGCCGCCCAGCAGGTAGCAAGGCTGGAATCCGCCATGTGCCTGCACCGCCGCTgg                                                                                      <  2:1452871/100‑1 (MQ=255)
       tGCGGTTAACCGCTACAAATTGCTGAAAAAAGATGCGCGTACCATCGCCGCCCAGCAGGTAGCAAGGCTGGAATCCGCCATGTGCCTGCACCGCCGCTgg                                                                                      <  2:2052195/100‑1 (MQ=255)
          ggTTAACCGCTACAAATTGCTGAAAAAAGATGCGCGTACCATCGCCGCCCAGCAGGTAGCAAGGCTGGAATCCGCCATGTGCCTGCACCGCCGCTGGTCg                                                                                   >  1:937873/1‑100 (MQ=255)
           gTTAACCGCTACAAATTGCTGAAAAAAGATGCGCGTACCATCGCCGCCCAGCAGGTAGCAAGGCTGGAATCCGCCATGTGCCTGCACCGCCGCTGGTCGc                                                                                  >  2:1654931/1‑100 (MQ=255)
            ttAACCGCTACAAATTGCTGAAAAAAGATGCGCGTACCATCGCCGCCCAGCAGGTAGCAAGGCTGGAATCCGCCATGTGCCTGCGCCGCCGCTGGTCGCt                                                                                 <  1:668410/100‑1 (MQ=255)
             tAACCGCTACAAATTGCTGAAAAAAGATGCGCGTACCATCGCCGCCCAGCAGGTAGCAAGGCTGGAATCCGCCATGTGCCTGCACCGCCGCTGGTCGCTg                                                                                <  1:720450/100‑1 (MQ=255)
             tAACCGCTACAAATTGCTGAAAAAAGATGCGCGTACCATCGCCGCCCAGCAGGTAGCAAGGCTGGAATCCGCCATGTGCCTGCACCGCCGCTGGTCGCTg                                                                                <  1:129417/100‑1 (MQ=255)
             tAACCGCTACAAATTGCTGAAAAAAGATGCGCGTACCATCGCCGCCCAGCAGGTAGCAAGGCTGGAATCCGCCATGTGCCTGCACCGCCGCTGGTCGCTg                                                                                <  1:1511676/100‑1 (MQ=255)
             tAACCGCTACAAATTGCTGAAAAAAGATGCGCGTACCATCGCCGCCCAGCAGGTAGCAAGGCTGGAATCCGCCATGTGCCTGCACCGCCGCTGGTCGCTg                                                                                <  1:2361210/100‑1 (MQ=255)
             tAACCGCTACAAATTGCTGAAAAAAGATGCGCGTACCATCGCCGCCCAGCAGGTAGCAAGGCTGGAATCCGCCATGTGCCTGCACCGCCGCTGGTCGCTg                                                                                >  1:62686/1‑100 (MQ=255)
                ccGCTACAAATTGCTGAAAAAAGATGCGCGTACCATCGCCGCCCAGCAGGTAGCAAGGCTGGAATCCGCCATGTGCCTGCACCGCCGCTGGTCGCTGGaa                                                                             >  1:223675/1‑100 (MQ=255)
                 cGCTACAAATTGCTGAAAAAAGATGCGCGTACCATCGCCGCCCAGCAGGTAGCAAGGCTGGAATCCGCCATGTGCCTGCACCGCCGCTGGTCGCTGGaaa                                                                            <  2:1546091/100‑1 (MQ=255)
                 cGCTACAAATTGCTGAAAAAAGATGCGCGTACCATCGCCGCCCAGCAGGTAGCAAGGCTGGAATCCGCCATGTGCCTGCACCGCCGCTGGTCGCTGGaaa                                                                            <  1:1440644/100‑1 (MQ=255)
                   cTACAAATTGCTGAAAAAAGATGCGCGTACCATCGCCGCCCAGCAGGTAGCAAGGCTGGAATCCGCCATGTGCCTGCACCGCCGCTGGTCGCTGGAAAAc                                                                          <  2:2275482/100‑1 (MQ=255)
                     aCAAATTGCTGAAAAAAGATGCGCGTACCATCGCCGCCCAGCAGGTAGCAAGGCTGGAATCCGCCATGTGCCTGCACCGCCGCTGGTCGCTGGAAAACtt                                                                        >  2:475221/1‑100 (MQ=255)
                     aCAAATTGCTGAAAAAAGATGCGCGTACCATCGCCGCCCAGCAGGTAGCAAGGCTGGAATCCGCCATGTGCCTGCACCGCCGCTGGTCGCTGGAAAACtt                                                                        >  2:676098/1‑100 (MQ=255)
                       aaaTTGCTGAAAAAAGATGCGCGTACCATCGCCGCCCAGCAGGTAGCAAGGCTGGAATCCGCCATGTGCCTGCACCGCCGCTGGTCGCTGGAAAACTTcc                                                                      <  2:1135015/100‑1 (MQ=255)
                       aaaTTGCTGAAAAAAGATGCGCGTACCATCGCCGCCCAGCAGGTAGCAAGGCTGGAATCCGCCATGTGCCTGCACCGCCGCTGGTCGCTGGAAAACTTcc                                                                      <  1:134790/100‑1 (MQ=255)
                        aaTTGCTGAAAAAAGATGCGCGTACCATCGCCGCCCAGCAGGTAGCAAGGCTGGAATCCGCCATGTGCCTGCACCGCCGCTGGTCGCTGGAAAACTTCCa                                                                     >  2:2566794/1‑100 (MQ=255)
                         aTTGCTGAAAAAAGATGCGCGTACCATCGCCGCCCAGCAGGTAGCAAGGCTGGAATCCGCCATGTGCCTGCACCGCCGCTGGTCGCTGGAAAACTTCCAg                                                                    <  2:2381419/100‑1 (MQ=255)
                           tGCTGAAAAAAGATGCGCGTACCATCGCCGCCCAGCAGGTAGCAAGGCTGGAATCCGCCATGTGCCTGCACCGCCGCTGGTCGCTGGAAGACTTCCAGc                                                                   <  1:2655466/99‑1 (MQ=255)
                           tGCTGAAAAAAGATGCGCGTACCATCGCCGCCCAGCAGGTAGCAAGGCTGGAATCCGCCATGTGCCTGCACCGCCGCTGGTCGCTGGAAAACTTCCAGct                                                                  <  2:1283826/100‑1 (MQ=255)
                           tGCTGAAAAAAGATGCGCGTACCATCGCCGCCCAGCAGGTAGCAAGGCTGGAATCCGCCATGTGCCTGCACCGCCGCTGGTCGCTGGAAAACTTCCAGct                                                                  >  1:1193216/1‑100 (MQ=255)
                           tGCTGAAAAAAGATGCGCGTACCATCGCCGCCCAGCAGGTAGCAAGGCTGGAATCCGCCATGTGCCTGCACCGCCGCTGGTCGCTGGAAAACTTCCAGct                                                                  <  2:2813577/100‑1 (MQ=255)
                            gCTGAAAAAAGATGCGCGTACCATCGCCGCCCAGCAGGTAGCAAGGCTGGAATCCGCCATGTGCCTGCACCGCCGCTGGTCGCTGGAAAACTTCCAGctc                                                                 >  2:2337379/1‑100 (MQ=255)
                               gAAAAAAGATGCGCGTACCATCGCCGCCCAGCAGGTAGCAAGGCTGGAATCCGCCATGTGCCTGCACCGCCGCTGGTCGCTGGAAAACTTCCAGCTCTTc                                                              >  1:1040653/1‑100 (MQ=255)
                               gAAAAAAGATGCGCGGAGGATTGCCACCCCGCCGAAAGCAAGGATGGAATCCTCGATCGTCCCGCACCGCCCCTGGTCGCTGGATCACTTGAAGCTCtt                                                               >  1:202873/1‑99 (MQ=255)
                                 aaaaaGATGCGCGTACCATCGCCGCCCAGCAGGTAGCAAGGCTGGAATCCGCCATGTGCCTGCACCGCCGCTGGTCGCTGGAAAACTTCCAGCTCTTCCt                                                            >  2:757231/1‑100 (MQ=255)
                                   aaaGATGCGCGTACCATCGCCGCCCAGCAGGTAGCAAGGCTGGAATCCGCCATGTGCCTGCACCGCCGCTGGTCGCTGGAAAACTTCTAGCTCTTCCTgg                                                          <  1:2647325/100‑1 (MQ=255)
                                    aaGATGCGCGTACCATCGCCGCCCAGCAGGTAGCAAGGCTGGAATCCGCCATGTGCCTGCACCGCCGCTGGTCGCTGGAAAACTTCCAGCTCTTCCTGGt                                                         <  1:55150/100‑1 (MQ=255)
                                    aaGATGCGCGTACCATCGCCGCCCAGCAGGTAGCAAGGCTGGAATCCGCCATGTGCCTGCACCGCCGCTGGTCGCTGGAAAACTTCCAGCTCTTCCTGGt                                                         <  2:3052172/100‑1 (MQ=255)
                                    aaGATGCGCGTACCATCGCCGCCCAGCAGGTAGCAAGGCTGGAATCCGCCATGTGCCTGCACCGCCGCTGGTCGCTGGAAAACTTCCAGCTCTTCCTGGt                                                         <  2:2787983/100‑1 (MQ=255)
                                    aaGATGCGCGTACCATCGCCGCCCAGCAGGTAGCAAGGCTGGAATCCGCCATGTGCCTGCACCGCCGCTGGTCGCTGGAAAACTTCCAGCTCTTCCTGGt                                                         <  1:1194206/100‑1 (MQ=255)
                                     aGATGCGCGTACCATCGCCGCCCAGCAGGTAGCAAGGCTGGAATCCGCCATGTGCCTGCACCGCCGCTGGTCGCTGGAAAACTTCCAGCTCTTCCTGGtt                                                        <  1:2260092/100‑1 (MQ=255)
                                      ggtGCGCGTACCAACGCCGCCCAGCAGGTAGCAAGGCTGGAATCCGCCATGTGCCTGCACCGCCGCTGGTCGCTGGAAAACTTCCAGCTCTTCCTGGTTg                                                       >  2:2077512/3‑100 (MQ=255)
                                      gATGCGCGTACCATCGCCGCCCAGCAGGTAGCAAGGCTGGAATCCGCCATGTGCCTGCACCGCCGCTGGTCGCTGGAAAACTTCCAGCTCTTCCTGGTTg                                                       <  1:3064388/100‑1 (MQ=255)
                                        tGCGCGTACCATCGCCGCCCAGCAGGTAGCAAGGCTGGAATCCGCCATGTGCCTGCACCGCCGCTGGTCGCTGGAAAACTTCCAGCTCTTCCTGGTTGAg                                                     <  2:62686/100‑1 (MQ=255)
                                        tGCGCGTACCATCGCCGCCCAGCAGGTAGCAAGGCTGGAATCCGCCATGTGCCTGCACCGCCGCTGGTCGCTGGAAAACTTCCAGCTCTTCCTGGTTGAg                                                     <  2:223675/100‑1 (MQ=255)
                                            cgTACCATCGCCGCCCAGCAGGTAGCAAGGCTGGAATCCGCCATGTGCCTGCACCGCCGCTGGTCGCTGGAAAACTTCCAGCTCTTCCTGGTTGAGCATc                                                 >  2:1839058/1‑100 (MQ=255)
                                             gTACCATCGCCGCCCAGCAGGCAGCAAGGCTGGAATCCGCCATGTGCCTGCACCGCCGCTGGTCGCTGGAAAACTTCCAGCTCTTCCTGGTTGAGCATcc                                                <  1:1829582/100‑1 (MQ=255)
                                                ccATCGCCGCCCAGCAGGTAGCAAGGCTGGAATCCGCCATGTGCCTGCACCGCCGCTGGTCGCTGGAAAACTTCCAGCTCTTCCTGGTTGAGCATCCGCt                                             >  2:668196/1‑100 (MQ=255)
                                                ccATCGCCGCCCAGCAGGTAGCAAGGCTGGAATCCGCCATGTGCCTGCACCGCCGCTGGTCGCTGGAAAACTTCCAGCTCTTCCTGGTTGAGCATCCGCt                                             <  2:1619285/100‑1 (MQ=255)
                                                ccATCGCCGCCCAGCAGGTAGCAAGGCTGGAATCCGCCATGTGCCTGCACCGCCGCTGGTCGCTGGAAAACTTCCAGCTCTTCCTGGTTGAGCATCCGCt                                             >  2:2963109/1‑100 (MQ=255)
                                                ccATCGCCGCCCAGCAGGTAGCAAGGCTGGAATCCGCCATGTGCCTGCACCGCCGCTGGTCGCTGGAAAACTTCCAGCTCTTCCTGGTTGAGCATCCGCt                                             >  2:1457279/1‑100 (MQ=255)
                                                 cATCGCCGCCCAGCAGGTAGCAAGGCTGGAATCCGCCATGTGCCTGCACCGCCGCTGGTCGCTGGAAAACTTCCAGCTCTTCCTGGTTGAGCATCCGCTg                                            <  1:2688372/100‑1 (MQ=255)
                                                   tCGCCGCCCAGCAGGTAGCAAGGCTGGAATCCGCCATGTGCCTGCACCGCCGCTGGTCGCTGGAAAACTTCCAGCTCTTCCTGGTTGAGCATCCGCTGGt                                          <  1:1077125/100‑1 (MQ=255)
                                                    cgccgcCCAGCAGGTAGCAAGGCTGGAATCCGTCATGTGCCTGCACCGCCGCTGGTCGCTGGAAAACTGCCAGCTCTTCCTGGTTGTGCATCCGCTGGtt                                         <  2:2057338/100‑1 (MQ=255)
                                                     gccgccCAGCAGGTAGCAAGGCTGGAATCCGCCATGTGCCTGCACCGCCGCTGGTCGCTGGAAAACTTCCAGCTCTTCCTGGTTGAGCATCCGCTGGTTc                                        >  1:1219176/1‑100 (MQ=255)
                                                     gccgccCAGCAGGTAGCAAGGCTGGAATCCGCCATGTGCCTGCACCGCCGCTGGTCGCTGGAAAACTTCCAGCTCTTCCTGGTTGAGCATCCGCTGGTTc                                        >  2:2484080/1‑100 (MQ=255)
                                                        gccCAGCAGGTAGCAAGGCTGGAATCCGCCATGTGCCTGCACCGCCGCTGGTCGCTGGAAAACTTCCAGCTCTTCCTGGTTGAGCATCCGCTGGTTCGTc                                     <  2:2132657/100‑1 (MQ=255)
                                                        gccCAGCAGGTAGCAAGGCTGGAATCCGCCATGTGCCTGCACCGCCGCTGGTCGCTGGAAAACTTCCAGCTCTTCCTGGTTGAGCATCCGCTGGTTCGTc                                     <  2:18772/100‑1 (MQ=255)
                                                         cccAGCAGGTAGCAAGGCTGGAATCCGCCATGTGCCTGCACCGCCGCTGGTCGCTGGAAAACTTCCAGCTCTTCCTGGTTGAGCATCCGCTGGTTCGTCa                                    <  2:1716116/100‑1 (MQ=255)
                                                           cagcagGTAGCAAGGCTGGAATCCGCCATGTGCCTGCACCGCCGCTGGTCGCTGGAAAACTTCCAGCTCTTCCTGGTTGAGCATCCGCTGGTTCGTCa                                    <  1:452732/98‑1 (MQ=255)
                                                           cagcagGTAGCAAGGCTGGAATCCGCCATGTGCCTGCACCGCCGCTGGTCGCTGGAAAACTTCCAGCTCTTCCTGGTTGAGCATCCGCTGGTTCGTCACt                                  <  1:1778624/100‑1 (MQ=255)
                                                           cagcagGTAGCAAGGCTGGAATCCGCCATGTGCCTGCACCGCCGCTGGTCGCTGGAAAACTTCCAGCTCTTCCTGGTTGAGCATCCGCTGGTTCGTCACt                                  >  2:2679926/1‑100 (MQ=255)
                                                             gcagGTAGCAAGGCTGGAATCCGCCATGTGCCTGCACCGCCGCTGGTCGCTGGAAAACTTCCAGCTCTTCCTGGTTGAGCATCCGCTGGTTCGTCACtt                                 <  2:1399408/99‑1 (MQ=255)
                                                              cagGTAGCAAGGCTGGAATCCGCCATGTGCCTGCACCGCCGCTGGTCGCTGGAAAACTTCCAGCTCTTCCTGGTTGAGCATCCGCTGGTTCGTCACTTaa                               >  1:806827/1‑100 (MQ=255)
                                                               agGTAGCAAGGCTGGAATCCGCCATGTGCCTGCACCGCCGCTGGTCGCTGGAAAACTTCCAGCTCTTCCTGGTTGAGCATCCGCTGGTTCGTCACTTAAc                              >  1:397412/1‑100 (MQ=255)
                                                                  tAGCAAGGCTGGAATCCGCCATGTGCCTGCACCGCCGCTGGTCGCTGGAAAACTTCCAGCTCTTCCTGGTTGAGCATCCGCTGGTTCGTCACTTAACccg                           >  1:1308458/1‑100 (MQ=255)
                                                                    gCAAGGCTGGAATCCGCCATGTGCCTGCACCGCCGCTGGTCGCTGGAAAACTTCCAGCTCTTCCTGGTTGAGCATCCGCTGGTTCGTCACTTAACccgcc                         >  2:1527803/1‑100 (MQ=255)
                                                                    gCAAGGCTGGAATCCGCCATGTGCCTGCACCGCCGCTGGTCGCTGGAAAACTTCCAGCTCTTCCTGGTTGAGCATCCGCTGGTTCGTCACTTAACccgcc                         >  2:2124355/1‑100 (MQ=255)
                                                                    gCAAGGCTGGAATCCGCCATGTGCCTGCACCGCCGCTGGTCGCTGGAAAACTTCCAGCTCTTCCTGGTTGAGCATCCGCTGGTTCGTCACTTAACccgcc                         >  1:907713/1‑100 (MQ=255)
                                                                      aaGGCTGGAATCCGCCATGTGCCTGCACCGCCGCTGGTCGCTGGAAAACTTCCAGCTCTTCCTGGTTGAGCATCCGCTGGTTCGTCACTTAACCCGCCGt                       >  1:1836549/1‑100 (MQ=255)
                                                                      aaGGCTGGAATCCGCCATGTGCCTGCACCGCCGCTGGTCGCTGGAAAACTTCCAGCTCTTCCTGGTTGAGCATCCGCTGGTTCGTCACTTAACCCGCCGt                       >  1:157275/1‑100 (MQ=255)
                                                                          cTGGAATCCGCCATGTGCCTGCACCGCCGCTGGTCGCTGGAAAACTTCCAGCTCTTCCTGGTTGAGCATCCGCTGGTTCGTCACTTAACCCGCCGTCTGa                   >  1:2657712/1‑100 (MQ=255)
                                                                             cgatctgCCATGTGCCTGCACCGCCGCTGGTCGCTGGAAAACTTCCAGCTCTTCCTGGTTGAGCATCCGCTGGTTCGTCACTTAACCCGCCGTCTGAttt                <  2:1054147/94‑1 (MQ=255)
                                                                               aTCCGCCATGTGCCTGCACCGCCGCTGGTCGCTGGAAAACTTCCAGCTCTTCCTGGTTGAGCATCCGCGGGTTCGTCACTTAACCCGCCGTCTGATTTgg              <  1:292299/100‑1 (MQ=255)
                                                                                tCCGCCATGTGCCTGCACCGCCGCTGGTCGCTGGAAAACTTCCAGCTCTTCCTGGTTGAGCATCCGCTGGTTCGTCACTTAACCCGCCGTCTGATTTggg             >  2:414193/1‑100 (MQ=255)
                                                                                tCCGCCATGTGCCTGCACCGCCGCTGGTCGCTGGAAAACTTCCAGCTCTTCCTGGTTGAGCATCCGCTGGTTCGTCACTTAACCCGCCGTCTGATTTggg             <  2:366539/100‑1 (MQ=255)
                                                                                 ccGCCATGTGCCTGCACCGCCGCTGGTCGCTGGAAAACTTCCAGCTCTTCCTGGTTGAGCATCCGCTGGTTCGTCACTTAACCCGCCGTCTGATTTgggg            <  2:1836556/100‑1 (MQ=255)
                                                                                  cGCCATGTGCCTGCACCGCCGCTGGTCGCTGGAAAACTTCCAGCTCTTCCTGGTTGAGCATCCGCTGGTTCGTCACTTAACCCGCCGTCTGATTTg               <  1:1784391/96‑1 (MQ=255)
                                                                                  cGCCATGTGCCTGCACCGCCGCTGGTCGCTGGAAAACTTCCAGCTCTTCCTGGTTGAGCATCCGCTGGTTCGTCACTTAACCCGCCGTCTGATTTGGGGc           >  2:921393/1‑100 (MQ=255)
                                                                                  cGCCATGTGCCTGCACCGCCGCTGGTCGCTGGAAAACTTCCAGCTCTTCCTGGTTGAGCATCCGCTGGTTCGTCACTTAACCCGCCGTCTGATTTGGGGc           <  1:757229/100‑1 (MQ=255)
                                                                                   gCCATGTGCCTGCACCGCCGCTGGTCGCTGGAAAACTTCCAGCTCTTCCTGGTTGAGCATCCGCTGGTTCGTCACTTAACCCGCCGTCTGAtttagatcg          >  1:1054143/1‑94 (MQ=255)
                                                                                     cATGTGCCTGCACCGCCGCTGGTCGCTGGAAAACTTCCAGCTCTTCCTGGTTGAGCATCCGCTGGTTCGTCACTTAACCCGCCGTCTGATTTGGGGCGtt        <  1:676097/100‑1 (MQ=255)
                                                                                     cATGTGCCTGCACCGCCGCTGGTCGCTGGAAAACTTCCAGCTCTTCCTGGTTGAGCATCCGCTGGTTCGTCACTTAACCCGCCGTCTGATTTGGGGCGtt        <  1:475219/100‑1 (MQ=255)
                                                                                      aTGTGCCTGCACCGCCGCTGGTCGCTGGAAAACTTCCAGCTCTTCCTGGTTGAGCATCCGCTGGTTCGTCACTTAACCCGCCGTCTGATTTGGGGCGttt       <  1:668195/100‑1 (MQ=255)
                                                                                       tgtgCCTGCACCGCCGCTGGTCGCTGGAAAACTTCCAGCTCTTCCTGGTTGAGCATCCGCTGGTTCGTCACTTAACCCGCCGTCTGATTTGGGGCGTTta      >  2:341849/1‑100 (MQ=255)
                                                                                       tgtgCCTGCACCGCCGCTGGTCGCTGGAAAACTTCCAGCTCTTCCTGGTTGAGCATCCGCTGGTTCGTCAATTAACCCGCCGTCTGATTTGGGGCGTTta      >  1:1498881/1‑100 (MQ=255)
                                                                                          gCCTGCACCGCCGCTGGTCGCTGGAAAACTTCCAGCTCTTCCTGGTTGAGCATCCGCTGGTTCGTCACTTAACCCGCCGTCTGATTTGGGGCGTTTATAg   >  2:136162/1‑100 (MQ=255)
                                                                                          gCCTGCACCGCCGCTGGTCGCTGGAAAACTTCCAGCTCTTCCTGGTTGAGCATCCGCTGGTTCGTCACTTAACCCGCCGTCTGATTTGGGGCGTTTATAg   <  1:2963090/100‑1 (MQ=255)
                                                                                           ccTGCACCGCCGCTGGTCGCTGGAAAACTTCCAGCTCTTCCTGGTTGAGCATCCGCTGGTTCGTCACTTAACCCGCCGTCTGATTTGGGGCGTTTATAgc  >  2:474876/1‑100 (MQ=255)
                                                                                                |                                                                                              
CGAACGATGCGGTTAACCGCTACAAATTGCTGAAAAAAGATGCGCGTACCATCGCCGCCCAGCAGGTAGCAAGGCTGGAATCCGCCATGTGCCTGCGCCGCCGCTGGTCGCTGGAAAACTTCCAGCTCTTCCTGGTTGAGCATCCGCTGGTTCGTCACTTAACCCGCCGTCTGATTTGGGGCGTTTATAGC  >  NC_000913/2199491‑2199681

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: