Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,200,173 C→T 5.3% pseudogene (852/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,200,1730CT 5.3% 241.8 / 2.4 75pseudogene (852/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 C (40/30);  new base T (2/2);  total (43/32)
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.06e-01

CGGATGGATGATCAAACCTTTGGGTCGCTGGTCGTTAATCATGGAAATCGATGAAGGCTTTGCGGTAGGCATGTCGCCAGCCGAACTCAGCGCTGAGCAGCTCTTAAGCAAGCTGTGGCTATGGGAAGGCAAAGCAGAAAGATATGGCTGGGGGAGTAATTCAACACAGGAAGCGCAGTTCTCCGTAATCGATGCCATCACCGCCAGCGAGCTAATTAACGATATTGAAGCCCTGTTTGAATAAGGAAAGAGCATGGACAAGGAATTACCGTGGCTGGCGGATAACGCCCAAC  >  NC_000913/2200026‑2200318
                                                                                                                                                   |                                                                                                                                                 
cGGATGGATGATCAAACATTTGGGTCGCTGGTCGTTAATCATGGAAATCGATGAAGGCTTTGCGGTATGCATGTCGCCAGCCGAACTCAGCGCTGAGCAGCTCTTAAGCAAGCTGTGGCTATGGGAAGGCAAAGCAGAAAGATATGGCt                                                                                                                                                  >  4:43519/1‑149 (MQ=255)
    tggatgATCAAACCTTTGGGTCGCTGGTCGTTAATCATGGAAATCGATGAAGGCTTTGCGGTAGGCATGTCGCCAGCCGAACTCAGCGCTGAGCAGCTCTTAAGCAAGCTGTGGCTATGGGAAGGCAAAGCAGAAAGATATGGCTgggg                                                                                                                                              <  7:283650/149‑1 (MQ=255)
      gatgatCAAACCTTTGGGTCGCTGGTCGTTAATCATGGAAATCGATGAAGGCTTTGCGGTAGGCATGTCGCCAGCCGAACTCAGCGCTGAGCAGCTCTTAAGCAAGCTGTGGCTATGGGAAGGCAAAGCAGAAAGATATGGCTgggg                                                                                                                                              >  2:206913/1‑147 (MQ=255)
           tCAAACCTTTGGGTCGCTGGTCGTTAATCATGGAAATCGATGAAGGCTTTGCGGTAGGCATGTCGCCAGCCGAACTCAGCGCTGAGCAGCTCTTAAGCAAGCTGTGGCTATGGGAAGGCAAAGCAGAAAGATATGGCTGGGGGAGTAAt                                                                                                                                       <  1:176595/149‑1 (MQ=255)
           tCAAACCTTTGGGTCGCTGGTCGTTAATCATGGAAATCGATGAAGGCTTTGCGGTAGGCATGTCGCCAGCCGAACTCAGCGCTGAGCAGCTCTTAAGCAAGCTGTGGCTATGGGAAGGCAAAGCAGAAAGATATGGCTGGGGGAGTAAt                                                                                                                                       >  3:223930/1‑149 (MQ=255)
           tCAAACCTTTGGGTCGCTGGTCGTTAATCATGGAAATCGATGAAGGCTTTGCGGTAGGCATGTCGCCAGCCGAACTCAGCGCTGAGCAGCTCTTAAGCAAGCTGTGGCTATGGGAAGGCAAAGCAGAAAGATATGGCTGGGGGAGTAAt                                                                                                                                       >  8:188857/1‑149 (MQ=255)
           tCAAACCTTTGGGTCGCTGGTCGTTAATCATGGAAATCGATGAAGGCTTTGCGGTAGGCATGTCGCCAGCCGAACTCAGCGCTGAGCAGCTCTTAAGCAAGCTGTGGCTATGGGAAGGCAAAGCAGAAAGATATGGCTGGGGGAGTAAt                                                                                                                                       >  7:248694/1‑149 (MQ=255)
           tCAAAACTTTCGGTCGCTGGTAGTTAATCATGGAAATCGATGAAGGCTTTGCTGTAGGCATGTCTCCAGCCGAACTCAGAGCTGAGCACCTCTTAAGCAAGGTGTGGCGATGGGAAGGCAAAGCAGAAAGATATGGCTGGGGGAGTAAt                                                                                                                                       >  6:314423/1‑149 (MQ=255)
             aaaCCTTTGGGTCGCTGGTCGTTAATCATGGAAATCGATGAAGGCTTTGCGGTAGGCATGTCGCCAGCCGAACTCAGCGCTGAGCAGCTCTTAAGCAAGCTGTGGCTATGGGAAGGCAAAGCAGAAAGATATGGCTGGGGGAGTAAtt                                                                                                                                      >  6:102614/1‑148 (MQ=255)
             aaaCCTTTGGGTCGCTGGTCGTTAATCATGGAAATCGATGAAGGCTTTGCGGTAGGCATGTCGCCAGCCGAACTCAGCGCTGAGCAGCTCTTAAGCAAGCTGTGGCTATGGGAAGGCAAAGCAGAAAGATACGGCTGGGGGAGTAAtt                                                                                                                                      >  8:259022/1‑148 (MQ=255)
               aCCTTTGGGTCGCTGGTCGGTAATCATGGAAATCGATGAAGGCTTTGCGGTAGGCATGTCGCCAGCCGAACTCAGCGCTGAGCAGCGCTTAAGCAAGCTGTGGCTATGGGAAGGCAAAGCAGAAAGATATGGCTGGGGGAGTAATTCaa                                                                                                                                   <  8:372364/149‑1 (MQ=255)
                    tGGGTCGCTGGTCGTTAATCATGGAAATCGATGAAGGCTTTGCGGTAGGCATGTCGCCAGCCGAACTCAGCGCTGAGCAGCTCTTAAGCAAGCTGTGGCTATGGGAAGGCAAAGCAGAAAGATATGGCTGGGGGAGTAATTCAAcaca                                                                                                                               <  1:239552/148‑1 (MQ=255)
                       gTCTCTGGTCGTTAATCATGGAAATCGATGAAGGCTTTGCGGTAGGCATGTCGCCAGCCGAACTCAGCGCTGAGCAGCTCTTAAGCAAGCTGTGGCTATGGGAAGGCAAAGCAGAAAGATATGGCTGGGGGTGTAATTCAACACAgg                                                                                                                             >  1:116821/1‑147 (MQ=255)
                              gTCGTTAATCATGGAAATCGATGAAGGCTTTGCGGTAGGCATGTCGCCAGCCGAACTCAGCGCTGAGCAGCTCTTAAGCAAGCTGTGGCTATGGGAAGGCAAAGCAGAAAGATATGGCTGGGGGAGTAATTCAACACAGGAAGCGCAg                                                                                                                     <  8:326142/148‑1 (MQ=255)
                                 gTTAATCATGGAAATCGATGAAGGCTTTGCGGTAGGCATGTCGCCAGCCGAACTCAGCGCTGAGCAGCTCTTAAGCAAGCTGTGGCTATGGGAAGGCAAAGCAGAAAGATATGGCTGGGGGAGTAATTCAACACAGGAAGCGCAGTtct                                                                                                                 >  6:157524/1‑149 (MQ=255)
                                 gTTAATCATGGAAATCGATGAAGGCTTTGCGGTAGGCATGTCGCCAGCCGAACTCAGCGCTGAGCAGCTCTTAAGCAAGCTGTGGCTATGGGAAGGCAAAGCAGAAAGATATGGCTGGGGGAGTAATTCAACACAGGAAGCGCAGTtct                                                                                                                 >  7:181352/1‑149 (MQ=255)
                                    aaTCATGGAAATCGATGAAGGCTTTGCGGTAGGCATGTCGCCAGCCGAACTCAGCGCTGAGCAGCTCTTAAGCAAGCTGTGGCTATGGGAAGGCAAAGCAGAAAGATATGGCTGGGGGAGTAATTCAACACAGGAAGCGCAGTTCTcc                                                                                                               <  7:6324/148‑1 (MQ=255)
                                    aaTCATGGAAATCGATGAAGGCTTTGCGGTAGGCATGTCGCCAGCCGAACTCAGCGCTGAGCAGCTCTTAAGCAAGCTGTGGCTATGGGAAGGCAAAGCAGAAAGATATGGCTGGGGGAGTAATTCAACACAGGAAGCGCAGTTCTcc                                                                                                               <  4:138094/148‑1 (MQ=255)
                                    aaTCATGGAAATCGATGAAGGCTTTGCGGTAGGCATGTCGCCAGCCGAACTCAGCGCTGAGCAGCTCTTAAGCAAGCTGTGGCTATGGGAAGGCAAAGCAGAAAGATATGGCTGGGGGAGTAATTCAACACAGGAAGCGCAGTTCTCCg                                                                                                              >  1:176668/1‑149 (MQ=255)
                                     aTCATGGAAATCGATGAAGGCTTTGCGGTAGGCATGTCGCCAGCCGAACTCAGCGCTGAGCAGCTCTTAAGCAAGCTGTGGCTATGGGAAGGCAAAGCAGAAAGATATTGCTGGGGGAGTAATTCAACACAGGAAGCGCAGTTCTCCg                                                                                                              >  7:277280/1‑148 (MQ=255)
                                     aTCATGGAAATCGATGAAGGCTTTGCGGTAGGCATGTCGCCAGCCGAACTCAGCGCTGAGCAGCTCTTAAGCAAGCTGTGGCTATGGGAAGGCAAAGCAGAAAGATATGGCTGGGGGAGTAATTCAACACAGGAAGCGCAGTTCTCCg                                                                                                              <  2:176668/148‑1 (MQ=255)
                                        aTGGAAATCGATGAAGGCTTTGCGGTAGGCATGTCGCCAGCCGAACTCAGCGCTGAGCAGCTCTTAAGCAAGCTGTGGCTATGGGAAGGCAAAGCAGAAAGATATGGCTGGGGGAGTAATTCAACACAGGAAGCGCAGTTATCCGTaa                                                                                                           <  5:359631/148‑1 (MQ=255)
                                                cGATGAAGGCTTTGCGGTAGGCATGTCGCCAGCCGAACTCAGCGCTGAGCAGCTCTTAAGCAAGCTGTGGCTATGGGAAGGCAAAGCAGAAAGATATGGCTGGGGGAGTAATTCAACACAGGAAGCGCAGTTCTCCGTAATCGATGcc                                                                                                   >  2:330009/1‑148 (MQ=255)
                                                      aGGCTTTGCGGTAGGCATGTCGCCAGCCGAACTCAGCGCTGAGCAGCTCTTAAGCAAGCTGTGGCTATGGGAAGGCAAAGCAGAAAGATATGGCTGGGGGAGTAATTCAACACAGGAAGCGCAGTTCTCCGTAATCGATGCCATCACCg                                                                                            >  8:16864/1‑149 (MQ=255)
                                                        gCTTTGCGGTAGGCATGTCGCCAGCCGAACTCAGCGCTGAGCAGCTCTTAAGCAAGCTGTGGCTATGGGAAGGCAAAGCAGAAAGATATGGCTGGGGGAGTAATTCAACACAGGAAGCGCAGTTCTCCGTAATCGATGCCATCACCGc                                                                                           >  1:234562/1‑148 (MQ=255)
                                                             gCGGTAGGCATGTCGCCAGCCGAACTCAGCGCTGAGCAGCTCTTAAGCAAGCTGTGGCTATGGGAAGGCAAAGCAGAAAGATATGGCTGGGGGAGTAATTCAACACAGGAAGCGCAGTTCTCCGTAATCGATGCCATCACCGCCAGc                                                                                       >  8:117940/1‑147 (MQ=255)
                                                                gTAGGCATGTCGCCAGCCGAACTCAGCGCTGAGCAGCTCTTAAGCAAGCTGTGGCTATGGGAAGGCAAAGCAGAAAGATATGGCTGGGGGAGTAATTCAACACAGGAAGCGCAGTTCTCCGTAATCGATGCCATCACCGCCAGCGAGc                                                                                   >  4:344324/1‑148 (MQ=255)
                                                                gTAGGCATGTCGCCAGCCGAACTCAGCGCTGAGCAGCTCTTAAGCAAGCTGTGGCTATGGGAAGGCAAAGCAGAAAGATATGGCTGGGGGAGTAATTCAACACAGGAAGCGCAGTTCTCCGTAATCGATGCCATCACCGCCAGCGAGc                                                                                   >  1:154553/1‑148 (MQ=255)
                                                                gTAGGCATGTCGCCAGCCGAACTCAGCGCTGAGCAGCTCTTAAGCAAGCTGTGGCTATGGGAAGGCAAAGCAGAAAGAGATGGCTGGGGGAGTATTTCAACCCAGGAAGCGCAGTTCTCCGTAATCGATGCCATCCCCGCCAGCGAGc                                                                                   >  1:278945/1‑148 (MQ=255)
                                                                 tAGGCATGTCGCCAGCCGAACTCAGCGCTGAGCAGCTCTTAAGCAAGCTGTGGCTATGGGAAGGCAAAGCAGAAAGATATGGTTGGGGGAGTAATTCAACACAGGAAGCGCAGTTCTCCGTAAt                                                                                                          >  6:367606/1‑124 (MQ=255)
                                                                 tAGGCATGTCGCCAGCCGAACTCAGCGCTGAGCAGCTCTTAAGCAAGCTGTGGCTATGGGAAGGCAAAGCAGAAAGATATGGTTGGGGGAGTAATTCAACACAGGAAGCGCAGTTCTCCGTAAt                                                                                                          <  5:367606/124‑1 (MQ=255)
                                                                    gCATGTCGCCAGCCGAACTCAGCGCTGAGCAGCTCTTAAGCAAGCTGTGGCTATGGGAAGGCAAAGCAGAAAGATATGGCTGGGGGAGTAATTCAACACAGGAAGCGCAGTTCTCCGTAATCGATGCCATCACCGCCAGCGAGCTAAt                                                                               >  4:116683/1‑148 (MQ=255)
                                                                    gCATGTCGCCAGCCGAACTCAGCGCTGAGCAGCTCTTAAGCAAGCTGTGGCTATGGGAAGGCAAAGCAGAAAGATATGGCGGGGGGAGTAATTCACCACAGGAAGCGCACTTCTCCGTAATCGATGCCATCACCGACCGCGAGCTAAt                                                                               >  6:372270/1‑148 (MQ=255)
                                                                     cATGTCGCCAGCCGAACTCAGCGCTGAGCAGCTCTTAAGCAAGCTGTGGCTATGGGAAGGCAAAGCAGAAAGATATGGCTGGGGGAGTAATTCAACACAGGAAGCGCAGTTCTCCGTAATCGATGCCATCACCGCCAGCGAGCTAAtt                                                                              >  2:46950/1‑148 (MQ=255)
                                                                        gTCGCCAGCCGAACTCAGCGCTGAGCAGCTCTTAAGCAAGCTGTGGCTATGGGAAGCCAAAGCAGAAAGATATGGCTGGGGGAGTAATTCAACACAGGAAGCGCAGTTCTCCGTAATCGATGCCATCACCGCCAGCGAGCTAATTAAc                                                                           >  7:24004/1‑148 (MQ=255)
                                                                             cAGCCGAACTCAGCGCTGAGCAGCTCTTAAGCAAGCTGTGGCTATGGGAAGGCAAAGCAGAAAGATATGGTTGGGGGAGTAATTCAACACAGGAAGCGCAGTTCTCCGTAATCGATGCCATCACCGCCAGCGAGCTAATTAACGATAtt                                                                     >  2:41071/1‑149 (MQ=255)
                                                                             cAGCCGAACTCAGCGCTGAGCAGCTCTTAAGCAAGCTGTGGCTATGGGAAGGCAAAGCAGAAAGATATGGCTGGGGGAGTAATTCAACACAGGAAGCGCAGTTCTCCGTAATCGATGCCATCACCGCCAGCGAGCTAATTAACGATAtt                                                                     <  6:130698/149‑1 (MQ=255)
                                                                               gCCGAACTCAGCGCTGAGCAGCTCTTAAGCAAGCTGTGGCTATGGGAAGGCAAAGCAGAAAGATATGGCTGGGGGAGTAATTCAACACAGGAAGCGCAGTTCTCCGTAATCGATGCCATCACCGCCAGCGAGCTAATTAACGATATTGa                                                                   <  1:46950/149‑1 (MQ=255)
                                                                                ccGAACTCAGCGCTGAGCAGCTCTTAAGCAAGCTGTGGCTATGGGAAGGCAAAGCAGAAAGATATGGCTGGGGGAGTAATTCAACACAGGAAGCGCAGTTCTCCGTAATCGATGCCATCACCGCCAGCGAGCTAATTAACGATATTGaa                                                                  <  3:43519/149‑1 (MQ=255)
                                                                                ccGAACTCAGCGCTGAGCAGCTCTTAAGCAAGCTGTGGCTATGGGAAGGCAAAGCAGAAAGATATGGCTGGGGGAGTAATTCAACACAGGAAGCGCAGTTCTCCGTAATCGATGCCATCACCGCCAGCGAGCTAATTAACGATATTGaa                                                                  <  7:117940/149‑1 (MQ=255)
                                                                                      tCAGCGCTGAGCAGCTCTTAAGCAAGCTGTGGCTATGGGAAGGCAAAGCAGAAAGATATGGCTGGGGGAGTAATTCAACACAGGAAGCGCAGTTCTCCGTAATCGATGCCATCACCGCCAGCGAGCTAATTAACGATATTGAAGCCCTg                                                            >  5:178390/1‑149 (MQ=255)
                                                                                         gcgcTGAGCAGCTCTTAAGCAAGCTGTGGCTATGGGAAGGCAAAGCAGAAAGATATGGCTGGGGGAGTAATTCAACACAGGAAGCGCAGTTCTCCGTAATCGATGCCATCACCGCCAGCGAGCTAATTAACGATATTGAAGCCCTGttt                                                         <  8:248694/149‑1 (MQ=255)
                                                                                             tGAGCAGCTCTTAAGCAAGCTGTGGCTATGGGAAGGCAAAGCAGAAAGATATGGCTGGGGGAGTAATTCAACACAGGAAGCGCAGTTCTCCGTAATCGATGCCATCACCGCCAGCGAGCTAATTAACGATATTGAAGCCCTGTTTGaa                                                      >  6:3872/1‑148 (MQ=255)
                                                                                                gcagcTCTTAAGCAAGCTGTGGCTATGGGAAGGCAAAGCAGAAAGATATGGCTGGGGGAGTAATTCAACACAGGAAGCGCAGTTCTCCGTAATCGATGCCATCACCGCCAGCGAGCTAATTAACGATATTGAAGCCCTGTTTGAATaa                                                   <  7:188857/148‑1 (MQ=255)
                                                                                                 cagcTCTTAAGCAAGCTGTGGCTATGGGAAGGCAAAGCAGAAAGATATGGCTGGGGGAGTAATTCAACACAGGAAGCGCAGTTCTCCGTAATCGATGCCATCACCGCCAGCGAGCTAATTAACGATATTGAAGCCCTGTTTGAATAAgg                                                 <  7:65395/149‑1 (MQ=255)
                                                                                                   gcTCTTAAGCAAGCTGTGGCTATGGGAAGGCAAAGCAGAAAGATATGGCTGGGGGTGTAATTCAACACAGGAAGCGCAGTTCTCCGTAATCGATGCCATCACCGCCAGCGAGCTAATTAACGATATTGAAGCCCTGTTTGAATAAGGaa                                               <  2:116821/149‑1 (MQ=255)
                                                                                                   gcTCTTAAGCAAGCTGTGGCTATGGGAAGACAAAGCAGAAAGATATGGCTGGGGGAGTAATTCAACACAGGAAGCGCAGTTCTCCGTAATCGATGCCATCACCGCCAGCGAGCTAATTAACGATATTGAAGCCCTGTTTGAATAAGGaa                                               <  8:24004/149‑1 (MQ=255)
                                                                                                    ctctTAAGCAAGCTGTGGCTATGGGAAGGCAAAGCAGAAAGATATGGCTGGGGGAGTAATTCAACACAGGAAGCGCAGTTCTCCGTAATCGATGCCATCACCGCCAGCGAGCTAATTAACGATATTGaa                                                                  <  1:67186/129‑1 (MQ=255)
                                                                                                    ctctTAAGCAAGCTGTGGCTATGGGAAGGCAAAGCAGAAAGATATGGCTGGGGGAGTAATTCAACACAGGAAGCGCAGTTCTCCGTAATCGATGCCATCACCGCCAGCGAGCTAATTAACGATATTGAAGCCCTGTTTGAATAAGGa                                                >  1:268743/1‑147 (MQ=255)
                                                                                                    ctctTAAGCAAGCTGTGGCTATGGGAAGGCAAAGCAGAAAGATATGGCTGGGGGAGTAATTCAACACAGGAAGCGCAGTTCTCCGTAATAGATGCCATCACCGCCAGCGAGCTAATTAACGATATTGaa                                                                  >  2:67186/1‑129 (MQ=255)
                                                                                                      ctTAAGCAAGCTGTGGCTATGGGAAGGCAAAGCAGAAAGATATGGCTGGGGGAGTAATTCAACACAGGAAGCGCAGTTCTCCGTAATCGATGCCATCACCGCCAGCGAGCTAATTAACGATATTGAAGCCCTGTTTGAATAAGGaaa                                              >  6:174736/1‑147 (MQ=255)
                                                                                                          agcaagcTGTGGCTATGGGAAGGCAAAGCAGAAAGATATGGCTGGGGGAGTAATTCAACACAGGAAGCGCAGTTCTCCGTAATCGATGCCATCACCGCCAGCGAGCTAATTAACGATATTGAAGCCCTGTTTGAATAAGGAAAGAGCAt                                        <  5:174736/149‑1 (MQ=255)
                                                                                                          agcaagcTGTGGCTATGGGAAGGCAAAGCAGAAAGATATGGCTGGGGGAGTAATTCAACACAGGAAGCGCAGTTCTCCGTAATCGATGCCATCACCGCCAGCGAGCTAATTAACGATATTGAAGCCCTGTTTGAATAAGGAAAGAGCAt                                        <  2:376575/149‑1 (MQ=255)
                                                                                                           gcaagcTGTGGCTATGGGAAGGCAAAGCAGAAAGATATGGCTGGGGGAGTAATTCAACACAGGAAGCGCAGTTCTCCGTAATCGATGCCATCACCGCCAGCGAGCTAATTAACGATATTGAAGCCCTGTTTGAATAAGGAAAGAGCATg                                       >  4:318216/1‑149 (MQ=255)
                                                                                                           gcaagcTGTGGCTATGGGAAGGCAAAGCAGAAAGATATGGCTGGGGGAGTAATTCAACACAGGAAGCGCAGTTCTCCGTAATCGATGCCATCACCGCCAGCGAGCTAATTAACGATATTGAAGCCCTGTTTGAATAAGGAAAGAGCATg                                       <  1:330009/149‑1 (MQ=255)
                                                                                                             aagcTGTGGCTATGGGAAGGCAAAGCAGAAAGATATGGCTGGGGGAGTAATTCACCACAGGAAGCGCAGTTCTCCGTAATCGATGCCATCACCGACAGCGAGCTAATTAACGATATTGAAGCACTGTTTGAATAAGGAAAGAGCAtcg                                      >  7:382267/1‑146 (MQ=255)
                                                                                                             aagcTGTGGCTATCGGAAGGCAAAGCAGAAAGATATGGCTGGGGGAGTAATTCAACACAGGAAGCGCAGATCTCCGTAATCGATGCCATCACCGCCAGCGAGCTAATTAACGATATTGAAGCCCTGTTTGAATAAGGAAAGAGCATgg                                      >  7:387639/1‑148 (MQ=255)
                                                                                                                   tgGCTATGGGAAGGCAAAGCAGAAAGATATGGCTGGGGGAGTAATTCAACACAGGAAGCGCAGTTCTCCGTAATCGATGCCATCACCGCCAGCGAGCTAATTAACGATATTGAAGCCCTGTTTGAATAAGGAAAGAGCATGGACAAgg                                <  1:206913/148‑1 (MQ=255)
                                                                                                                   tgGCGATGGGAAGGCAAAGCAGAACGATATGGCGGGGGGAGTAATTCACCACAGGAAGCGCAGTTCTCCGTAATCGATGCCATCACCGCCAGCGAGCTAATTAACGATATTGAAGCCCTGTTTGAATAAGGAAAGAGCATGGACAAGGa                               <  2:278945/149‑1 (MQ=255)
                                                                                                                    ggCTATGGGAAGGCAAAGCAGAAAGATACGGCTGGGGGAGTAATTCAACACAGGAAGCGCAGTTCTCCGTAATCGATGCCATCACCGCCAGCGAGCTAATTAACGATATTGAAGCCCTGTTTGAATAAGGAAAGAGCATGGACAAGGaa                              <  7:259022/149‑1 (MQ=255)
                                                                                                                      cTATGGGAAGGCAAAGCAGAAAGATATGGCTGGGGGAGTAATTCAACACAGGAAGCGCAGTTCTCCGTAATCGATGCCATCACCGCCAGCGAGCTAATTAACGa                                                                         >  8:31190/1‑104 (MQ=255)
                                                                                                                      cTATGGGAAGGCAAAGCAGAAAGATATGGCTGGGGGAGTAATTCAACACAGGAAGCGCAGTTCTCCGTAATCGATGCCATCACCGCCAGCGAGCTAATTAACGATATTGAAGCCCTGTTTGAATAAGGAAAGAGCATGGACAAGGAAt                             >  4:105191/1‑148 (MQ=255)
                                                                                                                         tGGGAAGGCAAAGCAGAAAGATATGGCTGGGGGAGTAATTCAAGACAGGAAGCGCAGTTCTCCGTAATCGATGCCATCACCGCCAGCGAGCTAATTAACGATATTGAAGCCCTGTTTGAATAAGGAAAGAGCATGGACAAGGAATTAc                          <  5:314423/148‑1 (MQ=255)
                                                                                                                           ggAAGGCAAAGCAGAAAGATATGGCTGGGGGAGTAATTCAACACAGGAAGCGCAGTTCTCCGTAATCGATGCCATCACCGCCAGCGAGCTAATTAAAGATATTGAAGCCCTGGTTGATTAAGGAAAGAGCAGGGACAAGGAATTCCCGt                       >  8:93017/1‑149 (MQ=255)
                                                                                                                           ggAAGGCAAAGCAGAAAGATATGGATGGGTGAGTAATTCAACACAGGAAGCGCAGTTCTCCGTAATCGATGCCATCACCGCCAGCGAGCTAATTAACGATATTGAAGCCCTGTTTGAATAAGGAAAGAGCAGGGACAAGGAATTAcctt                       >  3:307653/1‑147 (MQ=255)
                                                                                                                                  aaaGCAGAAAGATATGGCTGGGGGAGTAATTCAACACAGGAAGCGCAGTTCTCCGTAATCGATGCCATCACCGCCAGCGAGCTAATTAACGATATTGAAGCCCTGTTTGAATAAGGAAAGAGCATGGACAAGGAATTACCGtggctggc                >  6:339605/1‑149 (MQ=255)
                                                                                                                                  aaaGCAGAAAGATATGGCTGGGGGAGTAATTCAACACAGGAAGCGCAGTTCTCCGTAATCGATGCCATCACCGCCAGCGAGCTAATTAACGATATTGAAGCCCTGTTTGAATAAGGAAAGAGCATGGACAAGGAATTACCGtggctggc                <  3:105191/149‑1 (MQ=255)
                                                                                                                                  aaaGCAGAAAGATATGGCTGGGGGAGTAATTCAACACAGGAAGCGCAGTTCTCCGTAATCGATGCCATCACCGCCAGCGAGCTAATTAACGATATTGAAGCCCTGTTTGAATAAGGAAAGAGCATGGACAAGGAATTACCGtggctggc                <  3:116683/149‑1 (MQ=255)
                                                                                                                                     gCAGAAAGATATGGTTGGGGGAGTAATTCAACACAGGAAGCGCAGTTCTCCGTAATCGATGCCATCACCGCCAGCGAGCTAATTAACGATATTGAAGCCCTGTTTGAATAAGGAAAGAGCATGGACAAGGAATTACCGTGGCTGGCGGa             <  1:41071/149‑1 (MQ=255)
                                                                                                                                      cAGAAAGATATGGCTGGGGGAGTAATTCAACACAGGAAGCGCAGTTCTCCGTAATCGATGCCATCACCGCCAGCGAGCTAATTAACGATATTGAAGCCCTGTTTGAATAAGGAAAGAGCATGGACAAGGAATTACCGTGGCTGGCGGa             >  7:262725/1‑148 (MQ=255)
                                                                                                                                          aaGATATGGCTGGGGGAGTAATTCAACACAGGAAGCGCAGTTCTCCGTAATCGATGCCATCACCGCCAGCGAGCTAATTAACGATATTGAAGCCCTGTTTGAATAAGGAAAGAGCATGGACAAGGAATTACCGTGGCTGGCGGATAACg        <  5:3872/149‑1 (MQ=255)
                                                                                                                                               atGGCTGGGGGAGTAATTCAACACAGGAAGCGCAGTTCTCCGTAATCGATGCCATCACCGCCAGCGAGCTAATTAACGATATTGAAGCCCTGTTTGAATAAGGAAAGAGCATGGACAAGGa                               >  7:230969/1‑121 (MQ=255)
                                                                                                                                               atGGCTGGGGGAGTAATTCAACACAGGAAGCGCAGTTCTCCGTAATCGATGCCATCACCGCCAGCGAGCTAATTAACGATATTGAAGCCCTGTTTGAATAAGGAAAGAGCATGGACAAGGa                               <  8:230969/121‑1 (MQ=255)
                                                                                                                                               atGGCTGGGGGAGTAACTCAACACAGGAAGCGCAGTTCTCCGTAATCGATGCCATCACCGCCAGCGAGCTAATTAACGATATTGAAGCCCTGTTTGAATAAGGAAAGAGCATGGAc                                    >  6:331971/1‑116 (MQ=255)
                                                                                                                                                tGGCTGGGGGAGTAATTCAACACAGGAAGCGCAGTTCTCCGTAATCGATGCCATCACCGCCAGCGAGCTAATTAACGATATTGAAGCCCTGTTTGAATAAGGAAAGAGCATGGACAAGGAATTACCGTGGCTGGCGGATAACGCCCAAc  >  6:181204/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
CGGATGGATGATCAAACCTTTGGGTCGCTGGTCGTTAATCATGGAAATCGATGAAGGCTTTGCGGTAGGCATGTCGCCAGCCGAACTCAGCGCTGAGCAGCTCTTAAGCAAGCTGTGGCTATGGGAAGGCAAAGCAGAAAGATATGGCTGGGGGAGTAATTCAACACAGGAAGCGCAGTTCTCCGTAATCGATGCCATCACCGCCAGCGAGCTAATTAACGATATTGAAGCCCTGTTTGAATAAGGAAAGAGCATGGACAAGGAATTACCGTGGCTGGCGGATAACGCCCAAC  >  NC_000913/2200026‑2200318

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