Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,199,608 A→T 8.3% pseudogene (287/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,6080AT 8.3% 234.4 / 9.7 93pseudogene (287/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 A (42/42);  new base T (2/6);  total (45/48)
Fisher's exact test for biased strand distribution p-value = 2.71e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.96e-01

AAACCGAACAAAAGCGATGATGAAACGCGGGCGAACGATGCGGTTAACCGCTACAAATTGCTGAAAAAAGATGCGCGTACCATCGCCGCCCAGCAGGTAGCAAGGCTGGAATCCGCCATGTGCCTGCGCCGCCGCTGGTCGCTGGAAAACTTCCAGCTCTTCCTGGTTGAGCATCCGCTGGTTCGTCACTTAACCCGCCGTCTGATTTGGGGCGTTTATAGCGCCGAAAACCAGCTACTGGCTTGCTTTCGCGTAGCAGAAGATAACAGCTCCAGCACCGCTGACGATGATCTTTTCA  >  NC_000913/2199460‑2199757
                                                                                                                                                    |                                                                                                                                                     
aaaCCGAACAAAAGCGATGATTAAACGCGGGCGAACGATGCGGTTAACCGCTACAAATTGGTGAAAAAAGATGCGCGTACCATCGCCGCCCAGCAGGTAGCAAGGCTGGAATCCGCCATGTGCCTGCGCCGCCGCTGGTCTCTGGAAAAc                                                                                                                                                      <  2:361667/150‑1 (MQ=255)
aaaCCGAACAAAAGCGATGATGAAACGCGGGCGAACCATGCGGTTAACCGCTACAAATTGCTGCAAAAAGATGCGCGTACCATCGCCGCCCAGCAGGTAGCAAGGCTGGAATCCGCCATGTGCCTGCGCCGCCGCTGGTCGCTGGAAAAc                                                                                                                                                      <  2:799159/150‑1 (MQ=255)
aaaCCGAACAAAAGCGATGAGGAAACGCGGGCGAACGATGCGGTTAACCGCTACAAATTGCTGAAAAAAGATGCGCGTACCATCGCCGCCCAGCAGGTAGCAAGGCTGGAATCCGCCATGTGCCTGCGCCGCCGCTGGTCGCTGGAAAAc                                                                                                                                                      <  1:769912/150‑1 (MQ=255)
   ccGAACAAAAGCGATGATGAAACGCGGGCGAACGATGCGGTTAACCGCTACAAATTGCTGAAAAAAGATGCGCGTACCATCGCCGCCCAGCAGGTAGCAAGGCTGGAATCCGCCATGTGCCTGCGCCGCCGCTGGTCGCTGGAAAACTTc                                                                                                                                                   <  2:1118689/150‑1 (MQ=255)
     gAACAAAAGCGATGATGAAACGCGGGCGAACGATGCGGTTAACCGCTACAAATTGCTGAAAAAAGATGCGCGTACCATCGCCGCCCAGCAGGTAGCAAGGCTGGAATCCGCCATGTGCCTGCGCCGCCGCTGGTCGCTGGAAAACTTCCa                                                                                                                                                 >  1:794984/1‑150 (MQ=255)
       aCAAAAGCGATGATGAAACGCGGGCGAACGATGCGGTTAACCGCTACAAATTGCTGAAAAAAGATGCGCGTACCATCGCCGCCCAGCAGGTAGCAAGGCTGGAATCCGCCATGTGCCTGCGCCGCCGCTGGTCGCTGGAAAACTTCCAGc                                                                                                                                               >  1:1875986/1‑150 (MQ=255)
        cAAAAGCGATGATGAAACGCGGGCGAACGATGCGGTTAACCGCTACAAATTGCTGAAAAAAGATGCGCGTACCATCGCCGCCCAGCAGGTAGCAAGGCTGGAATCCGCCATGTGCCTGCGCCGCCGCTGGTCGCTGGAAGTTTTCCAGcg                                                                                                                                              <  1:842206/150‑2 (MQ=255)
        cAAAAGCGATGATGAAACGCGGGCGAACGATGCGGTTAACCGCTACAAATTGCTGAAAAAAGATGCGCGTACCATCGCCGCCCAGCAGGTAGCAAGGCTGGAATCCGCCATGTGCCTGCGCCGCCGCTGGTCGCTGGAAAACTTCCAGct                                                                                                                                              >  2:1853119/1‑150 (MQ=255)
         aacaGCGATGATGAAACGCGGGCGAACGATGCAGTTAACCGCTACAAATTGCTGAAACAAGATGCGCGTACCATCGCCGCCCAGCAGGTAGCAAGGCTGGAATCCGCCATGTGCCTGCGCCGCCGCTGGGCGCTGGAAGTTTTCCAGcga                                                                                                                                             <  2:1244581/147‑3 (MQ=255)
         aaaaGCGATGATGAAACGCGGGCGAACGATGCGGTTAACCGCTACAAATTGCTGAAAAAAGATGCGCGTACCATCGCCGCCCAGCAGGTAGCAAGGCTGGAATCCGCCATGTGCCTGCGCCGCCGCTGGTCGCTGGAAAACTTCCAGctc                                                                                                                                             >  2:787142/1‑150 (MQ=255)
          aacgCGATGAGGAAACGCGGGCGAACGATGCGGGTAACCGCTACAAATTGCTGAAAAAAGATGCGCGTACCATCGCCGCCCAGCAGGTAGCAAGGCTGGAATCCGCCATGTGCCTGCGCCGCCGCTGGTCGCTGGAAGTTTTCCAGcgac                                                                                                                                            <  1:546962/147‑4 (MQ=255)
          aaaGCGATGATGAAACGCGGGCGAACGATGCGGTTAACCGCTACAAATTGCTGAAAAAAGATGCGCGTACCATCGCCGCCCAGCAGGTAGCAAGGCTGGAATCCGCCATGTGCCTGCGCCGCCGCTGGTCGCTGGAAGTTTTCCAGcgac                                                                                                                                            <  1:1473872/150‑4 (MQ=255)
           aaGCGATGATGAAACGCGGGCGAACGATGCGGTTAACCGCTACAAATTGCTGAAAAAAGATGCGCGTACCATCGCCGCCCAGCAGGTAGCAAGGCTGGAATCCGCCATGTGCCTGCGCCGCCGCTGGTCGCTGGAAAACTTCCAGCTCtt                                                                                                                                           >  1:433638/1‑150 (MQ=255)
             gCGATGATGAAACGCGGGCGAACGATGCGGTTAACCGCTACAAATTGCTGAAAAAAGATGCGCGGACCATCGCCGCCCAGCAGGTAGCAAGGCTGGAATCCGCCATGGGCCTGCGCCGCCGCTGGTCGCTGGAAGTTTTCCAGcgaccag                                                                                                                                         <  2:1441440/150‑7 (MQ=255)
              cGATGATGAAACGCGGGCGAACGATGCGGGTAACCGCTACAAATTGCTGAAAAAAGATGCGCGGACCATCGCCGCCCAGCAGGTAGCAAGGCTGGAATCCGCCATGTGCCGGCGCCGCCGCTGGTCGCTGGAAGTTTTCCAGcgaccagc                                                                                                                                        <  1:1529171/150‑8 (MQ=255)
                  gatgaAACGCGGGCGAACGATGCGGTTAACCGCTACAAATTGCTGAAAAAAGATGCGCGTACCATCGCCGCCGAGCAGGTAGCAAGGCTGGAATCCGCCATGTGCCTGCGCCGCCGCTGGTCGCTGGAAAACTTCCAGCTCTTCCTGGtt                                                                                                                                    <  1:289572/150‑1 (MQ=255)
                   atgaAACGCGGGCGAACGATGCGGTTAACCGCTACAAATTGCTGAAAAAAGATGCGCGTACCATCGCCGCCCAGCAGGTAGCAAGGCTGGAATCCGCCATGTGCCTGCGCCGCCGCTGGTCGCTGGAAAACTTCCAGCTCTTCCTGGTTg                                                                                                                                   <  1:1233437/150‑1 (MQ=255)
                     gaAACGCGGGCGAACGATGCTGTTAACCGCTACAAATTGCTGAAAAAAGATGCGCGTACCATCGCCGCCCAGCAGGTAGCAAGGCTGGAATCCGCCATGTGCCTGCGCCGCCGCTGGTCGCTGGAAAACTTCCAGCTCTTCCTGGTTGAg                                                                                                                                 <  1:1006513/150‑1 (MQ=255)
                     gaAACGCGGGCGAACGATGCGGTTAACCGCTACAAATTGCTGAAAAAAGATGCGCGTACCATCGCCGCCCAGCAGGTAGCAAGGCTGGAATCCGCCATGTGCCTGCGCCGCCGCTGGTCGCTGGAAAACTTCCAGCTCTTCCTGGTTGAg                                                                                                                                 <  2:1990165/150‑1 (MQ=255)
                       aaCGCGGGCGAACGATGCGGTTAACCGCTACAAATTGCTGAAAAAAGATGCGCGTACCATCGCCGCCCAGCAGGTAGCAAGGCTGGAATCCGCCATGTGCCTGCGCCGCCGCTGGTCGCTGGAAAACTTCCAGCTCTTCCTGGTTGAGCa                                                                                                                               <  2:863019/150‑1 (MQ=255)
                        aCGCGGGCGAACGATGCGGTTAACCGCTACAAATTGCTGAAAAAAGATGCGCGTACCATCGCCGCCCAGCAGGTAGCAAGGCTGGAATCCGCCATGTGCCTGCGCCGCCGCTGGTCGCTGGAAAACTTCCAGCTCTTCCTGGTTGAGCAt                                                                                                                              >  2:1321665/1‑150 (MQ=255)
                             ggCGAACGATGCGGTTAACCGCTACAAATTGCTGAAAAAAGATGCGCGTACCATCGCCGCCCAGCAGGTAGCAAGGCTGGAATCCGCCATGTGCCTGCGCCGCCGCTGGTCGCTGGAAAACTTCCAGCTCTTCCTGGTTGAGCATCCGCt                                                                                                                         <  2:1613980/150‑1 (MQ=255)
                              gCGAACGATGCGGTTAACCGCTACAAATTGCTGAAAAAAGATGCGCGTACCATCGCCGCCCAGCAGGTAGCAAGGCTGGAATCCGCCATGTGCCTGCGCCGCCGCTGGTAGCTGGAAAACTTCCAGCTCGTCCTGGTTGAGCATCCGCTg                                                                                                                        >  2:1563892/1‑150 (MQ=255)
                               cGAACGATGCGGTTAACCGCTACAAATTGCTGAAAAAAGATGCGCGTACCATCGCCGCCCAGCTGGTAGCAAGGCTGGAATCCGCCATGTGCCTGCGCCGCCGCTGGTCGCTGGAAAACTTCCAGCTCTTCCTGGTTGAGCATCCGCTgg                                                                                                                       >  2:1858687/1‑150 (MQ=255)
                               cGAACGATGCGGTTAACCGCTACAAATTGCTGAAAAAAGATGCGCGTACCATCGCCGCCCAGCAGGTAGCAAGGCTGGAATCCGCCATGTGCCTGCGCCGCCGCTGGTCGCTGGAAAACTTCCAGCTCTTCCTGGTTGAGCATCCGCTgg                                                                                                                       >  2:1787612/1‑150 (MQ=255)
                               cGAACGATGCGGTTAACCGCTACAAATTGCTGAAAAAAGATGCGCGTACCATCGCCGCCCAGCAGGTAGCAAGGCTGGAATCCGCCATGTGCCTGCGCCGCCGCTGGTCGCTGGAAAACTTCCAGCTCTTCCTGGTTGAGCATCCGCTgg                                                                                                                       <  2:1950630/150‑1 (MQ=255)
                                     aTGCGGTTAACCGCTACAAATTGCTGAAAAAAGATGCGCGTACCATCGCCGCCCAGCAGGTAGCAAGGCTGGAATCCGCCATGTGCCTGCGCCGCCGCTGGTCGCTGGAAAACTTCCAGCTCTTCCTGGTTGAGCATCCGCTGGTTCGTc                                                                                                                 >  2:1647234/1‑150 (MQ=255)
                                     aTGCGGTTAACCGCTACAAATTGCTGAAAAAAGATGCGCGTACCATCGCCGCCCAGCAGGTAGCAAGGCTGGAATCCGCCATGTGCCTGCGCCGCCGCTGGTCGCTGGAAAACTTCCAGCTCTTCCTGGTTGAGCATCCGCTGGTTCGTc                                                                                                                 <  2:490162/150‑1 (MQ=255)
                                         ggTTAACCGCTTCAAATTGCTGAAAAAAGAAGCGCGTACCATCGCCGCCCAGCAGGTAGCAAGGCTGGAATCCGCGATGTGCCTGCGCCGCCGCTGGGCGCTGGAAAACTTCCAGCTCTTCCTGGTTGAGCATCCGCTGGTTCGTCACtt                                                                                                             <  2:1248507/150‑1 (MQ=255)
                                             aaCCGCTACAAATTGCTGAAAAAAGATGCGCGTACCATCGCCGCCCAGCAGGTAGCAAGGCTGGAATCCGCCATGTGCCTGCGCCGCCGCTGGTCGCTGGAAAACTTCCAGCTCTTCCTGGTTGAGCATCCGCTGGTTCGTCACTTAAcc                                                                                                         <  1:478837/150‑1 (MQ=255)
                                             aaCCGCTACAAATTGCTGAAAAAAGATGCGCGTACCATCGCCGCCCAGCAGGTAGCAAGGCTGGAATCCGCCATGTGCCTGCGCCGCCGCTGGTCGCTGGAAAACTTCCAGCTCTTCCTGGTTGAGCATCCGCTGGTTCGTCACTTAAcc                                                                                                         >  1:973736/1‑150 (MQ=255)
                                               ccGCTACAAATTGCTGAAAAAAGATGCGCGTACCATCGCCGCCCAGCAGGTAGCAAGGCTGGAATCCGCCATGTGCCTGCGCCGCCGCACGTCGCTGGAAAACTTCCAGCTCGTCCTGGTTGCGCATCCGCTGGTTCGTCACTTAAccgg                                                                                                       >  1:134113/1‑148 (MQ=255)
                                                   tACAAATTGCTGAAAAAAGATGCGCGTACCATCGCCGCCCAGCAGGTAGCAAGGCTGGAATCCGCCATGTGCCTGCGCCGCCGCTGGTCGCTGGAAAACTTCCAGCTCTTCCTGGTTGAGCATCCGCTGGTTCGTCACTTAACCCGCCGt                                                                                                   >  1:483015/1‑150 (MQ=255)
                                                           gCTGAAAAAAGATGCGCGTACCATCGCCGCCCAGCAGGTAGCAAGGCTGGAATCCGCCATGTGCCTGCGCCGCCGCTGGTCTCTGGAAAACTTCCAGCTCTTCCTGGTTGAGCATCCGCTGGTTCGTCACTTACCCCGCCGTCTGATTTg                                                                                           >  2:1511073/1‑150 (MQ=255)
                                                           gCTGAAAAAAGATGCGCGTACCATCGCCGCCCAGCAGGTAGCAAGGCTGGAATCCGCCATGTGCCTGCGCCGCCGCTGGTCGCTGGAAAACTTCCAGCTCTTCCTGGTTGAGCATCCGCTGGTTCGTCACTTAACCCGCCGTCTGATTTg                                                                                           >  1:1152384/1‑150 (MQ=255)
                                                           gCTGAAAAAAGATGCGCGTACCATCGCCGCCCAGCAGGTAGCAAGGCTGGAATCCGCCATGTGCCTGCGCCGCCGCTGGTCGCTGGAAAACTTCCAGCTCTTCCTGGTTGAGCATCCGCTGGTTCGTCACTTAACCCGCCGTCTGATTTg                                                                                           >  1:1815991/1‑150 (MQ=255)
                                                              gAAAAAAGATGCGCGTACCATCGCCGCCCAGCAGGTAGCAAGGCTGGAATCCGCCATGTGCCTGCGCCGCCGCTGGTCGCTGGAAAACTTCCAGCTCTTCCTGGTTGAGCATCCGCTGGTTCGGCCCTTAACCCGCCGTCGGCTTTgggg                                                                                        >  1:886398/1‑150 (MQ=255)
                                                              gAAAAAAGATGCGCGTACCATCGCCGCCCAGCAGGTAGCAAGGCTGGAATCCGCCATGTGCATGCGCCGCCGCTGGTCGCTGGAAAACTTCCAGCTCTTCCTGGTTGAGCATCCGCTGGTTCGTCACTTAACCCGCCGTCTGATTTgggg                                                                                        <  2:1300337/150‑1 (MQ=255)
                                                               aaaaaaGATGCGCGTACCATCGCCGCCCAGCAGGTAGCAAGGCTGGAATCCGCCATGTGCCTGCGCCGCCGCTGGTCGCTGGAAAACTTCCAGCTCTTCCTGGTTGAGCATCCGCTGGTTCGTCACTTAACCCGCCGTCTGATTTGGGGc                                                                                       >  2:1818622/1‑150 (MQ=255)
                                                                  aaaGATGCGCGTACCATCGCCGCCCAGCAGGTAGCAAGGCTGGAATCAGCCATGTGCCTGCTCCGCCGCTGGTCGCTGGAACACTTCCAGCTCTTGCTGGTTGAGCATCCGCTGGTGCGTCACTTAACCCGCCGGCGGAGTTgggggggt                                                                                    >  1:283950/1‑146 (MQ=255)
                                                                    aGATGCGCGTACCATCGCCGCCCAGCAGGTAGCAAGGCTGGAATCCGCCATGTGCCTGCGCCGCCGCTGGTCGCTGGAAAACTTCCAGCTCTTCCTGGTTGAGCATCCGCTGGTTCGTCACTTAACCCGCCGTCTGATTTGGGGCGTTta                                                                                  <  1:244374/150‑1 (MQ=255)
                                                                       tGCGCGTACCATCGCCGCCCAGCAGGTAGCAAGGCTGGAATCCGCCATGTGCCTGCGCCGCCGCTGGTCGCTGGAAAACTTCCAGCTCTTCCTGGTTGAGCATCCGCTGGTTCGTCACTTAACCCGCCGTCTGATTTGGGGCGTTTATAg                                                                               <  2:1160727/150‑1 (MQ=255)
                                                                            gTACCATCGCCGCCCAGCAGGTAGCAAGGCTGGAATCCGCCATGTGCCTGCGCCGCCGCTGGTCGCTGGAAAACTTCCAGCTCTTCCTGGTTGAGCATCCGCTGGTTCGTCACTTAACCCGCCGTCTGATTTGGGGCGTTTATAGCGCCg                                                                          >  1:1617128/1‑150 (MQ=255)
                                                                              aCCATCGCCGCCCAGCAGGTAGCAAGGCTGGAATCCGCCATGTGCCTGCGCCGCCGCTGGTCGCTGGAAAACTTCCAGCTCTTCCTGGTTGAGCATCCGCTGGTTCGTCACTTAACCCGCCGTCTGATTTGGGGCGTTTATAGCGCCGaa                                                                        >  1:1570549/1‑150 (MQ=255)
                                                                                    gccgccCAGCAGGTAGCAAGGCTGGAATCCGCCATGTGCCTGCGCCGCCGCTGGTCGCTGGAAACCTTCCAGCTCTTCCTGGTTGAGCATCCGCTGGTTCGTCACTTAACACGCCGTCTGATATGGGGCGTTTATAGCGCCGAAAACCAg                                                                  >  1:99895/1‑150 (MQ=255)
                                                                                    gccgccCAGCAGGTAGCAAGGCTGGAATCCGCCATGTGCCTGCGCCGCCGCTGGTCGCTGGAAAACTTCCAGCTCTTCCTGGTTGAGCATCCGCTGGTTCGTCACTTAACCCGCCGTCTGATTTGGGGCGTTTATAGCGCCGAAAACCAg                                                                  >  1:592191/1‑150 (MQ=255)
                                                                                    gccgccCAGCAGGTAGCAAGGCTGGAATCCGCCATGTGCCTGCGCCGCCGCTGGTCGCTGGAAAACTTCCAGCTCTTCCTGGTTGAGCATCCGCTGGTTCGTCACTTAACCCGCCGTCTGATTTGGGGCGTTTATAGCGCCGAAAACCAg                                                                  >  2:1806606/1‑150 (MQ=255)
                                                                                        cccAGCAGGTAGCAAGGCTGGAATCCGCCATGTGCCTGCGCCGCCGCTGGTCGCTGGAAAACTTCCAGCTCTTCCTGGTTGAGCATCCGCTGGTTCGTCACTTAACCCGCCGTCTGATTTGGGGCGTTTATAGCGCCGAAAACCAGCTAc                                                              >  1:1046079/1‑150 (MQ=255)
                                                                                        cccAGCAGGTAGCAAGGCTGGAATCCGCCATGTGCCTGCGCCGCCGCTGGTCGCTGGAAAACTTCCAGCTCTTCCTGGTTGAGCATCCGCTGGTTCGTCACTTAACCCGCCGTCTGATTTGGGGCGTTTATAGCGCCGAAAACCAGCTAc                                                              >  1:1856082/1‑150 (MQ=255)
                                                                                         ccAGCAGGTAGCAAGGCTGGAATCCGCCATGTGCCTGCGCCGCCGCTGGTCGCTGGAAAACTTCCCCCTCTTCCTGGTTCAGCATCCGCTGGATCGTCACTTAACCCGCCGTCTGAATTGGGGCGTTTATAGCGCCGAAAACCAGCTACt                                                             >  1:1958951/1‑150 (MQ=255)
                                                                                               ggTAGCAAGGCTGGAATCCGCCATGTGCCTGCGCCGCCGCTGGTCGCTGGAAAACTTCCAGCTCTTCCTGGTTGAGCATCCGCTGGTTCGTCACTTAACCCGCCGTCTGATTTGGGGCGTTTATAGCGCCGAACACCAGCTACTGgcttg                                                       >  2:599723/1‑150 (MQ=255)
                                                                                               ggTAGCAAGGCTGGAATCCGCCATGTGCCTGCGCCGCCGCTGGTCGCTGGAAAACTTCCAGCTCTTCCTGGTTGAGCATCCGCTGGTTCGTCACTTAACCCGCCGTCTGATTTGGGGCGTTTATAGCGCCGAAAACCAGCTACTGgcttg                                                       <  2:1186133/150‑1 (MQ=255)
                                                                                                gTAGCAAGTCTGGAATCCGCCATGTGCCTGCGCCGCCGCTGGTCGCTGGAAAACTTCCAGCTCTTCCTGGTTGAGCATCCGCTGGTTCGTCACTTAACCCGCCGTCTGATTTGGGGCGTTTATAGCGCCGAAAACCAGCTACTGgcttgc                                                      <  2:370847/150‑1 (MQ=255)
                                                                                                  tgCAAGGCTGGAATCCGCCATGTGCCTGCGCCGCCGCTGGTCGCTGGAAAACTTCCAGCTCTTCATGGTTGAGCATCCGCTGGTTCGTCACTTAACCCGCCGTCTGATTTGGGGCGTTTATAGCGCCGAAAACCAGCTACTGgcttgctt                                                    <  2:415328/149‑1 (MQ=255)
                                                                                                    cAAGGCTGGAATCCGCCATGTGCCTGCGCCGCCGCTGGTCGCTGGAAAACTTCCAGCTCTTCCTGGTTGAGCATCCGCTGGTTCGTCACTTAACCCGCCGTCTGATTTGGGGCGTTTATAGCGCCGAAAACCAGCTACTGGCTTGCTTTc                                                  <  1:1983994/150‑1 (MQ=255)
                                                                                                    cAAGGCTGGAATCCGCCATGTGCCTGCGCCGCCGCTGGTCGCTGGAAAACTTCCAGCTCTTCCTGGTTGAGCATCCGCTGGTTCGTCACTTAACCCGCCGTCTGATTTGGGGCGTTTATAGCGCCGAAAACCAGCTACTGGCTTGCTTTc                                                  >  2:17377/1‑150 (MQ=255)
                                                                                                    cAAGGCTGGAATCCGCCATGTGCCTGCGCCGCCGCTGGTCGCTGGAAAACTTCCAGCTCTTCCTGGTTGAGCATCCGCTGGTTCGTCACTTAACCCGCCGTCTGATTTGGGGCGTTTATAGCGCCGAAAACCAGCTACTGGCTTGCTTTc                                                  >  1:759795/1‑150 (MQ=255)
                                                                                                    cAAGGCGGGAATCCGCCATGTGCCTGCGCCGCCGCTGGTCTCTGGAAAACTTCAAGCTCTTGCGGGGTGAGCATCCGCTGGTTAGTCACTTCACCCGACGGCTGAGTTGGGGCTGTTATAGCGAAGCAAACCAGCTACGGGCTTGCTTTc                                                  >  2:119698/1‑150 (MQ=255)
                                                                                                     aaGGCTGGAATCCGCCATGTGCCTGCGCCGCCGCTGGTCGCTGGAAAACTTCCAGCTCTTCCTGGTTGAGCATCCGCTGGTTCGTCACTTAACCCGCCGTCTGATTTGGGGCGTTTATAGCGCCGAAAACCAGCTACTGGCTTGCTTTcg                                                 <  1:906099/150‑1 (MQ=255)
                                                                                                     aaGGCTGGAATCCGCCATGTGCCTGCGCCGCCGCTGGTCGCTGGAAAACTTCCAGCTCTTCCTGGTTGAGCATCCGCTGGTTCGTCACTTAACCCGCCGTCTGATTTGGGGCGTTTATAGCGCCGAAAACCAGCTACTGGCTTGCTTTcg                                                 <  2:1454523/150‑1 (MQ=255)
                                                                                                     aaGGCTGGAATCCGCCATGTGCCTGCGCCGCCGCTGGTCGCTGGAAAACTTCCAGCTCTTCCTGGTTGAGCATCCGCTGGTTCGTCACTTAACCCGCCGTCTGATTTGGGGCGTTTATAGCGCCGAAAACCAGCTACTGGCTTGCTTTcg                                                 <  1:1647385/150‑1 (MQ=255)
                                                                                                     aaGGCTGGAATCCGCCATGTGCCTGCGCCGCCGCGGGTCGCTGGAAAACTTCCAGCTCTTCCTGGTTGAGCATCCGCTGGTTCGTCACTTAACCCGCCGTCTGATTTGGGTCGTTTATAGCGCCGAAAACCAGCTACTGGCTTGCTTTcg                                                 <  2:1771433/150‑1 (MQ=255)
                                                                                                       ggCTGGAATCCGCCATGTGCCTGCGCCGCCGCTGGTCGCTGGAAAACTTCCAGCTCTTCCTTGTTGAGCATCCGCTGGTTCGTCACTTAACGCGCCGTCTGATTTGGGGCGTTTATAGCGCCGAAAACCAGCTACTGGCTTGCTTTcgcg                                               <  2:243320/150‑1 (MQ=255)
                                                                                                       ggCTGGAATCCGCCATGTGCCTGCGCCGCCGCTGGTCGCTGGAAAACTTCCAGCTCTTCCTGGTTGAGCATCCGCTGGTTCGTCACTTAACCCGCCGTCTGATTTGGGGCGTTTATAGCGCCGAAAACCAGCTACTGGCTTGCTTTcgcg                                               <  1:46208/150‑1 (MQ=255)
                                                                                                         cTGGAATCCGCCATGTGCCTGCGCCGCCGCTGGTCGCTGGACAACTTCCAGCTCTTCCTGGTTGAGCCTCCGCTGGTTCGTCACTTAACCCCCCGACTGATTTGGGGCGTTTATAGCGCCGAAAACCAGCTACTGGCTTGCTTTCGCGTa                                             >  1:1043286/1‑150 (MQ=255)
                                                                                                           ggAATCCGCCATGTGCCTGCGCCGCCGCTGGTCGCTGGAAAACTTCCAGCTCTTCCTGGTTGAGCATCCGCTGGTTCGTCACTTAACCCGCCGTCTGATTTGGGGCGTTTATAGCGCCGAAAACCAGCTACTGGCTTGCTTTCGCGTAGc                                           >  2:1901222/1‑150 (MQ=255)
                                                                                                           ggAATCCGCCATGTGCCTGCGCCGCCGCTGGTCGCTGGAAAACTTCCAGCTCTTCCTGGTTGAGCATCCGCTGGTTCGTCACTTAACCCGCCGTCTGATTTGGGGCGTTTATAGCGCCGAAAACCAGCTACTGGCTTGCTTTCGCGTAGc                                           <  1:1633050/150‑1 (MQ=255)
                                                                                                           ggAATCCGCCATGTGCCTGCGCCGCCGCTGGTCGCTGGAAAACTTCCAGCTCTTCCTGGTTGAGCATCCGCTGGTTCGTCACTTAACCCGCCGTCTGATTTGGGGCGTTTATAGCGCCGAAAACCAGCTACTGGCTTGCTTTCGCGGAGc                                           >  1:1846336/1‑150 (MQ=255)
                                                                                                            gAATCCGCCATGTGCCTGCGCCGCCGCTGGTCGCTGGAAAACTTCCAGCTCTTCCTGGTTGAGCATCCGCTGGTTCGTCACTTAACCCGCCGTCTGATTTGGGGCGTTTATAGCGCCGAAAACCAGCTACTGGCTTGCTTTCGCGTAGCa                                          <  1:899515/150‑1 (MQ=255)
                                                                                                               tCCGCCATGTGCCTGCGCCGCCGCTGGTCGCTGGAAAACTTCCAGCTCTTCCTGGTTGAGCATCCGCTGGTTCGTCACTTAACCCGCCGTCTGATTTGGGGCGTTTATAGCGCCGAAAACCAGCTACTGGCTTGCTTTCGCGTAGCagaa                                       <  1:137356/150‑1 (MQ=255)
                                                                                                                ccGGCATGTGCCTGCGCCGCCGCTGGGCGCTGGAAAACTTCCAGCTCTGCCTGGTTGAGCATCCGCTGGTTCGTCACTTAACCCGCCGTCTGATTTGGGGCGTTTATAGCGCCGAAAACCAGCTACTGGCTTGCTTTCGCGTAGCagaag                                      <  2:835473/150‑1 (MQ=255)
                                                                                                                ccGCCATGTGCCTGCGCCGCCGCTGGTCGCTGGAAAACTTCCAGCTCTTCCTGGTTGAGCATCCGCTGGTTCGTCACTTAACCCGCCGTCTGATTTGGGGCGTTTATAGCGCCGAAAACCAGCTACTGGCTTGCTTTCGCGTAGCagaag                                      <  2:418112/150‑1 (MQ=255)
                                                                                                                ccGCCATGTGCCTGCGCCGCCGCTGGTCGCTGGAAAACTTCCAGCTCTTCCTGGTTGAGCATCCGCTGGTTCGTCACTTAACCCGCCGTCTGATTTGGGGCGTTTATAGCGCCGAAAACCAGCTACTGGCTTGCTTTCGCGTAGCagaag                                      <  2:836769/150‑1 (MQ=255)
                                                                                                                ccGCCATGTGCCTGCGCCGCCGCTGGTCGCTGGAAAACTTCCAGCTCTTCCTGGTTGAGCATCCGCTGGTTCGTCACTTAACCCGCCGTCTGATTTGGGGCGTTTATAGCGCCGAAAACCAGCTACTGGCTTGCTTTCGCGTAGCagaag                                      >  1:2001527/1‑150 (MQ=255)
                                                                                                                 cGCCATGTGCCTGCGCCGCCGCTGGTCGCTGGAAAACTTCCAGCTCTTCCTGGTTGAGCATCCGCTGGTTCGTCACTTAACCCGCCGTCTGATTTGGGGCGTTTATAGCGCCGAAAACCAGCTACTGGCTTGCTTTCGCGTAGCagaaga                                     >  2:214491/1‑150 (MQ=255)
                                                                                                                  gCCATGTGCCTGCGCCGCCGCTGGTCGCTGGAAAACTTCCAGCTCTTCCTGGTTGAGCATCCGCTGGTTCGTCACTTAACCCGCCGTCTGATTTGGGGCGTTTATAGCGCCGAAAACCAGCTACTGGCTTGCTTTCGCGTAGCAGAAGAt                                    <  2:1555413/150‑1 (MQ=255)
                                                                                                                      tgtgCCTGCGCCGCCGCTGGTCGCTGGAAAACTTCCAGCTCTTCCTGGTTGAGCATCCGCTGGTTCGTCACTTAACCCGCCGTCTGATTTGGGGCGTTTATAGCGCCGAAAACCAGCTACTGGCTTGCTTTCGCGTAGCAGAAGATAACa                                <  2:861027/150‑1 (MQ=255)
                                                                                                                      tgtgCCTGCGCCGCCGCTGGTCGCTGGAAAACTTCCAGCTCTTCCTGGTTGAGCATCCGCTGGTTCGTCACTTAACCCGCCGTCTGATTTGGGGCGTTTATAGCGCCGAAAACCAGCTACTGGCTTGCTTTCGCGTAGCAGAAGATAACa                                <  1:1827169/150‑1 (MQ=255)
                                                                                                                       gtgCCTGCGCCGCCGCTGGTCGCTGGAAAACTTCCAGCTCTTCCTGGTTGAGCATCCGCTGGTTCGTCACTTAACCCGCCGTCTGATTTGGGGCGTTTATAGCGCCGAAAACCAGCTACTGGCTTGCTTTCGCGTAGCAGAAGATAACAg                               <  1:452592/150‑1 (MQ=255)
                                                                                                                       gtgCCTGCGCCGCCGCTGGTCGCTGGAAAACTTCCAGCTCTTCCTGGTTGAGCATCCGCTGGTTCGTCACTTAACCCGCCGTCTGATTTGGGGCGTTTATAGCGCCGAAAACCAGCTACTGGCTTGCTTTCGCGTAGCAGAAGATAACAg                               >  2:1188699/1‑150 (MQ=255)
                                                                                                                         gCCTGCGCCGCCGCTGGTCGCTGGAAAACTTCCAGCTCTTCCTGGTTGAGCATCCGCTGGTTCGTCACTTAACCCGCCGTCTGATTTGGGGCGTTTATAGCGCCGAAAACCAGCTACTGGCTTGCTTTCGCGTAGCAGAAGATAACAGCt                             <  2:1617267/150‑1 (MQ=255)
                                                                                                                         gCCTGCGCCGCCGCTGGTCGCTGGAAAACTTCCAGCTCTTCCTGGTTGAGCATCCGCTGGTTCGTCACTTAACCCGCCGTCTGATTTGGGGCGTTTATAGCGCCGAAAACCAGCTACTGGCTTGCTTTCGCGTAGCAGAAGATAACAGCt                             >  1:1066582/1‑150 (MQ=255)
                                                                                                                           cTGCGCCGCCGCTGGTCGCTGGAAAACTTCCAGCTCTTCCTGGTTGAGCATCCGCTGGTTCGTCACTTAACCCGCCGTCTGATTTGGGGCGTTTATAGCGCCGAAAACCAGCTACTGGCTTGCTTTCGCGTAGCAGAAGATAACAGCTcc                           >  1:235027/1‑150 (MQ=255)
                                                                                                                            tGCGCCGCCGCTGGTCGCTGGAAAACTTCCAGCTCTTCCTGGTTGAGCATCCGCTGGTTCGTCACTTAACCCGCCGTCTGATTTGGGGCGTTTATAGCGCCGAAAACCAGCTACTGGCTTGCTTTCGCGTAGCAGAAGATAACAGCTCCa                          <  1:1858687/150‑1 (MQ=255)
                                                                                                                                ccgccgcTGGTCGCTGGAAAACTTCCAGCTCTTCCTGGTTGAGCATCCGCTGGTTCGTCACTTAACCCGCCGTCTGATTTGGGGCGTTTATAGCGCCGAAAACCAGCTACTGGCTTGCTTTCGCGTAGCAGAAGATAACAGCTCCAGCAc                      >  1:563864/1‑150 (MQ=255)
                                                                                                                                 cgccgcTGGTCGCTGGAAAACTTCCAGCTCTTCCTGGTTGAGCATCCGCTGGTTCGTCACTTAACCCGCCGTCTGATTTGGGGCGTTTATAGCGCCGAAAACCAGCTACTGGCTTGCTTTCGCGTAGCAGAAGATAACAGCTCCAGCAcc                     >  2:1161188/1‑150 (MQ=255)
                                                                                                                                    cgcTGGTCGCTGGAAAACTTCCAGCTCTTCCTGGTTGAGCATCCGCTGGTTCGTCACTTAACCCGCCGTCTGATTTGGGGCGTTTATAGCGCCGAAAACCAGCTACTGGCTTGCTTTCGCGTAGCAGAAGATAACAGCTCCAGCACCGCt                  >  2:626354/1‑150 (MQ=255)
                                                                                                                                     aggaagagCTGGAAGTTTTCCAGCTCTTCCTGGTTGAGCATCCGCTGGTTCGTCACTTAACCCGCCGTCTGATTTGGGGCGTTTATAGCGCCGAAAACCAGCTACTGGCTTGCTTTCGCGTAGCAGAAGATAACAGCTCCAGCACCGCTg                 >  2:649832/8‑150 (MQ=255)
                                                                                                                                       tGGTCGCTGGAAAACTTCCAGCTCTTCCTGGTTGAGCATCCGCTGGTTCGTCACTTAACCCGCCGTCTGATTTGGGGCGTTTATAGCGCCGAAAACCAGCTACTGGCTTGCTTTCGCGTAGCAGAAGATAACAGCTCCAGCACCGCTGAc               <  1:658897/150‑1 (MQ=255)
                                                                                                                                        ggTCGCTGGAAAACTTCCAGCTCTTCCTGGTTGAGCATCCGCTGGTTCGTCACTTAACCCGCCGTCTGATTTGGGGCGTTTATAGCGCCAAAAACCAGCTACTGGCTTGCTTTCGCGTAGCAGAAGATAACAGCTCCAGCACCGCTGACg              <  1:139549/150‑1 (MQ=255)
                                                                                                                                         agagCTGGAAGTTTTCCAGCTCTTCCTGGTTGAGCATCCGCTGGTTCGTCACTTAACCCGCCGTCTGATTTGGGGCGTTTATAGCGCCGAAAACCAGCTACTGGCTTGCTTTCGCGTAGCAGAAGATAACAGCTCCAGCACCGCTGACga             >  2:1680220/4‑150 (MQ=255)
                                                                                                                                           cGCTGGAAAACTTCCAGCTCTTCCTGGTTGAGCATCCGCTGGTTCGTCACTTAACCCGCCGTCTGATTTGGGGCGTTTATAGCGCCGAAAACCAGCTACTGGCTTGCTTTCGCGTAGCAGAAGATAACAGCTCCAGCACCGCTGACgatg           <  1:707493/150‑1 (MQ=255)
                                                                                                                                             cTGGAAAACTTCCAGCTCTTCCTGGTTGAGCATCCGCTGGTTCGTCACTTAACCCGCCGTCTGATTTGGGGCGTTTATAGCGCCGAAAACCAGCTACTGGCTTGCTTTCGCGTAGCAGAAGATAACAGCTCCAGCACCGCTGCCgatgat         >  2:386562/1‑150 (MQ=255)
                                                                                                                                                 aaaaCTTCCCGCGCTGCCTGGTTGAGCATCCGCTGGTTAGTCACTTCACCCGCCGTCTGATTGTGGGAGTTTATAGCGCCGAAAACCAGCTACTGTCTTGCTTTCGCGTAGCAGAAGTTAACAGCTCCAGCACCGCTGACGATGATCttt     <  1:329307/150‑1 (MQ=255)
                                                                                                                                                 aaaaCTTCCAGCTCTTCCTGGTTGAGCATCCGCTGGTTCGTCACTTAACCCGCCGTCTGATTTGGGGCGTTTATAGCGCAGAAAACCACCTACTGGCTTGCTTTCGCGTAGCAGAAGATAACAGCTCCAGCACCGCTGACGATGATCttt     >  1:809694/1‑150 (MQ=255)
                                                                                                                                                 aaaaCTTCCAGCTCTTCCTGGTTGAGCATCCGCTGGTTCGGCACTTAACCCGCCGTCTGATTTGGGGCGTTTATAGCGCCGAAAACCAGCTACTGGCTTGCTTTCGCGTAGCAGAAGATAACAGCTCCAGCACCGCTGACGATGATCttt     <  1:214491/150‑1 (MQ=255)
                                                                                                                                                  aaaCTTCCAGCTCTTCCTGGTTGAGCATCCGCTGGTTCGTCACTTAACCCGCCGTCTGATTTCGGGCGTTTATAGCGCCGAAAACCAGCAACTGGCTTGCTTTCGCGTAGCAGAAGATAACCGCTACAGCACTGCTGACGATCATCtatt    >  1:1771830/1‑150 (MQ=255)
                                                                                                                                                   aaCTTCCAGCTCTTCCTGGTTGAGCATCCGCTGGTTCGTCACTTAACCCGCCGTCTGATTTGGGGCGTTTATAGCGCCGAAAACCAGCTACTGGCTTGCTTTCGCGTAGCAGAAGATAACAGCTCCAGCACCGCTGACGATGATCTTTTc   <  2:1535844/150‑1 (MQ=255)
                                                                                                                                                    aCTTCCAGCTCTTCCTGGTTGAGCATCCGCTGGTTCGTCACTTAACCCGCCGTCTGATTTGGGGCGTTTATAGCGCCGAAAACCAGCTACTGGCTTGCTTTCGCGTAGCAGAAGATAACAGCTCCAGCACCGCTGACGATGATCTTTTCa  <  2:192267/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
AAACCGAACAAAAGCGATGATGAAACGCGGGCGAACGATGCGGTTAACCGCTACAAATTGCTGAAAAAAGATGCGCGTACCATCGCCGCCCAGCAGGTAGCAAGGCTGGAATCCGCCATGTGCCTGCGCCGCCGCTGGTCGCTGGAAAACTTCCAGCTCTTCCTGGTTGAGCATCCGCTGGTTCGTCACTTAACCCGCCGTCTGATTTGGGGCGTTTATAGCGCCGAAAACCAGCTACTGGCTTGCTTTCGCGTAGCAGAAGATAACAGCTCCAGCACCGCTGACGATGATCTTTTCA  >  NC_000913/2199460‑2199757

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