Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,380,174 C→A 9.7% S458S (TCG→TCT)  aaeB ← p‑hydroxybenzoic acid efflux system component

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,380,1740CA 9.7% 578.3 / 23.2 164S458S (TCG→TCT) aaeBp‑hydroxybenzoic acid efflux system component
Reads supporting (aligned to +/- strand):  ref base C (63/84);  new base A (6/11);  total (69/95)
Fisher's exact test for biased strand distribution p-value = 6.13e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 165 total aligned reads displayed.

GATTTATCCCGCACCAGCAAAATAACGGTGAACGCGAGCACACAGCCGACGATTTGCCCTAATGCGCTGTCGAGAAACTGACTGAAATGGAAAGTCATCGGGTTATCCAGCACGATAATATTTATGGTGCTGGCCAGTGCCCCCATCGAGCCCAGTCGCCGTTTCTGTACTTCTATACCGAGGAAGAATCCCAGCACTGCCAGGCTAATGCACAGCAGCAACATGCTCTGTTGGGTATTAGGGATAATCACCAAAAAGTAGAGCAGCCCTAACGGCAGCGCGGCCAGCGTCCCG  >  CP009273/3380028‑3380321
                                                                                                                                                  |                                                                                                                                                   
 aTTTATCCCGCACCAGCAAAATAACGGTGAACGCGAGCACACAGCCGACGATTTGCCCTAATGCGCTGTCGAGAAACTGACTGAAATGGAAAGTCATCGGGTTATCCAGCACGATAATATTTATGGTGCTGGCCAGTGCCCCCATCGa                                                                                                                                                   >  1:965611/1‑148 (MQ=255)
 aTTTATCCCGCACCAGCAAAATAACGGTGAACGCGAGCACACAGCCGACGATTTGCCCTAATGCGCTGTCGAGAAACTGACTGAAATGGAAAGTCATCGGGTTATCCAGCACGATAATATTTATGGTGCTGGCCAGTGCCCCCATAGa                                                                                                                                                   >  2:134571/1‑148 (MQ=255)
     aTCCCGCACCAGCAAAATAACGGTGAACGCGAGCACACAGCCGACGATTTGCCCTAATGCGCTGTCGAGAAACTGACTGAAATGGAAAGTCATCGGGTTATCCAGCACGATAATATTTATGGTGCTGGCCAGTGCCCCCATCGAGccc                                                                                                                                               <  2:1787099/148‑1 (MQ=255)
      tCCCGCACCAGCAAAATAACGGTGAACGCGAGCACACAGCCGACGATTTGCCCTAATGCGCTGTCGAGAAACTGACTGAAATGGAAAGTCATCGGGTTATCCAGCACGATAATATTTATGGTGCTGGCCAGTGCCCCCATCGAGCCCa                                                                                                                                              >  1:3062242/1‑148 (MQ=255)
       cccGCACCAGCAAAATAACGGTGAACGCGAGCACACAGCCGACGATTTGCCCTAATGCGCTGTCGAGAAACTGACTGAAATGGAAAGTCATCGGGTTATCCAGCACGATAATATTTATGGTGCTGGCCAGTGCCCCCATCGAGCCCAg                                                                                                                                             <  2:244920/148‑1 (MQ=255)
       cccGCACCAGCAAAATAACGGTGAACGCGAGCACACAGCCGACGATTTGCCCTAATGCGCTGTCGAGAAACTGACTGAAATGGAAAGTCATCGGGTTATCCAGCACGATAATATTTATGGTGCTGGCCAGTGCCCCCATCGAGCCCAg                                                                                                                                             <  1:3141095/148‑1 (MQ=255)
        ccGCACCAGCAAAATAACGGTGAACGCGAGCACACAGCCGACGATTTGCCCTAATGCGCTGTCGAGAAACTGACTGAAATGGAAAGTCATCGGGTTATCCAGCACGATAATATTTATGGTGCTGGCCAGTGCCCCCATCGAGCCCAGt                                                                                                                                            <  1:2164342/148‑1 (MQ=255)
          gCACCAGCAAAATAACGGTGAACGCGAGCACACAGCCGACGATTTGCCCTAATGCGCTGTCGAGAAACTGACTGAAATGGAAAGTCATCGGGTTATCCAGCACGATAATATTTATGGTGCTGGCCAGTGCCCCCATCGAGCCCAGTCg                                                                                                                                          <  1:2689660/148‑1 (MQ=255)
              cAGCAAAATAACGGTGAACGCGAGCACACAGCCGACGATTTGCCCTAATGCGCTGTCGAGAAACTGACTGAAATGGAAAGTCATCGGGTTATCCAGCACGATAATATTTATGGTGCTGGCCAGTGCCCCCATCGAGCCCAGTCGCCGt                                                                                                                                      <  1:1925872/148‑1 (MQ=255)
                  aaaaTAACGGTGAACGCGAGCACACAGCCGACGATTTGCCCTAATGCGCTGTCGAGAAACTGACTGAAATGGAAAGTCATCGGGTTATCCAGCACGATAATATTTATGGTGCTGGCCAGTGCCCCCATAGAGCCCAGTCGCCGTTTCt                                                                                                                                  <  1:133440/148‑1 (MQ=255)
                  aaaaTAACGGTGAACGCGAGCACACAGCCGACGATTTGCCCTAATGCGCTGTCGAGAAACTGACTGAAATGGAAAGTCATCGGGTTATCCAGCACGATAATATTTATGGTGCTGGCCAGTGCCCCCATAGAGCCCAGTCGCCGTTTCt                                                                                                                                  <  1:133924/148‑1 (MQ=255)
                    aaTAACGGTGAACGCGAGCACACAGCCGACGATTTGCCCTAATGCGCTGTCGAGAAACTGACTGAAATGGAAAGTCATCGGGTTATCCAGCACGATAATATTTATGGTGCTGGCCAGTGCCCCCATCGAGCCCAGTCGCCGTTTCTGt                                                                                                                                >  2:1948654/1‑148 (MQ=255)
                       aaCGGTGAACGCGAGCACACAGCCGACGATTTGCCCTAATGCGCTGTCGAGAAACTGACTGAAATGGAAAGTCATCGGGTTATCCAGCACGATAATATTTATGGTGCTGGCCAGTGCCCCCATCGAGCCCAGTCGCCGTTTCTGTACt                                                                                                                             <  2:808105/148‑1 (MQ=255)
                        aCGGTGAACGCGAGCACACAGCCGACGATTTGCCCTAATGCGCTGTCGAGAAACTGACTGAAATGGAAAGTCATCGGGTTATCCAGCACGATAATATTTATGGTGCTGGCCAGTGCCCCCATCGAGCCCAGTCGCCGTTTCTGTACtt                                                                                                                            <  2:2881643/148‑1 (MQ=255)
                             gAACGCGAGCACACAGCCGACGATTTGCCCTAATGCGCTGTCGAGAAACTGACTGAAATGGAAAGTCATCGGGTTATCCAGCACGATAATATTTATGGTGCTGGCCAGTGCCCCCATCGAGCCCAGTCGCCGTTTCTGTACTTCtata                                                                                                                       <  1:686190/148‑1 (MQ=255)
                              aaCGCGAGCACACAGCCGACGATTTGCCCTAATGCGCTGTCGAGAAACTGACTGAAATGGAAAGTCATCGGGTTATCCAGCACGATAATATTTATGGTGCTGGCCAGTGCCCCCATCGAGCCCAGTCGCCGTTTCTGTACTTCTATAc                                                                                                                      >  1:690048/1‑148 (MQ=255)
                                  cgAGCACACAGCCGACGATTTGCCCTAATGCGCTGTCGAGAAACTGACTGAAATGGAAAGTCATCGGGTTATCCAGCACGATAATATTTATGGTGCTGGCCAGTGCCCCCATAGAGCCCAGTCGCCGTTTCTGTACTTCTATACCGAg                                                                                                                  <  1:1018565/148‑1 (MQ=255)
                                     gCACACAGCCGACGATTTGCCCTAATGCGCTGTCGAGAAACTGACTGAAATGGAAAGTCATCGGGTTATCCAGCACGATAATATTTATGGTGCTGGCCAGTGCCCCCATCGAGCCCAGTCGCCGTTTCTGTACTTCTATACCGAGgaa                                                                                                               >  1:424106/1‑148 (MQ=255)
                                          caGCCGACGATTTGCCCTAATGCGCTGTCGAGAAACTGACTGAAATGGAAAGTCATCGGGTTATCCAGCACGATAATATTTATGGTGCTGGCCAGTGCCCCCATCGAGCCCAGTCGCCGTTTCTGTACTTCTATACCGAGGAAGAATc                                                                                                          <  1:1085718/148‑1 (MQ=255)
                                          caGCCGACGATTTGCCCTAATGCGCTGTCGAGAAACTGACTGAAATGGAAAGTCATCGGGTTATCCAGCACGATAATATTTATGGTGCTGGCCAGTGCCCCCATCGAGCCCAGTCGCCGTTTCTGTACTTCTATACCGAGGAAGAATc                                                                                                          <  1:1890799/148‑1 (MQ=255)
                                          caGCCGACGATTTGCCCTAATGCGCTGTCGAGAAACTGACTGAAATGGAAAGTCATCGGGTTATCCAGCACGATAATATTTATGGTGCTGGCCAGTGCCCCCATCGAGCCCAGTCGCCGTTTCTGTACTTCTATACCGAGGAAGAATc                                                                                                          >  2:400625/1‑148 (MQ=255)
                                                acgaTTTGCCCTAATGCGCTGTCGAGAAACTGACTGAAATGGAAAGTCATCGGGTTATCCAGCACGATAATATTTATGGTGCTGGCCAGTGCCCCCATCGAGCCCAGTCGCCGTTTCTGTACTTCTATACCGAGGAAGAATCCCAGCa                                                                                                    <  1:2352616/148‑1 (MQ=255)
                                                acgaTTTGCCCTAATGCGCTGTCGAGAAACTGACTGAAATGGAAAGTCATCGGGTTATCCAGCACGATAATATTTATGGTGCTGGCCAGTGCCCCCATCGAGCCCAGTCGCCGTTTCTGTACTTCTATACCGAGGAAGAATCCCAGCa                                                                                                    <  2:2528026/148‑1 (MQ=255)
                                                   aTTTGCCCTAATGCGCTGTCGAGAAACTGACTGAAATGGAAAGTCATCGGGTTATCCAGCACGATAATATTTATGGTGCTGGCCAGTGCCCCCATCGAGCCCAGTCGCCGTTTCTGTACTTCTATACCGAGGAAGAATCCCAGCACTg                                                                                                 <  2:3297375/148‑1 (MQ=255)
                                                    tttGCCCTAATGCGCTGTCGAGAAACTGACTGAAATGGAAAGTCATCGGGTTATCCAGCACGATAATATTTATGGTGCTGGCCAGTGCCCCCATCGAGCCCAGTCGCCGTTTCTGTACTTCTATACCGAGGAAGAATCCCAGCACTGc                                                                                                <  2:1956868/148‑1 (MQ=255)
                                                    tttGCCCTAATGCGCTGTCGAGAAACTGACTGAAATGGAAAGTCATCGGGTTATCCAGCACGATAATATTTATGGTGCTGGCCAGTGCCCCCATCGAGCCCAGTCGCCGTTTCTGTACTTCTATACCGAGGAAGAATCCCAGCACTGc                                                                                                >  2:762813/1‑148 (MQ=255)
                                                     ttGCCCTAATGCGCTGTCGAGAAACTGACTGAAATGGAAAGTCATCGGGTTATCCAGCACGATAATATTTATGGTGCTGGCCAGTGCCCCCATCGAGCCCAGTCGCCGTTTCTGTACTTCTATACCGAGGACGAATCCCAGCACTGcc                                                                                               >  2:1760261/1‑148 (MQ=255)
                                                     ttGCCCTAATGCGCTGTCGAGAAACTGACTGAAATGGAAAGTCATCGGGTTATCCAGCACGATAATATTTATGGTGCTGGCCAGTGCCCCCATCGAGCCCAGTCGCCGTTTCTGTACTTCTATACCGAGGAAGAATCCCAGCACTGcc                                                                                               >  1:2884972/1‑148 (MQ=255)
                                                     ttGCCCTAATGCGCTGTCGAGAAACTGACTGAAATGGAAAGTCATCGGGTTATCCAGCACGATAATATTTATGGTGCTGGCCAGTGCCCCCATCGAGCCCAGTCGCCGTTTCTGTACTTCTATACCGAGGAAGAATCCCAGCACTGcc                                                                                               <  1:3153688/148‑1 (MQ=255)
                                                      tGCCCTAATGCGCTGTCGAGAAACTGACTGAAATGGAAAGTCATCGGGTTATCCAGCACGATAATATTTATGGTGCTGGCCAGTGCCCCCATCGAGCCCAGTCGCCGTTTCTGTACTTCTATACCGAGGAAGAATCCCAGCACTGCCa                                                                                              <  1:1432307/148‑1 (MQ=255)
                                                      tGCCCTAATGCGCTGTCGAGAAACTGACTGAAATGGAAAGTCATCGGGTTATCCAGCACGATAATATTTATGGTGCTGGCCAGTGCCCCCATCGAGCCCAGTCGCCGTTTCTGTACTTCTATACCGAGGAAGAATCCCAGCACTGCCa                                                                                              <  2:2313881/148‑1 (MQ=255)
                                                       gCCCTAATGCGCTGTCGAGAAACTGACTGAAATGGAAAGTCATCGGGTTATCCAGCACGATAATATTTATGGTGCTGGCCAGTGCCCCCATCGAGCCCAGTCGCCGTTTCTGTACTTCTATACCGAGGAAGAATCCCAGCACTGCCAg                                                                                             >  2:1733231/1‑148 (MQ=255)
                                                       gCCCTAATGCGCTGTCGAGAAACTGACTGAAATGGAAAGTCATCGGGTTATCCAGCACGATAATATTTATGGTGCTGGCCAGTGCCCCCATCGAGCCCAGTCGCCGTTTCTGTACTTCTATACCGAGGAAGAATCCCAGCACTGCCAg                                                                                             >  2:1732372/1‑148 (MQ=255)
                                                         ccTAATGCGCTGTCGAGAAACTGACTGAAATGGAAAGTCATCGGGTTATCCAGCACGATAATATTTATGGTGCTGGCCAGTGCCCCCATCGAGCCCAGTCGCCGTTTCTGTACTTCTATACCGAGGAAGAATCCCAGCACTGCCAGGc                                                                                           >  2:2552626/1‑148 (MQ=255)
                                                         ccTAATGCGCTGTCGAGAAACTGACTGAAATGGAAAGTCATCGGGTTATCCAGCACGATAATATTTATGGTGCTGGCCAGTGCCCCCATCGAGCCCAGTCGCCGTTTCTGTACTTCTATACCGAGGAAGAATCCCAGCACTGCCAGGc                                                                                           <  2:414043/148‑1 (MQ=255)
                                                         ccTAATGCGCTGTCGAGAAACTGACTGAAATGGAAAGTCATCGGGTTATCCAGCACGATAATATTTATGGTGCTGGCCAGTGCCCCCATCGAGCCCAGTCGCCGTTTCTGTACTTCTATACCGAGGAAGAATCCCAGCACTGCCAGGc                                                                                           >  2:121138/1‑148 (MQ=255)
                                                         ccTAATGCGCTGTCGAGAAACTGACTGAAATGGAAAGTCATCGGGTTATCCAGCACGATAATATTTATGGTGCTGGCCAGTGCCCCCATAGAGCCCAGTCGCCGTTTCTGTACTTCTATACCGAGGAAGAATCCCAGCACTGCCAGGc                                                                                           <  1:702189/148‑1 (MQ=255)
                                                          cTAATGCGCTGTCGAGAAACTGACTGAAATGGAAAGTCATCGGGTTATCCAGCACGATAATATTTATGGTGCTGGCCAGTGCCCCCATCGAGCCCAGTCGCCGTTTCTGTACTTCTATACCGAGGAAGAATCCCAGCACTGCCAGGCt                                                                                          >  2:2994494/1‑148 (MQ=255)
                                                              tGCGCTGTCGAGAAACTGACTGAAATGGAAAGTCATCGGGTTATCCAGCACGATAATATTTATGGTGCTGGCCAGTGCCCCCATCGAGCCCAGTCGCCGTTTCTGTACTTCTATACCGAGGAAGAATCCCAGCACTGCCAGGCTAATg                                                                                      >  1:2735423/1‑148 (MQ=255)
                                                               gcgcTGTCGAGAAACTGACTGAAATGGAAAGTCATCGGGTTATCCAGCACGATAATATTTATGGTGCTGGCCAGTGCCCCCATCGAGCCCAGTCGCCGTTTCTGTACTTCTATACCGAGGAAGAATCCCAGCACTGCCAGGCTAATGc                                                                                     <  2:343008/148‑1 (MQ=255)
                                                                cgcTGTCGAGATACTGACTGAAATGGAAAGTCATCGGGTTATCCAGCACGATAATATTTATGGTGCTGGCCAGTGCCCCCATCGAGCCCAGTCGCCGTTTCTGTACTTCTATACCGAGGAAGAATCCCAGCACTGCCAGGCTAATGca                                                                                    <  2:1385645/148‑1 (MQ=255)
                                                                        agaAACTGACTGAAATGGAAAGTCATCGGGTTATCCAGCACGATAATATTTATGGTGCTGGCCAGTGCCCCCATCGAGCCCAGTCGCCGTTTCTGTACTTCTATACCGAGGAAGAATCCCAGCACTGCCAGGCTAATGCACAgcagca                                                                            <  2:863617/148‑1 (MQ=255)
                                                                         gaAACTGACTGAAATGGAAAGTCATCGGGTTATCCAGCACGATAATATTTATGGTGCTGGCCAGTGCCCCCATCGAGCCCAGTCGCCGTTTCTGTACTTCTATACCGAGGAAGAATCCCAGCACTGCCAGGCTAATGCACAGCAGCaa                                                                           >  2:617637/1‑148 (MQ=255)
                                                                         gaAACTGACTGAAATGGAAAGTCATCGGGTTATCCAGCACGATAATATTTATGGTGCTGGCCAGTGCCCCCATCGAGCCCAGTCGCCGTTTCTGTACTTCTATACCGAGGAAGAATCCCAGCACTGCCAGGCTAATGCACAGCAGCaa                                                                           >  1:813225/1‑148 (MQ=255)
                                                                         gaAACTGACTGAAATGGAAAGTCATCGGGTTATCCAGCACGATAATATTTATGGTGCTGGCCAGTGCCCCCATCGAGCCCAGTCGCCGTTTCTGTACTTCTATACCGAGGAAGAATCCCAGCACTGCCAGGCTAATGCACAGCAGCaa                                                                           <  2:3062242/148‑1 (MQ=255)
                                                                          aaaCTGACTGAAATGGAAAGTCATCGGGTTATCCAGCACGATAATATTTATGGTGCTGGCCAGTGCCCCCATCGAGCCCAGTCGCCGTTTCTGTACTTCTATACCGAGGAAGAATCCCAGCACTGCCAGGCTAATGCACAGCAGCAAc                                                                          >  1:491974/1‑148 (MQ=255)
                                                                           aaCTGACTGAAATGGAAAGTCATCGGGTTATCCAGCACGATAATATTTATGGTGCTGGCCAGTGCCCCCATCGAGCCCAGTCGCCGTTTCTGTACTTCTATACCGAGGAAGAATCCCAGCACTGCCAGGCTAATGCACAGCAGCAACa                                                                         <  2:2055032/148‑1 (MQ=255)
                                                                               gactgaAATGGAAAGTCATCGGGTTATCCAGCACGATAATATTTATGGTGCTGGCCAGTGCCCCCATCGAGCCCAGTCGCCGTTTCTGTACTTCTATACCGAGGAAGAATCCCAGCACTGCCAGGCTAATGCACAGCAGCAACATGct                                                                     <  2:1219263/148‑1 (MQ=255)
                                                                                actgaAATGGAAAGTCATCGGGTTATCCAGCACGATAATATTTATGGTGCTGGCCAGTGCCCCCATCGAGCCCAGTCGCCGTTTCTGTACTTCTATACCGAGGAAGAATCCCAGCACTGCCAGGCTAATGCACAGCAGCAACATGctc                                                                    <  1:400625/148‑1 (MQ=255)
                                                                                actgaAATGGAAAGTCATCGGGTTATCCAGCACGATAATATTTATGGTGCTGGCCAGTGCCCCCATCGAGCCCAGTCGCCGTTTCTGTACTTCTATACCGAGGAAGAATCCCAGCACTGCCAGGCTAATGCACAGCAGCAACATGctc                                                                    <  2:2583994/148‑1 (MQ=255)
                                                                                  tgaAATGGAAAGTCATCGGGTTATCCTGCACGATAATATTTATGGTGCTGGCCAGTGCCCCCATCGAGCCCAGTCGCCGTTTCTGTACTTCTATACCGAGGAAGAATCCCAGCACTGCCAGGCTAATGCACAGCAGCAACATGTTCTg                                                                  <  1:2088861/148‑1 (MQ=255)
                                                                                  tgaAATGGAAAGTCATCGGGTTATCCAGCACGATAATATTTATGGTGCTGGCCAGTGCCCCCATCGAGCCCAGTCGCCGTTTCTGTACTTCTATACCGAGGAAGAATCCCAGCACTGCCAGGCTAATGCACAGCAGCAACCTGCTCTg                                                                  <  2:157738/148‑1 (MQ=255)
                                                                                  tgaAATGGAAAGTCATCGGGTTATCCAGCACGATAATATTTATGGTGCTGGCCAGTGCCCCCATCGAGCCCAGTCGCCGTTTCTGTACTTCTATACCGAGGAAGAATCCCAGCACTGCCAGGCTAATGCACAGCAGCAACATGCTCTg                                                                  <  2:157751/148‑1 (MQ=255)
                                                                                  tgaAATGGAAAGTCATCGGGTTATCCAGCACGATAATATTTATGGTGCTGGCCAGTGCCCCCATCGAGCCCAGTCGCCGTTTCTGTACTTCTATACCGAGGAAGAATCCCAGCACTGCCAGGCTAATGCACAGCAGCAACATGCTCTg                                                                  <  1:650337/148‑1 (MQ=255)
                                                                                   gaAATGGAAAGTCATCGGGTTATCCAGCACGATAATATTTATGGTGCTGGCCAGTGCCCCCATCGAGCCCAGTCGCCGTTTCTGTACTTCTATACCGAGGAAGAATCCCAGCACTGCCAGGCTAATGCACAGCAGCAACATGCTCTGt                                                                 >  2:1780422/1‑148 (MQ=255)
                                                                                   gaAATGGAAAGTCATCGGGTTATCCAGCACGATAATATTTATGGTGCTGGCCAGTGCCCCCATCGAGCCCAGTCGCCGTTTCTGTACTTCTATACCGAGGAAGAATCCCAGCACTGCCAGGCTAATGCACAGCAGCAACATGCTCTGt                                                                 <  2:2083918/148‑1 (MQ=255)
                                                                                   gaAATGGAAAGTCATCGGGTTATCCAGCACGATAATATTTATGGTGCTGGCCAGTGCCCCCATCGAGCCCAGTCGCCGTTTCTGTACTTCTATACCGAGGAAGAATCCCAGCACTGCCAGGCTAATGCACAGCAGCAACATGCTCTGt                                                                 <  2:2084296/148‑1 (MQ=255)
                                                                                    aaaTGGAAAGTCATCGGGTTATCCAGCACGATAATATTTATGGTGCTGGCCAGTGCCCCCATCGAGCCCAGTCGCCGTTTCTGTACTTATATACCGAGGAAGAATCCCAGCACTGCCAGGCTAATGCACAGCAGCAACATGCTCTGtt                                                                <  2:2924942/148‑1 (MQ=255)
                                                                                      aTGGAAAGTCATCGGGTTATCCAGCACGATAATATTTATGGTGCTGGCCAGTGCCCCCATAGAGCCCAGTCGCCGTTTCTGTACTTCTATACCGAGGAAGAATCCCAGCACTGCCAGGCTAATGCACAGCAGCAACATGCTCTGTTgg                                                              <  2:1079274/148‑1 (MQ=255)
                                                                                       tGGAAAGTCATCGGGTTATCCAGCACGATAATATTTATGGTGCTGGCCAGTGCCCCCATCGAGCCCAGTCGCCGTTTCTGTACTTCTATACCGAGGAAGAATCCCAGCACTGCCAGGCTAATGCACAGCAGCAACATGCTCTGTTggg                                                             >  2:1401778/1‑148 (MQ=255)
                                                                                        ggAAAGTCATCGGGTTATCCAGCACGATAATATTTATGGTGCTGGCCAGTGCCCCCATCGAGCCCAGTCGCCGTTTCTGTACTTCTATACCGAGGAAGAATCCCAGCACTGCCAGGCTAATGCACAGCAGCAACATGCTCTGTTGGGt                                                            >  1:2768418/1‑148 (MQ=255)
                                                                                           aaGTCATCGGGTTATCCAGCACGATAATATTTATGGTGCTGGCCAGTGCCCCCATCGAGCCCAGTCGCCGTTTCTGTACTTCTATACCGAGGAAGAATCCCAGCACTGCCAGGCTAATGCACAGCAGCAACATGCTCTGTTGGGTAtt                                                         <  1:588472/148‑1 (MQ=255)
                                                                                            aGTCATCGGGTTATCCAGCACGATAATATTTATGGTGCTGGCCAGTGCCCCCATCGAGCCCAGTCGCCGTTTCTGTACTTCTATACCGAGGAAGAATCCCAGCACTGCCAGGCTAATGCACAGCAGCAACATGCTCTGTTGGGTATTa                                                        >  1:2602720/1‑148 (MQ=255)
                                                                                             gTCATCGGGTTATCCAGCACGATAATATTTATGGTGCTGGCCAGTGCCCCCATCGAGCCCAGTCGCCGTTTCTGTACTTCTATACCGAGGAAGAATCCCAGCACTGCCAGGCTAATGCACAGCAGCAACATGCTCTGTTGGGTATTAg                                                       <  1:40928/148‑1 (MQ=255)
                                                                                               cATCGGGTTATCCAGCACGATAATATTTCTGGTGCTGGCCAGTGCCCCCATAGAGCCCAGTCGCCGTTTCTGTACTTCTATACCGAGGAAGAATCCCAGCACTGCCAGGCTAATGCACAGCAGCAACATGCTCTGTTGGGTATTAggg                                                     <  2:3214143/148‑1 (MQ=255)
                                                                                               cATCGGGTTATCCAGCACGATAATATTTATGGTGCTGGCCAGTGCCCCCATCGAGCCCAGTCGCCGTTTCTGTACTTCTATACCGAGGAAGAATCCCAGCACTGCCAGGCTAATGCACAGCAGCAACATGCTCTGTTGGGTATTAggg                                                     <  2:95296/148‑1 (MQ=255)
                                                                                                aTCGGGTTATCCAGCACGATAATATTTATGGTGCTGGCCAGTGCCCCCATCGAGCCCAGTCGCCGTTTCTGTACTTCTATACCGAGGAAGAATCCCAGCACTGCCAGGCTAATGCACAGCAGCAACATGCTCTGTTGGGTATTAGGGa                                                    <  2:1543009/148‑1 (MQ=255)
                                                                                                aTCGGGTTATCCAGCACGATAATATTTATGGTGCTGGCCAGTGCCCCCATCGAGCCCAGTCGCCGTTTCTGTACTTCTATACCGAGGAAGAATCCCAGCACTGCCAGGCTAATGCACAGCAGCAACATGCTCTGTTGGGTATTAGGGa                                                    <  1:1185267/148‑1 (MQ=255)
                                                                                                 tCGGGTTATCCAGCACGATAATATTTATGGTGCTGGCCAGTGCCCCCATCGAGCCCAGTCGCCGTTTCTGTACTTCTATACCGAGGAAGAATCCCAGCACTGCCAGGCTAATGCACAGCAGCAACATGCTCTGTTGGGTATTAGGGAt                                                   <  1:499580/148‑1 (MQ=255)
                                                                                                  cGGGTTATCCAGCACGATAATATTTATGGTGCTGGCCAGTGCCCCCATCGAGCCCAGTCGCCGTTTCTGTACTTCTATACCGAGGAAGAATCCCAGCACTGCCAGGCTAATGCACAGCAGCAACATGCTCTGTTGGGTATTAGGGATa                                                  <  1:1662746/148‑1 (MQ=255)
                                                                                                  cGGGTTATCCAGCACGATAATATTTATGGTGCTGGCCAGTGCCCCCATCGAGCCCAGTCGCCGTTTCTGTACTTCTATACCGAGGAAGAATCCCAGCACTGCCAGGCTAATGCACAGCAGCAACATGCTCTGTTGGGTATTAGGGATa                                                  >  2:3042527/1‑148 (MQ=255)
                                                                                                  cGGGTTATCCAGCACGATAATATTTATGGTGCTGGCCAGTGCCCCCATCGAGCCCAGTCGCCGTTTCTGTACTTCTATACCGAGGAAGAATCCCAGCACTGCCAGGCTAATGCACAGCAGCAACATGCTCTGTTGGGTATTAGGGATa                                                  >  2:3042197/1‑148 (MQ=255)
                                                                                                  cGGGTTATCCAGCACGATAATATTTATGGTGCTGGCCAGTGCCCCCATCGAGCCCAGTCGCCGTTTCTGTACTTCTATACCGAGGAAGAATCCCAGCACTGCCAGGCTAATGCACAGCAGCAACATGCTCTGTTGGGTATTAGGGATa                                                  >  2:2841452/1‑148 (MQ=255)
                                                                                                  cGGGTTATCCAGCACGATAATATTTATGGTGCTGGCCAGTGCCCCCATCGAGCCCAGTCGCCGTTTCTGTACTTCTATACCGAGGAAGAATCCCAGCACTGCCAGGCTAATGCACAGCAGCAACATGCTCTGTTGGGTATTAGGGATa                                                  >  2:2841373/1‑148 (MQ=255)
                                                                                                   gggTTATCCAGCACGATAATATTTATGGTGCTGGCCAGTGCCCCCATAGAGCCCAGTCGCCGTTTCTGTACTTCTATACCGAGGAAGAATCCCAGCACTGCCAGGCTAATGCACAGCAGCAACATGCTCTGTTGGGTATTAGGGATaa                                                 <  1:2590811/148‑1 (MQ=255)
                                                                                                   gggTTATCCAGCACGATAATATTTATGGTGCTGGCCAGTGCCCCCATAGAGCCCAGTCGCCGTTTCTGTACTTCTATACCGAGGAAGAATCCCAGCACTGCCAGGCTAATGCACAGCAGCAACATGCTCTGTTGGGTATTAGGGATaa                                                 <  1:2590223/148‑1 (MQ=255)
                                                                                                       tATCCAGCACGATAATATTTATGGTGCTGGCCAGTGCCCCCATCGAGCCCAGTCGCCGTTTCTGTACTTCTATACCGAGGAAGAATCCCAGCACTGCCAGGCTAATGCACAGCAGCAACATGCTCTGTTGGGTATTAGGGATAATCAc                                             >  2:3136978/1‑148 (MQ=255)
                                                                                                       tATCCAGCACGATAATATTTATGGTGCTGGCCAGTGCCCCCATCGAGCCCAGTCGCCGTTTCTGTACTTCTATACCGAGGAAGAATCCCAGCACTGCCAGGCTAATGCACAGCAGCAACATGCTCTGTTGGGTATTAGGGATAATCAc                                             >  2:3137045/1‑148 (MQ=255)
                                                                                                         tCCAGCACGATAATATTTATGGTGCTGGCCAGTGCCCCCATCGAGCCCAGTCGCCGTTTGTGTACTTCTATACCGAGGAAGAAAACCAGCACTGCCAGGCTAATGCACAGCAGCAACAAGCTCTGTTGGGTATTATGGCTAATCAccc                                           >  1:8592/1‑147 (MQ=255)
                                                                                                          ccAGCACGATAATATTTATGGTGCTGGCCAGTGCCCCCATAGAGCCCAGTCGCCGTTTCTGTACTTCTATACCGAGGAAGAATCCCAGCACTGCCAGGCTAATGCACAGCAGCAACATGCTCTGTTGGGTATTAGGGATAATCACCaa                                          >  2:2476975/1‑148 (MQ=255)
                                                                                                              cACGATAATATTTATGGTGCTGGCCAGTGCCCCCATCGAGCCCAGTCGCCGTTTCTGTACTTCTATACCGAGGAAGAATCCCAGCACTGCCAGGCTAATGCACAGCAGCAACATGCTCTGTTGGGTATTAGGGATAATCACCAAAAAg                                      <  1:1966788/148‑1 (MQ=255)
                                                                                                                 gATAATATTTATGGTGCTGGCCAGTGCCCCCATCGAGCCCAGTCGCCGTGTCTGTACTTATATACCGAGGAAGAATCCCAGCCCAGCCAGGCTAATGCACAGCAGCAACATGCTCTGTTGGGTATTAGGGATAATCACCAAAAAGTag                                   >  1:1588191/1‑148 (MQ=255)
                                                                                                                  ataataTTTATGGTGCTGGCCAGTGCCCCCATAGAGCCCAGTCGCCGTTTCTGTACTTCTATACCGAGGAAGAATCCCACCACTGCCAGGCTAATGCACAGCAGCCACATGCTCTGTTTGGTATTAGGGATAATCACCCAAAAGTaga                                  >  1:635744/1‑148 (MQ=255)
                                                                                                                   taataTTTATGGTGCTGGCCAGTGCCCCCATCGAGCCCAGTCGCCGTTTCTGTACTTCTATACCGAGGAAGAATCCCAGCACTGCCAGGCTAATGCACAGCAGCAACATGCTCTGTTGGGTATTAGGGATAATCACCAAAAAGTagag                                 >  2:1029252/1‑148 (MQ=255)
                                                                                                                   taataTTTATGGTGCTGGCCAGTGCCCCCATCGAGCCCAGTCGCCGTTTCTGTACTTCTATACCGAGGAAGAATCCCAGCACTGCCAGGCTAATGCACAGCAGCAACATGCTCTGTTGGGTATTAGGGATAATCACCAAAAAGTagag                                 >  2:1649151/1‑148 (MQ=255)
                                                                                                                    aataTTTATGGTGCTGGCCAGTGCCCCCATCGAGCCCAGTCGCCGTTTCTGTACTTCTATACCGAGGAAGAATCCCAGCACTGCCAGGCTAATGCACAGCAGCAACATGCTCTGTTGGGTATTAGGGATAATCACCAAAAAGTAGagc                                >  2:2812048/1‑148 (MQ=255)
                                                                                                                    aataTTTATGGTGCTGGCCAGTGCCCCCATCGAGCCCAGTCGCCGTTTCTGTACTTCTATACCGAGGAAGAATCCCAGCACTGCCAGGCTAATGCACAGCAGCAACATGCTCTGTTGGGTATTAGGGATAATCACCAAAAAGTAGagc                                >  1:2511852/1‑148 (MQ=255)
                                                                                                                     atatTTATGGTGCTGGCCAGTGCCCCCATCGAGCCCAGTCGCCGTTTCTGTACTTCTATACCGAGGAAGAATCCCAGCACTGCCAGGCTAATGCACAGCAGCAACATGCTCTGTTGGGTATTAGGGATAATCACCAAAAAGTAGagca                               >  1:3169506/1‑148 (MQ=255)
                                                                                                                     atatTTATGGTGCTGGCCAGTGCCCCCATCGAGCCCAGTCGCCGTTTCTGTACTTCTATACCGAGGAAGAATCCCAGCACTGCCAGGCTAATGCACAGCAGCAACATGCTCTGTTGGGTATTAGGGATAATCACCAAAAAGTAGagca                               <  1:1780422/148‑1 (MQ=255)
                                                                                                                       atTTATGGTGCTGGCCAGTGCCCCCATCGAGCCCAGTCGCCGTTTCTGTACTTCTATACCGAGGAAGAATCCCAGCACTGCCAGGCTAATGCACAGCAGCAACATGCTCTGTTGGGTATTAGGGATAATCACCAAAAAGTAGagcagc                             <  1:2329265/148‑1 (MQ=255)
                                                                                                                        tttATGGTGCTGGCCAGTGCCCCCATAGAGCCCAGTCGCCGTTTCTGTACTTCTATACCGAGGAAGAATCCCAGCACTGCCAGGCTAATGCACAGCAGCAACATGCTCTGTTGGGTATTAGGGATAATCACCAAAAAGTAGAGCAGcc                            <  2:3292809/148‑1 (MQ=255)
                                                                                                                         ttATGGTGCTGGCCAGTGCCCCCATCGAGCCCAGTCGCCGTTTCTGTACTTCTATACCGAGGAAGAATCCCAGCACTGCCAGGCTAATGCACAGCAGCAACATGCTCTGTTGGGTATTAGGGATAATCACCAAAAAGTAGAGCAGccc                           <  1:988904/148‑1 (MQ=255)
                                                                                                                                gCTGGCCAGTGCCCCCATCGAGCCCAGTCGCCGTTTCTGTACTTCTATACCGAGGAAGAATCCCAGCACTGCCAGGCTAATGCACAGCAGCAACATGCTCTGTTGGGTATTAGGGATAATCACCAAAAAGTAGAGCAGCCCTAACGGc                    <  1:333794/148‑1 (MQ=255)
                                                                                                                                gCTGGCCAGTGCCCCCATCGAGCCCAGTCGCCGTTTCTGTACTTCTATACCGAGGAAGAATCCCAGCACTGCCAGGCTAATGCACAGCAGCAACATGCTCTGTTGGGTATTAGGGATAATCACCAAAAAGTAGAGCAGCCCTAACGGc                    >  2:449376/1‑148 (MQ=255)
                                                                                                                                  tGGCCAGTGCCCCCATCGAGCCCAGTCGCCGTTTCTGTACTTCTATACCGAGGAAGAATCCCAGCACTGCCAGGCTAATGCACAGCAGCAACATGCTCTGTTGGGTATTAGGGATAATCACCAAAAAGTAGAGCAGCCCTAACGGCAg                  <  1:1393828/148‑1 (MQ=255)
                                                                                                                                    gCCAGTGCCCCCATCGAGCCCAGTCGCCGTTTCTGTACTTCTATACCGAGGAAGAATCCCAGCACTGCCAGGCTAATGCACAGCAGCAACATGCTCTGTTGGGTATTAGGGATAATCACCAAAAAGTAGAGCAGCCCTAACGGCAgcg                <  1:2778939/148‑1 (MQ=255)
                                                                                                                                        gTGCCCCCATCGAGCCCTGTCGCCGTTTCTGTACTTCTATACCGAGGAAGAATCCCAGCACTGCCAGGCTAATGCACAGCAGCAACATGCTCTGTTGGGTATTAGGGATAATCACCAAAAAGTAGAGCAGCCCTAACGGCAGCGCGGc            <  1:3050169/148‑1 (MQ=255)
                                                                                                                                               cATAGAGCCCAGTCGCCGTTTCTGTACTTCTATACCGAGGAAGAATCCCAGCACTGCCAGGCTAATGCACAGCAGCAACATGCTCTGTTGGGTATTAGGGATAATCACCAAAAAGTAGAGCAGCCCTAACGGCAGCGCGGCCAGCGTc     >  2:779411/1‑148 (MQ=255)
                                                                                                                                                aTCGAGCCCAGTCGCCGTTTCTGTACTTCTATACCGAGGAAGAATCCCAGCACTGCCAGGCTAATGCACAGCAGCAACATGCTCTGTTGGGTATTAGGGATAATCACCAAAAAGTAGAGCAGCCCTAACGGCAGCGCGGCCAGCGTcc    >  1:234634/1‑148 (MQ=255)
                                                                                                                                                  cGAGCCCAGTCGCCGTTTCTGTACTTCTATAACGAGGAAGAATACCAGCACTGCCACGCTAATGCACAGCAGCAACAAGCTCTGTTGGGTATTATGGATAATCACCAAAAAGTAGAGCAGCCCAAACGGCAGCGCGGCCAGCGTCCCg  >  1:1192770/1‑148 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
GATTTATCCCGCACCAGCAAAATAACGGTGAACGCGAGCACACAGCCGACGATTTGCCCTAATGCGCTGTCGAGAAACTGACTGAAATGGAAAGTCATCGGGTTATCCAGCACGATAATATTTATGGTGCTGGCCAGTGCCCCCATCGAGCCCAGTCGCCGTTTCTGTACTTCTATACCGAGGAAGAATCCCAGCACTGCCAGGCTAATGCACAGCAGCAACATGCTCTGTTGGGTATTAGGGATAATCACCAAAAAGTAGAGCAGCCCTAACGGCAGCGCGGCCAGCGTCCCG  >  CP009273/3380028‑3380321

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: ‑