Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,346,015 G→T 41.9% L72I (CTC→ATC)  dtpC ← dipeptide and tripeptide permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,346,0150GT41.9% 260.7 / 280.0 238L72I (CTC→ATC) dtpCdipeptide and tripeptide permease
Reads supporting (aligned to +/- strand):  ref base G (73/65);  new base T (46/54);  total (119/119)
Fisher's exact test for biased strand distribution p-value = 3.58e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 239 total aligned reads displayed.

GCTGATGTTTGATTTGAATAAACCGTAGCCACAAATAATGATTGCCAGCGCCAGATACAGGCTAAAGGTTGAATTTGTATCAATACCCAGCACCACATGGCCAAGGGTCATTAACAGCGCGCCGGCAATCACTGCAGTGCGGTTGCCGAGCAGGCGGTCGGCAAGCCAGCCGCCGAGAATAGGGGTAACGTAAACCAGAGAAGCATATGCGCTGAACAGGCTGATGGCATGGTTATCATCAAAACCAAGCTGATGGGTGAGATAGAGAATGAGTAAGGCACGCATGCCGTAAAAACTG  >  CP009273/4345866‑4346163
                                                                                                                                                     |                                                                                                                                                    
  tGATGTTTGATTTGAATAAACCGTAGCCACAAATAATGATTGCCAGCGCCAGATACAGGCTAAAGGTTGAATTTGTATCAATACCCAGCACCACATGGCCAAGGGTCATTAACAGCGCGCCGGCAATCACTGCAGTGCGGTTGCCGATCa                                                                                                                                                    >  2:4147900/1‑150 (MQ=255)
  tGATGTTTGATTTGAATAAACCGTAGCCACAAATAATGATTGCCAGCGCCAGATACAGGCTAAAGGTTGAATTTGTATCAATACCCAGCACCACATGGCCAAGGGTCATTAACAGCGCGCCGGCAATCACTGCAGTGCGGTTGCCGATCa                                                                                                                                                    >  1:749531/1‑150 (MQ=255)
  tGATGTTTGATTTGAATAAACCGTAGCCACAAATAATGATTGCCAGCGCCAGATACAGGCTAAAGGTTGAATTTGTATCAATACCCAGCACCACATGGCCAAGGGTCATTAACAGCGCGCCGGCAATCACTGCAGTGCGGTTGCCGATCa                                                                                                                                                    >  2:3221170/1‑150 (MQ=255)
  tGATGTTTGATTTGAATAAACCGTAGCCACAAATAATGATTGCCAGCGCCAGATACAGGCTAAAGGTTGAATTTGTATCAATACCCAGCACCACATGGCCAAGGGTCATTAACAGCGCGCCGGCAATAACTGCAGTGCGGTTGCCGAGCa                                                                                                                                                    >  2:3839895/1‑150 (MQ=255)
   gATGTTTGATTTGAATAAACCGTAGCCACAAATAATGATTGCCAGCGCCAGATACAGGCTAAAGGTTGAATTTGTATCAATACCCAGCACCACATGGCCAAGGGTCATTAACAGCGCGCCGGCAATAACTGCAGTGCGGTTGCCGAGCAg                                                                                                                                                   <  1:4085693/150‑1 (MQ=255)
   gATGTTTGATTTGAATAAACCGTAGCCACAAATAATGATTGCCAGCGCCAGATACAGGCTAAAGGTTGAATTTGTATCAATACCCAGCACCACATGGCCAAGGGTCATTAACAGCGCGCCGGCAATAACTGCAGTGCGGTTGCCGAGCAg                                                                                                                                                   >  2:3006512/1‑150 (MQ=255)
   gATGTTTGATTTGAATAAACCGTAGCCACAAATAATGATTGCCAGCGCCAGATACAGGCTAAAGGTTGAATTTGTATCAATACCCAGCACCACATGGCCAAGGGTCATTAACAGCGCGCCGGCAATAACTGCAGTGCGGTTGCCGAGCAg                                                                                                                                                   <  1:2288690/150‑1 (MQ=255)
      gTTTGATTTGAATAAACCGTAGCCACAAATAATGATTGCCAGCGCCAGATACCGGCTAAAGGTTGAATTTGTATCAATACCCAGCACCACATGGCCAAGGGTCATTAACAGCGCGCCGGCAATCACTGCAGTGCGGTTGCCGAGCAGGCg                                                                                                                                                >  1:1366306/1‑150 (MQ=255)
      gTTTGATTTGAATAAACCGTAGCCACAAATAATGATTGCCAGCGCCAGATACAGGCTAAAGGTTGAATTTGTATCAATACCCAGCACCACATGGCCAAGGGTCATTAACAGCGCGCCGGCAATCACTGCAGTGCGGTTGCCGATCAGGCg                                                                                                                                                >  1:4156857/1‑150 (MQ=255)
        ttgatttgaATAAACCGTAGCCACAAATAATGATTGCCAGCGCCAGATACAGGCTAAAGGTTGAATTTGTATCAATACCCAGCACCACATGGCCAAGGGTCATTAACAGCGCGCCGGCAATCACTGCAGTGCGGTTGCCGAGCAGGCGGt                                                                                                                                              >  2:656106/1‑150 (MQ=255)
           atttgaATAAACCGTAGCCACAAATAATGATTGCCAGCGCCAGATACAGGCTAAAGGTTGAATTTGTATCAATACCCAGCACCACATGGCCAAGGGTCATTAACAGCGCGCCGGCAATCACTGCAGTGCGGTTGCCGAGCAGGCGGTCgg                                                                                                                                           <  1:4143726/150‑1 (MQ=255)
               gaATAAACCGTAGCCACAAATAATGATTGCCAGCGCCAGATACAGGCTAAAGGTTGAATTTGTATCAATACCCAGCACCACATGGCCAAGGGTCATTAACAGCGCGCCGGCAATAACTGCAGTGCGGTTGCCGAGCAGGCGGTCGGCAag                                                                                                                                       >  1:3303074/1‑150 (MQ=255)
                 aTAAACCGTAGCCACAACTAATGATTGCCAGCGCCAGATACAGGCTAAAGGTTGAATTTGTATCAATACCCAGCACCACATGGCCAAGGGTCATTAACAGCGCGCCGGCAATCACTGCAGTGCGGTTGCCGATCAGGCGGTCGGCAagcc                                                                                                                                     <  2:2660017/150‑1 (MQ=255)
                 aTAAACCGTAGCCACAAATAATGATTGCCAGCGCCAGATACAGGCTAAAGGTTGAATTTGTATCAATACCCAGCACCACATGGCCAAGGGTCATTAACAGCGCGCCGGCACTAACTGCAGTGCGGTTGCCGAGCAGGCGGTCGGCAagcc                                                                                                                                     >  2:919351/1‑150 (MQ=255)
                 aTAAACCGTAGCCACAAATAATGATTGCCAGCGCCAGATACAGGCTAAAGGTTGAATTTGTATCAATACCCAGCACCACATGGCCAAGGGTCATTAACAGCGCGCCGGCAATAACTGCAGTGCGGTTGCCGAGCAGGCGGTCGGCAagcc                                                                                                                                     >  1:2184596/1‑150 (MQ=255)
                   aaaCCGTAGCCACAAATAATGATTGCCAGCGCCAGATACAGGCTAAAGGTTGAATTTGTATCAATACCCAGCACCACATGGCCAAGGGTCATTAACAGCGCGCCGGCAATAACTGCAGTGCGGTTGCCGAGCAGGCGGTCGGCAagccag                                                                                                                                   >  2:2773494/1‑150 (MQ=255)
                     aCCGTAGCCACAAATAATGATTGCCAGCGCCAGATACAGGCTAAAGGTTGAATTTGTATCAATACCCAGCACCACATGGCCAAGGGTCATTAACAGCGCGCCGGCAATAACTGCAGTGCGGTTGCCGAGCAGGCGGTCGGCAagccagcc                                                                                                                                 >  1:1931571/1‑150 (MQ=255)
                      ccGGAGCCACAAATAATGATTGCCAGCGCCAGATACAGGCTAAAGGTTGAATTTGTATCAATACCCAGCACCACATGGCCAAGGGTCATTAACAGCGCGCCGGCAATAACTGCAGTGCGGTTGCCGAGCAGGCGGTCGGCAAGCCAgccg                                                                                                                                <  2:1795490/150‑1 (MQ=255)
                          aGCCACAAATAATGATTGCCAGCGCCAGATACAGGCTAAAGGTTGAATTTGTATCAATACCCAGCACCACATGGCCAAGGGTCATTAACAGCGCGCCGGCAATAACTGCAGTGCGGTTGCCGAGCAGGCGGTCGGCAAGCCAGCCGCCga                                                                                                                            <  2:2954684/150‑1 (MQ=255)
                            ccACAAATAATGATTGCCAGCGCCAGATACAGGCTAAAGGTTGAATTTGTATCAATACCCAGCACCACATGGCCAAGGGTCATTAACAGCGCGCCGGCAATAACTGCAGTGCGGTTGCCGAGCAGGCGGTCGGCAAGCCAGCCGCCgaga                                                                                                                          <  1:3212123/150‑1 (MQ=255)
                              acaAATAATGATTGCCAGCGCCAGATACAGGCTAAAGGTTGAATTTGTATCAATACCCAGCACCACATGGCCAAGGGTCATTAACAGCGCGCCGGCAATCACTGCAGTGCGGTTGCCGAGCAGGCGGTCGGCAAGCCAGCCGCCGAGAAt                                                                                                                        <  1:274945/150‑1 (MQ=255)
                                aaaTAATGATTGCCAGCGCCAGATACAGGCTAAAGGTTGAATTTGTATCAATACCCAGCACCACATGGCCAAGGGTCATTAACAGCGCGCCGGCAATAACTGCAGTGCGGTTGCCGAGCAGGCGGTCGGCAAGCCAGCCGCCGAGAATAg                                                                                                                      >  2:3683023/1‑150 (MQ=255)
                                 aataatGATTGCCAGCGCCAGATACAGGCTAAAGGTTGAATTTGTATCAATACCCAGCACCACATGGCCAAGGGTCATTAACAGCGCGCCGGCAATCACTGCAGTGCGGTTGCCGAGCAGGCGGTCGGCAAGCCAGCCGCCGAGCATAgg                                                                                                                     >  2:3448950/1‑150 (MQ=255)
                                  ataatGATTGCCAGCGCCAGATACAGGCTAAAGGTTGAATTTGTATCAATACCCAGCACCACATGGCCAAGGGTCATTAACCGCGCGCCGGCAATCACTGCAGTGCGGTTGCCGATCAGGCGGTCGGCAAGCCAGCCGCCGAGAATAggg                                                                                                                    >  2:1926031/1‑150 (MQ=255)
                                   taatGATTGCCAGCGCCAGATACAGGCTAAAGGTTGAATTTGTATCAATACCCAGCACCACATGGCCAAGGGTCATTAACCGCGCGCCGGCAATCACTGCAGTGCGGTTGCCGATCAGGCGGTCGGCAAGCCAGCCGCCGAGAATAgggg                                                                                                                   >  2:2990936/1‑150 (MQ=255)
                                   taatGATTGCCAGCGCCAGATACAGGCTAAAGGTTGAATTTGTATCAATACCCAGCACCACATGGCCAAGGGTCATTAACAGCGCGCCGGCAATCACTGCAGTGCGGTTGCCGATCAGGCGGTCGGCAAGCCAGCCGCCGAGAATAgggg                                                                                                                   <  2:3240448/150‑1 (MQ=255)
                                     atGATTGCCAGCGCCAGATACAGGCTAAAGGTTGAATTTGTATCAATACCCAGCACCACATGGCCAAGGGTCATTAACAGCGCGCCGGCAATCACTGCAGTGCGGTTGCCGATCAGGCGGTCGGCAAGCCAGCCGCCGAGAATAGGGGTa                                                                                                                 <  1:2481772/150‑1 (MQ=255)
                                      tGATTGCCAGCGCCAGATACAGGCTAAAGGTTGAATTTGTATCAATACCCAGCACCACATGGCCAAGGGTCATTAACAGCGCGCCGGCAATCACTGCAGTGCGGTTGCCGAGCAGGCGGTCGGCAAGCCAGCCGCCGAGAATAGGGGTaa                                                                                                                >  1:711045/1‑150 (MQ=255)
                                       gATTGCCAGCGCCAGATACAGGCTAAAGGTTGAATTTGTATCAATACCCAGCACCACATGGCCAAGGGTCATTAACCGCGCGCCGGCAATCACTGCAGTGCGGTTGCCGAGCAGGCGGTCGGCAAGCCAGCCGCCGATAATAGGGGTAAc                                                                                                               <  2:914075/150‑1 (MQ=255)
                                          tGCCAGCGCCAGATACAGGCTAAAGGTTGAATTTGTATCAATACCCAGCACCACATGGCCAAGGGTCATTAACAGCGCGCCGGCAATCACTGCAGTGCGGTTGCCGATCAGGCGGTCGGCAAGCCAGCCGCCGAGAATAGGGGTAACGTa                                                                                                            >  1:1984313/1‑150 (MQ=255)
                                          tGCCAGCGCCAGATACAGGCTAAAGGTTGAATTTGTATCAATACCCAGCACCACATGGCCAAGGGTCATTAACAGCGCGCCGGCAATCACTGCAGTGCGGTTGCCGAGCAGGCGGTCGGCAAGCCAGCCGCCGAGAATAGGGGTAACGTa                                                                                                            >  1:2331023/1‑150 (MQ=255)
                                          tGCCAGCGCCAGATACAGGCTAAAGGTTGAATTTGTATCAATACCCAGCACCACATGGCCAAGGGTCATTAACAGCGCGCCGGCAATAACTGCAGTGCGGTTGCCGAGCAGGCGGTCGGCAAGCCAGCCGCCGAGAATAGGGGTAACGTa                                                                                                            <  2:2309623/150‑1 (MQ=255)
                                           gCCAGCGCCAGATACAGGCTAAAGGTTGAATTTGTATCAATACCCAGCACCACCTGGCCAAGGGTCATTAACAGCGCGCCGGCAATAACTGCAGTGCGGTTGCCGAGCAGGCGGTCGGCAAGCCAGCCGCCGAGAATAGGGGTAACGTaa                                                                                                           <  2:3956414/150‑1 (MQ=255)
                                             cAGCGCCAGATACAGGCTAAAGGTTGAATTTGTATCAATACCCAGCACCACATGGCCAAGGGTCATTAACAGCGCGCCGGCAATCACTGCAGTGCGGTTGCCGATCAGGCGGTCGGCAAGCCAGCCGCCGAGAATAGGGGTAACGTAAAc                                                                                                         <  1:118438/150‑1 (MQ=255)
                                             cAGCGCCAGATACAGGCTAAAGGTTGAATTTGTATCAATACCCAGCACCACATGGCCAAGGGTCATTAACAGCGCGCCGGCAATCACTGCAGTGCGGTTGCCGAGCAGGCGGTCGGCAAGCCAGCCGCCGAGAATAGGGGTAACGTAAAc                                                                                                         >  2:1781232/1‑150 (MQ=255)
                                             cAGCGCCAGATACAGGCTAAAGGTTGAATTTGTATCAATACCCAGCACCACATGGCCAAGGGTCATTAACAGCGCGCCGGCAATCACTGCAGTGCGGTTGCCGAGCAGGCGGTCGGCAAGCCAGCCGCCGAGAATAGGGGTAACGTAAAc                                                                                                         >  2:1206446/1‑150 (MQ=255)
                                             cAGCGCCAGATACAGGCTAAAGGTTGAATTTGTATCAATACCCAGCACCACATGGCCAAGGGTCATTAACAGCGCGCCGGCAATAACTGCAGTGCGGTTGCCGAGCAGGCGGTCGGCAAGCCCGCCGCCGAGAATAGGGGTAACGTAAAc                                                                                                         >  2:3839211/1‑150 (MQ=255)
                                              aGCGCCAGATACAGGGTAAAGGTTGAATTTGTATCAATACCCAGCACCACATGGCCAAGGGTCATTAACAGCGCGCCGGCAATAACTGCAGTGCGGTTGCCGAGCAGGCGGTCGGCAAGCCAGCCGCCGAGAATAGGGGTAACGTAAAcc                                                                                                        >  1:2647695/1‑150 (MQ=255)
                                                cgcCAGATACAGGCTAAAGGTTGAATTTGTATCAATACCCAGCACCACATGGCCAAGGGTCATTAACAGCGCGCCGGCAATCACTGCAGTGCGGTTGCCGATCAGGCGGTCGGCAAGCCAGCCGCCGAGAATAGGGGTAACGTAAACCag                                                                                                      <  1:3363600/150‑1 (MQ=255)
                                                cgcCAGATACAGGCTAAAGGTTGAATTTGTATCAATACCCAGCACCACATGGCCAAGGGTCATTAACAGCGCGCCGGCAATAACTGCAGTGCGGTTGCCGAGCAGGCGGTCGGCAAGCCAGCCGCCGAGAATAGGGGTAACGTAAACCag                                                                                                      <  2:2390737/150‑1 (MQ=255)
                                                 gcCAGATACAGGCTAAAGGTTGAATTTGTATCAATACCCAGCACCACATGGCCAAGGGTCATTAACAGCGCGCCGGCAATAACTGCAGTGCGGTTGCCGAGCAGGCGGTCGGCAAGCCAGCCGCCGAGAATAGGGGTAACGTAAACCaga                                                                                                     <  1:2240132/150‑1 (MQ=255)
                                                     gATACAGGCTAAAGGTTGAATTTGTATCAATACCCAGCACCACATGGCCAAGGGTCATTAACAGCGCGCCGGCAATCACTGCAGTGCGGTTGCCGATCAGGCGGTCGGCAAGCCAGCCGCCGAGAATAGGGGTAACGTAAACCAgagacg                                                                                                 >  2:1305035/1‑148 (MQ=255)
                                                     gATACAGGCTAAAGGTTGAATTTGTATCAATACCCAGCACCACATGGCCAAGGGTCATTAACAGCGCGCCGGCAATAACTGCAGTGCGGTTGCCGAGCAGGCGGTCGGCAAGCCAGCCGCCGAGAATAGGGGTAACGTAAACCAGAGAAg                                                                                                 <  2:2456631/150‑1 (MQ=255)
                                                      aTACAGGCTAAAGGTTGAATTTGTATCAATACCCAGCACCACATGGCCAAGGGTCATTAACAGCGCGCCGGCAATCACTGCAGTGCGGTTGCCGATCAGGCGGTCGGCAAGCCAGCCGCCGAGACTAGGGGTAACGTAAACCAGAGAAGc                                                                                                >  2:2680101/1‑150 (MQ=255)
                                                          aGGCTAAAGGTTGAATTTGTATCAATACCCAGCACCACATGGCCAAGGGTCATTAACAGCGCGCCGGCAATCACTGCAGTGCGGTTGCCGATCAGGCGGTCGGCAAGCCAGCCGCCGAGAATAGGGGTAACGTAAACCAGAGAAGCatat                                                                                            >  1:3787592/1‑150 (MQ=255)
                                                          aGGCTAAAGGTTGAATTTGTATCAATACCCAGCACCACATGGCCAAGGGTCATTAACAGCGCGCCGGCAATCACTGCAGTGCGGTTGCCGATCAGGCGGTCGGCAAGCCAGCCGCCGAGAATAGGGGTAACGTAAACCAGAGAAGCatat                                                                                            >  1:2544149/1‑150 (MQ=255)
                                                            gCTAAAGGTTGAATTTGTATCAATACCCAGCACCACATGGCCAAGGGTCATTAACAGCGCGCCGGCAATCACTGCAGTGCGGTTGCCGATCATGCGGTCGGCAAGCCAGCCGCCGAGAATAGGGGTAACGTAAACCAGAGAAGCATATgc                                                                                          >  2:1638891/1‑150 (MQ=255)
                                                                aaGGTTGAATTTGTATCAATACCCAGCACCACATGGCCAAGGGTCATTAACAGCGCGCCGGCAATAACTGCAGTGCGGTTGCCGAGCAGGCGGTCGGCAAGCCAGCCGCCGAGAATAGGGGTAACGTAAACCAGAGACGCATATGCGCTg                                                                                      >  2:3741389/1‑150 (MQ=255)
                                                                 aGGTTGAATTTGTATCAATACCCAGCACCACATGGCCAAGGGTCATTAACAGCGCGCCGGCAATCACTGCAGTGCGGTTGCCGATCAGGCGGTCGGCAAGCCAGCCGCCGAGAATAGGGGTAACGTAAACCAGAGAAGCATATGCGCTGa                                                                                     <  2:1217622/150‑1 (MQ=255)
                                                                   gTTGAATTTGTATCAATACCCAGCACCACATGGCCAAGGGTCATTAACAGCGCGCCGGCAATCACTGCAGTGCGGTTGCCGATCAGGCGGTCGGCAAGCCAGCCGCCGAGAATAGGGGTAACGTAAACCAGAGAAGCATATGCGCTGAAc                                                                                   >  2:684237/1‑150 (MQ=255)
                                                                   gTTGAATTTGTATCAATACCCAGCACCACATGGCCAAGGGTCAGTAACAGCGCGCCGGCACTAACTGCAGTGCGGTTGCCGAGCAGGCGGTCGGCAAGCCAGCCGCCGAGACGAGGGGTAACGTAAACCAGAGACGCATATGCGCTGAAc                                                                                   >  1:3998459/1‑150 (MQ=255)
                                                                     tGAATTTGTATCAATACCCAGCACCACATGGCCAAGGGTCATTAACAGCGCGCCGGCAATCACTGCAGTGCGGTTGCCGATCAGGCGGTCGGCAAGCCAGCCGCCGAGAATAGGGGTAACGTAAACCAGAGAAGCATATGCGCTGAACAg                                                                                 <  1:3046922/150‑1 (MQ=255)
                                                                        aTTTGTATCAATACCCAGCACCACATGGCCAAGGGTCATTAACAGCGCGCCGGCAATCACTGCAGTGCGGTTGCCGATCAGGCGGTCGGCAAGCCAGCCGCCGAGAATAGGGGTAACGTAAACCAGAGAAGCATATGCGCTGAACAGGCt                                                                              >  1:754190/1‑150 (MQ=255)
                                                                             tATCAATACCCAGCACCACATGGCCAAGGGTCATTAACAGCGCGCCGGCAATAACTGCAGTGCGGTTGCCGAGCAGGCGGTCGGCAAGCCAGCCGCCGAGAATAGGGGTAACGTAAACCAGAGAAGCATATGCGCTGAACAGGCTGATgg                                                                         >  1:3609535/1‑150 (MQ=255)
                                                                                     cccAGCACCACATGGCCAAGGGTCATTAACAGCGCGCCGGCAATCACTGCAGTGCGGTTGCCGATCAGGCGGTCGGCAAGCCAGCCGCCGAGAATAGGGGTAACGTAAACCAGAGACGCATATGCGCTGAACAGGCTGATGGCATGGTTa                                                                 >  2:2003240/1‑150 (MQ=255)
                                                                                      ccAGCACCACATGTCCAAGGGTCATTAACAGCGCGCCGGCAATCACTGCAGTGCGGTTGCCGATCAGGCGGTCGGCAAGCCAGCCGCCGAGAATAGGGGTAACGTAAACCAGAGAAGCATATGCGCTGAACAGGCTGATGGCATGGTTat                                                                >  1:1504553/1‑150 (MQ=255)
                                                                                      ccAGCACCACATGGCCAAGGGTCATTAACAGCGCGCCGGCAATCACTGCAGTGCGGTTGCCGATCAGGCGGTCGGCAAGCCAGCCGCCGAGAATAGGGGTAACGTAAACCAGAGAAGCATATGCGCTGAACAGGCTGATGGCATGGTTat                                                                <  1:45432/150‑1 (MQ=255)
                                                                                      ccAGCACCACATGGCCAAGGGTCATTAACAGCGCGCCGGCAATCAATGCAGTGCGGTTGCCGATCAGGCGTTCGGCAAGCCAGCCGCCGAGAATAGGGGTACCGAAAACCAGAGAATCAAATGCGCTGAACAGGCTGATGGCATGGTTat                                                                >  1:3760979/1‑150 (MQ=255)
                                                                                      ccAGCACCACATGGCCAAGGGTCATTAACAGCGCGCCGGCAATAACTGCAGTGCGGTTGCCGAGCAGGCGGTCGGCAAGCCAGCCGCCGAGAATAGGGGTAACGTAAACCAGAGAAGCATATGCGCTGAACAGGCTGATGGCATGGTTat                                                                <  1:3847091/150‑1 (MQ=255)
                                                                                       cAGCACCACATGGCCAAGGGTCATTAACAGCGCGCCGGCAATCACTGCAGTGCGGTTGCCGATCAGGCGGTCGGCAAGCCAGCCGCCGAGAATAGGGGTAACGTAAACCAGAGAAGCATATGCGCTGAACAGGCTGATGGCATGGTTatc                                                               <  2:2897818/150‑1 (MQ=255)
                                                                                       cAGCACCACATGGCCAAGGGTCATTAACAGCGCGCCGGCAATAACTGCAGTGCGGTTGCCGAGCAGGCGGTCGGCAAGCCAGCCGCCGAGAATAGGGGTAACGTAAACCAGAGAAGCATATGCGCTGAACAGGCTGATGGCATGGTTatc                                                               >  2:3698391/1‑150 (MQ=255)
                                                                                            ccacCTGGCCACGGGTCATTCACCGCGCGCCGGCAATCACTGCAGTGCGGTTGCCGATCAGGCGGTCGGCAAGCCAGCCGCCGAGAATAGGGGTAACGTAAACCAGAGAAGCATATGCGCTGAACAGGCTGATGGCATGGTTATCATCaa                                                          <  2:3811963/150‑1 (MQ=255)
                                                                                             cacaTGGCCAAGGGTCATTAACAGCGCGCCGGCAATCACTGCAGTGCGGTTGCCGAGCAGGCGGTCGGCAAGCCAGCCGCCGAGAATAGGGGTAACGTAAACCAGAGAAGCATATGCGCTGAACAGGCTGATGGCATGGTTATCATCaaa                                                         >  1:2470999/1‑150 (MQ=255)
                                                                                                   gCCAAGGGTCATTAACAGCGCGCCGGCAATCACTGCAGTGCGGTTGCCGATCAGGCGGTCGGCAAGCCAGCCGCCGAGAATAGGGGTAACGTAAACCAGAGAAGCATATGCGCTGAACAGGCTGATGGCATGGTTATCATCAAAACCAAg                                                   <  2:1066801/150‑1 (MQ=255)
                                                                                                   gCCAAGGGTCATTAACAGCGCGCCGGCAATAACTGCAGTGCGGTTGCCGAGCAGGCGGTCGGCAAGCCAGCCGCCGAGAATAGGGGTAACGTAAACCAGAGAAGCATATGCGCTGAACAGGCTGATGGCATGGTTATCATCAAAACCAAg                                                   <  2:15101/150‑1 (MQ=255)
                                                                                                    ccAAGGGTCATTAACAGCGCGCCGGCAATAACTGCAGTGCGGTTGCCGAGCAGGCGGTCGGCAAGCCAGCCGCCGAGAATAGGGGTAACGTAAACCAGAGAAGCATATGCGCTGAACAGGCTGATGGCATGGTTATCATCAAAACCAAGc                                                  >  1:627357/1‑150 (MQ=255)
                                                                                                     cAAGGGTCATTAACAGCGCGCCGGCAATAACTGCAGTGCGGTTGCCGAGCAGGCGGTCGGCAAGCCAGCCGCCGAGAATAGGGGTAACGTAAACCAGAGAAGCATATGCGCTGAACAGGCTGATGGCATGGTTATCATCAAAACCAAGCt                                                 >  1:2832911/1‑150 (MQ=255)
                                                                                                      aaGGGTCATTAACAGCGCGCCGGCAATAACTGCAGTGCGGTTGCCGAGCAGGCGGTCGGCAAGCCAGCCGCCGAGAATAGGGGTAACGTAAACCAGAGAAGCATATGCGCTGAACAGGCTGATGGCATGGTTATCATCAAAACCAAGCTg                                                >  2:492390/1‑150 (MQ=255)
                                                                                                         ggTCATTAACAGCGCGCCGGCAATAACTGCAGTGCGGTTGCCGAGCAGGCGGTCGGCAAGCCAGCCGCCGAGAATAGGGGTAACGTAAACCAGAGAAGCATATGCGCTGAACAGGCTGATGGCATGGTTATCATCAAAACCAAGCTGATg                                             >  1:3696200/1‑150 (MQ=255)
                                                                                                           taaTTAACAGCGCGCCGGCAATAACTGCAGTGCGGTTGCCGAGCAGGCGGTCGGCAAGCCAGCCGCCGAGAATAGGGGTAACGTAAACCAGAGAAGCATATGCGCTGAACAGGCTGATGGCATGGTTATCATCAAAACCAAGCTGATggg                                           <  2:1555084/148‑1 (MQ=255)
                                                                                                           tCATTAACAGCGCGCCGGCAATCACTGCAGTGCGGTTGCCGATCAGGCGGTCGGCAAGCCAGCCGCCGAGAATAGGGGTAACGTAAACCAGAGAAGCATATGCGCTGAACAGGCTGATGGCATGGTTATCATCAAAACCAAGCTGATggg                                           <  2:2665722/150‑1 (MQ=255)
                                                                                                           tCATTAACAGCGCGCCGGCAATCACTGCAGTGCGGTTGCCGAGCAGGCGGTCGGCAAGCCAGCCGCCGAGAATAGGGGTAACGTAAACCAGAGAAGCATATGCGCTGAACAGGCTGATGGCATGGTTATCATCAAAACCAAGCTGATggg                                           <  1:3982319/150‑1 (MQ=255)
                                                                                                           tCATTAACAGCGCGCCGGCAAGAACTGCAGTGCGGTTGCCGAGCAGGCGGTCGGCAAGCCAGCCGCCGAGAATAGGGGTAACGTAAACCAGAGAAGCATATGCGCTGAACAGGCTGATGGCATGGTTATCATCAAAACCAAGCTGATggg                                           <  2:1555097/150‑1 (MQ=255)
                                                                                                            cATTAACAGCGCGCCGGCAATCACTGCAGTGCGGTTGCCGATCAGGCGGTCGGCAAGCCAGCCGCCGAGAATAGGGGTAACGTAAACCAGAGAAGCATATGCGCTGAACAGGCTGATGGCATGGTTATCATCAAAACCAAGCTGATGGGt                                          <  1:2116508/150‑1 (MQ=255)
                                                                                                               tAACAGCGCGCCGGCAATCGCTGCAGTGCGGTTGCCGATCAGGCGGTCGGCAAGCCAGCCGCCGAGAATAGGGGTAACGTAAACCAGAGAAGCATATGCGCTGAACAGGCTGATGGCATGGTTATCATCAAAACCAAGCTGATGGGTgag                                       <  2:4023460/150‑1 (MQ=255)
                                                                                                               tAACAGCGCGCCGGCAATCACTGCAGTGCGGTTGCCGATCAGGCGGTCGGCAAGCCAGCCGCCGAGAATAGGGGTAACGTAAACCAGAGAAGCATATGCGCTGAACAGGCTGATGGCATGGTTATCATCAAAACCAAGCTGATGGGTgag                                       >  1:1724054/1‑150 (MQ=255)
                                                                                                                aaCAGCGCGCCGGCAATCACTGCAGTGCGGTTGCCGATCAGGCGGTCGGCAAGCCAGCCGCCGAGAATAGGGGTAACGTAAACCAGAGAAGCATATGCGCTGAACAGGCTGATGGCATGGTTATCATCAAAACCAAGCTGATGGGTgaga                                      >  1:2650327/1‑150 (MQ=255)
                                                                                                                aaCAGCGCGCCGGCAATAACTGCAGTGCGGTTGCCGAGCAGGCGGTCGGCAAGCCAGCCGCCGAGAATAGGGGTAACGTAAACCAGAGAAGCATATGCGCTGAACAGGCTGATGGCATGGTTATCATCAAAACCAAGCTGATGGGTgaga                                      >  2:2513587/1‑150 (MQ=255)
                                                                                                                aaCAGCGCGCCGGCAATAACTGCAGTGCGGTTGCCGAGCAGGCGGTCGGCAAGCCAGCCGCCGAGAATAGGGGTAACGTAAACCAGAGAAGCATATGCGCTGAACAGGCTGATGGCATGGTTATCATCAAAACCAAGCTGATGGGTgaga                                      <  2:1170478/150‑1 (MQ=255)
                                                                                                                  cAGCGCGCCGGCAATCACTGCAGTGCGGTTGCCGATCAGGCGGTCGGCAAGCCAGCCGCCGAGAATAGGGGTAACGTAAACCAGAGAAGCATATGCGCTGAACAGGCTGATGGCATGGTTATCATCAAAACCAAGCTGATGGGTGAGATa                                    <  2:21143/150‑1 (MQ=255)
                                                                                                                  cAGCGCGCCGGCAATCACTGCAGTGCGGTTGCCGATCAGGCGGTCGGCAAGCCAGCCGCCGAGAATAGGGGTAACGTAAACCAGAGAAGCATATGCGCTGAACAGGCTGATGGCATGGTTATCATCAAAACCAAGCTGATGGGTGAGATa                                    <  2:3545674/150‑1 (MQ=255)
                                                                                                                  cAGCGCGCCGGCAATAACTGCAGTGCGGTTGCCGAGCAGGCGGTCGGCAAGCCAGCCGCCGAGAATAGGGGTAACGTAAACCAGAGAAGCATATGCGCTGAACAGGCTGATGGCATGGTTATCATCAAAACCAAGCTGATGGGTGAGATa                                    >  1:100162/1‑150 (MQ=255)
                                                                                                                   aGCGCGCCGGCAATCACTGCAGTGCGGTTGCCGATCAGGCGGTCGGCAAGCCAGCCGCCGAGAATAGGGGTAACGTAAACCAGAGAAGCATATGCGCTGAACAGGCTGATGGCATGGTTATCATCATAACCAAGCTGATGGGTGAGATag                                   >  2:3912858/1‑150 (MQ=255)
                                                                                                                     cgcgcCGGCACTAACTGCAGTGCGGTTGCCGAGCAGGCGGTCGGCAAGCCAGCCGCCGAGAATAGGGGTAACGTAAACCAGAGAAGCATATGCGCTGAACAGGCTGATGGCATGGTTATCATCAAAACCAAGCTGATGGGTGAGATagag                                 >  2:747510/1‑150 (MQ=255)
                                                                                                                     cgcgcCGGCAATCACTGCAGTGCGGTTGCCGAGCAGGCGGTCGGCAAGCCAGCCGCCGAGAATAGGGGTAACGTAAACCAGAGAAGCATATGCGCTGAACAGGCTGATGGCATGGTTATCATCAAAACCAAGCTGATGGGTGAGATagag                                 <  1:3879798/150‑1 (MQ=255)
                                                                                                                     cgcgcCGGCAATAACTGCAGTGCGGTTGCCGAGCAGGCGGTCGGCAAGCCAGCCGCCGAGACTAGGGGTAACGTAAACCAGAGAAGCATATGCGCTGAACAGGCTGATGGCATGGTTATCATCAAAACCAAGCTGATGGGTGAGATagag                                 >  1:142950/1‑150 (MQ=255)
                                                                                                                      gcgcCGGCAATAACTGCAGTGCGGTTGCCGAGCAGGCGGTCGGCAAGCCAGCCGCCGAGAATAGGGGTAACGTAAACCAGAGAAGCATATGCGCTGAACAGGCTGATGGCATGGTTATCATCAAAACCAAGCTGATGGGTGAGATAgaga                                <  2:3568243/150‑1 (MQ=255)
                                                                                                                           ggCAATCACTGCAGTGCGGTTGCCGAGCAGGCGGTCGGCAAGCCAGCCGCCGAGAATAGGGGTAACGTAAACCAGAGAAGCATATGCGCTGAACAGGCTGATGGCATGGTTATCATCAAAACCAAGCTGATGGGTGAGATAGAGAATGAg                           >  2:2275411/1‑150 (MQ=255)
                                                                                                                            gCAATCACTGCAGTGCGGTTGCCGATCAGGCGGTCGGCAAGCCAGCCGCCGAGAATAGGGGTAACGTAAACCAGAGAAGCATATGCGCTGAACAGGCTGATGGCATGGTTATCATCAAAACCAAGCTGATGGGTGAGATAGAGAATGAGt                          >  1:1463333/1‑150 (MQ=255)
                                                                                                                              aaTCACTGCAGTGCGGTTGCCGATCAGGCGGTCGGCAAGCCAGCCGCCGAGAATAGGGGTAACGTAAACCAGAGAAGCATATGCGCTGAACAGGCTGATGGCATGGTTATCATCAAAACCAAGCTGATGGGTGAGATAGAGAATGAGTaa                        >  2:3686497/1‑150 (MQ=255)
                                                                                                                              aaTAACTGCAGTGCGGTTGCCGAGCAGGCGGTCGGCAAGCCAGTCGCCGAGAATAGGGGTAACGTAAACCAGAGAAGCATATGCGCTGAACAGGCTGATGGCATGGTTATCATCAAAACCAAGCTGATGGGTGAGATAGAGAATGAGTaa                        <  2:3775200/150‑1 (MQ=255)
                                                                                                                               ataaCTGCAGTGCGGTTGCCGAGCAGGCGGTCGGCAAGCCAGCCGCCGAGAATAGGGGTAACGTAAACCAGAGAAGCATATGCGCTGAACAGGCTGATGGCATGGTTATCATCAAAACCAAGCTGATGGGTGAGATAGAGAATGAGTAAg                       >  1:2863215/4‑150 (MQ=255)
                                                                                                                               aTCACTGCAGTGCGGTTGCCGATCAGGCGGTCGGCAAGCCAGCCGCCGAGAATAGGGGTAACGTAAACCAGAGAAGCATATGCGCTGAACAGGCTGATGGCATGGTTATCATCAAAACCAAGCTGATGGGTGAGATAGAGAATGAGTAAg                       >  1:2689188/1‑150 (MQ=255)
                                                                                                                                     gCAGTGCGGTTGCCGAGCAGGCGGTCGGCAAGCCAGCCGCCGAGAATAGGGGTAACGTAAACCAGAGAAGCATATGCGCTGAACAGGCTGATGGCATGGTTATCATCAAAACCAAGCTGATGGGTGAGATAGAGAATGAGTAAGGCACGc                 >  1:881509/1‑150 (MQ=255)
                                                                                                                                        gTGCGGTTGCCGAGCAGGCGGTCGGCAAGCCAGCCGCCGAGAATAGGGGTAACGTAAACCAGAGAAGCATATGCGCTGAACAGGCTGATGGCATGGTTATCATCAAAACCAAGCTGATGGGTGAGATAGAGAATGAGTAAGGCACGCATg              <  2:4052229/150‑1 (MQ=255)
                                                                                                                                         tGCGGTTGCCGATCAGGCGGTCGGCAAGCCAGCCGCCGAGAATAGGGGTAACGTAAACCAGAGAAGCATATGCGCTGAACAGGCTGATGGCATGGTTATCATCAAAACCAAGCTGATGGGTGAGATAGAGAATGAGTAAGGCACGCATGc             <  1:243949/150‑1 (MQ=255)
                                                                                                                                            ggTTGCCGAGCAGGCGGTCGGCAAGCCAGCCGCCGAGAATAGGGGTAACGTAAACCAGAGAAGCATATGCGCTGAACAGGCTGATGGCATGGTTATCATCAAAACCAAGCTGATGGGTGAGATAGAGAATGAGTAAGGCACGCATGCCGt          <  1:391879/150‑1 (MQ=255)
                                                                                                                                             gTTGCCGAGCAGGCGGTCGGCAAGCCAGCCGCCGAGAATAGGGGTAACGTAAACCAGAGAAGCATATGCGCTGAACAGGCTGATGGCATGGTTATCATCAAAACCAAGCTGATGGGTGAGATAGAGAATGAGTAAGGCACGCATGCCGTa         <  1:1649185/150‑1 (MQ=255)
                                                                                                                                                 ccGATCAGTCGGTCGGCAAGCCAGCCGCCGAGAATAGGGGTAACGTAAACCAGAGAAGCATATGCGCTGAACAGGCTGATGGCATGGTTATCATCAAAACAAAGCTGATGGGTGAGATAGAGAATGAGTAAGGAACGCATGCCGTaaaaa     >  2:3882363/1‑150 (MQ=255)
                                                                                                                                                  cGATCAGGCGGTCGGCAAGCCAGCCGCCGAGAATAGGGGTAACGTAAACCAGAGAAGCATATGCGCTGAACAGGCTGATGGCATGGTTATCATCAAAACCAAGCTGATGGGTGAGATAGAGAATGAGTAAGGCACGCATGCCGTAAAAAc    <  1:3618443/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
GCTGATGTTTGATTTGAATAAACCGTAGCCACAAATAATGATTGCCAGCGCCAGATACAGGCTAAAGGTTGAATTTGTATCAATACCCAGCACCACATGGCCAAGGGTCATTAACAGCGCGCCGGCAATCACTGCAGTGCGGTTGCCGAGCAGGCGGTCGGCAAGCCAGCCGCCGAGAATAGGGGTAACGTAAACCAGAGAAGCATATGCGCTGAACAGGCTGATGGCATGGTTATCATCAAAACCAAGCTGATGGGTGAGATAGAGAATGAGTAAGGCACGCATGCCGTAAAAACTG  >  CP009273/4345866‑4346163

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: