Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,192,571 C→A 6.8% pseudogene (2619/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
*CP0092732,192,5710CA 6.8% 760.3 / 20.3 204pseudogene (2619/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 (55/135);  new base A (6/8);  total (61/143)
Fisher's exact test for biased strand distribution p-value = 3.63e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.89e-01

Only 100 of 222 total aligned reads displayed.

GGGATTATTGCAGGTGAAAGACGTCTGTTCCGCCGACTCACTGGCGGGATTTGCCTGGGATCTGTTTACCGCCTGGCAGACCGCTGGCGCGCCGTCGAAAGAGAGTTGGGCGTTCACTGCGTTAGGCGTTCTCGGTAACGATGACACCGCCCGCAAACTGACGCCATTAATACGCGCCTGGCCTGGTGAATCCCAGCATAAACGCGCCACCGTTGGGTTGGATATTCTCGCTGCTATCGGTAGTGATATCGCCCTTATGCAGCTTAACGGCATCGCCCAGAAACTGAAATTCAAA  >  CP009273/2192424‑2192718
                                                                                                                                                   |                                                                                                                                                   
gggATTATTGCAGGTGAAAGACGTCTGTTCCGCCGACTCACTGGCGGGATTTGCCTGGGATCTGTTTACCGCCTGGCAGACCGCTGGCGCGCCGTCGAAAGAGAGTTGGGCGTTCACTGCGTTAGGCGTTCTCGGTAACGATGACAcc                                                                                                                                                     >  2:2435574/1‑148 (MQ=255)
  gATTATTGCAGGTGAAAGACGTCTGTTCCGCCGACTCACTGGCGGGATTTGCCTGGGATCTGTTTACCGCCTGGCAGACCGCTGGCGCGCCGTCGAAAGAGAGTTGGGCGTTCACTGCGTTAGGCGTTCTCGGTAACGATGACAccgc                                                                                                                                                   >  2:149270/1‑148 (MQ=255)
  gATTATTGCAGGTGAAAGACGTCTGTTCCGCCGACTCACTGGCGGGATTTGCCTGGGATCTGTTTACCGCCTGGCAGACCGCTGGCGCGCCGTCGAAAGAGAGTTGGGCGTTCACTGCGTTAGGCGTTCTCGGTAACGATGACAccgc                                                                                                                                                   >  2:3110402/1‑148 (MQ=255)
    ttattGCAGGTGAAAGACGTCTGTTCCGCCGACTCACTGGCGGGATTTGCCTGGGATCTGTTTACCGCCTGGCAGACCGCTGGCGCGCCGTCGAAAGAGAGTTGGGCGTTCACTGCGTTAGGCGTTCTCGGTAACGATGACAccgccc                                                                                                                                                 >  1:2010523/1‑148 (MQ=255)
        tGCAGGTGAAAGACGTCTGTTCCGCCGACTCACTGGCGGGATTTGCCTGGGATCTGTTTACCGCCTGGCAGACCGCTGGCGCGCCGTCGAAAGAGAGTTGGGCGTTCACTGCGTTAGGCGTTCTCGGTAACGATGACACCGCCCGCaa                                                                                                                                             >  1:1909221/1‑148 (MQ=255)
         gCAGGTGAAAGACGTCTGTTCCGCCGACTCACTGGCGGGATTTGCCTGGGATCTGTTTACCGCCTGGCAGACCGCTGGCGCGCCGTCGAAAGAGAGTTGGGCGTTCACTGCGTTAGGCGTTCTCGGTAACGATGACACCGCCCGCaaa                                                                                                                                            <  2:1760850/148‑1 (MQ=255)
         gCAGGTGAAAGACGTCTGTTCCGCCGACTCACTGGCGGGATTTGCCTGGGATCTGTTTACCGCCTGGCAGACCGCTGGCGCGCCGTCGAAAGAGAGTTGGGCGTTCACTGCGTTAGGCGTTCTCGGTAACGATGACACCGCCCGCaaa                                                                                                                                            <  1:4406/148‑1 (MQ=255)
         gCAGGTGAAAGACGTCTGTTCCGCCGACTCACTGGCGGGATTTGCCTGGGATCTGTTTACCGCCTGGCAGACCGCTGGCGCGCCGTCGAAAGAGAGTTGGGCGTTCACTGCGTTAGGCGTTCTCGGTAACGATGACACCGCCCGCaaa                                                                                                                                            <  1:602464/148‑1 (MQ=255)
           aGGTGAAAGACGTCTGTTCCGCCGACTCACTGGCGGGATTTGCCTGGGATCTGTTTACCGCCTGGCAGACCGCTGGCGCGCCGTCGAAAGAGAGTTGGGCGTTCACTGCGTTAGGCGTTCTCGGTAACGATGACACCGCCCGCAAACt                                                                                                                                          <  2:1976984/148‑1 (MQ=255)
           aGGTGAAAGACGTCTGTTCCGCCGACTCACTGGCGGGATTTGCCTGGGATCTGTTTACCGCCTGGCAGACCGCTGGCGCGCCGTCGAAAGAGAGTTGGGCGTTCACTGCGTTAGGCGTTCTCGGTAACGATGACACCGCCCGCAAACt                                                                                                                                          <  2:3536800/148‑1 (MQ=255)
             gTGAAAGACGTCTGTTCCGCCGACTCACTGGCGGGATTTGCCTGGGATCTGTTTACCGCCTGGCAGACCGCTGGCGCGCCGTCGAAAGAGAGTTGGGCGTTCACTGCGTTAGGCGTTCTCGGTAACGATGACACCGCCCGCAAACTGa                                                                                                                                        <  2:218617/148‑1 (MQ=255)
             gTGAAAGACGTCTGTTCCGCCGACTCACTGGCGGGATTTGCCTGGGATCTGTTTACCGCCTGGCAGACCGCTGGCGCGCCGTCGAAAGAGAGTTGGGCGTTCACTGCGTTAGGCGTTCTCGGTAACGATGACACCGCCCGCAAACTGa                                                                                                                                        <  2:218635/148‑1 (MQ=255)
                aaaGACGTCTGTTCCGCCGACTCACTGGCGGGATTTGCCTGGGATCTGTTTACCGCCTGGCAGACCGCTGGCGCGCCGTCGAAAGAGAGTTGGGCGTTCACTGCGTTAGGCGTTCTCGGTAACGATGACACCGCCCGCAAACTGACGc                                                                                                                                     <  1:1460414/148‑1 (MQ=255)
                aaaGACGTCTGTTCCGCCGACTCACTGGCGGGATTTGCCTGGGATCTGTTTACCGCCTGGCAGACCGCTGGCGCGCCGTCGAAAGAGAGTTGGGCGTTCACTGCGTTAGGCGTTCTCGGTAACGATGACACCGCCCGCAAACTGACGc                                                                                                                                     <  1:2943509/148‑1 (MQ=255)
                  aGACGTCTGTTCCGCCGACTCACTGGCGGGATTTGCCTGGGATCTGTTTACCGCCTGGCAGACCGCTGGCGCGCCGTCGAAAGAGAGTTGGGCGTTCACTGCGTTAGGCGTTCTCGGTAACGATGACACCGCCCGCAAACTGACGCCa                                                                                                                                   <  1:1661886/148‑1 (MQ=255)
                     cGTCTGTTCCGCCGACTCACTGGCGGGATTTGCCTGGGATCTGTTTACCGCCTGGCAGACCGCTGGCGCGCCGTCGAAAGAGAGTTGGGCGTTCACTGCGTTAGGCGTTCTCGGTAACGATGACACCGCCCGCAAACTGACGCCATTa                                                                                                                                <  1:1872662/148‑1 (MQ=255)
                     cGTCTGTTCCGCCGACTCACTGGCGGGATTTGCCTGGGATCTGTTTACCGCCTGGCAGACCGCTGGCGCGCCGTCGAAAGAGAGTTGGGCGTTCACTGCGTTAGGCGTTCTCGGTAACGATGACACCGCCCGCAAACTGACGCCATTa                                                                                                                                <  1:2450596/148‑1 (MQ=255)
                      gTCTGTTCCGCCGACTCACTGGCGGGATTTGCCTGGGATCTGTTTACCGCCTGGCAGACCGCTGGCGCGCCGTCGAAAGAGAGTTGGGCGTTCACTGCGTTAGGCGTTCTCGGTAACGATGACACCGCCCGCAAACTGACGCCATTaa                                                                                                                               <  1:3383876/148‑1 (MQ=255)
                       tCTGTTCCGCCGACTCACTGGCGGGATTTGCCTGGGATCTGTTTACCGCCTGGCAGACCGCTGGCGCGCCGTCGAAAGAGAGTTGGGCGTTCACTGCGTTAGGCGTTCTCGGTAACGATGACACCGCCCGCAAACTGACGCCATTAAt                                                                                                                              >  2:1559817/1‑148 (MQ=255)
                       tCTGTTCCGCCGACTCACTGGCGGGATTTGCCTGGGATCTGTTTACCGCCTGGCAGACCGCTGGCGCGCCGTCGAAAGAGAGTTGGGCGTTCACTGCGTTAGGCGTTCTCGGTAACGATGACACAGCCCGCAAACTGACGCCATTAAt                                                                                                                              >  2:499237/1‑148 (MQ=255)
                        cTGTTCCGCCGACTCACTGGCGGGATTTGCCTGGGATCTGTTTACCGCCTGGCAGACCGCTGGCGCGCCGTCGAAAGAGAGTTGGGCGTTCACTGCGTTAGGCGTTCTCGGTAACGATGACACCGCCCGCAAACTGACGCCATTAATa                                                                                                                             <  2:2028100/148‑1 (MQ=255)
                        cTGTTCCGCCGACTCACTGGCGGGATTTGCCTGGGATCTGTTTACCGCCTGGCAGACCGCTGGCGCGCCGTCGAAAGAGAGTTGGGCGTTCACTGCGTTAGGCGTTCTCGGTAACGATGACACCGCCCGCAAACTGACGCCATTAATa                                                                                                                             <  1:181843/148‑1 (MQ=255)
                        cTGTTCCGCCGACTCACTGGCGGGATTTGCCTGGGATCTGTTTACCGCCTGGCAGACCGCTGGCGCGCCGTCGAAAGAGAGTTGGGCGTTCACTGCGTTAGGCGTTCTCGGTAACGATGACACCGCCCGCAAACTGACGCCATTAATa                                                                                                                             <  1:1507869/148‑1 (MQ=255)
                          gTTCCGCCGACTCACTGGCGGGATTTGCCTGGGATCTGTTTACCGCCTGGCAGACCGCTGGCGCGCCGTCGAAAGAGAGTTGGGCGTTCACTGCGTTAGGCGTTCTCGGTAACGATGACACCGCCCGCAAACTGACGCCATTAATAcg                                                                                                                           <  1:2166732/148‑1 (MQ=255)
                              cgccgACTCACTGGCGGGATTTGCCTGGGATCTGTTTACCGCCTGGCAGACCGCTGGCGCGCCGTCGAAAGAGAGTTGGGCGTTCACTGCGTTAGGCGTTCTCGGTAACGATGACACCGCCCGCAAACTGACGCCATTAATACGCgcc                                                                                                                       <  1:2758646/148‑1 (MQ=255)
                              cgccgACTCACTGGCGGGATTTGCCTGGGATCTGTTTACCGCCTGGCAGACCGCTGGCGCGCCGTCGAAAGAGAGTTGGGCGTTCACTGCGTTAGGCGTTCTCGGTAACGATGACACCGCCCGCAAACTGACGCCATTAATACGCgcc                                                                                                                       <  1:1704646/148‑1 (MQ=255)
                                   aCTCACTGGCGGGATTTGCCTGGGATCTGTTTACCGCCTGGCAGACCGCTGGCGCGCCGTCGAAAGAGAGTTGGGCGTTCACTGCGTTAGGCGTTCTCGGTAACGATGACACCGCCCGCAAACTGACGCCATTAATACGCgcctggcc                                                                                                                  <  2:439969/148‑1 (MQ=255)
                                   aCTCACTGGCGGGATTTGCCTGGGATCTGTTTACCGCCTGGCAGACCGCTGGCGCGCCGTCGAAAGAGAGTTGGGCGTTCACTGCGTTAGGCGTTCTCGGTAACGATGACACCGCCCGCAAACTGACGCCATTAATACGCgcctggcc                                                                                                                  >  1:2607416/1‑148 (MQ=255)
                                     tCACTGGCTGGATTTGCCTGGGATCTGTTTACCGCCTGGCAGACCGCTGGCGCGCCGTCGAAAGAGAGTTGGGCGTTCACTGCGTTAGGCGTTCTCGGTAACGATGACACCGCCCGCAAACTGACGCCATTAATACGCgcctggcctg                                                                                                                <  2:1467415/148‑1 (MQ=255)
                                     tCACTGGCGGGATTTGCCTGGGATCTGTTTACCGCCTGGCAGACCGCTGGCGCGCCGTCGAAAGAGAGTTGGGCGTTCACTGCGTTAGGCGTTCTCGGTAACGATGACACCGCCCGCAAACTGACGCCATTAATACGCgcctggcctg                                                                                                                >  1:2936816/1‑148 (MQ=255)
                                      cACTGGCGGGATTTGCCTGGGATCTGTTTACCGCCTGGCAGACCGCTGGCGCGCCGTCGAAAGAGAGTTGGGCGTTCACTGCGTTAGGCGTTCTCGGTAACGATGACACCGCCCGCAAACTGACGCCATTAATACGCGcctggcctgg                                                                                                               >  2:1097524/1‑148 (MQ=255)
                                      cACTGGCGGGATTTGCCTGGGATCTGTTTACCGCCTGGCAGACCGCTGGCGCGCCGTCGAAAGAGAGTTGGGCGTTCACTGCGTTAGGCGTTCTCGGTAACGATGACACCGCCCGCAAACTGACGCCATTAATACGCGcctggcctgg                                                                                                               <  2:3329340/148‑1 (MQ=255)
                                         tGGCGGGATTTGCCTGGGATCTGTTTACCGCCTGGCAGACCGCTGGCGCGCCGTCGAAAGAGAGTTGGGCGTTCACTGCGTTAGGCGTTCTCGGTAACGATGACACCGCCCGCAAACTGACGCCATTAATACGCGCCTGGCCTGGTGa                                                                                                            >  2:492553/1‑148 (MQ=255)
                                           gCGGGATTTGCCTGGGATCTGTTTACCGCCTGGCAGACCGCTGGCGCGCCGTCGAAAGAGAGTTGGGCGTTCACTGCGTTAGGCGTTCTCGGTAACGATGACACCGCCCGCAAACTGACGCCATTAATACGCGCCTGGCCTGGTGAAt                                                                                                          <  2:310741/148‑1 (MQ=255)
                                           gCGGGATTTGCCTGGGATCTGTTTACCGCCTGGCAGACCGCTGGCGCGCCGTCGAAAGAGAGTTGGGCGTTCACTGCGTTAGGCGTTCTCGGTAACGATGACACCGCCCGCAAACTGACGCCATTAATACGCGCCTGGCCTGGTGAAt                                                                                                          <  1:3496054/148‑1 (MQ=255)
                                            cGGGATTTGCCTGGGATCTGTTTACCGCCTGGCAGACCGCTGGCGCGCCGTCGAAAGAGAGTTGGGCGTTCACTGCGTTAGGCGTTCTCGGTAACGATGACACCGCCCGCAAACTGACGCCATTAATACGCGCCTGGCCTGGTGAATc                                                                                                         <  1:1052337/148‑1 (MQ=255)
                                            cGGGATTTGCCTGGGATCTGTTTACCGCCTGGCAGACCGCTGGCGCGCCGTCGAAAGAGAGTTGGGCGTTCACTGCGTTAGGCGTTCTCGGTAACGATGACACCGCCCGCAAACTGACGCCATTAATACGCGCCTGGCCTGGTGAATc                                                                                                         <  2:2423520/148‑1 (MQ=255)
                                              ggATTTGCCTGGGATCTGTTTACCGACTGGCAGACCGCTGGCGCGCCGTCGAAAGAGAGTTGGGCGTTCACTGCGTTAGGCGTTCTCGGTAACGATGACACCGCCCGCAAACTGACGCCATTAATACGCGCCTGGCCTGGTGAATccc                                                                                                       >  2:1634561/1‑148 (MQ=255)
                                               gATTTGCCTGGGATCTGTTTACCGCCTGGCAGACCGCTGGCGCGCCGTCGAAAGAGAGTTGGGCGTTCACTGCGTTAGGCGTTCTCGGTAACGATGACACCGCCCGCAAACTGACGCCATTAATACGCGCCTGGCCTGGTGAATCCCa                                                                                                      >  1:2837396/1‑148 (MQ=255)
                                                  ttGCCTGGGATCTGTTTACCGCCTGGCAGACCGCTGGCGCGCCGTCGAAAGAGAGTTGGGCGTTCACTGCGTTAGGCGTTCTCGGTAACGATGACACCGCCCGCAAACTGACGCCATTAATACGCGCCTGGCCTGGTGAATCCCAGCa                                                                                                   >  2:2131117/1‑148 (MQ=255)
                                                  ttGCCTGGGATCTGTTTACCGCCTGGCAGACCGCTGGCGCGCCGTCGAAAGAGAGTTGGGCGTTCACTGCGTTAGGCGTTCTCGGTAACGATGACACCGCCCGCAAACTGACGCCATTAATACGCGCCTGGCCTGGTGAATCCCAGCa                                                                                                   >  1:2125731/1‑148 (MQ=255)
                                                     ccTGGGATCTGTTTACCGCCTGGCAGACCGCTGGCGCGCCGTCGAAAGAGAGTTGGGCGTTCACTGCGTTAGGCGTTCTCGGTAACGATGACACCGCCCGCAAACTGACGCCATTAATACGCGCCTGGCCTGGTGAATCCCAGCATaa                                                                                                <  2:220158/148‑1 (MQ=255)
                                                     ccTGGGATCTGTTTACCGCCTGGCAGACCGCTGGCGCGCCGTCGAAAGAGAGTTGGGCGTTCACTGCGTTAGGCGTTCTCGGTAACGATGACACCGCCCGCAAACTGACGCCATTAATACGCGCCTGGCCTGGTGAATCCCAGCATaa                                                                                                <  2:803667/148‑1 (MQ=255)
                                                     ccTGGGATCTGTTTACCGCCTGGCAGACCGCTGGCGCGCCGTCGAAAGAGAGTTGGGCGTTCACTGCGTTAGGCGTTCTCGGTAACGATGACACAGCCCGCAAACTGACGCCATTAATACGCGCCTGGCCTGGTGAATCCCAGCATaa                                                                                                <  1:499237/148‑1 (MQ=255)
                                                          gATCTGTTTACCGCCTGGCAGACCGCTGGCGCGCCGTCGAAAGAGAGTTGGGCGTTCACTGCGTTAGGCGTTCTCGGTAACGATGACACCGCCCGCAAACTGACGCCATTAATACGCGCCTGGCCTGGTGAATCCCAGCATAAAcgcg                                                                                           <  1:1341295/148‑1 (MQ=255)
                                                           aTCTGTTTACCGCCTGGCAGACCGCTGGCGCGCCGTCGAAAGAGAGTTGGGCGTTCACTGCGTTAGGCGTTCTCGGTAACGATGACACCGCCCGCAAACTGACGCCATTAATACGCGCCTGGCCTGGTGAATCCCAGCATAAACgcgc                                                                                          >  1:2634764/1‑148 (MQ=255)
                                                           aTCTGTTTACCGCCTGGCAGACCGCTGGCGCGCCGTCGAAAGAGAGTTGGGCGTTCACTGCGTTAGGCGTTCTCGGTAACGATGACACCGCCCGCAAACTGACGCCATTAATACGCGCCTGGCCTGGTGAATCCCAGCATAAACgcgc                                                                                          >  1:799112/1‑148 (MQ=255)
                                                            tCTGTTTACCGCCTGGCAGACCGCTGGCGCGCCGTCGAAAGAGAGTTGGGCGTTCACTGCGTTAGGCGTTCTCGGTAACGATGACACCGCCCGCAAACTGACGCCATTAATACGCGCCTGGCCTGGTGAATCCCAGCATAAACGCGcc                                                                                         >  1:1962761/1‑148 (MQ=255)
                                                             cTGTTCACCGCCTGGCAGACCGCTGGCGCGCCGTCGAAAGAGAGTTGGGCGTTCACTGCGTTAGGCGTTCTCGGTAACGATGACACCGCCCGCAAACTGACGCCATTAATACGCGCCTGGCCTGGTGAATCCCAGCATAAACGCGCCa                                                                                        <  1:858550/148‑1 (MQ=255)
                                                                    ccGCCTGGCAGACCGCTGGCGCGCCGTCGAAAGAGAGTTGGGCGTTCACTGCGTTAGGCGTTCTCGGTAACGATGACACAGCCCGCAAACTGACGCCATTAATACGCGCCTGGCCTGGTGAATCCCAGCATAAACGCGCCACCgttgg                                                                                 <  2:1880393/148‑1 (MQ=255)
                                                                      gCCTGGCAGTCCGCTGGCGCGCCGTCGAAAGAGATTTGGGCGTTCACTGCGTTAGGCGTTCTCGGTAACGATGACACCGCCCGCAAACTGACGCCATTAATACGCGCCTGGCCTTGTGAATCCCAGCATAAACGCGCCACCgttgggt                                                                               <  2:1592431/148‑1 (MQ=255)
                                                                      gCCTGGCAGACCGCTGGCGCGCCGTCGAAAGAGAGTTGGGCGTTCACTGCGTTAGGCGTTCTCGGTAACGATGACACCGCCCGCAAACTGACGCCATTAATACGCGCCTGGCCTGGTGAATCCCAGCATAAACGCGCCACCgttgggt                                                                               <  2:1049339/148‑1 (MQ=255)
                                                                      gCCTGGCAGACCGCTGGCGCGCCGTCGAAAGAGAGTTGGGAGTTCACTGCGTTAGGCGTTCTCGGTAACGATGACACCGCCCGCAAACTGACGCCATTAATACGCGCCTGGCCTGGTGAATCCCAGCATAAACGCGCCACCgttgggt                                                                               <  2:3547686/148‑1 (MQ=255)
                                                                         tGGCAGACCGCTGGCGCGCCGTCGAAAGAGAGTTGGGCGTTCACTGCGTTAGGCGTTCTCGGTAACGATGACACCGCCCGCAAACTGACGCCATTAATACGCGCCTGGCCTGGTGAATCCCAGCATAAACGCGCCACCgttgggttgg                                                                            <  2:1909645/148‑1 (MQ=255)
                                                                          ggCAGACCGCTGGCGCGCCGTCGAAAGAGAGTTGGGCGTTCACTGCGTTAGGCGTTCTCGGTAACGATGACACCGCCCGCAAACTGACGCCATTAATACGCGCCTGGCCTGGTGAATCCCAGCATAAACGCGCCACCGTTGGGTTGGa                                                                           >  1:534077/1‑148 (MQ=255)
                                                                           gCAGACCGCTGGCGCGCCGTCGAAAGAGAGTTGGGCGTTCACTGCGTTAGGCGTTCTCGGTAACGATGACACCGCCCGCAAACTGACGCCATTAATACGCGCCTGGCCTGGTGAATCCCAGCATAAACGCGCCACCGTTGGGTTGGat                                                                          <  2:1357889/148‑1 (MQ=255)
                                                                             aGACCGCTGGCGCGCCGTCGAAAGAGAGTTGGGCGTTCACTGCGTTAGGCGTTCTCGGTAACGATGACACCGCCCGCAAACTGACGCCATTAATACGCGCCTGGCCTGGTGAATCCCAGCATAAACGCGCCACCGTTGGGTTGGatat                                                                        <  1:1981495/148‑1 (MQ=255)
                                                                               aCCGCTGGCGCGCCGTCGAAAGAGAGTTGGGCGTTCACTGCGTTAGGCGTTCTCGGTAACGATGACACCGCCCGCAAACTGACGCCATTAATACGCGCCTGGCCTGGTGAATCCCAGCATAAACGCGCCACCGTTGGGTTGGATATtc                                                                      <  2:886734/148‑1 (MQ=255)
                                                                               aCCGCTGGCGCGCCGTCGAAAGAGAGTTGGGCGTTCACTGCGTTAGGCGTTCTCGGTAACGATGACACCGCCCGCAAACTGACGCCATTAATACGCGCCTGGCCTGGTGAATCCCAGCATAAACGCGCCACCGTTGGGTTGGATATtc                                                                      <  2:2607416/148‑1 (MQ=255)
                                                                                ccGCTGGCGCGCCGTCGAAAGAGAGTTGGGCGTTCACTGCGTTAGGCGTTCTCGGTAACGATGACACCGCCCGCAAACTGACGCCATTAATACGCGCCTGGCCTGGTGAATCCCAGCATAAACGCGCCACCGTTGGGTTGGATATtct                                                                     >  2:289535/1‑148 (MQ=255)
                                                                                ccGCTGGCGCGCCGTCGAAAGAGAGTTGGGCGTTCACTGCGTTAGGCGTTCTCGGTAACGATGACACCGCCCGCAAACTGACGCCATTAATACGCGCCTGGCCTGGTGAATCCCAGCATAAACGCGCCACCGTTGGGTTGGATATtct                                                                     <  1:112203/148‑1 (MQ=255)
                                                                                   cTGGCGCGCCGTCGAAAGAGAGTTGGGCGTTCACTGCGTTAGGCGTTCTCGGTAACGATGACACCGCCCGCAAACTGACGCCATTAATACGCGCCTGGCCTGGTGAATCCCAGCATAAACGCGCCACCGTTGGGTTGGATATTCTCgc                                                                  <  2:2978012/148‑1 (MQ=255)
                                                                                     ggCGCGCCGTCGAAAGAGAGTTGGGCGTTCACTGCGTTAGGCGTTCTCGGTAACGATGACACCGCCCGCAAACTGACGCCATTAATACGCGCCTGGCCTGGTGAATCCCAGCATAAACGCGCCACCGTTGGGTTGGATATTCTCgctg                                                                <  1:37529/148‑1 (MQ=255)
                                                                                        gcgcCGTCGAAAGAGAGTTGGGCGTTCACTGCGTTAGGCGTTCTCGGTAACGATGACACCGCCCGCAAACTGACGCCATTAATACGCGCCTGGCCTGGTGAATCCCAGCATAAACGCGCCACCGTTGGGTTGGATATTCTCGCTGCTa                                                             <  2:3225486/148‑1 (MQ=255)
                                                                                          gcCGTCGAAAGAGAGTTGGGCGTTCACTGCGTTAGGCGTTCTCGGTAACGATGACACCGCCCGCAAACTGACGCCATTAATACGCGCCTGGCCTGGTGAATCCCAGCATAAACGCGCCACCGTTGGGTTGGATATTCTCGCTGCTATc                                                           <  2:3500549/148‑1 (MQ=255)
                                                                                             gTCGAAAGAGAGTTGGGCGTTCACTGCGTTAGGCGTTCTCGGTAACGATGACACCGCCCGCAAACTGACGCCATTAATACGCGCCTGGCCTGGTGAATCCCAGCATAAACGCGCCACCGTTGGGTTGGATATTCTCGCTGCTATCGGt                                                        <  1:3110402/148‑1 (MQ=255)
                                                                                              tCGAAAGAGAGTTGGGCGTTCACTGCGTTAGGCGTTCTCGGTAACGATGACACCGCCCGCAAACTGACGCCATTAATACGCGCCTGGCCTGGTGAATCCCAGCATAAACGCGCCACCGTTGGGTTGGATATTCTCGCTGCTATCGGTa                                                       <  2:401829/148‑1 (MQ=255)
                                                                                              tCGAAAGAGAGTTGGGCGTTCACTGCGTTAGGCGTTCTCGGTAACGATGACACCGCCCGCAAACTGACGCCATTAATACGCGCCTGGCCTGGTGAATCCCAGCATAAACGCGCCACCGTTGGGTTGGATATTCTCGCTGCTATCGGTa                                                       <  2:3126587/148‑1 (MQ=255)
                                                                                              tCGAAAGAGAGTTGGGCGTTCACTGCGTTAGGCGTTCTCGGTAACGATGACACCGCCCGCAAACTGACGCCATTAATACGCGCCTGGCCTGGTGAATCCCAGCATAAACGCGCCACCGTTGGGTTGGATATTCTCGCTGCTATCGGTa                                                       <  1:2501468/148‑1 (MQ=255)
                                                                                               cGAAAGAGAGTTGGGCGTTCACTGCGTTAGGCGTTCTCGGTAACGATGACACCGCCCGCAAACTGACGCCATTAATACGCGCCTGGCCTGGTGAATCCCAGCATAAACGCGCCACCGTTGGGTTGGATATTCTCGCTGCTATCGGTAg                                                      >  2:1190098/1‑148 (MQ=255)
                                                                                               cGAAAGAGAGTTGGGCGTTCACTGCGTTAGGCGTTCTCGGTAACGATGACACCGCCCGCAAACTGACGCCATTAATACGCGCCTGGCCTGGTGAATCCCAGCATAAACGCGCCACCGTTGGGTTGGATATTCTCGCTGCTATCGGTAg                                                      <  2:650092/148‑1 (MQ=255)
                                                                                                gAAAGAGAGTTGGGCGTTCACTGCGTTAGGCGTTCTCGGTAACGATGACACCGCCCGCAAACTGACGCCATTAATACGCGCCTGGCCTGGTGAATCCCAGCATAAACGCGCCACCGTTGGGTTGGATATTCTCGCTGCTATCGGTAGt                                                     <  2:285907/148‑1 (MQ=255)
                                                                                                 aaaGAGAGTTGGGCGTTCACTGCGTTAGGCGTTCTCGGTAACGATGACACCGCCCGCAAACTGACGCCATTAATACGCGCCTGGCCTGGTGAATCCCAGCATAAACGCGCCACCGTTGGGTTGGATATTCTCGCTGCTATCGGTAGTg                                                    <  2:3194296/148‑1 (MQ=255)
                                                                                                  aaGAGAGTTGGGCGTTCACTGCGTTAGGCGTTCTCGGTAACGATGACACCGCCCGCAAACTGACGCCATTACTACGCGCCTGGCCTGGTGAATCCCAGCATAAACGCGCCACCGTTGGGTTGGATATTCTCGCTGCTATCGGTAGTGa                                                   <  2:3273177/148‑1 (MQ=255)
                                                                                                  aaGAGAGTTGGGCGTTCACTGCGTTAGGCGTTCTCGGTAACGATGACACCGCCCGCAAACTGACGCCATTAATACGCGCCTGGCCTGGTGAATCCCAGCATAAACGCGCCACCGTTGGGTTGGATATTCTCGCTGCTATCGGTAGTGa                                                   >  1:3346371/1‑148 (MQ=255)
                                                                                                   agagagTTGGGCGTTCACTGCGTTAGGCGTTCTCGGTAACGATGACACCGCCCGCAAACTGACGCCATTAATACGCGCCTGGCCTGGTGAATCCCAGCATAAACGCGCCACCGTTGGGTTGGATATTCTCGCTGCTATCGGTAGTGat                                                  <  1:1482052/148‑1 (MQ=255)
                                                                                                   agagagTTGGGCGTTCACTGCGTTAGGCGTTCTCGGTAACGATGACACCGCCCGCAAACTGACGCCATTAATACGCGCCTGGCCTGGTGAATCCCAGCATAAACGCGCCACCGTTGGGTTGGATATTCTCGCTGCTATCGGTAGTGat                                                  <  1:2564261/148‑1 (MQ=255)
                                                                                                         ttGGGCGTTCACTGCGTTAGGCGTTCTCGGTAACGATGACACCGCCCGCAAACTGACGCCATTAATACGCGCCTGGCCTGGTGAATCCCAGCATAAACGCGCCACCGTTGGGTTGGATATTCTCGCTGCTATCGGTAGTGATATCGcc                                            >  1:968918/1‑148 (MQ=255)
                                                                                                          tGGGCGTTCACTGCGTTAGGCGTTCTCGGTAACGATGACACCGCCCGCAAACTGACGCCATTAATACGCGCCTGGCCTGGTGAATCCCAGCATAAACGCGCCACCGTTGGGTTGGATATTCTCGCTGCTATCGGTAGTGATATCGccc                                           <  1:3236109/148‑1 (MQ=255)
                                                                                                          tGGGCGTTCACTGCGTTAGGCGTTCTCGGTAACGATGACACCGCCCGCAAACTGACGCCATTAATACGCGCCTGGCCTGGTGAATCCCAGCATAAACGCGCCACCGTTGGGTTGGATATTCTCGCTGCTATCGGTAGTGATATCGccc                                           <  2:234628/148‑1 (MQ=255)
                                                                                                             gCGTTCACTGCGTTAGGCGTTCTCGGTAACGATGACACCGCCCGCAAACTGACGCCATTAATACGCGCCTGGCCTGGTGAATCCCAGCATAAACGCGCCACCGTTGGGTTGGATATTCTCGCTGCTATCGGTAGTGATATCGCCCTTa                                        <  1:3105882/148‑1 (MQ=255)
                                                                                                             gCGTTCACTGCGTTAGGCGTTCTCGGTAACGATGACACCGCCCGCAAACTGACGCCATTAATACGCGCCTGGCCTGGTGAATCCCAGCATAAACGCGCCACCGTTGGGTTGGATATTCTCGCTGCTATCGGTAGTGATATCGCCCTTa                                        <  1:307136/148‑1 (MQ=255)
                                                                                                               gTTCAATGCGTTAGGCGTTCTCGGTAACGATGACACCGCCCGCAAACTGACGCCATTAATACGCGCCTGGCCTGGTGAATCCCAGCATAAACGCGCCACCGTTGGGTTGGATATTCTCGCTGCTATCGGTAGTGATATCGCCCTTATg                                      <  1:924171/148‑1 (MQ=255)
                                                                                                                 tCACTGCGTTAGGCGTTCTCGGTAACGATGACACCGCCCGCAAACTGACGCCATTAATACGCGCCTGGCCTGGTGAATCCCAGCATAAACGCGCCACCGTTGGGTTGGATATTCTCGCTGCTATCGGTAGTGATATCGCCCTTATGCa                                    >  2:2445864/1‑148 (MQ=255)
                                                                                                                  cACTGCGTTAGGCGTTCTCGGTAACGATGACACCGCCCGCAAACTGACGCCATTAATACGCGCCTGGCCTGGTGAATCCCAGCATAAACGCGCCACCGTTGGGTTGGATATTCTCGCTGCTATCGGTAGTGATATCGCCCTTATGCAg                                   <  2:1601585/148‑1 (MQ=255)
                                                                                                                   aCTGCGTTAGGCGTTCTCGGTAACGATGACACCGCCCGCAAACTGACGCCATTAATACGCGCCTGGCCTGGTGAATCCCAGCATAAACGCGCCACCGTTGGGTTGGATATTCTCGCTGCTATCGGTAGTGATATCGCCCTTATGCAGc                                  >  1:1399134/1‑148 (MQ=255)
                                                                                                                   aCTGCGTTAGGCGTTCTCGGTAACGATGACACCGCCCGCAAACTGACGCCATTAATACGCGCCTGGCCTGGTGAATCCCAGCATAAACGCGCCACCGTTGGGTTGGATATTCTCGCTGCTATCGGTAGTGATATCGCCCTTATGCAGc                                  >  1:3337390/1‑148 (MQ=255)
                                                                                                                           aGGCGTTCTCGGTAACGATGACACAGCCCGCAAACTGACGCCATTAATACGCGCCTGGCCTGGTGAATCCCAGCATAAACGCGCCACCGTTGGGTTGGATATTCTCGCTGCTATCGGTAGTGATATCGCCCTTATGCAGCTTAACGGc                          >  2:3251567/1‑148 (MQ=255)
                                                                                                                            ggCGTTCTCGGTAACGATGACACCGCCCGCAAACTGACGCCATTAATACGCGCCTGGCCTGGTGAATCCCAGCATAAACGCGCCACCGTTGGGTTGGATATTCTCGCTGCTATCGGTAGTGATATCGCCCTTATGCAGCTTAACGGCa                         >  2:1512339/1‑148 (MQ=255)
                                                                                                                            ggCGTTCTCGGTAACGATGACACCGCCCGCAAACTGACGCCATTAATACGCGCCTGGCCTGGTGAATCCCAGCATAAACGCGCCACCGTTGGGTTGGATATTCTCGCTGCTATCGGTAGTGATATCGCCCTTATGCAGCTTAACGGCa                         >  2:3417699/1‑148 (MQ=255)
                                                                                                                              cGTTCTCGGTAACGATGACACCGCCCGCAAACTGACGCCATTAATACGCGCCTGGCCTGGTGAATCCCAGCATAAACGCGCCACCGTTGGGTTGGATATTCTCGCTGCTATCGGTAGTGATATCGCCCTTATGCAGCTTAACGGCATc                       >  2:85336/1‑148 (MQ=255)
                                                                                                                                 tctcGGTAACGATGACACCGCCCGCAAACTGACGCCATTAATACGCGCCTGGCCTGGTGAATCCCAGCATAAACGCGCCACCGTTGGGTTGGATATTCTCGCTGCTATCGGTAGTGATATCGCCCTTATGCAGCTTAACGGCATCGcc                    <  1:1167849/148‑1 (MQ=255)
                                                                                                                                  ctcGGTAACGATGACACCGCCCGCAAACTGACGCCATTAATACGCGCCTGGCCTGGTGAATCCCAGCATAAACGCGCCACCGTTGGGTTGGATATTCTCGCTGCTATCGGTAGTGATATCGCCCTTATGCAGCTTAACGGCATCGccc                   <  2:2182639/148‑1 (MQ=255)
                                                                                                                                       tAACGATGACACCGCCCGCAAACTGACGCCATTAATACGCGCCTGGCCTGGTGAATCCCAGCATAAACGCGCCACCGTTGGGTTGGATATTCTCGCTGCTATCGGTAGTGATATCGCCCTTATGCAGCTTAACGGCATCGCCCAGaaa              <  1:1324054/148‑1 (MQ=255)
                                                                                                                                         aCGATGACACCGCCCGCAAACTGACGCCATTAATACGCGCCTGGCCTGGTGAATCCCAGCATAAACGCGCCACCGTTGGGTTGGATATTCTCGCTGCTATCGGTAGTGATATCGCCCTTATGCAGCTTAACGGCATCGCCCAGAAACt            >  2:2769/1‑148 (MQ=255)
                                                                                                                                         aCGATGACACCGCCCGCAAACTGACGCCATTAATACGCGCCTGGCCTGGTGAATCCCAGCATAAACGCGCCACCGTTGGGTTGGATATTCTCGCTGCTATCGGTAGTGATATCGCCCTTATGCAGCTTAACGGCATCGCCCAGAAACt            >  2:56739/1‑148 (MQ=255)
                                                                                                                                            aTGACACCGCCCGCAAACTGACGCCATTAATACGCGCCTGGCCTGGTGAATCCCAGCATAAACGCGCCACCGTTGGGTTGGATATTCTCGCTGCTATCGGTAGTGATATCGCCCTTATGCAGCTTAACGGCATCGCCCAGAAACTGaa         <  1:6970/148‑1 (MQ=255)
                                                                                                                                                  ccgcccgcAAACTGACGCCATTAATACGCGCCTGGCCTGGTGAATCCCAGCATAAACGCGCCACCGTTGGGTTGGATATTCTCGCTGCTATCGGTAGTGATATCGCCCTTATGCAGCTTAACGGCATCGCCCAGAAACTGAAATTCaa   <  1:2908668/148‑1 (MQ=255)
                                                                                                                                                  ccgcccgcAAACTGACGCCATTAATACGCGCCTGGCCTGGTGAATCCCAGCATAAACGCGCCACCGTTGGGTTGGATATTCTCGCTGCTATCGGTAGTGATATCGCCCTTATGCAGCTTAACGGCATCGCCCAGAAACTGAAATTCaa   >  2:3476824/1‑148 (MQ=255)
                                                                                                                                                  ccgcccgcAAACTGACGCCATTAATACGCGCCTGGCCTGGTGAATCCCAGCATAAACGCGCCACCGTTGGGTTGGATATTCTCGCTGCTATCGGTAGTGATATCGCCCTTATGCAGCTTAACGGCATCGCCCAGAAACTGAAATTCaa   >  1:3226608/1‑148 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
GGGATTATTGCAGGTGAAAGACGTCTGTTCCGCCGACTCACTGGCGGGATTTGCCTGGGATCTGTTTACCGCCTGGCAGACCGCTGGCGCGCCGTCGAAAGAGAGTTGGGCGTTCACTGCGTTAGGCGTTCTCGGTAACGATGACACCGCCCGCAAACTGACGCCATTAATACGCGCCTGGCCTGGTGAATCCCAGCATAAACGCGCCACCGTTGGGTTGGATATTCTCGCTGCTATCGGTAGTGATATCGCCCTTATGCAGCTTAACGGCATCGCCCAGAAACTGAAATTCAAA  >  CP009273/2192424‑2192718

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