Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,198,999 C→T 5.5% pseudogene (2526/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,198,9990CT 5.5% 231.3 / 2.5 72pseudogene (2526/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 C (35/33);  new base T (2/2);  total (37/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.49e-01

GTATTAAAAGCAAATGCCCAATCACTGCCAGATAGCGCCCTCCTCCACCTCGGTGAAATGCTCCGCTTCCCTCAGGAAGAGGCTCTGTATCCGGGATTATTGCAGGTGAAAGACGTCTGTTCCGCCGACTCACTGGCGGGATTTGCCTGGGATCTGTTTACCGCCTGGCAGACCGCTGGCGCGCCGTCGAAAGAGAGTTGGGCGTTCACTGCGTTAGGCGTTCTCGGTAACGATGACACCGCCCGCAAACTGACGCCATTAATACGCGCCTGGCCTGGTGAATCCCAGC  >  NC_000913/2198853‑2199141
                                                                                                                                                  |                                                                                                                                              
gTATTAAAAGCAAATGCCCAATCACTGCCAGATAGCGCCCTCCTCCACCTCGGTGAAATGCTCCGCTTCCCTCAGGAAGAGGCTCTGTATCCGGGATTATTGCAGGTGAAAGACGTCTGTTCCGCCGACTCACTGGCGGGATTTGCCTg                                                                                                                                              >  4:163160/1‑149 (MQ=255)
 tATTAAAAGCAAATGCCCAATCACTGCCAGATAGCGCCCTCCTCCACCTCGGTGAAATGCTCCGCTTCCCTCAGGAAGAGGCTCTGTATCCGGGATTATTGCAGGTGAAAGACGTCTGTTCCGCCGACTCACTGGCGGGATTTGCCTgg                                                                                                                                             <  2:192043/149‑1 (MQ=255)
   ttAAAAGCAAATGCCCAATCACTGCCAGATAGCGCCCTCCTCCACCTCGGTGAAATGCTCCGCTTCCCTCAGGAAGAGGCTCTGTATCCGGGATTATTGCAGGTGAAAGACGTCTGTTCCGCCGACTCACTGGCGGGATTTGCCTGGGa                                                                                                                                           <  3:242356/149‑1 (MQ=255)
   ttAAAAGCAAATGCCCAATCACTGCCAGATAGCGCCCTCCTCCACCTCGGTGAAATGCTCCGCTTCCCTCAGGAAGAGGCTCTGTATCCGGGATTATTGCAGGTGAAAGACGTCTGTTCCGCCGACTCACTGGCGGGATTTGCCTGGGa                                                                                                                                           >  7:255995/1‑149 (MQ=255)
   ttAAAAGCAAATGCCCAATCACTGCCAGATAGCGCCCTCCTCCACCTCGGTGAAATGCGCCGCTTCCCTCAGGAAGAGGCTCTGTAGCCGGGATTATTGCAGGTGAAAGACGTCTGTTCCGCCGACTCACTGGCGGGATTTGCCTggg                                                                                                                                            <  4:98392/148‑1 (MQ=255)
     aacaGCAAATGCCCAATCACTGCCAGATAGCGCCCTCCTCCACCTCGGGGAAATTCTCCGCTTCCCTCAGGAAGAGGCTCTGTATCCGGGAGTATTGCAGGTGAAAGACGTCTGTTCCGCCGACTCACTGGCGGGATTTGCCTGGGAt                                                                                                                                          <  6:94212/145‑1 (MQ=255)
       aaGCAAATGCCCAATCACTGCCAGATAGCGCCCTCCTCCACCTCGGTGAAATGCTCCGCTTCCCTCAGGAAGAGGCTCTGTATCCGGGATTATTGCAGGTGAAAGACGTCTGTTCCGCCGACTCACTGGCGGGATTTGCCTGGGATCTg                                                                                                                                       >  2:40517/1‑149 (MQ=255)
       aaGCAAATGCCCAATCACTGCCAGATAGCACCCTCCTCCACCTCGGTGAAATGCTCCGCTTCCCTCAGGAAGAGGCTCTGTATCCGGGATTATTGCAGGTGAAAGACGTCTGTTCCGCCGACTCACTGGCGGGATTTGCCTGGGATCTg                                                                                                                                       <  2:35358/149‑1 (MQ=255)
         gCAAATGCCCAATCACTGCCAGATAGCGCCCTCCTCCACCTCGGTGAAATGCTCCGCTTCCCTCAGGAAGAGGCTCTGTATCCGGGATTATTGCAGGTGAAAGACGTCTGTTCCGCCGACTCACTGGCGGGATTTGCCTGGGATCTGtt                                                                                                                                     >  7:310217/1‑149 (MQ=255)
         gCAAATGCCCAATCACTGCCAGATAGCGCCCTCCTCCACCTCGGTGAAATGCTCCGCTTCCCTCAGGAAGAGGCTCTGTATCCGGGATTATTGCAGGTGAAAGACGTCTGTTCCGCCGACTCACTGGCGGGATTTGCCTGGGATCTGtt                                                                                                                                     >  7:319199/1‑149 (MQ=255)
             aTGCCCAATCACTGCCAGATAGCGCCCTCCTCCACCTCGGTGAAATGCTCCGCTTCCCTCAGGAAGAGGCTCTGTATCCGGGATTATTGCAGGTGAAAGACGTCTGTTCCGCCGACTCACTGGCGGGATTTGCCTGGGATCTGTTTAcc                                                                                                                                 >  7:272538/1‑149 (MQ=255)
             aTGCCCAATCACTGCCAGATAGCGCCCTCCTCCACCTCGGTGAAATGCTCCGCTTCCCTCAGGAAGAGGCTCTGTATCCGGGATTATTGCAGGTGAAAGACGTCTGTTCCGCCGACTCACTGGCGGGATTTGCCTGGGATCTGTTTAcc                                                                                                                                 <  5:259335/149‑1 (MQ=255)
                cccAATCACTGCCAGATAGCGCCCTCCTCCACCTCGGTGAAATGCTCCGCTTCCCTCAGGAAGAGGCTCTGTATCCGGGATTATTGCAGGTGAAAGACGTCTGTTCCGCCGACTCACTGGCGGGATTTGCCTGGGATCTGTTTACCGcc                                                                                                                              >  2:168165/1‑149 (MQ=255)
                 ccAATCACTGCCAGATAGCGCCCTCCTCCACCTCGGTGAAATGCTCCGCTTCCCTCAGGAAGAGGCTCTGTATCCGGGATTATTGCAGGTGAAAGACGTCTGTTCCGCCGACTCACTGGCGGGATTTGCCTGGGATCTGTTTACCGCCt                                                                                                                             >  8:319590/1‑149 (MQ=255)
                  cAATCACTGCCAGATAGCGCCCTCCTCCACCTCGGTGAAATGCTCCGCTTCCCTCAGGAAGAGGCTCTGTATCCGGGATTATTGCAGGTGAAAGACGTCTGTTCCGCCGACTCACTGGCGGGATTTGCCTGGGATCTGTTTACCGCCTg                                                                                                                            <  7:319590/149‑1 (MQ=255)
                   aaTCACTGCCAGATAGCGCCCTCCTCCACCTCGGTGAAATGCTCCGCTTCCCTCAGGAAGAGGCTCTGTATCCGGGATTATTGCAGGTGAAAGACGTCTGTTCCGCCGACTCACTGGCGGGATTTGCCTGGGATCTGTTTACCGCCTgg                                                                                                                           >  5:340804/1‑149 (MQ=255)
                   aaTCACTGCCAGATAGCGCCCTCCTCCACCTCGGTGAAATGCTCCGCTTCCCTCAGGAAGAGGCTCTGTATCCGGGATTATTGCAGGTGAAAGACGTCTGTTCCGCCGACTCACTGGCGGGATTTGCCTGGGATCTGTTTACCGCCTgg                                                                                                                           >  2:220570/1‑149 (MQ=255)
                     tCACTGCCAGATAGCGCCCTCCTCCACCTCGGTGAAATGCTCCGCTTCCCTCAGGAAGAGGCTCTGTATCCGGGATTATTGCAGGTGAAAGACGTCTGTTCCGCCGACTCACTGGCGGGATTTGCCTGGGATCTGTTTACCGCCTgg                                                                                                                           >  3:147831/1‑147 (MQ=255)
                     tCACTGCCAGATAGCGCCCTCCTCCACCTCGGTGAAATGCTCCGCTTCCCTCAGGAAGAGGCTCTGTATCCGGGATTATTGCAGGTGAAAGACGTCTGTTCCGCCGACTCACTGGCGGGATTTGCCTGGGATCTGTTTACCGCCTGGCa                                                                                                                         <  8:125157/149‑1 (MQ=255)
                         tGCCAGATAGCGCCCTCCTCCACCTCGGTGAAATGCTCCGCTTCCCTCAGGAAGAGGCTCTGTATCCGGGATTATTGCAGGTGAAAGACGTCTGTTCCGCCGACTCACTGGCGGGATTTGCCTGGGATCTGTTTACCGCCTGGCAGAcc                                                                                                                     >  8:86582/1‑149 (MQ=255)
                               aTAGCGCCCTCCTCCACCTCGGTGAAATGCTCCGCTTCCCTCAGGAAGAGGCTCTGTATCCGGGATTATTGCAGGTGAAAGACGTCTGTTCCGCCGACTCACTGGCGGGATTTGCCTGGGATCTGTTTACCGCCTGGCAGACCGCTGgc                                                                                                               <  5:138230/149‑1 (MQ=255)
                                  gcgcCCTCCTCCACCTCGGTGAAATGCTCCGCTTCCCTCAGGAAGAGGCTCTGTATCCGGGATTATTGCAGGTGAAAGACGTCTGTTCCGCCGACTCACTGGCGGGATTTGCCTGGGATCTGTTTACCGCCTGGCAGACCGCTGgcgcg                                                                                                            >  8:1011/1‑149 (MQ=255)
                                         cctccACCTCGGTGAAATGCTCCGCTTCCCTCAGGAAGAGGCTCTGTATCCGGGATTATTGCAGGTGAAAGACGTCTGTTCCGCCGACTCACTGGCGGGATTTGCCTGGGATCTGTTTACCGCCTGGCAGACCGCTGGCGCGCCGTCg                                                                                                      >  2:112471/1‑148 (MQ=255)
                                          ctccACCTCGGTGAAATGCTCCGCTTCCCTCAGGAAGAGGCTCTGTATCCGGGATTATTGCAGGTGAAAGACGTCTGTTCCGCCGACTCACTGGCGGGATTTGCCTGGGATCTGTTTACCGCCTGGCAGACCGCTGGCGCGCCGTCGa                                                                                                     >  1:57194/1‑148 (MQ=255)
                                             cACCTCGGTGAAATGCTCCGCTTCCCTCAGGAAGAGGCTCTGTATCCGGGATTATTGCAGGTGAAAGACGTCTGTTCCGCCGACTCACTGGCGGGATTTGCCTGGGATCTGTTTACCGCCTGGCAGACCGCTGGCGCGCCGTCGaaa                                                                                                   <  1:220570/147‑1 (MQ=255)
                                             cACCTCGGTGAAATGCTCCGCTTCCCTCAGGAAGAGGCTCTGTATCCGGGATTATTGCAGGTGAAAGACGTCTGTTCCGCCGACTCACTGGCGGGATTTGCCTGGGATCTGTTTACCGCCTGGCAGACCGCTGGCGCGCCGTCGAAaga                                                                                                 <  8:310217/149‑1 (MQ=255)
                                               ccTCGGTGAAATGCTCCGCTTCCCTCAGGAAGAGGCTCTGTATCCGGGATTATTGCAGGTGAAAGACGTCTGTTCCGCCGACTCACTGGCGGGATTTGCCTGGGATCTGTTTACCGCCTGGCAGACCGCTGGCGCGCCGTCGAAagaga                                                                                               <  6:301428/149‑1 (MQ=255)
                                                cTCGGTGAAATGCTCCGCTTCCCTCAGGAAGAGGCTCTGTATCCGGGATTATTGCAGGTGAAAGACGTCTGTTCCGCCGACTCACTGGCGGGATTTGCCTGGGATCTGt                                                                                                                                      >  2:74688/1‑109 (MQ=255)
                                                cTCGGTGAAATGCTCCGCTTCCCTCAGGAAGAGGCTCTGTATCCGGGATTATTGCAGGTGAAAGACGTCTGTTCCGCCGACTCACTGGCGGGATTTGCCTGGGATCTGt                                                                                                                                      <  1:74688/109‑1 (MQ=255)
                                                 tCGGTGAAATGCTCCGCTTCCCTCAGGAAGAGGCTCTGTATCCGGGATTATTGCAGGTGAAAGACGTCTGTTCCGCCGACTCACTGGCGGGATTTGCCTGGGATCTGTTTACCGCCTGGCAGACCGCTGGCGCGCCGTCGAAAGAGAGt                                                                                             >  1:326297/1‑149 (MQ=255)
                                                 tCGGTGAAATGCTCCGCTTCCCTCAGGAAGAGGCTCTGTATCCGGGATTATTGCAGGTGAAAGACGTCTGTTCCGCCGACTCACTGGCGGGATTTGCCTGGGATCTGTTTACCGCCTGGCAGACCGCTGGCGCGCCGTCGAAAGAGAGt                                                                                             <  2:326297/149‑1 (MQ=255)
                                                   ggTGAAATGCTCCGCTTCCCTCAGGAAGAGGCTCTGTATCCGGGATTATTGCAGGTGAAAGACGTCTGTTCCGCCGACTCACTGGCGGGATTTGCCTGGGATCTGTTTACCGCCTGGCAGACCGCTGGCGCGCCGTCGAAAGAGAGTTg                                                                                           <  5:143973/149‑1 (MQ=255)
                                                      gAAATGCTCCGCTTCCCTCAGGAAGAGGCTCTGTATCCGGGATTATTGCAGGTGAAAGACGTCTGTTCCGCCGACTCACTGGCGGGATTTGCCTGGGATCTGTTTACCGCCTGGCAGACCGCTGGCGCGCCGTCGAAAGAGAGTTGGGc                                                                                        >  7:276561/1‑149 (MQ=255)
                                                                 cTTCCCTCAGGAAGAGGCTCTGTATCCGGGATTATTGCAGGTGAAAGACGTCTGTTCCGCCGACTCACTGGCGGGATTTGCCTGGGATCTGTTTACCGCCTGGCAGACCGCTGGCGCGCCGTCGAAAGAGAGTTGGGCGTTCACTGCGt                                                                             <  8:319199/149‑1 (MQ=255)
                                                                 cTTCCCTCAGGAAGAGGCTCTGTATCCGGGATTATTGCAGGTGAAAGACGTCTGTTCCGCCGACTCACTGGCGGGATTTGCCTGGGATCTGTTTACCGCCTGGCAGACCGCTGGCGCGCCGTCGAAAGAGAGTTGGGCGTTCACTGCGt                                                                             <  5:275448/149‑1 (MQ=255)
                                                                 cTTCCCTCAGGAAGAGGCTCTGTATCCGGGATTATTGCAGGTGAAAGACGTCTGTTCCGCCGACTCACTGGCGGGATTTGCCTGGGATCTGTTTACCGCCTGGCAGACCGCTGGCGCGCCGTCGAAAGAGAGTTGGGCGTTCACTGCGt                                                                             <  7:86582/149‑1 (MQ=255)
                                                                     ccTCAGGAAGAGGCTCTGTATCCGGGATTATTGCAGGTGAAAGACGTCTGTTCCGCCGACTCACTGGCGGGATTTGCCTGGGATCTGTTTACCGCCTGGCAGACCGCTGGCGCGCCGTCGAAAGAGAGTTGGGCGTTCACTGCGTTAg                                                                          <  1:1414/148‑1 (MQ=255)
                                                                        cAGGAAGAGGCTCTGTATCCGGGATTATTGCAGGTGAAAGACGTCTGTTCCGCCGACTCACTGGCGGGATTTGCCTGGGATCTGTTTACCGCCTGGCAGACCGCTGGCGCGCCGTCGAAAGAGAGTTGGGCGTTCACTGCGTTAGGCGt                                                                      >  7:64565/1‑149 (MQ=255)
                                                                         aGGAAGAGGCTCTGTATCCGGGATTATTGCAGGTGAAAGACGTCTGTTCCGCCGACTCACTGGCGGGATTTGCCTGGGATCTGTTTACCGCCTGGCAGACCGCTGGCGCGCCGTCGAAAGAGAGTTGGGCGTTCACTGCGTTAGGCGtt                                                                     >  5:37768/1‑149 (MQ=255)
                                                                              gagGCTCTGTATCCGGGATTATTGCAGGTGAAAGCCGTCTGTTCCGCCGACTCACTGGCGGGATTTGCCTGGGATCTGTTTACCGCCTGGCAGACCGCTGGCGCGCCGTCGAAAGAGAGTTGGGCGTTCACTGCGTTAGGCGTTCTCgg                                                                <  1:111080/149‑1 (MQ=255)
                                                                                 gCTCTGTATCCGGGATTATTGCAGGTGAAAGACGTCTGTTCCGCCGACTCACTGGCGGGATTTGCCTGGGATCTGTTTACCGCCTGGCAGACCGCTGGCGCGCCGTCGAAAGAGAGTTGGGCGTTCACTGCGTTAGGCGTTCTCGGTaa                                                             >  4:370114/1‑149 (MQ=255)
                                                                                 gCTCTGTATCCGGGATTATTGCAGGTGAAAGACGTCTGTTCCGCCGACTCACTGGCGGGATTTGCCTGGGATCTGTTTACCGCCTGGCAGACCGCTGGCGCGCCGTCGAAAGAGAGTTGGGCGTTCACTGCGTTAGGCGTTCTCGGTaa                                                             >  3:141971/1‑149 (MQ=255)
                                                                                   tctGTATCCGGGATTATTGCAGGTGAAAGACGTCTGTTCCGCCGACTCACTGGCGGGATTTGCCTCGGATCTGTTTACCGCCTGGCAGACCGCTGGCGCGCCGTCGAAAGAGAGTTGGGCGTTCACTGCGTTAGGCGTTCTCGGTAAc                                                            >  3:349831/1‑148 (MQ=255)
                                                                                       tATCCGGGATTATTGCAGGTGAAAGACGTCTGTTCCGCCGACTCACTGGCGGGATTTGCCTGGGATCTGTTTACCGCCTGGCAGACCGCTGGCGCGCCGTCGAAAGAGAGTTGGGCGTTCACTGCGTTAGGCGTTCTCGGTAACGAt                                                         <  1:198940/147‑1 (MQ=255)
                                                                                       tATCCGGGATTATTGCAGGTGAAAGACGTCTGTTCCGCCGACTCACTGGCGGGATTTGCCTGGGATCTGTTTACCGCCTGGCAGACCGCTGGCGCGCCGTCGAAAGAGAGTTGGGCGTTCACTGCGTTAGGCGTTCTCGGTAACGAt                                                         <  8:255995/147‑1 (MQ=255)
                                                                                       tATCCGGGATTATTGCAGGTGAAAGACGTCTGTTCCGCCGACTCACTGGCGGGATTTGCCTGGGATCTGTTTACCGCCTGGCAGACCGCTGGCGCGCCGTCGAAAGAGAGTTGGGCGTTCACTGCGTTAGGCGTTCTCGGTAACGAt                                                         <  6:37768/147‑1 (MQ=255)
                                                                                       tATCCGGGATTATTGCAGGTGAAAGACGTCTGTTCCGCCGACTCACTGGCGGGATTTGCCTGGGATCTGTTTACCGCCTGGCAGACCGCTGGCGCGCCGTCGAAAGAGAGTTGGGCGTTCACTGCGTTAGGCGTTCTCGGTAACGATGa                                                       <  1:112471/149‑1 (MQ=255)
                                                                                         tCCGGGATTATTGCAGGTGAAAGACGTCTGTTCCGCCGACTCACTGGCGGGATTTGCTTGGGATCTGTTTACCGCCTGGCAGACCGCTGGCGCGCCGTCGAAAGAGAGTTGGGc                                                                                        >  3:211801/1‑114 (MQ=255)
                                                                                         tCCGGGATTATTGCAGGTGAAAGACGTCTGTTCCGCCGACTCACTGGCGGGATTTGCTTGGGATCTGTTTACCGCCTGGCAGACCGCTGGCGCGCCGTCGAAAGAGAGTTGGGc                                                                                        >  7:216564/1‑114 (MQ=255)
                                                                                         tCCGGGATTATTGCAGGTGAAAGACGTCTGTTCCGCCGACTCACTGGCGGGATTTGCTTGGGATCTGTTTACCGCCTGGCAGACCGCTGGCGCGCCGTCGAAAGAGAGTTGGGc                                                                                        <  4:211801/114‑1 (MQ=255)
                                                                                         tCCGGGATTATTGCAGGTGAAAGACGTCTGTTCCGCCGACTCACTGGCGGGATTTGCTTGGGATCTGTTTACCGCCTGGCAGACCGCTGGCGCGCCGTCGAAAGAGAGTTGGGc                                                                                        <  8:216564/114‑1 (MQ=255)
                                                                                           cGGGATTATTGCAGGTGAAAGACGTCTGTTCCGCCGACTCACTGGCGGGATTTGCCTGGGATCTGTTTACCGCCTGGCAGACCGCTGGCGCGCCGTCGAAAGAGAGTTGGGCGTTCACTGCGTTAGGCGTTCTCGGTAACGATGACAcc                                                   >  4:58773/1‑149 (MQ=255)
                                                                                            gggATTATTGCAGGTGAAAGACGTCTGTTCCGCCGACTCACTGGCGGGATTTGCCTGGGATCTGTTTACCGCCTGGCAGACCGCTGGCGCGCCGTCGAAAGAGAGTTGGGCGTTCACTGCGTTAGGCGTTCTCg                                                                 >  3:365880/1‑134 (MQ=255)
                                                                                            gggATTATTGCAGGTGAAAGACGTCTGTTCCGCCGACTCACTGGCGGGATTTGCCTGGGATCTGTTTACCGCCTGGCAGACCGCTGGCGCGCCGTCGAAAGAGAGTTGGGCGTTCACTGCGTTAGGCGTTCTCg                                                                 <  4:365880/134‑1 (MQ=255)
                                                                                             ggATTATTGCAGGTGAAAGACGTCTGTTCCGCCGACTCACTGGCGGGATTTGCCTGGGATCTGTTTACCGCCTGGCAGACCGCTGGCGCGCCGTCGAAAGAGAGTTGGGCGTTCACTGCGTTAGGCGTTCTCGGTAACGATGACAccgc                                                 >  4:359197/1‑149 (MQ=255)
                                                                                               attattGCAGGTGAAAGACGTCTGTTCCGCCGACTCACTGGCGGGATTTGCCTGGGATCTGTTTACCGCCTGGCAGACCGCTGGCGCGCCGTCGAAAGAGAGTTGGGCGTTCACTGCGTTAGGCGTTCTCGGTAACGATGACAccgccc                                               <  7:1011/149‑1 (MQ=255)
                                                                                                      cAGGTGAAAGACGTCTGTTCCGCCGACTCACTGGCGGGATTTGCCTGGGATCTGTTTACCGCCTGGCAGACCGCTGGCGCGCCGTCGAAAGAGAGTTGGGCGTTCACTGCGTTAGGCGTTCTCGGTAACGATGACACCGCCCGCAAACt                                        >  7:73436/1‑149 (MQ=255)
                                                                                                      cAGGTGAAAGACGTCTGTTCCGCCGACTCACTGGCGGGATTTGCCTGGGATCTGTTTACCGCCTGGCAGACCGCTGGCGCGCCGTCGAAAGAGAGTTGGGCGTTCACTGCGTTAGGCGTTCGCGGTAACGATGACACCGCCCGCAAACt                                        >  5:383605/1‑149 (MQ=255)
                                                                                                            aaaGACGTCTGTTCCGCCGACTCACTGGCGGGATTTGCCTGGGATCTGTTTACCGCCTGGCAGACCGCTGGCGCGCCGTCGAAAGAGAGTTGGGCGTTCACTGCGTTAGGCGGTCTCGGTAACGATGACACCGCCCGCAAACTGACGc                                   <  1:28910/148‑1 (MQ=255)
                                                                                                              aGACGTCTGTTCCGCCGACTCACTGGCGGGATTTGCCTGGGATCTGTTTACCGCCTGGCAGACCGCTGGCGCGCCGTCGAAAGAGAGTTGGGCGTTCACTGCGTTAGGCGTTCTCGGTAACGATGACACCGCCCGCAAACTGACGCCAt                                <  1:168165/149‑1 (MQ=255)
                                                                                                               gACGTCTGTTCCGCCGACTCACTGGCGGGATTTGCCTGGGATCTGTTTACCGCCTGGCAGACCGCTGGCGCGCCGTCGAAAGAGAGTTGGGCGTTCACTGCGTTAGGCGTTCTCGGTAACGATGACACCGCCCGCAAACTGACGCCAtt                               <  6:340804/149‑1 (MQ=255)
                                                                                                                aCGTCTGTTCCGCCGACTCACTGGCGGGATTTGCCTGGGATCTGTTTACCGCCTGGCAGCCCGCTGGCGCGCCGTCGAAAGAGAGTTGGGCGTTCACTGCGTTAGGCGTTCTCGGTAACGATGACCCCGCCCGCCAACTGACGCCAtt                               >  6:343545/1‑148 (MQ=255)
                                                                                                                  gTCTGTTCCGCCGACTCACTGGCGGGATTTGCCTGGGATCTGTTTACCGCCTGGCAGACCGCTGGCGCGCCGTCGAAAGAGAGTTGGGCGTTCACTGCGTTAGGCGTTCTCGGTAACGATGa                                                       >  5:160811/1‑122 (MQ=255)
                                                                                                                  gTCTGTTCCGCCGACTCACTGGCGGGATTTGCCTGGGATCTGTTTACCGCCTGGCAGACCGCTGGCGCGCCGTCGAAAGAGAGTTGGGCGTTCACTGCGTTAGGCGTTCTCGGTAACGATGa                                                       <  6:160811/122‑1 (MQ=255)
                                                                                                                  gTCTGTTCCGCCGACTCACTGGCGGGATTTGCCTGGGATCTGTTTACCGCCTGGCAGACCGCTGGCGCGCCGTCGAAAGAGAGTTGGGCGTTCACTGCGTTAGGCGTTCTCGGTAACGATGACACCGCCCGCAAACTGACGCCATTaa                             >  1:37589/1‑148 (MQ=255)
                                                                                                                   tCTGTTCCGCCGACTCACTGGCGGGATTTGCCTGGGATCTGTTTACCGCCTGGCAGACCGCTGGCGCGCCGTCGAAAGAGAGTTGGGCGTTCACTGCGTTAGGCGTTCTCGGTAACGATGACACCGCCCGCAAACTGACGCCATTAATa                           <  3:163160/149‑1 (MQ=255)
                                                                                                                             cgACTCACTGGCGGGATTTGCCTGGGATCTGTTTACCGCCTGGCAGACCGCTGGCGCGCCGTCGAAAGAGAGTTGGGCGTTCACTGCGTTAGGCGTTCTCGGTAACGATGACACCGCCCGCAAACTGACGCCATTAATACGCgcctggc                 >  7:387341/1‑149 (MQ=255)
                                                                                                                                cTCACTGGCGGGATTTGCCTGGGATCTGTTTACCGCCTGGCAGACCGCTGGCGCGCCGTCGAAAGAGAGTTTGGCGTTCACTGCGTTAGGCGTTCTCTGTAACGATGACACCGCCCGCAAACTGACGCCATTAATACGCgcctggcctg              <  8:272538/149‑1 (MQ=255)
                                                                                                                                 tCACTGGCGGGATTTGCCTGGGATCTGTTTACCGCCTGGCAGACCGCTGGCGCGCCGTCGAAAGAGAGTTGGGCGTTCACTGCGTTAGGCGTTCTCGGTAACGATGACACCGCCCGCAAACTGACGCCATTAATACGCGcctggcctgg             >  4:99192/1‑149 (MQ=255)
                                                                                                                                   aCTGGCGGGATTTGCCTGGGATCTGTTTACCGCCTGGCAGACCGCTGGCGCGCCGTCGAAAGAGAGTTGGGCGTTCACTGCGTTAGGCGTTCTCGGTAACGATGACACCGCCCGCAAACTGACGCCATTAATACGCGCCTGGCCTGGt            >  2:197825/1‑148 (MQ=255)
                                                                                                                                     tGGCGGGATTTGCCTGGGATCTGTTTACCGCCTGGCAGACCGCTGGCGCGCCGTCGAAAGAGAGTTGGGCGTTCACTGCGTTAGGCGTTCTCGGTAACGATGACACCGCCCGCAAACTGACGCCATTAATACGCGCCTGGCCTGGTGaa         >  1:83293/1‑149 (MQ=255)
                                                                                                                                            aTTTGCCTGGGATCTGTTTACCGCCTGGCAGACCGCTGGCGCGCCGTCGAAAGAGAGTTGGGCGTTCACTGCGTTAGGCGTTCTCGGTAACGATGACACCGCCCGCAAACTGACGCCATTAATACGCGCCTGGCCTGGTGAATCCCAGc  <  8:64565/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                              
GTATTAAAAGCAAATGCCCAATCACTGCCAGATAGCGCCCTCCTCCACCTCGGTGAAATGCTCCGCTTCCCTCAGGAAGAGGCTCTGTATCCGGGATTATTGCAGGTGAAAGACGTCTGTTCCGCCGACTCACTGGCGGGATTTGCCTGGGATCTGTTTACCGCCTGGCAGACCGCTGGCGCGCCGTCGAAAGAGAGTTGGGCGTTCACTGCGTTAGGCGTTCTCGGTAACGATGACACCGCCCGCAAACTGACGCCATTAATACGCGCCTGGCCTGGTGAATCCCAGC  >  NC_000913/2198853‑2199141

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