Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,189,888 C→A 7.4% 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,8880CA 7.4% 798.7 / 30.9 216intergenic (+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 (97/103);  new base A (8/8);  total (105/111)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.46e-01

Only 100 of 229 total aligned reads displayed.

CGGTCCTGGCGGGAAAGATATTTTCCTGCTAAGCCCGGATGCCGGTGCTAAACCGGGTCATCAGGTGAAATAATCCCCCTTCAAGGCGCTGCATCGACAGCGCCTTTTCTTTATAAATTCCTAAAGTTGTTTTCTTGCGATTTTGTCTCTCTCTAACCCGCATAAATACTGGTAGCATCTGCATTCAACTGGATAAAATTACAGGGATGCAGAATGAGACACTTTATCTATCAGGACGAAAAATCACATAAATTCAGGGCAGTTGAGCAACAGGGAAACGAGTTGCATATCAGTTG  >  CP009273/2189740‑2190035
                                                                                                                                                    |                                                                                                                                                   
   tCCTGGCGGGAAAGATATTTTCCTGCTAAGCCCGGATGCCGGTGCTAAACCGGGTCATCAGGTGAAACAATCCCCCTTCAAGGCGCTGCATCGACAGCGCCTTTTCTTTATAAATTCCTAAAGTTGTTTTCTTGCGATTTTGtctctctc                                                                                                                                                 <  2:4256005/150‑1 (MQ=255)
    ccTGGCGGGAAAGATATTTTCCTGCTAAGCCCGGATGCCGGTGCTAAACCGGGTCATCAGGTGAAATAATCCCCCTTCAAGGCGCTGCATCGACAGCGCCTTTTCTTTATAAATTCCTAAAGTTGTTTTCTTGCGATTTTGTctatctct                                                                                                                                                <  2:472040/150‑1 (MQ=255)
       ggCGGGAAAGATATTTTCCTGCTAAGCCCGGATGCCGGTGCTAAACCGGGTCATCAGGTGAAATAATCCCCCTTCAAGGCGCTGCATCGACAGCGCCTTTTCTTTATAAATTCCTAAAGTTGTTTTCTTGCGATTTTGTCTCTCTCTAAc                                                                                                                                             >  2:396809/1‑150 (MQ=255)
        gCGGGAAAGATATTTTCCTGCTAAGCCCGGATGCCGGTGCTAAACCGGGTCATCAGGTGAAATAATCCCCCTTCAAGGCGCTGCATCGACAGCGCCTTTTCTTTATAAATTCCTAAAGTTGTTTTCTTGCGATTTTGTCTCTCTCTAAcc                                                                                                                                            >  2:321911/1‑150 (MQ=255)
         cGGGAAAGATATTTTCCTGCTAAGCCCGGATGCCGGTGCTAAACCGGGTCATCAGGTGAAATAATCCCCCTTCAAGGCGCTGCATCGACAGCGCCTTTTCTTTATAAATTCCTAAAGTTGTTTTCTTGCGATTTTGTCTCTCTCTAAccc                                                                                                                                           >  2:83121/1‑150 (MQ=255)
             aaaGATATTTTCCTGCTAAGCCCGGATGCCGGTGCTAAACCGGGTCATCAGGTGAAATAATCCCCCTTCAAGGCGCTGCATCGACAGCGCCTTTTCTTTATAAATTCCTAAAGTTGTTTTCTTGCGATTTTGTCTCTCTCTAACCCGCAt                                                                                                                                       >  1:3603942/1‑150 (MQ=255)
               aGATATTTTCCTGCTAAGCCCGGATGCCGGTGCTAAACCGGGTCATCAGGTGAAATAATCCCCCTTCAAGGCGCTGCATCGACAGCGCCTTTTCTTTATAAATTCCTAAAGTTGTTTTCTTGCGATTTTGTCTCTCTCTAACCCGCATaa                                                                                                                                     >  1:2160909/1‑150 (MQ=255)
                gATATTTTCCTGCTAAGCCCGGATGCCGGTGCTAAACCGGGTCATCAGGTGAAATAATCCCCCTTCAAGGCGCTGCATCGACAGCGCCTTTTCTTTATAAATTCCTAAAGTTGTTTTCTTGCGATTTTGTCTCTCTCTAACCCGCATaaa                                                                                                                                    >  2:4382867/1‑150 (MQ=255)
                  tatTTTCCTGCTAAGCCCGGATGCCGGTGCTAAACCGGGTCATCAGGTGAAATAATCCCCCTTCAAGGCGCTGCATCGACAGCGCCTTTTCTTTATAAATTCCTAAAGTTGTTTTCTTGCGATTTTGTCTCTCTCTAACCCGCATAAATa                                                                                                                                  <  2:499965/150‑1 (MQ=255)
                  tatTTTCCTGCTAAGCCCGGATGCCGGTGCTAAACCGGGTCATCAGGTGAAATAATCCCCCTTCAAGGCGCTGCATCGACAGCGCCTTTTCTTTATAAATTCCTAAAGTTGTTTTCTTGCGATTTTGTCTCTCTCAAACCCGCATAAATa                                                                                                                                  <  2:1372231/150‑1 (MQ=255)
                   atTTTCCTGCTAAGCCCGGATGCCGGTGCTAAACCGGGTCATCAGGTGAAATAATCCCCCTTCAAGGCGCTGCATCGACAGCGCCTTTTCTTTATAAATTCCTAAAGTTGTTTTCTTGCGATTTTGTCTCTCTCTAACCCGCATAAATAc                                                                                                                                 >  1:2692125/1‑150 (MQ=255)
                       tCCTGCTAAGCCCGGATGCCGGTGCTAAACCGGGTCATCAGGTGAAATAATCCCCCTTCAAGGCGCTGCATCGACAGCGCCTTTTCTTTATAAATTCCTAAAGTTGTTTTCTTGCGATTTTGTCTCTCTCTAACCCGCATAAATACTGGt                                                                                                                             <  1:3358870/150‑1 (MQ=255)
                       tCCTGCTAAGCCCGGATGCCGGTGCTAAACCGGGTCATCAGGTGAAATAATCCCCCTTCAAGGCGCTGCATCGACAGCGCCTTTTCTTTATAAATTCCTAAAGTTGTTTTCTTGCGATTTTGTCTCTCTCTAACCCGCATAAATACTGGt                                                                                                                             <  1:3433585/150‑1 (MQ=255)
                        ccTGCTAAGCCCGGATGCCGGTGCTAAACCGGGTCATCAGGTGAAATAATCCCCCTTCAAGGCGCTGCATCGACAGCGCCTTTTCTTTATAAATTCCTAAAGTTGTTTTCTTGCGATTTTGTCTCTCTCTAACCCGCATAAATACTGGTa                                                                                                                            <  2:4334199/150‑1 (MQ=255)
                             tAAGCCCGGATGCCGGTGCTAAACCGGGTCATCAGGTGAAATAATCCCCCTTCAAGGCGCTGCATCGACAGCGCCTTTTCTTTATAAATTCCTAAAGTTGTTTTCTTGCGATTTTGTCTCTCTCTAACCCGCATAAATACTGGTAGCATc                                                                                                                       <  2:2440556/150‑1 (MQ=255)
                             tAAGCCCGGATGCCGGTGCTAAACCGGGTCATCAGGTGAAATAATCCCCCTTCAAGGCGCTGCATCGACAGCGCCTTTTCTTTATAAATTCCTAAAGTTGTTTTCTTGCGATTTTGTCTCTCTCTAACCCGCATAAATACTGGTAGCATc                                                                                                                       <  2:2440524/150‑1 (MQ=255)
                             tAAGCCCGGATGCCGGTGCTAAACCGGGTCATCAGGTGAAATAATCCCCCTTCAAGGCGCTGCATCGACAGCGCCTTTTCTTTATAAATTCCTAAAGTTGTTTTCTTGCGATTTTGTCTCTCTCTAACCCGCATAAATACTGGTAGCATc                                                                                                                       <  1:1021776/150‑1 (MQ=255)
                              aaGCCCGGATGCCGGTGCTAAACCGGGTCATCAGGTGAAATAATCCCCCTTCAAGGCGCTGCATCGACAGCGCCTTTTCTTTATAAATTCCTAAAGTTGTTTTCTTGCGATTTTGTCTCTCTCTAACCCGCATAAATACTGGTAGCATCt                                                                                                                      >  2:787978/1‑150 (MQ=255)
                               aGCCCGGATGCCGGTGCTAAACCGGGTCATCAGGTGAAATAATCCCCCTTCAAGGCGCTGCATCGACAGCGCCTTTTCTTTATAAATTCCTAAAGTTGTTTTCTTGCGATTTTGTCTCTCTCTAACCCGCATAAATACTGGTAGCATCTg                                                                                                                     >  2:2592604/1‑150 (MQ=255)
                                gCCCGGATGCCGGTGCTAAACCGGGTCATCAGGTGAAATAATCCCCCTTCAAGGCGCTGCATCGACAGCGCCTTTTCTTTATAAATTCCTAAAGTTGTTTTCTTGCGATTTTGTCTATCTCTAACCCGCATAAATACTGGTAGCATCTGc                                                                                                                    <  2:3527085/150‑1 (MQ=255)
                                    ggATGCCGGTGCTAAACCGGGTCATCAGGTGAAATAATCCCCCTTCAAGGCGCTGCATCGACAGCGCCTTTTCTTTATAAATTCCTCAAGTTGTTTTCTTGCGATTTTGTCTCTCTCTAACCCGCATAAATACTGGTAGCATCTGCATTc                                                                                                                >  2:211146/1‑150 (MQ=255)
                                    ggATGCCGGTGCTAAACCGGGTCATCAGGTGAAATAATCCCCCTTCAAGGCGCTGCATCGACAGCGCCTTTTCTTTATAAATTCCTAAAGTTGTTTTCTTGCGATTTTGTCTCTCTCTAACCCGCATAAATACTGGTAGCATCTGCATTc                                                                                                                >  1:1000173/1‑150 (MQ=255)
                                    ggATGCCGGTGCTAAACCGGGTCATCAGGTGAAATAATCCCCCTTCAAGGCGCTGCATCGACAGCGCCTTTTCTTTATAAATTCCTAAAGTTGTTTTCTTGCGATTTTGTCTCTCTCTAACCCGCATAAATACTGGTAGCATCTGCATTc                                                                                                                <  1:425748/150‑1 (MQ=255)
                                    ggATGCCGGTGCTAAACCGGGTCATCAGGTGAAATAATCCCCCTTCAAGGCGCTGCATCGACAGCGCCTTTTCTTTATAAATTCCTAAAGTTGTTTTCTTGCGATTTTGTCTATCTCTAACCCGCATAAATACTGGTAGCATCTGCATTc                                                                                                                <  2:382965/150‑1 (MQ=255)
                                       tGCCGGTGCTAAACCGGGTCATCAGGTGAAATAATCCCCCTTCAAGGCGCTGCATCGACAGCGCCTTTTCTTTATAAATTCCTAAAGTTGTTTTCTTGCGATTTTGTCTCTCTCTAACCCGCATAAATACTGGTAGCATCTGCATTCAAc                                                                                                             <  2:3305246/150‑1 (MQ=255)
                                         ccGGTGCTAAACCGGGTCATCAGGTGAAATAATCCCCCTTCAAGGCGCTGCATCGACAGCGCCTTTTCTTTATAAATTCCTAAAGTTGTTTTCTTGCGATTTTGTCTCTCTCTAACCCGCATAAATACTGGTAGCATCTGCATTCAACTg                                                                                                           <  2:1233048/150‑1 (MQ=255)
                                              gCTAAACCGGGTCATCAGGTGAAATAATCCCCCTTCAAGGCGCTGCATCGACAGCGCCTTTTCTTTATAAATTCCTAAAGTTGTTTTCTTGCGATTTTGTCTCTCTCTAACCCGCATAAATACTGGTAGCATCTGCATTCAACTGGATaa                                                                                                      >  1:1298746/1‑150 (MQ=255)
                                               cTAAACCGGGTCATCAGGTGAAATAATCCCCCTTCAAGGCGCTGCATCGACAGCGCCTTTTCTTTATAAATTCCTAAAGTTGTTTTCTTGCGATTTTGTCTCTCTCTAACCCGCATAAATACTGGTAGCATCTGCATTCAACTGGATaaa                                                                                                     >  2:1709895/1‑150 (MQ=255)
                                                tAAACCGGGTCATCAGGTGAAATAATCCCCCTTCAAGGCGCTGCATCGACAGCGCCTTTTCTTTATAAATTCCTAAAGTTGTTTTCTTGCGATTTTGTCTCTCTCTAACCCGCATAAATACTGGTAGCATCTGCGTTCAACTGGATaaaa                                                                                                    >  2:3412017/1‑150 (MQ=255)
                                                tAAACCGGGTCATCAGGTGAAATAATCCCCCTTCAAGGCGCTGCATCGACAGCGCCTTTTCTTTATAAATTCCTAAAGTTGTTTTCTTGCGATTTTGTCTCTCTCTAACCCGCATAAATACTGGTAGCATCTGCATTCAACTGGATaaaa                                                                                                    >  2:1011779/1‑150 (MQ=255)
                                                tAAACCGGGTCATCAGGTGAAATAATCCCCCTTCAAGGCGCTGCATCGACAGCGCCTTTTCTTTATAAATTCCTAAAGTTGTTTTCTTGCGATTTTGTCTCTCTCTAACCCGCATAAATACTGGTAGCATCTGCATTCAACTGGATaaaa                                                                                                    <  1:715705/150‑1 (MQ=255)
                                                tAAACCGGGTCATCAGGTGAAATAATCCCCCTTCAAGGCGCTGCATCGACAGCGCCTTTTCTTTATAAATTCCTAAAGTTGTTTTCTTGCGATTTTGTCTCTCTCTAACCCGCATAAATACTGGTAGCATCTGCATTCAACTGGATaaaa                                                                                                    >  2:3567864/1‑150 (MQ=255)
                                                tAAACCGGGTCATCAGGTGAAATAATCCCCCTTCAAGGCGCTGCATCGACAGCGCCTTTTATTTATAAATTCCTAAAGTTGTTTTCTTGCGATTTTGTCTCTCTCTAACCCGCATAAATACTGGTAGCATCTGCATTCAACTGGATaaaa                                                                                                    <  1:3106113/150‑1 (MQ=255)
                                                 aaaCCGGGTCATCAGGTGAAATAATCCCCCTTCAAGGCGCTGCATCGACAGCGCCTTTTCTTTATAAATTCCTAAAGTTGTTTTCTTGCGATTTTGTCTCTCTCTAACCCGCATAAATACTGGTAGCATCTGCATTCAACTGGATAAAAt                                                                                                   <  1:3997674/150‑1 (MQ=255)
                                                 aaaCCGGGTCATCAGGTGAAATAATCCCCCTTCAAGGCGCTGCATCGACAGCGCCTTTTCTTTATAAATTCCTAAAGTTGTTTTCTTGCGATTTTGTCTCTCTCTAACCCGCATAAATACTGGTAGCATCTGCATTCAACTGGATAAAAt                                                                                                   <  2:3953916/150‑1 (MQ=255)
                                                  aaCCGGGTCATCAGGTGAAATAATCCCCCTTCAAGGCGCTGCATCGACAGCGCCTTTTCTTTATAAATTCCTAAAGTTGTTTTCTTGCGATTTTGTCTCTCTCTAACCCGCATAAATACTGGTAGCATCTGCATTCAACTGGATAAAAtt                                                                                                  <  1:1962208/150‑1 (MQ=255)
                                                      gggTCATCAGGTGAAATAATCCCCCTTCAAGGCGCTGCATCGACAGCGCCTTTTCTTTATAAATTCCTAAAGTTGTTTTCTTGCGATTTTGTCTCTCTCTAACCCGCATAAATACTGGTAGCATCTGCATTCAACTGGATAAAATTACAg                                                                                              <  1:3791824/150‑1 (MQ=255)
                                                        gTCATCAGGTGAAATAATCCCCCTTCAAGGCGCTGCATCGACAGCGCCTTTTCTTTATAAATTCCTAAAGTTGTTTTCTTGCGATTTTGTCTCTCTCTAACCCGCATAAATACTGGTAGCATCTGCATTCAACTGGATAAAATTACAggg                                                                                            <  2:2538418/150‑1 (MQ=255)
                                                        gTCATCAGGTGAAATAATCCCCCTTCAAGGCGCTGCATCGACAGCGCCTTTTCTTTATAAATTCCTAAAGTTGTTTTCTTGCGATTTTGTCTCTCTCTAACCCGCATAAATACTGGTAGCATCTGCATTCAACTGGATAAAATTACAggg                                                                                            <  2:1526840/150‑1 (MQ=255)
                                                            tcaGGTGAAATAATCCCCCTTCAAGGCGCTGCATCGACAGCGCCTTTTCTTTATAAATTCCTAAAGTTGTTTTCTTGCGATTTTGTCTCTCTCTAACCCGCATAAATACTGGTAGCATCTGCATTCAACTGGATAAAATTACAGGGATGc                                                                                        >  2:4020391/1‑150 (MQ=255)
                                                              aGGTGAAATAATCCCCCTTCAAGGCGCTGCATCGACAGCGCCTTTTCTTTATAAATTCCTAAAGTTGTTTTCTTGCGATTTTGTCTCTCTCTAACCCGCATAAATACTGGTAGCATCTGCATTCAACTGGATAAAATTACAGGGATGCAg                                                                                      <  2:4446013/150‑1 (MQ=255)
                                                                gTGAAATAATCCCCCTTCAAGGCGCTGCATCGACAGCGCCTTTTCTTTATAAATTCCTAAAGTTGTTTTCTTGCGATTTTGTCTCTCTCTAACCCGCATAAATACTGGTAGCATCTGCATTCAACTGGATAAAATTACAGGGATGCAGaa                                                                                    >  2:939878/1‑150 (MQ=255)
                                                                 tGAAATAATCCCCCTTCAAGGCGCTGCATCGACAGCGCCTTTTCGTTATAAATTCCTAAAGTTGTTTTCTTGCGATTTTGTCTCTCTCTAACCCGCATAAATACTGGTAGCATCTGCATTCAACTGGATAAAATTACAGGGATGCAGAAt                                                                                   >  2:3604668/1‑150 (MQ=255)
                                                                 tGAAATAATCCCCATTCATGGCTCTGCATCGACAGCGCCTTTTCTTTATAAATTCCTAAAGTTGTTTTCTTGCGATTTTGTCTCTCTCTAACCCGCATAAATACTGGTAGCTTCTTCATTCAACTGGATAAAATTACAGGGATGCAGAAt                                                                                   <  1:4444941/150‑1 (MQ=255)
                                                                   aaaTAATCCCCCGTCAAGGCGCTGCATCGACAGCGCCTTTTCTTTATAAATTCCTAAAGTTGTTTTCTTGCGATTTTGTCTCTCTCTAACCCGCATAAATACTGGTAGCATCTGCATTCAACTGGATAAAATTACAGGGATGCAGAATga                                                                                 >  2:889292/1‑150 (MQ=255)
                                                                     ataatCCCCCTTCAAGGCGCTGCATCGACAGCGCCTTTTCTTTATAAATTCCTAAAGTTGTTTTCTTGCGATTTTGTCTCTCTCTAACCCGCATAAATACTGGTAGCATCTGCATTCAACTGGATAAAATTACAGGGATGCAGAATgaga                                                                               >  2:3865757/1‑150 (MQ=255)
                                                                     ataatCCCCCTTCAAGGCGCTGCATCGACAGCGCCTTTTCTTTATAAATTCCTAAAGTTGTTTTCTTGCGATTTTGTCTCTCTCTAACCCGCATAAATACTGGTAGCATCTGCATTCAACTGGATAAAATTACAGGGATGCAGAATgaga                                                                               >  1:631058/1‑150 (MQ=255)
                                                                     ataatCCCCCTTCAAGGCGCTGCATCGACAGCGCCTTTTCTTTATAAATTCCTAAAGTTGTTTTCTTGCGATTTTGTCTCTCTCTAACCCGCATAAATACTGGTAGCATCTGCATTCAACTGGATAAAATTACAGGGATGCAGAATgaga                                                                               <  2:1920741/150‑1 (MQ=255)
                                                                        atCCCCCTTCAAGGCGCTGCATCGACAGCGCCTTTTCTTTATAAATTCCTAAAGTTGTTTTCTTGCGATTTTGTCTCTCTCTAACCCGCATAAATACTGGTAGCATCTGCATTCAACTGGATAAAATTACAGGGATGCAGAATGAGacac                                                                            >  1:1769838/1‑150 (MQ=255)
                                                                        atCCCCCTTCAAGGCGCTGCATCGACAGCGCCTTTTCTTTATAAATTCCTAAAGTTGTTTTCTTGCGATTTTGTCTCTCTCTAACCCGCATAAATACTGGTAGCATCTGCATTCAACTGGATAAAATTACAGGGATGCAGAATGAGacac                                                                            <  1:4359644/150‑1 (MQ=255)
                                                                         tCCCCCTTCAAGGCGCTGCATCGACAGCGCCTTTTCTTTATAAATTCCCATAGTTGTTTTCTTGCGATTTTGTCTCTCTCTAACCCGCATAAATACTGGTAGCATCTGCATTCAACTGGATAAAATTACAGGGATGCAGAATGAGACACt                                                                           <  1:2631213/150‑1 (MQ=255)
                                                                          cccccTTCAAGGCGCTGCATCGACAGCGCCTTTTCTTTATAAATTCCTAAAGTTGTTTTCTTGCGATTTTGTCTCTCTCTAACCCGCATAAATACTGGTAGCATCTGCATTCAACTGGATAAAATTACAGGGATGCAGAATGAGACACtt                                                                          <  1:396796/150‑1 (MQ=255)
                                                                          cccccTTCAAGGCGCTGCATCGACAGCGCCTTTTCTTTATAAATTCCTAAAGTTGTTTTCTTGCGATTTTGTCTCTCTCTAACCCGCATAAATACTGGTAGCATCTGCATTCAACTGGATAAAATTACAGGGATGCAGAATGAGACACtt                                                                          <  2:3383280/150‑1 (MQ=255)
                                                                            cccTTCAAGGCGCTGCATCGACAGCGCCTTTTCTTTATAAATTCCTAAAGTTGTTTTCTTGCGATTTTGTCTCTCTCTAACCCGCATAAATACTGGTAGCATCTGCATTCAACTGGATAAAATTACAGGGATGCAGAATGAGACACTTta                                                                        <  1:1723251/150‑1 (MQ=255)
                                                                             ccTTCAAGGCGCTGCATCGACAGCGCCTTTTCTTTATAAATTCCTAAAGTTGTTTTCTTGCGATTTTGTCTCTCTCTAACCCGCATAAATACTGGTAGCATCTGCATTCAACTGGATAAAATTACAGGGATGCAGAATGAGACACTTtat                                                                       <  2:2793383/150‑1 (MQ=255)
                                                                               ttCAAGGCGCTGCATCGACAGCGCCTTTTCTTTATAAATTCCTAAAGTTGTTTTCTTGCGATTTTGTCTCTCTCTAACCCGCATAAATACTGGTAGCATCTGCATTCAACTGGATAAAATTACAGGGATGCAGAATGAGACACTTtatct                                                                     <  2:33699/150‑1 (MQ=255)
                                                                               ttCAAGGCGCTGCATCGACAGCGCCTTTTCTTTATAAAGTCCTAAAGTTGTTTTCTTGCGATTTTGTCTCTCTCTAACCCGCATAAATACTGGTAGCATCTGCATTCAACTGGATAAAATTACAGGGATGCAGAATGAGACACTTtatct                                                                     <  2:32935/150‑1 (MQ=255)
                                                                                tCAAGGCGCTGCATCGACAGCGCCTTTTCTTTATAAATTCCTAAAGTTGTTTTCTTGCGATTTTGTCTCTCTCTAACCCGCATAAATACTGGTAGCATCTGCATTCAACTGGATAAAATTACAGGGATGCAGAATGAGACACTTtatcta                                                                    >  1:3789866/1‑150 (MQ=255)
                                                                                tCAAGGCGCTGCATCGACAGCGCCTTTTCTTTATAAATTCCTAAAGTTGTTTTCTTGCGATTTTGTCTCTCTCTAACCCGCATAAATACTGGTAGCATCTGCATTCAAATGGATAAAATTACAGGGATGCAGAATGAGACACTTtatcta                                                                    >  1:2138372/1‑150 (MQ=255)
                                                                                 cAAGGCGCTGCATCGACAGCGCCTTTTCTTTATAAATTCCTAAAGTTGTTTTCTTGCGATTTTGTCTCTCTCTAACCCGCATAAATACTGGTAGCATCTGCATTCAACTGGATAAAATTACAGGGATGCAGAATGAGACACTTtatctat                                                                   >  1:497078/1‑150 (MQ=255)
                                                                                   aGGCGCTGCATCGACAGCGCCTTTTCTTTATAAATTCCTAAAGTTGTTTTCTTGCGATTTTGTCTATCTCTAACCCGCATAAATACTGGTAGCATCTGCATTCAACTGGATAAAATTACAGGGATGCAGAATGAGACACTTTATCTATCa                                                                 >  1:1647513/1‑150 (MQ=255)
                                                                                    ggCGCTGCATCGACAGCGCCTTTTCTTTATAAATTCCTAAAGTTGTTTTCTTGCGATTTTGTCTATCTCTAACCCGCATAAATACTGGTAGCATCTGCATTCAACTGGATAAAATTACAGGGATGCAGAATGAGACACTTTATCTATCAg                                                                <  1:4566283/150‑1 (MQ=255)
                                                                                    ggCGCTGCATCGACAGCGCCTTTTCTTTATAAATTCCTAAAGTTGTTTTCTTGCGATTTTGTCTATCTCTAACCCGCATAAATACTGGTAGCATCTGCATTCAACTGGATAAAATTACAGGGATGCAGAATGAGACACTTTATCTATCAg                                                                >  2:4134884/1‑150 (MQ=255)
                                                                                    ggCGCTGCATCGACAGCGCCTTTTCTTTATAAATTCCTAAAGTTGTTTTCTTGCGATTTTGTCTATCTCTAACCCGCATAAATACTGGTAGCATCTGCATTCAACTGGATAAAATTACAGGGATGCAGAATGAGACACTTTATCTATCAg                                                                >  2:4134816/1‑150 (MQ=255)
                                                                                      cgcTGCATCGACAGCGCCTTTTCTTTATAAATTCCTAAAGTTGTTTTCTTGCGATTTTGTCTCTCTCTAACCCGCATAAATACTGGTAGCATCTGCATTCAACTGGATAAAATTACAGGGATGCAGAATGAGACACTTTATCTATCAGGa                                                              >  1:4682303/1‑150 (MQ=255)
                                                                                      cgcTGCATCGACAGCGCCTTTTCTTTATAAATTCCTAAAGTTGTTTTCTTGCGATTTTGTCTCTCTCTAACCCGCATAAATACTGGTAGCATCTGCATTCAACTGGATAAAATTACAGGGATGCAGAATGAGACACTTTATCTATCAGGa                                                              <  1:2771176/150‑1 (MQ=255)
                                                                                       gcTGCATCGACAGCGCCTTTTCTTTATAAATTCCTAAAGTTGTTTTCTTGCGATTTTGTCTCTCTCTAACCCGCATAAATACTGGTAGCATCTGCATTCAACTGGATAAAATTACAGGGATGCAGAATGAGACACTTTATCTATCAGGAc                                                             <  1:1716799/150‑1 (MQ=255)
                                                                                       gcTGCATCGACAGCGCCTTTTCTTTATAAATTCCTAAAGTTGTTTTCTTGCGATTTTGTCTCTCTCTAACCCGCATAAATACTGGTAGCATCTGCATTCAACTGGATAAAATTACAGGGATGCAGAATGAGACACTTTATCTATCAGGAc                                                             <  2:2132476/150‑1 (MQ=255)
                                                                                       gcTGCATCGACAGCGCCTTTTCTTTATAAATTCCTAAAGTTGTTTTCTTGCGATTTTGTCTCTCTCTAACCCGCATAAATACTGGTAGCATCTGCATTCAACTGGATAAAATTACAGGGATGCAGAATGAGACACTTTATCTATCAGGAc                                                             <  2:1952445/150‑1 (MQ=255)
                                                                                          gCATCGACAGCGCAATTACTTTATAAATTCCTAAAGTTGTTTTCTTGCGATTTTGTCTCTCTCTAACCCGCATAAATACTGGTAGCATCTGCATTCAACTGGATAAAATTACAGGGATGCAGAATGAGACACTTTATCTATCAGGACGaa                                                          <  2:3743203/150‑1 (MQ=255)
                                                                                           cATCGACAGCGCCTTTTCTTTATAAATTCCTAAAGTTGTTTTCTTGCGATTTTGTCTCTCTCTAACCCGCATAAATACTGGTAGCATCTGCATTCAACTGGATAAAATTACAGGGATGCAGAATGAGACACTTTATCTATCAGGACGaaa                                                         <  1:2038214/150‑1 (MQ=255)
                                                                                              cGACAGCGCCTTTTCTTTATAAATTCCTAAAGTTGTTTTCTTGCGATTTTGTCTCTCTCTAACCCGCATAAATACTGGTAGCATCTGCATTCAACTGGATAAAATTACAGGGATGCAGAATGAGACACTTTATCTATCAGGACGAAAAAt                                                      >  1:717670/1‑150 (MQ=255)
                                                                                               gACAGCGCCTTTTCTTTATAAATTCCTAAAGTTGTTTTCTTGCGATTTTGTCTCTCTCTAACCCGCATAAATACTGGTAGCATCTGCATTCAACTGGATAAAATTACAGGGATGCAGAATGAGACACTTTATCTATCAGGACGAAAAATc                                                     >  2:3020168/1‑150 (MQ=255)
                                                                                                  aGCGCCTTTTCTTTATAAATTCCTAAAGTTGTTTTCTTGCGATTTTGTCTCTCTCTAACCCGCATAAATACTGGTAGCATCTGCATTCAACTGGATAAAATTACAGGGATGCAGAATGAGACACTTTATCTATCAGGACGAAAAATcaca                                                  >  2:2717525/1‑150 (MQ=255)
                                                                                                    cgcCTTTTCTTTATAAATTCCTAAAGTTGTTTTCTTGCGATTTTGTCTCTCTCTAACCCGCATAAATACTGGTAGCATCTGCATTCAACTGGATAAAATTACAGGGATGCAGAATGAGACACTTTATCTATCAGGACGAAAAATCACATa                                                <  2:1001796/150‑1 (MQ=255)
                                                                                                    cgcCTTTTCTTTATAAATTCCTAAAGTTGTTTTCTTGCGATTTTGTCTATCTCTAACCCGCATAAATACTGGTAGCATCTGCATTCAACTGGATAAAATTACAGGGATGCAGAATGAGACACTTTATCTATCAGGACGAAAAATCACATa                                                <  1:3365618/150‑1 (MQ=255)
                                                                                                       cTTTTCTTTATAAATTCCTAAAGTTGTTTTCTTGCGATTTTGTCTCTCTCTAACCCGCATAAATACTGGTAGCATCTGCATTCAACTGGATAAAATTACAGGGATGCAGAATGAGACACTTTATCTATCAGGACGAAAAATCACATAAAt                                             >  2:1409490/1‑150 (MQ=255)
                                                                                                            cTTTATAAATTCCTAAAGTTGTTTTCTTGCGATTTTGTCTCTCTCTAACCCGCATAAATACTGGTAGCATCTGCATTCAACTGGATAAAATTACAGGGATGCAGAATGAGACACTTTATCTATCAGGACGAAAAATCATATAAATTCAgg                                        >  2:2002477/1‑150 (MQ=255)
                                                                                                              ttATAAATTCCTAAAGTTGTTTTCTTGCGATTTTGTCTCTCTCTAACCCGCATAAATACTGGTAGCATCTGCATTCAACTGGATAAAATTACAGGGATGCAGAATGAGACACTTTATCTATCAGGACGAAAAATCACATAAATTCAGGGc                                      >  1:2211533/1‑150 (MQ=255)
                                                                                                              ttATAAATTCCTAAAGTTGTTTTCTTGCGATTTTGTCTCTCTCTAACCCGCATAAATACTGGTAGCATCTGCATTCAACTGGATAAAATTACAGGGATGCAGAATGAGACACTTTATCTATCAGGACGAAAAATCACATAAATTCAGGGc                                      >  2:861418/1‑150 (MQ=255)
                                                                                                               tataAATTCCTAAAGTTGTTTTCTTGCGATTTTGTCTCTCTCTAACCCGCATAAATACTGGTAGCATCTGCATTCAACTGGATAAAATTACAGGGATGCAGAATGAGACACTTTATCTATCAGGACGAAAAATCACATAAATTCAGGGCa                                     >  1:780538/1‑150 (MQ=255)
                                                                                                                 taAATTCCTAAAGTTGTTTTCTTGCGATTTTGTCTCTCTCTAACCCGCATAAATACTGGTAGCATCTGCATTCAACTGGATAAAATTACAGGGATGCAGAATGAGACACTTTATCTATCAGGACGAAAAATCACATAAATTCAGGGCAGt                                   <  2:2596120/150‑1 (MQ=255)
                                                                                                                 taAATTCCTAAAGTTGTTTTCTTGCGATTTTGTCTCTCTCTAACCCGCATAAATACTGGTAGCATCTGCATTCAACTGGATAAAATTACAGGGATGCAGAATGAGACACTTTATCTATCAGGACGAAAAATCACATAAATTCAGGGCAGt                                   <  2:4279514/150‑1 (MQ=255)
                                                                                                                   aagtCCTAAAGTTGTTTTCTTGCGATTTTGTCTCTCTCTAACCCGCATAAATACTGGTAGCATCTGCATTCAACTGGATAAAATTACAGGGATGCAGAATGAGACACTTTATCTATCAGGACGAAAAATCACATAAATTCAGGGCAGTTg                                 <  2:4536540/147‑1 (MQ=255)
                                                                                                                          aaaGTTGTTTTCTTGCGATTTTGTCTCTCTCTAACCCGCATAAAAACTGGTAGCATCTGCATTCAACTGGATAAAATTAAAGGGATGCAGAATGAGACACTTTATCTATCAGGACGAAAAATCACATAAATTAAGGGCAGTTGAGCAACa                          >  1:3284388/1‑150 (MQ=255)
                                                                                                                              ttgttTTCTTGCGATTTTGTCTCTCTCTAACCCGCATAAATACTGGTAGCATCTGCATTCAACTGGATAAAATTACAGGGATGCAGAATGAGACACTTTATCTATCAGGACGAAAAATCACATAAATTCAGGGCAGTTGAGCAACAGGGa                      <  1:3210803/150‑1 (MQ=255)
                                                                                                                               tgttTTCTTGCGATTTTGTCTCTCTCTAACCCGCATAAATACTGGTAGCATCTGCATTCAACTGGATAAAATTACAGGGATGCAGAATGAGACACTTTATCTATCAGGACGAAAAATCACATAAATTCAGGGCAGTTGAGCAACAGGGaa                     <  1:1296309/150‑1 (MQ=255)
                                                                                                                                gttTTCTTGCGATTTTGTCTCTCTCTAACCCGCATAAATACTGGTAGCATCTGCATTCAACTGGATAAAATTACAGGGATGCAGAATGAGACACTTTATCTATCAGGACGAAAAATCACATAAATTCAGGGCAGTTGAGCAACAGGGaaa                    >  2:3460441/1‑150 (MQ=255)
                                                                                                                                gttTTCTTGCGATTTTGTCTCTCTCTAACCCGCATAAATACTGGTAGCATCTGCATTCAACTGGATAAAATTACAGGGATGCAGAATGAGACACTTTATCTATCAGGACGAAAAATCACATAAATTCAGGGCAGTTGAGCAACAGGGaaa                    >  1:3784670/1‑150 (MQ=255)
                                                                                                                                gttTTCTTGCGATTTTGTCTCTCTCTAACCCGCATAAATACTGGTAGCATCTGCATTCAACTGGATAAAATTACAGGGATGCAGAATGAGACACTTTATCTATCAGGACGAAAAATCACATAAATTCAGGGCAGTTGAGCAACAGGGaaa                    <  1:1513800/150‑1 (MQ=255)
                                                                                                                                gttTTCTTGCGATTTTGTCTCTCTCTAACCCGCATAAATACTGGTAGCATCTGCATTCAACTGGATAAAATTACAGGGATGCAGAATGAGACACTTTATCTATCAGGACGAAAAATCACATAAATTCAGGGCAGTTGAGCAACAGGGaaa                    >  2:2031097/1‑150 (MQ=255)
                                                                                                                                 ttttCTTGCGATTTTGTCTCTCTCTAACCCGCATAAATACTGGTAGCATCTGCATTCAACTGGATAAAATTACAGGGATGCAGAATGAGACACTTTATCTATCAGGACGAAAAATCACATAAATTCAGGGCAGTTGAGCAACAGGGAAAc                   >  2:687052/1‑150 (MQ=255)
                                                                                                                                   ttCTTGCGATTTTGTCTCTCTCTAACCCGCATAAATACTGGTAGCATCTGCATTCAACTGGATAAAATTACAGGGATGCAGAATGAGACACTTTATCTATCAGGACGAAAAATCACATAAATTCAGGGCAGTTGAGCAACAGGGAAACGa                 >  1:3220345/1‑150 (MQ=255)
                                                                                                                                   ttCTTGCGATTTTGTCTCTCTCTAACCCGCATAAATACTGGTAGCATCTGCATTCAACTGGATAAAATTACAGGGATGCAGAATGAGACACTTTATCTATCAGGACGAAAAATCACATAAATTCAGGGCAGTTGAGCAACAGGGAAACGa                 <  2:232985/150‑1 (MQ=255)
                                                                                                                                     cTTGCGATTTTGTCTCTCTCTAACCCGCATAAATACTGGTAGCATCTGCATTCAACTGGATAAAATTACAGGGATGCAGAATGAGACACTTTATCTATCAGGACGAAAAATCACATAAATTCAGGGCAGTTGAGCAACAGGGAAACGAGt               >  2:935892/1‑150 (MQ=255)
                                                                                                                                       tGCGATTTTGTCTCTCTCTAACCCGCATAAATACTGGTAGCATCTGCATTCAACTGGATAAAATTACAGGGATGCAGAATGAGACACTTTATCTATCAGGACGAAAAATCACATAAATTCAGGGCAGTTGAGCAACAGGGAAACGAGTTg             <  2:3206529/150‑1 (MQ=255)
                                                                                                                                          gATTTTGTCTCTCTCTAACCCGCATAAATACTGGTAGCATCTGCATTCAACTGGATAAAATTACAGGGATGCAGAATGAGACACTTTATCTATCAGGACGAAAAATCACATAAATTCAGGGCAGTTGAGCAACAGGGAAACGAGTTGCat          <  1:3736268/150‑1 (MQ=255)
                                                                                                                                          gATTTTGTCTCTCTCTAACCCGCATAAATACTGGTAGCATCTGCATTCAACTGGATAAAATTACAGGGATGCAGAATGAGACACTTTATCTATCAGGACGAAAAATCACATAAATTCAGGGCAGTTGAGCAACAGGGAAACGAGTTGCat          >  2:3715269/1‑150 (MQ=255)
                                                                                                                                            ttttGTCTCTCTCTAACCCGCATAAATACTGGTAGCATCTGCATTCAACTGGATAAAATTACAGGGATGCAGAATGAGACACTTTATCTATCAGGACGAAAAATCACATAAATTCAGGGCAGTTGAGCAACAGGGAAACGAGTTGCatat        >  1:2549437/1‑150 (MQ=255)
                                                                                                                                                 tctctctcTAACCCGCATAAATACTGGTAGCATCTGCATTCAACTGGATAAAATTACAGGGATGCAGAATGAGACACTTTATCTATCAGGACGAAAAATCACATAAATTCAGGGCAGTTGAGCAACAGGGAAACGAGTTGCATATCAGtt   <  2:1164890/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
CGGTCCTGGCGGGAAAGATATTTTCCTGCTAAGCCCGGATGCCGGTGCTAAACCGGGTCATCAGGTGAAATAATCCCCCTTCAAGGCGCTGCATCGACAGCGCCTTTTCTTTATAAATTCCTAAAGTTGTTTTCTTGCGATTTTGTCTCTCTCTAACCCGCATAAATACTGGTAGCATCTGCATTCAACTGGATAAAATTACAGGGATGCAGAATGAGACACTTTATCTATCAGGACGAAAAATCACATAAATTCAGGGCAGTTGAGCAACAGGGAAACGAGTTGCATATCAGTTG  >  CP009273/2189740‑2190035

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: