Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,191,496 C→A 9.6% pseudogene (1544/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,191,4960CA 9.6% 733.7 / 37.7 225pseudogene (1544/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 (112/91);  new base A (11/11);  total (123/102)
Fisher's exact test for biased strand distribution p-value = 6.59e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.26e-01

Only 100 of 227 total aligned reads displayed.

AACCAGCACTGGCAGAACAGGTCATTCCCTGGCTCTCGACACCCGCAGTTGCCGTACTGAAATCATGCCAGCAACAACTGACACAGCCCTCAAACCATGCCAGCGCCGATCTACTGCCAGCCGTAGTAGTCTCCCCTCCCTGGCTTTCGAAAAAGAAAAAATCGCCGATTCCGGTGCTGGATTTAGCGCCATTAGGCATTGAGCCAATCTGTTATCTGACAGAAGAAATCAGTAATCAACTTTTGGCGAAATATATCTGGTATTCAAAACACATCACGGTTAGCCATGAAGAA  >  CP009273/2191349‑2191641
                                                                                                                                                   |                                                                                                                                                 
aaCCAGCACTGGCAGAACAGGTCATTCCCTGGCTCTCGACACCCGCAGTTGCCGTACTGAAATCATGCCAGCAACAACTGACACAGCCCTCAAACCATGCCAGCGCCGATCTACTGCCAGCCGTAGTAGTCTCCCCTCCCTGGCTTTc                                                                                                                                                   >  2:780448/1‑148 (MQ=255)
 aCCAGCACTGGCAGAACAGGTCATTCCCTGGCTCTCGACACCCGCAGTTGCCGTACTGAAATCATGCCAGCAACAACTGACACAGCCCTCAAACCATGCCAGCGCCGATCTACTGCCAGCCGTAGTAGTCTCCCCTCCCTGGCTTTCg                                                                                                                                                  >  1:1068661/1‑148 (MQ=255)
  ccAGCACTGGCAGAACAGGTCATTCCCTGGCTCTCGACACCCGCAGTTGCCGTACTGAAATCATGCCAGCAACAACTGACACAGCCCTCAAACCATGCCAGCGCCGATCTACTGCCAGCCGTAGTAGTCTCCCCTCCCTGGCTTTCga                                                                                                                                                 >  2:2228007/1‑148 (MQ=255)
  ccAGCACTGGCAGAACAGGTCATTCCCTGGCTCTCGACACCCGCAGTTGCCGTACTGAAATCATGCCAGCAACAACTGACACAGCCCTCAAACCATGCCAGCGCCGATCTACTGCCAGCCGTAGTAGTCTCCCCTCCCTGGCTTTCga                                                                                                                                                 >  2:2427953/1‑148 (MQ=255)
      cACTGGCAGAACAGGTCATTCCCTGGCTCTCGACACCCGCAGTTGCCGTACTGAAATCATGCCAGCAACAACTGACACAGCCCTCAAACCATGCCAGCGCCGATCTACTGCCAGCCGTAGTAGTCTCCCCTCCCTGGCTTTCgaaaaa                                                                                                                                             >  1:3282791/1‑148 (MQ=255)
      cACTGGCAGAACAGGTCATTCCCTGGCTCTCGACACCCGCAGTTGCCGTACTGAAATCATGCCAGCAACAACTGACACAGCCCTCAAACCATGCCAGCGCCGATCTACTGCCAGCCGTAGTAGTCTCCCCTCCCTGGCTTTCgaaaaa                                                                                                                                             <  1:95267/148‑1 (MQ=255)
       aCTGGCAGAACAGGTCATTCCCTGGCTCTCGACACCCGCAGTTGCCGTACTGAAATCATGCCAGCAACAACTGACACAGCCCTCAAACCATGCCAGCGCCGATCTACTGCCAGCCGTAGTAGTCTCCCCTCCCTGGCTTTCgaaaaag                                                                                                                                            <  2:658599/148‑1 (MQ=255)
        cTGGCAGAACAGGTCATTCCCTGGCTCTCGACACCCGCAGTTGCCGTACTGAAATCATGCCAGCAACAACTGACACAGCCCTCAAACCATGCCAGCGCCGATCTACTGCCAGCCGTAGTAGTCTCCCCTCCCTGGCTTTAgaaaaaga                                                                                                                                           >  1:2142458/1‑148 (MQ=255)
        cTGGCAGAACAGGTCATTCCCTGGCTCTCGACACCCGCAGTTGCCGTACTGAAATCATGCCAGCAACAACTGACACAGCCCTCAAACCATGCCAGCGCCGATCTACTGCCAGCCGTAGTAGTCTCCCCTCCCTGGCTTTAgaaaaaga                                                                                                                                           >  2:1078674/1‑148 (MQ=255)
         tGGCAGAACAGGTCATTCCCTGGCTCTCGACACCCGCAGTTGCCGTACTGAAATCATGCCAGCAACAACTGACACAGCCCTCAAACCATGCCAGCGCCGATCTACTGCCAGCCGTAGTAGTCTCCCCTCCCTGGCTTTCgaaaaagaa                                                                                                                                          >  2:279961/1‑148 (MQ=255)
           gCAGAACAGGTCATTCCCTGGCTCTCGACACCCGCAGTTGCCGTACTGAAATCATGCCAGCAACAACTGACACAGCCCTCAAACCATGCCAGCGCCGATCTACTGCCAGCCGTAGTAGTCTCCCCTCCCTGGCTTTagaaaaagaaaa                                                                                                                                        <  2:3198696/148‑1 (MQ=255)
            cataaCAGGTCATTCCCTGGCTCTCGACACCCGCAGTTGCCGTACTGAAATCATGCCAGCAACAACTGACACAGCCCTCAAACCATGCCAGCGCCGATCTACTGCCAGCCGTAGTAGTCTCCCCTCCCTGGCTTTCgaaaaagaaaaa                                                                                                                                       >  2:2368376/4‑148 (MQ=255)
            cAGAACAGGTCATTCCCTGGCTCTCGACACCCGCAGTTGCCGTACTGAAATCATGCCAGCAACAACTGACACAGCCCTCAAACCCTGCCAGCGCCGATCTACTGCCAGCCGTAGTAGTCTCCCCTCCCTGGCTTTCgaaaaagaaaaa                                                                                                                                       <  2:519752/148‑1 (MQ=255)
              gAACAGGTCATTCCCTGGCTCTCGACACCCGCAGTTGCCGTACTGAAATCATGCCAGCAACAACTGACACAGCCCTCAAACCATGCCAGCGCCGATCTACTGCCAGCCGTAGTAGTCTCCCCTCCCTGGCTTTCGAAAAAGAAAAAAt                                                                                                                                     >  1:2399666/1‑148 (MQ=255)
                aCAGGTCATTCCCTGGCTCTCGACACCCGCAGTTGCCGTACTTAAATCATGCCAGCAACAACTGACACAGCCCTCAAACCATGCCAGCGCCGATCTACTGCCAGCCGTAGTAGTCTCCCCTCCCTGGCTTTCGAAAAAGAAAAAATCg                                                                                                                                   <  1:53261/148‑1 (MQ=255)
                aCAGGTCATTCCCTGGCTCTCGACACCCGCAGTTGCCGTACTGAAATCATGCCAGCAACAACTGACACAGCCCTCAAACCATGCCAGCGCCGATCTACTGCCAGCCGTAGTAGTCTCCCCTCCCTGGCTTTCGAAAAAGAAAAAATCg                                                                                                                                   >  1:2597225/1‑148 (MQ=255)
                        ttCCCTGGCTCTCGACACCCGCAGTTGCCGTACTGAAATCATGCCAGCAACAACTGACACAGCCCTCAAACCATGCCAGCGCCGATCTACTGCCAGCCGTAGTAGTCTCCCCTCCCTGGCTTTCGAAAAAGAAAAAATCGCCGATTcc                                                                                                                           >  2:450597/1‑148 (MQ=255)
                         tCCCTGGCTCTCGACACCCGCAGTTGCCGTACTGAAATCATGCCAGCAACAACTGACACAGCCCTCAAACCATGCCAGCGCCGATCTACTGCCAGCCGTAGTAGTCTCCCCTCCCTGGCTTTCGAAAAAGAAAAAATCGCCGATTCCg                                                                                                                          >  1:2780942/1‑148 (MQ=255)
                         tCCCTGGCTCTCGACACCCGCAGTTGCCGTACTGAAATCATGCCAGCAACAACTGACACAGCCCTCAAACCATGCCAGCGCCGATCTACTGCCAGCCGTAGTAGTCTCCCCTCCCTGGCTTTCGAAAAAGAAAAAATCGCCGATTCCg                                                                                                                          >  1:2289002/1‑148 (MQ=255)
                                ctctCGACACCCGCAGTTGCCGTACTGAAATCATGCCAGCAACAACTGACACAGCCCTCAAACCATGCCAGCGCCGATCTACTGCCAGCCGTAGTAGTCTCCCCTCCCTGGCTTTCGAAAAAGAAAAACTCGCCGAGTCCGGTGCTgg                                                                                                                   >  2:2870352/1‑148 (MQ=255)
                                ctctCGACACCCGCAGTTGCCGTACTGAAATCATGCCAGCAACAACTGACACAGCCCTCAAACCATGCCAGCGCCGATCTACTGCCAGCCGTAGTAGTCTCCCCTCCCTGGCTTTCGAAAAAGAAAAAATCGCCGATTCCGGTGCTgg                                                                                                                   >  1:1249323/1‑148 (MQ=255)
                                ctctCGACACCCGCAGTTGCCGTACTGAAATCATGCCAGCAACAACTGACACAGCCCTCAAACCATGCCAGCGCCGATCTACTGCCAGCCGTAGTAGTCTCCCCTCCCTGGCTTTCGAAAAAGAAAAAATCGCCGATTCCGGTGCTgg                                                                                                                   >  2:3065615/1‑148 (MQ=255)
                                  ctcGACACCCGCAGTTGCCGTACTGAAATCATGCCAGCAACAACTGACACAGCCCAAAAACCATGCCAGCGCCGATCTACTGCCAGCCGTAGTAGTCTCCCCTCCCTGGCTTTCGAAAAAGAAAAAATCGCCGATTCCGGTGCTGGAt                                                                                                                 <  1:1367148/148‑1 (MQ=255)
                                    ccacacCCGCTGTGGCCGTCCTTAAATCATGCCAGCAACAACTCACACACCCCACAACCCCGGCCAGCGCCGATATCCTGCCCGCCGTTTTCTTCCCCCCCCCCGGGCTTTCCAAAAAAAAAAAATCGCCGATTCCGGTGCTGGAttt                                                                                                               <  1:335504/146‑1 (MQ=255)
                                           cGCAGTTGCCGTACTGAAATCATGCCAGCAACAACTGACACAGCCCTCAAACCATGCCAGCGCCGATCTACTGCCAGCCGTAGTAGTCTCCCCTCCCTGGCTTTCGAAAAAGAAAAAATCGCCGATTCCGGTGCTGGATTTAGCGCCa                                                                                                        >  2:1306934/1‑148 (MQ=255)
                                             cAGTTGCCGTACTGAAATCATGCCAGCAACAACTGACACAGCCCTCAAACCATGCCAGCGCCGATCTACTGCCAGCCGTAGTAGTCTCCCCTCCCTGGCTTTCGAAAAAGAAAAAATCGCCGATTCCGGTGCTGGATTTAGCGCCAtt                                                                                                      >  1:3051825/1‑148 (MQ=255)
                                                ttGCCGTACTGAAATCATGCCAGCAACAACTGACACAGCCCTCAAACCATGCCAGCGCCGATCTACTGCCAGCCGTAGTAGTCTCCCCTCCCTGGCTTTCGAAAAAGAAAAAATCGCCGATTCCGGTGCTGGATTTAGCGCCATTAgg                                                                                                   >  2:2500060/1‑148 (MQ=255)
                                                ttGCCGTACTGAAATCATGCCAGCAACAACTGACACAGCCCTCAAACCATGCCAGCGCCGATCTACTGCCAGCCGTAGTAGTCTCCCCTCCCTGGCTTTCGAAAAAGAAAAAATCGCCGATTCCGGTGCGGGATTTAGCGCCATTAgg                                                                                                   >  2:2509028/1‑148 (MQ=255)
                                                 tGCCGTACTTAAATCATGCCACCAAAAACTCACACAGCCCTCAAACCATGCCCGCGCCGATCTACTGCCAGCCGTAGTAGTCTCCCCTCCCTGGCTTTCTAAAAAGAAAAAATCGCCGATTCCGGTGCTGGATTTAGCGCCATTAGGc                                                                                                  <  1:196573/148‑1 (MQ=255)
                                                 tGCCGTACTGAAATCATGCCAGCAACAACTGACACAGCCCTCAAACCATGCCAGCGCCGATCTACTGCCAGCCGTAGTAGTCTCCCCTCCCTGGCTTTCGAAAAAGAAAAAATCGCCGATTCCGGTGCTGGATTTAGCGCCATTAGGc                                                                                                  >  1:2771372/1‑148 (MQ=255)
                                                 ggCCTTACTCAAAACAGGCCAGCAACAACTGACACCGCCCTCAAACCATGCCAGCGCCGATCTACTGCCAGCCGTAGTAGTCTCCCCTCCCTGGCTTTCGAAAAACAAAAAATCGCCGATTCCGGTGCTGGATTTAGCGCCATTAGGc                                                                                                  <  1:2392122/147‑1 (MQ=255)
                                                   ccGTACTGAAATCATGCCAGCAACAACTGACACAGCCCTCAAACCATGCCAGCGCCGATCTACTGCCAGCCGTAGTAGTCTCCCCTCCCTGGCTTTCGAAAAAGAAAAAATCGCCGATTCCGGTGCTGGATTTAGCGCCATTAGGCAt                                                                                                <  2:172166/148‑1 (MQ=255)
                                                      tACTGAAATCATGCCAGCAACAACTGACACAGCCCTCAAACCATGCCAGCGCCGATCTACTGCCAGCCGTAGTAGTCTCCCCTCCCTGGCTTTCGAAAAAGAAAAAATCGCCGATTCCGGTGCTGGATTTAGCGCCATTAGGCATTGa                                                                                             >  2:761573/1‑148 (MQ=255)
                                                       aCTGAAATCATGCCAGCAACAACTGACACAGCCCTCAAACCATGCCCGCGCCGATCTACTGCCAGCCGTAGTAGTCTCCCCTCCCGGGCTTTCTAAAAATAAAAACTCGCCGATTCCGGTGCTGGATTTAGCGCCATTAGGCATTGAg                                                                                            <  1:1055540/148‑1 (MQ=255)
                                                        cTGAAATCATGCCAGCAACAACTGACACAGCCCTCAAACCATGCCAGCGCCGATCTACTGCCAGCCGTAGTAGTCTCCCCTCCCTGGCTTTAGAAAAAGAAAACATCGCCGATTCCGGTGCTGGATTTAGCGCCATTAGGCATTGAGc                                                                                           <  2:3005611/148‑1 (MQ=255)
                                                          taaaTCATGCCAGCAACACCTGACACACCCCTCAAACCCGGCCAGCGCCGATCTACTGCCAGCCGTAGTAGTCTCCCCCCCCGGGCTTTCGAAAAACAAAAAATCGCCGATTCCGGTGCTGGATTTAGCGCCATTAGGCATTGAGCCa                                                                                         <  1:983765/147‑1 (MQ=255)
                                                          gAAATCATGCCAGCAACAACTGACACACCCCTCAAACCCTCCCAGCGCCTATCTACTGCCAGCCGTAGTAGTCTCCCCTCCCTGGCTTTCCAAAAATAAAAAATCGCCGATTCCGGTGCTGGATTTAGCGCCATTAGGCATTGAGCCa                                                                                         <  1:640754/148‑1 (MQ=255)
                                                          caaaTCATCCCAGCAACAACTCACACACCCCTCAACCCATGCCAGCCCCGATCAACTCCCAGCCGTAGTAGTCCCCCCTCCCTGGCTTTCGAAAAACAAAAAATCGCCGATTCCGGTGCTGGATTTAGCGCCATTAGGCATTGAGCCa                                                                                         <  1:2637105/147‑1 (MQ=255)
                                                           aaaTCATGCCAGCAACAACTGACACAGCCCTCAAACCATGCCAGCGCCGATCTACTGCCAGCCGTAGTAGTCTCCCCTCCCTGGCTTTCGAAAAATAAAAAATCGCCGATTCCGGTGCTGGATTTAGCGCCATTAGGCATTGAGCCaa                                                                                        <  1:2202473/148‑1 (MQ=255)
                                                           aaaTCATGCCAGCAACAACTGACACAGCCCTCAAACCATGCCAGCGCCGATCTACTGCCAGCCGTAGTAGTCTCCCCTCCCTGGCTTTCGAAAAAGAAAAAATCGCCGATTCCGGTGCTGGATTTAGCGCCATTAGGCATTGAGCCaa                                                                                        >  2:1534900/1‑148 (MQ=255)
                                                           aaaTCATGCCAGCAACAACTGACACAGCCCTCAAACCATGCCAGCGCCGATCTACTGCCAGCCGTAGTAGTCTCCCCTCCCTGGCTTTCGAAAAAGAAAAAATCGCCGATTCCGGTGCTGGATTTAGCGCCATTAGGCATTGAGCCaa                                                                                        >  1:1814763/1‑148 (MQ=255)
                                                            aaTCATGCCAGCAACAACTGACACAGCCCTCAAACCATGCCAGCGCCGATCTACTGCCAGCCGTAGTAGTCTCCCCTCCCTGGCTTTCGAAAAAGAAAAAATCGCCGATTCCGGTGCTGGATTTAGCGCCATTAGGCATTGAGCCAAt                                                                                       >  1:1610486/1‑148 (MQ=255)
                                                            aaTCATGCCAGCAACAACTGACACAGCCCTCAAACCATGCCAGCGCCGATCTACTGCCAGCCGTAGTAGTCTCCCCTCCCTGGCTTTCGAAAAAGAAAAAATCGCCGATTCCGGTGCTGGATTTAGCGCCATTAGGCATTGAGCCAAt                                                                                       <  2:2714769/148‑1 (MQ=255)
                                                             aacagcccaccaACACCTAACACCGCCCACAAACCATCCCAGCGCCGATCTACTGCCAGCCGTAGTTGTCTCCCCTCCCGGTCTTTCTAAAAATAAAAAATCGCCGATTCCGGTGCTGGATTTAGCGCCATTAGGCATTGAGCCAATc                                                                                      <  1:917187/142‑1 (MQ=255)
                                                             aTCATGCCAGCAACAACTGACACAGCCCTCAAACCATGCCAGCGCCGATCTACTGCCAGCCGTAGTAGTCTCCCCTCCCTGGCTTTCGAAAAAGAAAAAATCGCCGATTCCGGTGCTGGATTTAGCGCCATTAGGCATTGAGCCAATc                                                                                      >  2:2191610/1‑148 (MQ=255)
                                                               cATGCCAGCAACAACTGACACAGCCCTCAAACCATGCCAGCGCCGATCTACTGCCAGCCGTAGTAGTCTCCCCTCCCTGGCTTTCGAAAAAGAAAAAATCGCCGATTCCGGTGCTGGATTTAGCGCCATTAGGCATTGAGCCAATCTg                                                                                    <  2:1094182/148‑1 (MQ=255)
                                                                aTGCCAGCAACAACTGACACAGCCCTCAAACCATGCCAGCGCCGATCTACTGCCAGCCGTAGTAGTCTCCCCTCCCTGGCTTTCGAAAAAGAAAAAATCGCCGATTCCGGTGCTGGATTTAGCGCCATTAGGCATTGAGCCAATCTGt                                                                                   >  1:1881352/1‑148 (MQ=255)
                                                                 tGCCAGCAACAACTGACACAGCCCTCAAACCATGCCAGCGCCGATCTACTGCCAGCCGTAGTAGTCTCCCCTCCCTGGCTTTAGAAAAAGAAAAAATCGCCGATTCCGGTGCTGGATTTAGCGCCATTAGGCATTGAGCCAATCTGtt                                                                                  <  2:228002/148‑1 (MQ=255)
                                                                      gCAACACCTGACACAGCCCTCAAACCATGCCAGCGCCGATCTACTGCCAGCCGTAGTAGTCTCCCCTCCCTGGCTTTCGAAAAAGAAAAAATCGCCGATTCCGGTGCTGGATTTAGCGCCATTAGGCATTGAGCCAATCTGTTATCTg                                                                             <  2:308489/148‑1 (MQ=255)
                                                                      gCAACACCTCACACCGCCCTAAACCCATGCCAGCGCCGATCTCCTGCCCGCCGTCGTATTCGCCCCCCCGTGGATTTCGAAAAATAAAAAATCGCCGATTCCGGTGCTGGATTTAGCGCCATTAGGCATTGAGCCAATCTGTTATCTg                                                                             <  1:867553/148‑1 (MQ=255)
                                                                       caacaaCTGACACAGCCCTCAAACCATGCCAGCGCCGATCTACTGCCAGCCGTAGTAGTCTCCCCTCCCTGGCTTTCGAAAAAGAAAAAATCGCCGATTCCGGTGCTGGATTTAGCGCCATTAGGCATTGAGCCAATCTGTTATCTGa                                                                            >  2:1146924/1‑148 (MQ=255)
                                                                       caacaaCTGACACAGCCCTCAAACCATGCCAGCGCCGATCTACTGCCAGCCGTAGTAGTCTCCCCACCCTGGCTTTAGAAAAATAAAAAATCGCCGATTCCGGTGCTGGATTTAGCGCCATTAGGCATTGAGCCAATCTGTTATCTGa                                                                            <  2:598186/148‑1 (MQ=255)
                                                                        aacaccTCACACACCCCCCAACCCCGCCCCGCGCCGATCTACGGCCCGCCGTAGTAGTCTCCCCCCCCTGGCTTTCGAAAAAAAAAAAATCGCCGATTCCGGTGCTGGATTTAGCGCCATTAGGCATTGAGCCAATCTGTTATCTGAc                                                                           <  1:2507330/148‑1 (MQ=255)
                                                                           aacTGACACAGCCCTCAAACCATGCCAGCGCCGATCTACTGCCAGCCGTAGTAGTCTCCCCTCCCTGGCTTTCGAAAAAGAAAAAATCGCCGATTCCGGTGCTGGATTTAGCGCCATTAGGCATTGAGCCAATCTGTTATCTGACaga                                                                        >  2:3295086/1‑148 (MQ=255)
                                                                                 cacaGCCATCAAAACTGCCCAGCGCCTATCTACTGCCAGGCTTATTTGTCTCCCCTCCCTGGCTTTCTAAAAATAAAAAATCGCCGATTCCGGTGCTGGATTTAGCGCCATTAGGCATTGAGCCAATCTGTTATCTGACAGAAGAAAt                                                                  <  1:1026795/148‑1 (MQ=255)
                                                                                 cacaCCCCACAAACCATGCCCGCGCCGATCTACTGCCAGCCTTAGTAGTCTCCCCTCCCTGGATTTCTAAAAATAAAAAATCGCCGATTCCGGTGCTGGATTTAGCGCCATTAGGCATTGAGCCAATCTGTTATCTGACAGAAGAAAt                                                                  <  1:603583/148‑1 (MQ=255)
                                                                                  acaGCCCTCAAACCATCCCAGCGCCGATCTACTGCCAGCCGTAGTAGTCTCCCCTCCCTGGCTTTAGAAAAAGAAAAAATCGCCGATTCCGGTGCTGGATTTAGCGCCATTAGGCATTGAGCCAATCTGTTATCTGACAGAAGAAATc                                                                 <  2:332878/148‑1 (MQ=255)
                                                                                  acaGCCCTCAAACCAGCCCAGCGCCTATCTCCGGCCAGCCGTAGTAGTCTCCCCTCCCTGGCTTTCGAAAAATAAAAAATCGCCGATTCCGGTGCTGGATTTAGCGCCATTAGGCATTGAGCCAATCTGTTATCTGACAGAAGAAATc                                                                 <  1:2243464/148‑1 (MQ=255)
                                                                                         tCAAACCATGCCAGCGCCGATCTACTGCCAGCCGTAGTAGTCTCCCCTCCCTGGCTTTCGAAAAAGAAAAAATCGCCGATTCCGGTGCTGGATTTAGCGCCATTAGGCATTGAGCCAATCTGTTATCTGACAGAAGAAATCAGTAATc                                                          >  2:901667/1‑148 (MQ=255)
                                                                                          cAAACCATGCCAGCGCCGATCTACTGCCAGCCGTAGTAGTCTCCCCTCCCTGGCTTTCGAAAAAGAAAAAATCGCCGATTCCGGTGCTGGATTTAGCGCCATTAGGCATTGAGCCAATCTGTTATCTGACAGAAGAAATCAGTAATCa                                                         >  1:2805399/1‑148 (MQ=255)
                                                                                          aaaaCCATGCCAGCGCCTATCTACTGCCAGCCGTAGTAGTCTCCCCTCCCTGGATTTCTAAAAAGAAAAAATCGCCGATTCCGGTGCTGGATTTAGCGCCATTAGGCATTGAGCCAATCTGTTATCTGACAGAAGAAATCAGTAATCa                                                         <  1:2043220/147‑1 (MQ=255)
                                                                                           ataCCATGCCAGCGCCTATCTACGGCTCGGCGGATTAGTCTCCCCTCCCTGGATTTCGAAAAAGAAATAATCGCCGATTCCTGTGCTGGATTTAGCGCCTTTAGGCATTGAGCCAATCTGTTATCTGACAGAAGTAATCCGTAATCaa                                                        <  2:2061039/146‑1 (MQ=255)
                                                                                           aaaCCATGCCAGCGCCGATCTACTGCCAGCCGTAGTAGTCTCCCCTCCCGGGCTTTCGAAAACGAAAAAATCGCCGATTCCGGTGCTGGATTTAGCGCCATTAGGCATTGAGCCAATCTGTTATCTGACAGAAGAAATCAGTAATCaa                                                        <  2:2460159/148‑1 (MQ=255)
                                                                                              ccATGCCAGCGCCGATCTACTGCCAGCCGTAGTAGTCTCCCCTCCCTGGCTTTCGAAAAAGAAAAAATCGCCGATTCCGGTGCTGGATTTAGCGCCATTAGGCATTGAGCCAATCTGTTATCTGACAGAAGAAATCAGTAATCAACtt                                                     >  1:2673374/1‑148 (MQ=255)
                                                                                              ccATGCCAGCGCCGATCTACTGCCAGCCGTAGTAGTCTCCCCTCCCTGGCTTTCGAAAAAGAAAAAATCGCCGATTCCGGTGCTGGATTTAGCGCCATTAGGCATTGAGCCAATCTGTTATCTGACAGAAGAAATCAGTAATCAACtt                                                     >  2:1765852/1‑148 (MQ=255)
                                                                                               cATGCCAGCGCCGATCTACTGCCAGCCGTAGTAGTCTCCCCTCCCTGGCTTTCGAAAAATAAAAAATCGCCGATTCCGGTGCTGGATTTAGCGCCATTAGGCATTGAGCCAATCTGTTATCTGACAGAAGAAATCAGTAATCAACttt                                                    <  1:3008028/148‑1 (MQ=255)
                                                                                                 tGCCAGCGCCGATCTACTGCCAGCCGTAGTAGTCTCCCCTCCCTGGCTTTCGAAAAAGAAAAAATCGCCGATTCCGGTGCTGGATTTAGCGCCATTAGGCATTGAGCCAATCTGTTATCTGACAGAAGAAATCAGTAATCAACTTTTg                                                  >  2:2561468/1‑148 (MQ=255)
                                                                                                     aGCGCCGATCTACTGCCAGCCGTAGTAGTCTCCCCTCCCTGGCTTTCGAAAAAGAAAAAATCGCCGATTCCGGTGCTGGATTTAGCGCCATTAGGCATTGAGCCAATCTGTTATCTGACAGAAGAAATCAGTAATCAACTTTTGGCGa                                              >  1:1140266/1‑148 (MQ=255)
                                                                                                       cgcCGATCTACTGCCAGCCGTAGTATTCTCCCCTCCCTGGATTTCGAAAAACAAAAAATCGCCGATTCCGGTGCTGGATTTAGCGCCATTAGGCATTGAGCCAATCTGTTATCTGACAGAAGAAATCAGTAATCAACTTTTGGCGaaa                                            <  1:411920/148‑1 (MQ=255)
                                                                                                        gcCGATCTACTGCCAGCCGTAGTAGTCTCCCCTCCCTGGCTTTCGAAAAAGAAAAAATCGCCGATTCCGGTGCTGGATTTAGCGCCATTAGGCATTGAGCCAATCTGTTATCTGACAGAAGAAATCAGTAATCAACTTTTGGCGAAat                                           <  1:2285408/148‑1 (MQ=255)
                                                                                                         ccGATCTACTGCCAGCCGTAGTAGTCTCCCCTCCCTGGCTTTCGAAAAAGAAAAAATCGCCGATTCCGGTGCTGGATTTAGCGCCATTAGGCATTGAGCCAATCTGTTATCTGACAGAAGAAATCAGTAATCAACTTTTGGCGAAata                                          <  2:3133386/148‑1 (MQ=255)
                                                                                                          cGATCTACTGCCAGCCGTAGTAGTCTCCCCTCCCTGGCTTTCGAAAAAGAAAAAATCGCCGATTCCGGTGCTGGATTTAGCGCCATTAGGCATTGAGCCAATCTGTTATCTGACAGAAGAAATCAGTAATCAACTTTTGGCGAAatat                                         >  2:2883051/1‑148 (MQ=255)
                                                                                                          cGATCTACTGCCAGCCGTAGTAGTCTCCCCTCCCTGGCTTTCGAAAAAGAAAAAATCGCCGATTCCGGTGCTGGATTTAGCGCCATTAGGCATTGAGCCAATCTGTTATCTGACAGAAGAAATCAGTAATCAACTTTTGGCGAAatat                                         >  1:790044/1‑148 (MQ=255)
                                                                                                          cGATCTACTGCCAGCCGTAGTAGTCTCCCCTCCCTGGCATTCGAAAAAGAAAAAATCGCCGATTCCGGTGCTGGATTTAGCGCCATTAGGCATTGAGCCAATCTGTTATCTGACAGAAGAAATCAGTAATCAACTTTTGGCGAAatat                                         >  2:2882865/1‑148 (MQ=255)
                                                                                                          cGATCTACTGCCAGCCGTAGTAGTCTCCCCTCCCTGGATTTCTAAAAATAAAAAATCGCCGATTCCGGTGCTGGATTTAGCGCCATTAGGCATTGAGCCAATCTGTTATCTGACAGAAGAAATCAGTAATCAACTTTTGGCGAAatat                                         <  1:617656/148‑1 (MQ=255)
                                                                                                            aTCTACGGCCAGCCGTAGTAGTCGCCCCTCCCTGGCTTTCGAAAAAGAAAAAATCGCCGATTCCGGTGCTGGATTTAGCGCCATTAGGCATTGAGCCAATCTGTTATCTGACAGAAGAAATCAGTAATCAACTTTTGGCGAAatatat                                       <  2:887536/148‑1 (MQ=255)
                                                                                                             tCTACTGCCCGCCGTAGTAGTCCCCCCCCCCTGGCTTTCCAAAAACAAAAAATCGCCGATTCCGGTGCTGGATTTAGCGCCATTAGGCATTGAGCCAATCTGTTATCTGACAGAAGAAATCAGTAATCAACTTTTGGCGAAATATATc                                      <  1:2350108/148‑1 (MQ=255)
                                                                                                             tCTACTGCCAGCCGTAGTAGTCTCCCCTCCCTGGCTTTCGAAAAAGAAAAAATCGCCGATTCCGGTGCTGGATTTAGCGCCATTAGGCATTGAGCCAATCTGTTATCTGACAGAAGAAATCAGTAATCAACTTTTGGCGAAATATATc                                      >  1:2569505/1‑148 (MQ=255)
                                                                                                              cTACTGCCAGCCGTAGTAGTCTCCCCTCCCTGGCTTTAGAAAAAGAAAAAATCGCCGATTCCGGTGCTGGATTTAGCGCCATTAGGCATTGAGCCAATCTGTTATCTGACAGAAGAAATCAGTAATCAACTTTTGGCGAAATATATCt                                     <  2:2944250/148‑1 (MQ=255)
                                                                                                              cTACTGCCAGCCGTAGTAGTCTCCCCTCCCTGGCTTTAGAAAAAGAAAAAATCGCCGATTCCGGTGCTGGATTTAGCGCCATTAGGCATTGAGCCAATCTGTTATCTGACAGAAGAAATCAGTAATCAACTTTTGGCGAAATATATCt                                     <  2:374660/148‑1 (MQ=255)
                                                                                                               tACTGCCCGCCTTATTAGTCTCCCCTCCCTGGCTTTCGAAAAATAAAAAATCGCCGATTCCGGTGCTGGATTTAGCGCCATTAGGCATTGAGCCAATCTGTTATCTGACAGAAGAAATCAGTAATCAACTTTTGGCGAAATATATCTg                                    <  1:3191198/148‑1 (MQ=255)
                                                                                                                aCTGCCAGCCGTAGTAGTCTCCCCTCCCTGGCTTTCGAAAAAGAAAAAATCGCCGATTCCGGTGCTGGATTTAGCGCCATTAGGCATTGAGCCAATCTGTTATCTGACAGAAGAAATCAGTAATCAACTTTTGGCGAAATATATCTgg                                   >  2:371115/1‑148 (MQ=255)
                                                                                                                  tGCCAGCCGTAGTAGTCTCCCCTCCCTGGCTTTCGAAAAAGAAAAAATCGCCGATTCCGGTGCTGGATTTAGCGCCATTAGGCATTGAGCCAATCTGTTATCTGACAGAAGAAATCAGTAATCAACTTTTGGCGAAATATATCTGGTa                                 <  2:1796329/148‑1 (MQ=255)
                                                                                                                   gCCAGCCGTAGTAGTCTCCCCTCCCTGGCTTTCGAAAAAGAAAAAATCGCCGATTCCGGTGCTGGATTTAGCGCCATTAGGCATTGAGCCAATCTGTTATCTGACAGAAGAAATCAGTAATCAACTTTTGGCGAAATATATCTGGTAt                                >  1:638019/1‑148 (MQ=255)
                                                                                                                    ccAGCCGTAGTAGTCTCCCCTCCCTGGCTTTCGAAAAAGAAAAAATCGCCGATTCCGGTGCTGGATTTAGCGCCATTAGGCATTGAGCCAATCTGTTATCTGACAGAAGAAATCAGTAATCAACTTTTGGCGAAATATATCTGGTAtt                               <  2:3230675/148‑1 (MQ=255)
                                                                                                                              tagtCTCCCCTCCCTGGCTTTCGAAAAAGAAAAAATCGCCGATTCCGGTGCTGGATTTAGCGCCATTAGGCATTGAGCCAATCTGTTATCTGACAGAAGAAATCAGTAATCAACTTTTGGCGAAATATATCTGGTATTCAAAACACAt                     <  2:1127690/148‑1 (MQ=255)
                                                                                                                               agtCTCCCCTCCCTGGCTTTAGAAAAAGAAAAAATCGCCGATTCCGGTGCTGGATTTAGCGCCATTAGGCATTGAGCCAATCTGTTATCTGACAGAAGAAATCAGTAATCAACTTTTGGCGAAATATATCTGGTATTCAAAACACATc                    <  2:1232619/148‑1 (MQ=255)
                                                                                                                                 tctcCCCTCCCTGGCTTTCGAAAAAGAAAAAATCGCCGATTCCGGTGCTGGATTTAGCGCCATTAGGCAAAGATCCAATCTGTTATCTGACAGAAGAAATCAGTAATCAAAATTTTGCGAAATATATCTGGTATTCAAAAAACATCAc                  >  2:1276938/1‑148 (MQ=255)
                                                                                                                                  ctcccctcccTGGCTTTCGAAAAATAAAAAATCGCCGATTCCGGTGCTGGATTTAGCGCCATTAGGCATTGAGCCAATCTGTTATCTGACAGAAGAAATCAGTAATCAACTTTTGGCGAAATATATCTGGTATTCAAAACACATCACg                 <  1:2312226/148‑1 (MQ=255)
                                                                                                                                      cctccctGGCTTTCGAAAAAGAAAAAATCGCCGATTCCGGTGCTGGATTTAGCGCCATTAGGCATTGAGCCAATCTGTTATCTGACAGAAGAAATCAGTAATCAACTTTTGGCGAAATATATCTGGTATTCAAAACACATCACGGTTa             <  2:2007997/148‑1 (MQ=255)
                                                                                                                                        tccctGGCTTTCTAAAAAGAAAAAATCGCCGATTCCGGTGCTGGATTTAGCGCCATTAGGCATTGAGCCAATCTGTTATCTGACAGAAGAAATCAGTAATCAACTTTTGGCGAAATATATCTGGTATTCAAAACACATCACGGTTAGc           <  1:2932960/148‑1 (MQ=255)
                                                                                                                                         ccctGGCTTTCGAAAAAGAAAAAATCGCCGATTCCGGTGCTGGATTTAGCGCCATTAGGCATTGAGCCAATCTGTTATCTGACAGAAGAAATCAGTAATCAACTTTTGGCGAAATATATCTGGTATTCAAAACACATCACGGTTAGcc          >  2:2114768/1‑148 (MQ=255)
                                                                                                                                         ccctGGCTTTCGAAAAAGAAAAAATCGCCGATTCCGGTGCTGGATTTAGCGCCATTAGGCATTGAGCCAATCTGTTATCTGACAGAAGAAATCAGTAATCAACTTTTGGCGAAATATATCTGGTATTCAAAACACATCACGGTTAGcc          <  2:1964497/148‑1 (MQ=255)
                                                                                                                                         ccctGGCTTTCGAAAAAGAAAAAATCGCCGATTCCGGTGCTGGATTTAGCGCCATTAGGCATTGAGCCAATCTGTTATCTGACAGAAGAAATCAGTAATCAACTTTTGGCGAAATATATCTGGTATTCAAAACACATCACGGTTAGcc          >  1:2890287/1‑148 (MQ=255)
                                                                                                                                          cctGGCTTTCGAAAAAGAAAAAATCGCCGATTCCGGTGCTGGATTTAGCGCCATTAGGCATTGAGCCAATCTGTTATCTGACAGAAGAAATCAGTAATCAACTTTTGGCGAAATATATCTGGTATTCAAAACACATCACGGTTAGCCa         <  2:2597225/148‑1 (MQ=255)
                                                                                                                                              gCTTTCGAAAAAGAAAAAATCGCCGATTCCGGTGCTGGATTTAGCGCCATTAGGCATTGAGCCAATCTGTTATCTGACAGAAGAAATCAGTAATCAACTTTTGGCGAAATATATCTGGTATTCAAAACACATCACGGTTAGCCATgaa     >  2:2474750/1‑148 (MQ=255)
                                                                                                                                              gCTTTCGAAAAAGAAAAAATCGCCGATTCCGGTGCTGGATTTAGCGCCATTAGGCATTGAGCCAATCTGTTATCTGACAGAAGAAATCAGTAATCAACTTTTGGCGAAATATATCTGGTATTCAAAACACATCACGGTTAGCCATgaa     >  1:1246243/1‑148 (MQ=255)
                                                                                                                                              gCTTTCGAAAAAGAAAAAATCGCCGATTCCGGTGCTGGATTTAGCGCCATTAGGCATTGAGCCAATCTGTTATCTGACAGAAGAAATCAGTAATCAACTTTTGGCGAAATATATCTGGTATTCAAAACACATCACGGTTAGCCATgaa     >  2:2626102/1‑148 (MQ=255)
                                                                                                                                                tttCGAAAAAGAAAAAATCGCCGATTCCGGTGCTGGATTTAGCGCCATTAGGCATTGAGCCAATCTGTTATCTGACAGAAGAAATCAGTAATCAACTTTTGGCGAAATATATCTGGTATTCAAAACACATCACGGTTAGCCATgaaga   >  2:1808149/1‑148 (MQ=255)
                                                                                                                                                 ttCGAAAAAGAAAAAATCGCCGATTCCGGTGCTGGATTTAGCGCCATTAGGCATTGAGCCAATCTGTTATCTGACAGAAGAAATCAGTAATCAACTTTTGGCGAAATATATCTGGTATTCAAAACACATCACGGTTAGCCATgaagaa  >  1:2108864/1‑148 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
AACCAGCACTGGCAGAACAGGTCATTCCCTGGCTCTCGACACCCGCAGTTGCCGTACTGAAATCATGCCAGCAACAACTGACACAGCCCTCAAACCATGCCAGCGCCGATCTACTGCCAGCCGTAGTAGTCTCCCCTCCCTGGCTTTCGAAAAAGAAAAAATCGCCGATTCCGGTGCTGGATTTAGCGCCATTAGGCATTGAGCCAATCTGTTATCTGACAGAAGAAATCAGTAATCAACTTTTGGCGAAATATATCTGGTATTCAAAACACATCACGGTTAGCCATGAAGAA  >  CP009273/2191349‑2191641

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