Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,189,888 C→A 23.3% intergenic (+76/‑65) metG → / → yehH methionyl‑tRNA synthetase/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,189,8880CA23.3% 505.2 / 126.7 210intergenic (+76/‑65)metG/yehHmethionyl‑tRNA synthetase/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
Reads supporting (aligned to +/- strand):  ref base C (78/83);  new base A (23/26);  total (101/109)
Fisher's exact test for biased strand distribution p-value = 8.72e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

Only 100 of 224 total aligned reads displayed.

CCGGTCCTGGCGGGAAAGATATTTTCCTGCTAAGCCCGGATGCCGGTGCTAAACCGGGTCATCAGGTGAAATAATCCCCCTTCAAGGCGCTGCATCGACAGCGCCTTTTCTTTATAAATTCCTAAAGTTGTTTTCTTGCGATTTTGTCTCTCTCTAACCCGCATAAATACTGGTAGCATCTGCATTCAACTGGATAAAATTACAGGGATGCAGAATGAGACACTTTATCTATCAGGACGAAAAATCACATAAATTCAGGGCAGTTGAGCAACAGGGAAACGAGTTGCATATCAGTTGGG  >  CP009273/2189739‑2190037
                                                                                                                                                     |                                                                                                                                                     
ccGGTCCTGGCGGGAAAGATATTTTCCTGCTAAGCCCGGATGCCGGTGCTAAACCGGGTCATCAGGTGAAATAATCCCCCTTCAAGGCGCTGCATCGACAGCGCCTTTTCTTTATAAATTCCTAAAGTTGTTTTCTTGCGATTTTGtctc                                                                                                                                                       <  1:3952504/150‑1 (MQ=255)
   gTCCTGGCGGGAAAGATATTTTCCTGCTAAGCCCGGATGCCGGTGCTAAACCGGGTCATCAGGTGAAATAATCCCCCTTCAAGGCGCTGCATCGACAGCGCCTTTTCTTTATAAATTCCTAAAGTTGTTTTCTTGCGATTTTGtctctct                                                                                                                                                    >  2:3026826/1‑150 (MQ=255)
       ttgCGGGAAAGATATTTTCCTGCTAAGCCCGGATGCCGGTGCTAAACCGGGTCATCAGGTGAAATAATCCCCCTTCAAGGCGCTGCATCGACAGCGCCTTTTCTTTATACATTCCGAAAGTTGTATTCTTGCGATGTTGTCTATCTCTaa                                                                                                                                                >  2:4050066/3‑150 (MQ=255)
        ggCGGGAAAGATATTGTCCTGCTAAGCCCGGATGCCGGTGCTAAACCGGGTCATCAGGTGAAATAATCCCCCTTCACGGCGCTGCATCGACAGCGCCTTTTCTTTATAAATTCCTAAAGTTGTTTTCTTGCGATTTTGTCTATCTCTAAc                                                                                                                                               <  2:504049/150‑1 (MQ=255)
         gCGGGAAAGATATTTTCCTGCTAAGCCCGGATGCCGGTGCTAAACCGGGTCATCAGGTGAAATAATCCCCCTTCAAGGCGCTGCATCGACAGCGCCTTTTCTTTATAAATTCCTAAAGTTGTTTTCTTGCGATTTTGTCTCTCTCTAAcc                                                                                                                                              <  1:4070702/150‑1 (MQ=255)
          cGGGAAAGATATTTTCCTGCTAAGCCCGGATGCCGGTGCTAAACCGGGTCATCAGGTGAAATAATCCCCCATCAAGGCGCTGCATCGACAGCGCCTTTTCTTTATAAATTCCTAAAGTTGTTTTCTTGAGATTTAGTCTCTCTCTAAccc                                                                                                                                             >  2:3441481/1‑150 (MQ=255)
           tggAAAGATATTTTCCTGCTAAGCCCGGATGCCGGTGCTAACCCGGGTCATCAGGTGAAATAATCCCCCTTCAAGGCGCTGCATCGACAGCGCCTTTTCTTTATAAATTCCTAAAGTTGTTTTCTTGCGATTTTGTCTCGCTCTAACCCg                                                                                                                                            >  1:245748/2‑150 (MQ=255)
           gggAAAGATATTTTCCTGCTAAGCCCGGATGCCGGTGCTAAACCGGGTCATCAGGTGAAATAATCCCCCTTCAAGGCGCTGCATCGACAGCGCCTTTTCTTTATACATTCCTAAAGTTGTTTTCTTGCGATTTTGTGTCTCTCTAACCCg                                                                                                                                            >  1:3271446/1‑150 (MQ=255)
                aGATATTTTCCTGCTAAGCCCGGATGCCGGTGCTAAACCGGGTCATCAGGTGAAATAATCCCCCTTCAAGGCGCTGCATCGACAGCGCCTTTTCTTTATAAATTCCTAAAGTTGTTTTCTTGCGATTTTGTCTCTCTCTAACCCGCATaa                                                                                                                                       >  1:1758038/1‑150 (MQ=255)
                   tatTTTCCTGCTAAGCCCGGATGCCGGTGCTAAACCGGGTCATCAGGTGAAATAATCCCCCTTCAAGGCGCTGCATCGACAGCGCCTTTTCTTTATAAATTCCTAAAGTTGTTTTCTTGCGATTTTGTCTCTCTCTAACCCGCATAAATa                                                                                                                                    >  2:2284926/1‑150 (MQ=255)
                      tttCCTGCTAAGCCCGGATGCCGGTGCTAAACCGGGTCATCAGGTGAAATAATCCCCCTTCAAGGCGCTGCATCGACAGCGCCTTTTCTTTATAAATTCCTAAAGTTGTTTTCTTGCGATTTTGTCTATCTCTAACCCGCATAAATACTg                                                                                                                                 <  1:728013/150‑1 (MQ=255)
                       ttCCTGCTAAGCCCGGATGCCGGTGCTAAACCGGGTCATCAGGTGAAATAATCCCCCTTCAAGGCGCTGCATCGACAGCGCCTTTTCTTTATAAATTCCTAAAGTTGTTTTCTTGCGATTTTGTCTCTCTCTAACCCGCATAAATACTgg                                                                                                                                <  2:436826/150‑1 (MQ=255)
                       ttCCTGCTAAGCCCGGATGCCGGTGCTAAACCGGGTCATCAGGTGAAATAATCCCCCTTCAAGGCGCTGCATCGACAGCGCCTTTTCTTTATAAATTCCTAAAGTTGTTTTCTTGCGATTTTGTCTCTCTCTAACCCGCATAAATACTgg                                                                                                                                >  1:3178417/1‑150 (MQ=255)
                          cTGCTAAGCCCGGATGCCGGTGCTAAACCGGGTCATCAGGTGAAATAATCCCCCTTCAAGGCGCTGCATCGACAGCGCCTTTTCTTTATAAATTCCTAAAGTTGTTTTCTTGCGATTTTGTCTATCTCTAACCCGCATAAATACTGGTAg                                                                                                                             <  2:3222546/150‑1 (MQ=255)
                             cTAAGCCCGGATGCCGGTGCTAAACCGGGTCATCAGGTGAAATAATCCCCCTTCAAGGCGCTGCATCGACAGCGCCTTTTCTTTATAAATTCCTAAAGTTGTTTTCTTGCGATTTTGTCTCTCTCTAACCCGCATAAATACTGGTAGCAt                                                                                                                          <  2:3478234/150‑1 (MQ=255)
                             cTAAGCCCGGATGCCGGTGCTAAACCGGGTCATCAGGTGAAATAATCCCCCTTCAAGGCGCTGCATCGACAGCGCCTTTTCTTTATAAATTCCTAAAGTTGTTTTCTTGCGATTTTGTCTATCTCTAACCCGCATAAATACTGGTAGCAt                                                                                                                          >  2:598524/1‑150 (MQ=255)
                                aGCCCGGATGCCGGTGCTAAACCGGGTCATCAGGTGAAATAATCCCCCTTCAAGGCGCTGCATCGACAGCGCCTTTTCTTTATAAATTCCTAAAGTTGTTTTCTTGCGATTTTGTCTCTCTCTAACCCGCATAAATACTGGTAGCATCTg                                                                                                                       <  1:2769044/150‑1 (MQ=255)
                                aGCCCGGATGCCGGTGCTAAACCGGGTCATCAGGTGAAATAATCCCCCTTCAAGGCGCTGCATCGACAGCGCCTTTTCTTTATAAATTCCTAAAGTTGTTTTCTTGCGATTTTGTCTCTCTCTAACCCGCATAAATACTGGTAGCATCTg                                                                                                                       >  1:3067253/1‑150 (MQ=255)
                                aGCCCGGATGCCGGTGCTAAACCGGGTCATCAGGTGAAATAATCCCCCTTCAAGGCGCTGCATCGACAGCGCCTTTTCTTTATAAATTCCTAAAGTTGTTTTCTTGCGATTTTGTCTCTCTCTAACCCGCATAAATACTGGTAGCATCTg                                                                                                                       <  2:1077736/150‑1 (MQ=255)
                                    cGGATGCCGGTGCTAAACCGGGTCATCAGGTGAAATAATCCCCCTTCAAGGCGCTGCATCGACAGCGCCTTTTCTTTATAAATTCCTAAAGTTGTTTTCTTGCGATTTTGTCTCTCTCTAACCCGCATAAATACTGGTAGCATCTGCAtt                                                                                                                   <  1:3828074/150‑1 (MQ=255)
                                    cGGATGCCGGTGCTAAACCGGGTCATCAGGTGAAATAATCCCCCTTCAAGGCGCTGCATCGACAGCGCCTTTTCTTTATAAATTCCTAAAGTTGTTTTCTTGCGATTTTGTCTCTCTCTAACCCACATAAATACTGGTAGCATCTGCAtt                                                                                                                   <  1:523581/150‑1 (MQ=255)
                                    cGGATGCCGGTGCTAAACCGGGTCATCAGGTGAAATAATCCCCCTTCAAGGCGCTGCATCGACAGCGCCTTTTCTTTATAAATTCCTAAAGTTGTTTTCTTGCGATTTTGTCTCTCTCTAACCCACATAAATACTGGTAGCATCTGCAtt                                                                                                                   <  1:1493354/150‑1 (MQ=255)
                                     ggATGCCGGTGCTAAACCGGGTCATCAGGTGAAATAATCCCCCTTCAAGGCGCTGCATCGACAGCGCCTTTTCTTTATAAATTCCTAAAGTTGTTTTCTTGCGATTTTGTCTATCTCTAACCCGCATAAATACTGGTAGCATCTGCATTc                                                                                                                  >  1:1695572/1‑150 (MQ=255)
                                      gATGCCGGTGCTAAACCGGGTCATCAGGTGAAATAATCCCCCTTCAAGGCGCTGCATCGACAGCGCCTTTTCTTTATAAATTCCTAAAGTTGTTTTCTTGCGATTTTGTCTCTCTCTAACCCGCATAAATACTGGTAGCATCTGCATTCa                                                                                                                 <  2:1279219/150‑1 (MQ=255)
                                         gCCGGTGCTAAACCGGGTCATCAGGTGAAATAATCCCCCTTCAAGGCGCTGCATCGACAGCGCCTTTTCTTTATAAATTCCTAAAGTTGTTTTCTTGCGATTTTGTCTCTCTCTAACCCGCATAAATACTGGTAGCATCTGCATTCAACt                                                                                                              <  1:1214290/150‑1 (MQ=255)
                                         gCCGGTGCTAAACCGGGTCATCAGGTGAAATAATCCCCCTTCAAGGCGCTGCATCGACAGCGCCTTTTCTTTATAAATTCCTAAAGTTGTTTTCTTGCGATTTTGTCTCTCTCTAACCCGCATAAATACTGGTAGCATCTGCATTCAACt                                                                                                              <  2:863731/150‑1 (MQ=255)
                                           cGGTGCTAAACCGGGTCATCAGGTGAAATAATCCCCCTTCAAGGCGCTGCATCGACAGCGCCTTTTCTTTATAAATTCCTAAAGTTGTTTTCTTGCGATTTTGTCTCTCTCTAACCCGCATAAATACTGGTAGCATCTGCATTCAACTgg                                                                                                            <  2:3075231/150‑1 (MQ=255)
                                              tGCTAAACCGGGTCATCAGGTGAAATAATCCCCCTTCAAGGCGCTGCATCGACAGCGCCATTTCTTTATAAATTCCTAAAGTTGTTTTCTTGCGATTTTGTCTATCTCTAACCCGCATAAATACTGGTAGCATCTGCATTCAACTGGATa                                                                                                         <  1:1330494/150‑1 (MQ=255)
                                               gCTAAACCGGGTCATCAGGTGAAATAATCCCCCTTCAAGGCGCTGCATCGACAGCGCCTTTTCTTTATAAATTCCTAAAGTTGTTTTCTTGCGATTTTGTCTATCTCTAACCCGCATAAATACTGGTAGCATCTGCATTCAACTGGATaa                                                                                                        <  2:1001919/150‑1 (MQ=255)
                                               gCTAAACCGGGTCATCAGGTGAAATAATCCCCCTTCAAGGCGCTGCATCGACAGCGCCTTTTCTTTATAAATTCCTAAAGTTGTTTTCTTGCGATTTTGTCTATCTCTAACCCGCATAAATACTGGTAGCATCTGCATTCAACTGGATaa                                                                                                        <  2:3272306/150‑1 (MQ=255)
                                                   aaCCGGGTCATCAGGTGAAATAATCCCCCTTCAAGGCGCTGCATCGACAGCGCCTTTTCTTTATAAATTCCTAAAGTTGTTTTCTTGCGATTTTGTCTCTCTCTAACCCGCATAAATACTGGTAGCATCTGCATTCAACTGGATAAAAtt                                                                                                    >  2:2258248/1‑150 (MQ=255)
                                                      cGGGTCATCAGGTGAAATAATCCCCCTTCAAGGCGCTGCATCGACAGCGCCTTTTCTTTATAAATTCCTAAAGTTGTTTTCTTGCGATTTTGTCTCTCTCTAACCCGCATAAATACTGGTAGCATCTGCATTCAACTGGATAAAATTACa                                                                                                 <  1:1342944/150‑1 (MQ=255)
                                                      cGGGTCATCAGGTGAAATAATCCCCCTTCAAGGCGCTGCATCGACAGCGCCTTTTCTTTATAAATTCCTAAAGTTGTTTTCTTGCGATTTTGTCTCTCTCTAACCCGCATAAATACTGGTAGCATCTGCATTCAACTGGATAAAATTACa                                                                                                 <  2:3008453/150‑1 (MQ=255)
                                                        ggTCATCAGGTGAAATAATCCCCCTTCAAGGCGCTGCATCGACAGCGCCTTTTCTTTATAAATTCCTAAAGTTGTTTTCTTGCGATTTTGTCTCTCTCTAACCCGCATAAATACTGGTAGCATCTGCATTCAACTGGATAAAATTACAgg                                                                                               >  1:3119531/1‑150 (MQ=255)
                                                           catcaGGTGAAATAATCCCCCTTCAAGGCGCTGCATCGACAGCGCCTTTTCTTTATAAATTCCTAAAGTTGTTTTCTTGCGATTTTGTCTCTCTCTAACCCGCATAAATACTGGTAGCATCTGCATTCAACTGGATAAAATTACAGGGAt                                                                                            >  1:327173/1‑150 (MQ=255)
                                                           catcaGGTGAAATAATCCCCCTTCAAGGCGCTGCATCGACAGCGCCTTTTCTTTATAAATTCCTAAAGTTGTTTTCTTGCGATTTTGTCTCTCTCTAACCCGCATAAATACTGGTAGCATCTGCATTCAACTGGATAAAATTACAGGGAt                                                                                            >  1:3430300/1‑150 (MQ=255)
                                                              caGGTGAAATAATCCCCCTTCAAGGCGCTGCATCGACAGCGCCTTTTCTTTATAAATTCCTAAAGTTGTTTTCTTGCGATTTTGTCTCTCTCTAACCCGCATAAATACTGGTAGCATCTGCATTCAACTGGATAAAATTACAGGGATGCa                                                                                         >  1:2973313/1‑150 (MQ=255)
                                                                 gTGAAATAATCCCCCTTCAAGGCGCTGCATCGACAGCGCCTTTTCTTTATAAATTCCTAAAGTTGTTTTCTTGCGATTTTGTCTCTCTCTAACCCGCATAAATACTGGTAGCATCTGCATTCAACTGGATAAAATTACAGGGATGCAGaa                                                                                      <  1:833454/150‑1 (MQ=255)
                                                                  tGAAATAATCCCCCTTCAAGGCGCTGCATCGACAGCGCCTTTTCTTTATAAATTCCTAAAGTTGTTTTCTTGCGATTTTGTCTCTCTCTAACCCGCATAAATACTGGTAGCATCTGCATTCAACTGGATAAAATTACAGGGATGCAGAAt                                                                                     <  2:4162432/150‑1 (MQ=255)
                                                                    aaaTAATCCCCCTTCAAGGCGCTGCATCGACAGCGCCTTTTCTTTATAAATTCCTAAAGTTGTTTTCTTGCGATTTTGTCTCTCTCTAACCCGCATAAATACTGGTAGCATCTGCATTCAACTGGATAAAATTACAGGGATGCAGAATga                                                                                   <  1:3053658/150‑1 (MQ=255)
                                                                      ataatCCCCCTTCAAGGCGCTGCATCGACAGCGCCTGTTCTTTCTAAATTCCTAAAGTTGTTTTCTTGCGATTTTGTCTCTCTCTAACCCGCATAAATACTGGTAGCATCTGCATTCAACTGGCTAAAATTACAGGGATGCAGAATgaga                                                                                 >  2:1489459/1‑150 (MQ=255)
                                                                       taatCCCCCTTCAAGGCGCTGCATCGACAGCGCCTTTTCTTTATAAATTCCTAAAGTTGTTTTCTTGCGATTTTGTCTCTCTCTAACCCGCATAAATACTGGTAGCATCTGCATTCAACTGGATAAAATTACAGGGATGCAGAATGAGac                                                                                >  1:581597/1‑150 (MQ=255)
                                                                         atCCCCCTTCAAGGCGCTGCATCGACAGCGCCTTTTCTTTATAAATTCCTAAAGTTGTTTTCTTGCGATTTTGTCTCTCTCTAACCCGCATAAATACTGGTAGCATCTGCATTCAACTGGATAAAATTACAGGGATGCAGAATGAGacac                                                                              >  2:3168048/1‑150 (MQ=255)
                                                                           cccccTTCAAGGCGCTGCATCGACAGCGCCTTTTCTTTATAAATTCCTAAAGTTGTTTTCTTGCGATTTTGTCTCTCTCTAACCCGCATAAATACTGGTAGCATCTGCATTCAACTGGATAAAATTACAGGGATGCAGAATGAGACACtt                                                                            <  2:4233616/150‑1 (MQ=255)
                                                                             cccTTCAAGGCGCTGCATCGACAGCGCCTTTTCTTTATAAATTCCTAAAGTTGTTTTCTTGCGATTTTGTCTCTCTCTAACCCGCATAAATACTGGTAGCATCTGCATTCAACTGGATAAAATTACAGGGATGCAGAATGAGACACTTta                                                                          >  1:3472903/1‑150 (MQ=255)
                                                                              ccTTCAAGGCGCTGCATCGACAGCGCCTTTTCTTTATAAATTCCTAAAGTTGTTTTCTTGCGATTTTGTCTCTCTCTAACCCGCATAAATACTGGTAGCATCTGCATTCAACTGGATAAAATTACAGGGATGCAGAATGAGACACTTtat                                                                         <  1:2947901/150‑1 (MQ=255)
                                                                               cTTCAAGGCGCTGCATCGACAGCGCCTTTTCTTTATAAATTCCTAAAGTTGTTTTCTTGCGATTTTGTCTCTCTCTAACCCGCATAAATACTGGTAGCATCTGCATTCAACTGGATAAAATTACAGGGATGCAGAATGAGACACTTtatc                                                                        <  1:3642336/150‑1 (MQ=255)
                                                                                   aaGGCGCTGCATCGACAGCGCCTTTTCTTTATAAATTCCTAAAGTTGTTTTCTTGCGATTTTGTCTCTCTCTAACCCGCATAAATACTTGTATCATCTGCATTCAACTGTATAAAATTAAATTTATGAAGAATTAGACAATTtatatata                                                                    >  1:4166149/1‑149 (MQ=255)
                                                                                   aaGGCGCTGCATCGACAGCGCCTTTTCTTTATAAATTCCTAAAGTTGTTTTCTTGCGATTTTGTCTCTCTCTAACCCGCATAAATACTGGTAGCATCTGCATTCAACTGGATAAAATTACAGGGATGCAGAATGAGACACTTtatctatc                                                                    >  1:1650152/1‑150 (MQ=255)
                                                                                     ggCGCTGCATCGACAGCGCCTTTTCTTTATAAATTCCTAAAGTTGTTTTCTTGCGATTTTGTCTCTCTCTAACCCGCATAAATACTGGTAGCATCTGCATTCAACTGGATAAAATTACAGGGATGCAGAATGAGACACTTTATCTATCAg                                                                  >  1:2748583/1‑150 (MQ=255)
                                                                                       cgcTGCATCGACAGCGCCTTTTCTTTATAAATTCCTAAAGTTGTTTTCTTGCGATTTTGTCTCTCTCTAACCCGCATAAATACTGGTAGCATCTGCATTCAACTGGATAAAATTACAGGGATGCAGAATGAGACACTTTATCTATCAGGa                                                                >  2:1575777/1‑150 (MQ=255)
                                                                                       cgcTGCATCGACAGCGCCTTTTCTTTATAAATTCCTAAAGTTGTTTTCTTGCGATTTTGTCTCTCTCTAACCCGCATAAATACTGGTAGCATCTGCATTCAACTGGATAAAATTACAGGGATGCAGAATGAGACACTTTATCTATCAGGa                                                                <  2:3851995/150‑1 (MQ=255)
                                                                                           gCATCGACAGCGCCTTTTCTTTATAAATTCCTAAAGTTGTTTTCTTGCGATTTTGTCTCTCTCTAACCCGCATAAATACTGGTAGCATCTGCATTCAACTGGATAAAATTACAGGGATGCAGAATGAGACACTTTATCTATCAGGACGaa                                                            >  2:1193731/1‑150 (MQ=255)
                                                                                           gCAGCGACAGCGCCTTTTCTTTATAAATTCCTAAAGTTGTTTTCTTGCGATTTTGTCTATCTCTAACCCGCATAAATACTGGTAGCATCTGCATTCAACTGGATAAAATTACAGGGATGCAGAATGAGACACTTTATCTATCAGGACGaa                                                            <  1:2277302/150‑1 (MQ=255)
                                                                                            cATCGACAGCGCCTTTTCTTTATAAATTCCTAAAGTTGTTTTCTTGCGATTTTGTCTCTCTCTAACCCGCATAAATACTGGTAGCATCTGCATTCAACTGGATAAAATTACAGGGATGCAGAATGAGACACTTTATCTATCAGGACGaaa                                                           <  2:3635475/150‑1 (MQ=255)
                                                                                            cATCGACAGCGCCTTTTCTTTATAAATTCCTAAAGTTGTTTTCTTGCGATTTTGTCTCTCTCTAACCCGCATAAATACTGGTAGCATCTGCATTCAACTGGATAAAATTACAGGGATGCAGAATGAGACACTTTATCTATCAGGACGaaa                                                           >  1:3957136/1‑150 (MQ=255)
                                                                                             aTCGACAGCGCCTTTTCTTTATAAATTCCTAAAGTTGTTTTCTTGCGATTTTGTCTCTCTCTAACCCGCATAAATACTGGTAGCATCTGCATTCAACTGGATAAAATTACAGGGATGCAGAATGAGACACTTTATCTATCAGGACGaaaa                                                          <  1:3077319/150‑1 (MQ=255)
                                                                                               cGACAGCGCCTTTTCTTTATAAATTCCTAAAGTTGTTTTCTTGCGATTTTGTCTATCTCTAACCCGCATAAATACTGGTAGCATCTGCATTCAACTGGATAAAATTACAGGGATGCAGAATGAGACACTTTATCTATCAGGACGAAAAAt                                                        <  1:4281902/150‑1 (MQ=255)
                                                                                                  cAGCGCCTTTTCTTTATAAATTCCTAAAGTTGTTTTCTTGCGATTTTGTCTCTCTCTAACCCGCATAAATACTGGTAGCATCTGCATTCAACTGGATAAAATTACAGGGATGCAGAATGAGACACTTTATCTATCAGGACGAAAAATcac                                                     <  2:859317/150‑1 (MQ=255)
                                                                                                  cAGCGCCTTTTCTTTATAAATTCCTAAAGTTGTTTTCTTGCGATTTTGTCTCTCTCTAACCCGCATAAATACTGGTAGCATCTGCATTCAACTGGATAAAATTACAGGGATGCAGAATGAGACACTTTATCTATCAGGACGAAAAATcac                                                     >  2:1795420/1‑150 (MQ=255)
                                                                                                    gcgcCTTTTCTTTATAAATTCCTAAAGTTGTTTTCTTGCGATTTTGTCTATCTCTAACCCGCATAAATACTGGTAGCATCTGCATTCAACTGGATAAAATTACAGGGATGCAGAATGAGACACTTTATCTATCAGGACGAAAAATCACAt                                                   <  1:186951/150‑1 (MQ=255)
                                                                                                     cgcCTTTTCTTTATAAATTCCTAAAGTTGTTTTCTTGCGATTTTGTCTCTCTCTAACCCGCATAAATACTGGTAGCATCTGCATTCAACTGGATAAAATTACAGGGATGCAGAATGAGACACTTTATCTATCAGGACGAAAAATCACATa                                                  >  2:2925383/1‑150 (MQ=255)
                                                                                                      gcCTTTTCTTTATAAATTCCTAAAGTTGTTTTCTTGCGATTTTGTCTCTCTCTAACCCGCATAAATACTGGTAGCATCTGCATTCAACTGGATAAAATTACAGGGATGCAGAATGAGACACTTTATCTATCAGGACGAAAAATCACATaa                                                 >  1:1596815/1‑150 (MQ=255)
                                                                                                       ccTTTTCTTTATAAATTCCTAAAGTTGTTTTCTTGCGATTTTGTCTCTCTCTAACCCGCATAAATACTGGTAGCATCTGCATTCAACTGGATAAAATTACAGGGATGCAGAATGAGACACTTTATCTATCAGGACGAAAAATCACATaaa                                                >  1:2513743/1‑150 (MQ=255)
                                                                                                         ttttCTTTATAAATTCCTAAAGTTGTTTTCTTGCGATTTTGTCTATCTCTAACCCGCATAAATACTGGTAGCATCTGCATTCAACTGGATAAAATTACAGGGATGCAGAATGAGACACTTTATCTATCAGGACGAAAAATCACATAAAtt                                              <  1:709685/150‑1 (MQ=255)
                                                                                                          tttCTTTATAAATTCCTAAAGTTGTTTTCTTGCGATTTTGTCTCTCTCTAACCCGCATAAATACTGGTAGCATCTGCATTCAACTGGATAAAATTACAGGGATGCAGAATGAGACACTTTATCTATCAGGACGAAAAATCACATAAATTc                                             <  1:1098214/150‑1 (MQ=255)
                                                                                                           ttCTTTATAAATTCCTAAAGTTGTTTTCTTGCGATTTTGTCTCTCTCTAACCCGCATAAATACTGGTAGCATCTGCATTCAACTGGATAAAATTACAGGGATGCAGAATGAGACACTTTATCTATCAGGACGAAAAATCACATAAATTCa                                            >  1:3004814/1‑150 (MQ=255)
                                                                                                              tttATAAATTCCTAAAGTTGTTTTCTTGCGATTTTGTCTCTCTCTAACCCGCATAAATACTGGTAGCATCTGCATTCAACTGGATAAAATTACAGGGATGCAGAATGAGACACTTTATCTATCAGGACGAAAAATCACATAAATTCAggg                                         <  2:3697220/150‑1 (MQ=255)
                                                                                                                tataAAATCCTAAAGTTGTTTTCTTGCGATTTTGTCTCTCTCTAACCCGCATAAATACTGGTAGCATCTGCATTCAACTGGATAAAATTACAGGGATGCAGAATGAGACACTTTATCTATCAGGACGAAAAATCACATAAATTCAGGGCa                                       >  2:4248626/1‑150 (MQ=255)
                                                                                                                 ataAATTCCTAAAGTTGTTTTCTTGCGATTTTGTCTCTCTCTAACCCGCATAAATACTGGTAGCATCTGCATTCAACTGGATAAAATTACAGGGATGCAGAATGAGACACTTTATCTATCAGGACGAAAAATCACATAAATTAAGGGCAg                                      >  1:108425/1‑150 (MQ=255)
                                                                                                                  taAATTCCTAAAGTTGTTTTCTTGCGATTTTGTCTCTCTCTAACCCGCATAAATACTGGTAGCATCTGCATTCAACTGGATAAAATTACAGGGATGCAGAATGAGACACTTTATCTATCAGGACGAAAAATCACATAAATTCAGGGCAGt                                     <  2:2076509/150‑1 (MQ=255)
                                                                                                                  taAATTCCTAAAGTTGTTTTCTTGCGATTTTGTCTCTCTCTAACCCGCATAAATACTGGTAGCATCTGCATTCAACTGGATAAAATTACAGGGATGCAGAATGAGACACTTTATCTATCAGGACGAAAAATCACATAAATTCAGGGCAGt                                     >  1:4291573/1‑150 (MQ=255)
                                                                                                                   aaaTTCCTAAAGTTGTTTTCTTGCGATTTTGTCTCTCTCTAACCCGCATAAATACTGGTAGCATCTGCATTCAACTGGATAAAATTACAGGGATGCAGAATGAGACACTTTATCTATCAGGACGAAAAATCACATAAATTCAGGGCAGtt                                    <  1:2378931/150‑1 (MQ=255)
                                                                                                                      ttCCTAAAGTTGTTTTCTTGCGATTTTGTCTCTCTCTAACCCGCATAAATACTGGTAGCATCTGCATTCAACTGGATAAAATTACAGGGATGCAGAATGAGACACTTTATCTATCAGGACGAAAAATCACATAAATTCAGGGCAGTTGAg                                 <  2:46565/150‑1 (MQ=255)
                                                                                                                       tCCTAAAGTTGTTTTCTTGCGATTTTGTCTATCTCTAACCCGCATAAAAACTGGTAGCATCTGCATTCAACTGGATAAAATTACATGGATGCAGAATGAGACACATTATCTATCAGGACGAAAAATCACAAAAATTCAGGGCAGTTGAGc                                >  1:757744/1‑150 (MQ=255)
                                                                                                                         cTAAAGTTGTTTTCTTGCGATTTTGTCTCTCTCTAACCCGCATAAAAAATGGTAGCATCAGCATTCAACTGGATAAAATTACAGGGATGCAGAATGAGACACATTATCTATCAGGACGAAAAATCACATAAATTCAGGGCCGTTGAGCaa                              >  1:3354244/1‑150 (MQ=255)
                                                                                                                           aaaGTTGTTTTCTTGCGATTTTGTCTCTCTCTAACCCGCATAAATACTGGTAGCATCTGCATTCAACTGGATAAAATTACAGGGATGCAGAATGAGACACTTTATCTATCAGGACGAAAAATCACATAAATTCAGGGCAGTTGAGCAACa                            <  2:813285/150‑1 (MQ=255)
                                                                                                                            aaGTTGTTTTCTTGCGATTTTGTCTCTCTCTAACCCGCATAAATACTGGTAGCATCTGCATTCAACTGGATAAAATTACAGGGATGCAGAATGAGACACTTTATCTATCAGGACGAAAAATCACATAAATTCAGGGCAGTTGAGCAACAg                           >  1:2068508/1‑150 (MQ=255)
                                                                                                                            aaGTTGTTTTCTTGCGATTTTGTCTATCTCTAACCCGCATAAATACTGGTAGCATCTGCATTCAACTGGATAAAATTACAGGGATGCAGAATGAGACACTTTATCTATCAGGACGAAAAATCACATAAATTCAGGGCAGTTGAGCAACAg                           <  2:652722/150‑1 (MQ=255)
                                                                                                                             aGTTGTTTTCTTGCGATTTTGTCTCTCTCTAACCCGCATAAATACTGGTAGCATCTGCATTCAACTGGATAAAATTACAGGGATGCAGAATGAGACACTTTATCTATCAGGACGAAAAATCACATAAATTCAGGGCAGTTGAGCAACAgg                          >  2:4121389/1‑150 (MQ=255)
                                                                                                                              gttgttTTCTTGCGATTTTGTCTCTCTCTAACCCGCATAAATACTGGTAGCATCTGCATTCAACTGGATAAAATTACAGGGATGCAGAATGAGACACTTTATCTATCAGGACGAAAAATCACATAAATTCAGGGCAGTTGAGCAACAggg                         <  2:853359/150‑1 (MQ=255)
                                                                                                                                tgttTTCTTGCGATTTTGTCTCTCTCTAACCCGCATAAATACTGGTAGCATCTGCATTCAACTGGATAAAATTACAGGGATGCAGAATGAGACACTTTATCTATCAGGACGAAAAATCACATAAATTCAGGGCAGTTGAGCAACAGGGaa                       <  2:1419748/150‑1 (MQ=255)
                                                                                                                                tgttTTCTTGCGATTTTGTCTCTCTCTAACCCGCATAAATACTGGTAGCATCTGCATTCAACTGGATAAAATTACAGGGATGCAGAATGAGACACTTTATCTATCAGGACGAAAAATCACATAAATTCAGGGCAGTTGAGCAACAGGGaa                       <  2:2366321/150‑1 (MQ=255)
                                                                                                                                  ttttCTTGCGATTTTGTCTATCTCTAACCCGCATAAATACTGGTAGCATCTGCATTCAACTGGATAAAATTACAGGGATGCAGAATGAGACACTTTATCTATCAGGACGAAAAATCACATAAATTCAGGGCAGTTGAGCAACAGGGAAAc                     >  1:1425949/1‑150 (MQ=255)
                                                                                                                                   tttCTTGCGATTTTGTCTCTCTCTAACCCGCATAAATACTGGTAGCATCTGCATTCAACTGGATAAAATTACAGGGATGCAGAATGAGACACTTTATCTATCAGGACGAAAAATCACATAAATTCAGGGCAGTTGAGCAACAGGGAAACg                    >  2:2025250/1‑150 (MQ=255)
                                                                                                                                    ttCTTGCGATTTTGTCTCTCTCTAACCCGCATAAATACTGGTAGCATCTGCATTCAACTGGATAAAATTACAGGGATGCAGAATGAGACACTTTATCTATCAGGACGAAAAATCACATAAATTCAGGGCAGTTGAGCAACAGGGAAACGa                   <  2:2513743/150‑1 (MQ=255)
                                                                                                                                     tCTTGCGATTTTGTCTCTCTCTAACCCGCATAAATACTGGTAGCATCTGCATTCAACTGGATAAAATTACAGGGATGCAGAATGAGACACTTTATCTATCAGGACGAAAAATCACATAAATTCAGGGCAGTTGAGCAACAGGGAAACGAg                  <  1:1567628/150‑1 (MQ=255)
                                                                                                                                     tCTTGCGATTTTGTCTCTCTCTAACCCGCATAAATACTGGTAGCATCTGCATTCAACTGGATAAAATTACAGGGATGCAGAATGAGACACTTTATCTATCAGGACGAAAAATCACATAAATTCAGGGCAGTTGAGCAACAGGGAAACGAg                  <  1:1567602/150‑1 (MQ=255)
                                                                                                                                     tCTTGCGATTTTGTCTATCTCTAACCCGCATAAATACTGGTAGCATCTGCATTCAACTGGATAAAATTACAGGGATGCAGAATGAGACACTTTATCTATCAGGACGAAAAATCACATAAATTCAGGGCAGTTGAGCAACAGGGACACGAg                  >  1:1192924/1‑150 (MQ=255)
                                                                                                                                      cTTGCGATTTTGTCTATCTCTAACCCGCATAAATACTGGTAGCATCTGCATTCAACTGGATAAAATTACAGGGATGCAGAATGAGACACTTTATCTATCAGGACGAAAAATCACATAAATTCAGGGCAGTTGAGCAACAGGGAAACGAGt                 >  2:2132383/1‑150 (MQ=255)
                                                                                                                                           gATTTTGTCTATCTCTAACCCGCATAAATACTGGTAGCATCTGCATTCAACTGGATAAAATTACAGGGATGCAGAATGAGACACTTTATCTATCAGGACGAAAAATCACATAAATTCAGGGCAGTTGAGCAACAGGGAAACGAGTTGCat            <  2:4358326/150‑1 (MQ=255)
                                                                                                                                           gATTTTGTCTATCTCTAACCCGCATAAATACTGGTAGCATCTGCATTCAACTGGATAAAATTACAGGGATGCAGAATGAGACACTTTATCTATCAGGACGAAAAATCACATAAATTCAGGGCAGTTGAGCAACAGGGAAACGAGTTGCat            <  2:3931867/150‑1 (MQ=255)
                                                                                                                                            aTTTTGTCTCTCTCTAACCCGCATAAATACTGGTAGCATCTGCATTCAACTGGATAAAATTACAGGGATGCAGAATGAGACACTTTATCTATCAGGACGAAAAATCACATAAATTCAGGGCAGTTGAGCAACAGGGAAACGAGTTGCata           >  2:3615059/1‑150 (MQ=255)
                                                                                                                                               ttGTCTATCTCTAACCCGCATAAATACTGGTAGCATCTGCATTCAACTGGATAAAATTACAGGGATGCAGAATGAGACACTTTATCTATCAGGACGAAAAATCACATAAATTCAGGGCAGTTGAGCAACAGGGAAACGAGTTGCATATCa        <  2:1638238/150‑1 (MQ=255)
                                                                                                                                                tGTCTCTCTCTAACCCGCATAAATACTGGTAGCATCTGCATTCAACTGGATAAAATTACAGGGATGCAGAATGAGACACTTTATCTATCAGGACGAAAAATCACATAAATTCAGGGCAGTTGAGCAACAGGGAAACGAGTTGCATATCAg       >  2:1281578/1‑150 (MQ=255)
                                                                                                                                                 gTCTCTCTCTAACCCGCATAAATACTGGTAGCATCTGCATTCAACTGGATAAAATTACAGGGATGCAGAATGAGACACTTTATCTATCAGGACGAAAAATCACATAAATTCAGGGCAGTTGAGCAACAGGGAAACGAGTTGCATATCAGt      >  1:2469890/1‑150 (MQ=255)
                                                                                                                                                   ctctctctAACCCGCATAAATACTGGTAGCATCTGCATTCAACTGGATAAAATTACAGGGATGCAGAATGAGACACTTTATCTATCAGGACGAAAAATCACATAAATTCAGGGCAGTTGAGCAACAGGGAAACGAGTTGCATATCAGTTg    <  1:200684/150‑1 (MQ=255)
                                                                                                                                                   ctctctctAACCCGCATAAATACTGGTAGCATCTGCATTCAACTGGATAAAATTACAGGGATGCAGAATGAGACACTTTATCTATCAGGACGAAAAATCACATAAATTCAGGGCAGTTGAGCAACAGGGAAACGAGTTGCATATCAGTTg    <  1:2178161/150‑1 (MQ=255)
                                                                                                                                                   ctctctctAACCCGCATAAATACTGGTAGCATCTGCATTCAACTGGATAAAATTACAGGGATGCAGAATGAGACACTTTATCTATCAGGACGAAAAATCACATAAATTCAGGGCAGTTGAGCAACAGGGAAACGAGTTGCATATCAGTTg    <  1:1018751/150‑1 (MQ=255)
                                                                                                                                                     ctctctAACCCGCATAAATACTGGTAGCATCTGCATTCAACTGGATAAAATTACAGGGATGCAGAATGAGACACTTTATCTATCAGGACGAAAAATCACATAAATTCAGGGCAGTTGAGCAACAGGGAAACGAGTTGCATATCAGTTggg  >  1:1319584/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CCGGTCCTGGCGGGAAAGATATTTTCCTGCTAAGCCCGGATGCCGGTGCTAAACCGGGTCATCAGGTGAAATAATCCCCCTTCAAGGCGCTGCATCGACAGCGCCTTTTCTTTATAAATTCCTAAAGTTGTTTTCTTGCGATTTTGTCTCTCTCTAACCCGCATAAATACTGGTAGCATCTGCATTCAACTGGATAAAATTACAGGGATGCAGAATGAGACACTTTATCTATCAGGACGAAAAATCACATAAATTCAGGGCAGTTGAGCAACAGGGAAACGAGTTGCATATCAGTTGGG  >  CP009273/2189739‑2190037

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: