Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,380,174 C→A 30.8% 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,1740CA30.8% 630.4 / 267.0 325S458S (TCG→TCTaaeBp‑hydroxybenzoic acid efflux system component
Reads supporting (aligned to +/- strand):  ref base C (115/108);  new base A (56/45);  total (171/154)
Fisher's exact test for biased strand distribution p-value = 5.49e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 329 total aligned reads displayed.

CGCGATTTATCCCGCACCAGCAAAATAACGGTGAACGCGAGCACACAGCCGACGATTTGCCCTAATGCGCTGTCGAGAAACTGACTGAAATGGAAAGTCATCGGGTTATCCAGCACGATAATATTTATGGTGCTGGCCAGTGCCCCCATCGAGCCCAGTCGCCGTTTCTGTACTTCTATACCGAGGAAGAATCCCAGCACTGCCAGGCTAATGCACAGCAGCAACATGCTCTGTTGGGTATTAGGGATAATCACCAAAAAGTAGAGCAGCCCTAACGGCAGCGCGGCCAGCGTCCCGT  >  CP009273/3380025‑3380322
                                                                                                                                                     |                                                                                                                                                    
cgcgATTTATCCCGCACCAGCAAAATAACGGTGAACGCGAGCACACAGCCGACGATTTGCCCTAATGCGCTGTCGAGAAACTGACTGAAATGGAAAGTCATCGGGTTATCCAGCACGATAATATTTATGGTGCTGGCCAGTGCCCCCATc                                                                                                                                                      >  1:934858/1‑150 (MQ=255)
     tttATCCCGCACCAGCAAAATAACGGTGAACGCGAGCACACAGCCGACGATTTGCCCTAATGCGCTGTCGAGAAACTGACTGAAATGGAAAGTCATCGGGTTATCCAGCACGATAATATTTATGGTGCTGGCCAGTGCCCCCATCGAGcc                                                                                                                                                 >  1:586955/1‑150 (MQ=255)
        aTCCCGCACCAGCAAAATAACGGTGAACGCGAGCACACAGCCGACGATTTGCCCTAATGCGCTGTCGAGAAACTGACTGAAATGGAAAGTCATCGGGTTATCCAGCACGATAATATTTATGGTGCTGGCCAGTGCCCCCATCGAGCCCAg                                                                                                                                              >  2:89279/1‑150 (MQ=255)
         tCCCGCACCAGCAAAATAACGGTGAACGCGAGCACACAGCCGACGATTTGCCCTAATGCGCTGTCGAGAAACTGACTGAAATGGAAAGTCATCGGGTTATCCAGCACGATAATATTTATGGTGCTGGCCAGTGCCCCCATCGAGCCCAGt                                                                                                                                             <  1:1748953/150‑1 (MQ=255)
         tCCCGCACCAGCAAAATAACGGTGAACGCGAGCACACAGCCGACGATTTGCCCTAATGCGCTGTCGAGAAACTGACTGAAATGGAAAGTCATCGGGTTATCCAGCACGATAATATTTATGGTGCTGGCCAGTGCCCCCATAGAGCCCAGt                                                                                                                                             <  1:3254010/150‑1 (MQ=255)
          cccGCACCAGCAAAATAACGGTGAACGCGAGCACACAGCCGACGATTTGCCCTAATGCGCTGTCGAGAAACTGACTGAAATGGAAAGTCATCGGGTTATCCAGCACGATAATATTTATGGTGCTGGCCAGTGCCCCCATCGAGCCCAGTc                                                                                                                                            <  2:4339733/150‑1 (MQ=255)
            cGCACCAGCAAAATAACGGTGAACGCGAGCACACAGCCGACGATTTGCCCTAATGCGCTGTCGAGAAACTGACTGAAATGGAAAGTCATCGGGTTATCCAGCACGATAATATTTATGGTGCTGGCCAGTGCCCCCATCGAGCCCAGTCGc                                                                                                                                          <  2:3370474/150‑1 (MQ=255)
            cGCACCAGCAAAATAACGGTGAACGCGAGCACACAGCCGACGATTTGCCCTAATGCGCTGTCGAGAAACTGACTGAAATGGAAAGTCATCGGGTTATCCAGCACGATAATATTTATGGTGCTGGCCAGTGCCCCCATAGAGCCCAGTCGc                                                                                                                                          >  1:2998110/1‑150 (MQ=255)
               aCCAGCAAAATAACGGTGAACGCGAGCACACAGCCGACGATTTGCCCTAATGCGCTGTCGAGAAACTGACTGAAATGGAAAGTCATCGGGTTATCCAGCACGATAATATTTATGGTGCTGGCCAGTGCCCCCATCGAGCCCAGTCGCCGt                                                                                                                                       >  1:3178327/1‑150 (MQ=255)
                ccAGCAAAATAACGGTGAACGCGAGCACACAGCCGACGATTTGCCCTAATGCGCTGTCGAGAAACTGACTGAAATGGAAAGTCATCGGGTTATCCAGCACGATAATATTTATGGTGCTGGCCAGTGCCCCCATCGAGCCCAGTCGCCGtt                                                                                                                                      >  2:1792451/1‑150 (MQ=255)
                ccAGCAAAATAACGGTGAACGCGAGCACACAGCCGACGATTTGCCCTAATGCGCTGTCGAGAAACTGACTGAAATGGAAAGTCATCGGGTTATCCAGCACGATAATATTTATGGTGCTGGCCAGTGCCCCCATCGAGCCCAGTCGCCGtt                                                                                                                                      >  2:1733778/1‑150 (MQ=255)
                 cAGCAAAATAACGGTGAACGCGAGCACACAGCCGACGATTTGCCCTAATGCGCTGTCGAGAAACTGACTGAAATGGAAAGTCATCGGGTTATCCAGCACGATAATATTTATGGTGCTGGCCAGTGCCCCCATAGAGCCCAGTCGCCGttt                                                                                                                                     <  1:1257308/150‑1 (MQ=255)
                  aGCAAAATAACGGTGAACGCGAGCACACAGCCGACGATTTGCCCTAATGCGCTGTCGAGAAACTGACTGAAATGGAAAGTCATCGGGTTATCCAGCACGATAATATTTATGGTGCTGGCCAGTGCCCCCATAGAGCCCAGTCGCCGTTTc                                                                                                                                    >  1:4362540/1‑150 (MQ=255)
                  aGCAAAATAACGGTGAACGCGAGCACACAGCCGACGATTTGCCCTAATGCGCTGTCGAGAAACTGACTGAAATGGAAAGTCATCGGGTTATCCAGCACGATAATATTTATGGTGCTGGCCAGTGCCCCCATAGAGCCCAGTCGCCGTTTc                                                                                                                                    >  2:2440089/1‑150 (MQ=255)
                     aaaaTAACGGTGAACGCGAGCACACAGCCGACGATTTGCCCTAATGCGCTGTCGAGAAACTGACTGAAATGGAAAGTCATCGGGTTATCCAGCACGATAATATTTATGGTGCTGGCCAGTGCCCCCATCGAGCCCAGTCGCCGTTTCTGt                                                                                                                                 >  1:1533099/1‑150 (MQ=255)
                     aaaaTAACGGTGAACGCGAGCACACAGCCGACGATTTGCCCTAATGCGCTGTCGAGAAACTGACTGAAATGGAAAGTCATCGGGTTATCCAGCACGATAATATTTATGGTGCTGGCCAGTGCCCCCATCGAGCCCAGTCGCCGTTTCTGt                                                                                                                                 >  2:4270016/1‑150 (MQ=255)
                     aaaaTAACGGTGAACGCGAGCACACAGCCGACGATTTGCCCTAATGCGCTGTCGAGAAACTGACTGAAATGGAAAGTCATCGGGTTATCCAGCACGATAATATTTATGGTGCTGGCCAGTGCCCCCATCGAGCCCAGTCGCCGTTTCTGt                                                                                                                                 >  2:3926103/1‑150 (MQ=255)
                     aaaaTAACGGTGAACGCGAGCACACAGCCGACGATTTGCCCTAATGCGCTGTCGAGAAACTGACTGAAATGGAAAGTCATCGGGTTATCCAGCACGATAATATTTATGGTGCTGGCCAGTGCCCCCATCGAGCCCAGTCGCCGTTTCTGt                                                                                                                                 >  1:536653/1‑150 (MQ=255)
                     aaaaTAACGGTGAACGCGAGCACACAGCCGACGATTTGCCCTAATGCGCTGTCGAGAAACTGACTGAAATGGAAAGTCATCGGGTTATCCAGCACGATAATATTTATGGTGCTGGCCAGTGACCCAATCGAGCCCAGTCTCCGTTTCTGt                                                                                                                                 >  2:4551890/1‑150 (MQ=255)
                      aaaTAACGGTGAACGCGAGCACACAGCCGACGATTTGCCCTAATGCGCTGTCGAGAAACTGACTGAAATGGAAAGTCATCGGGTTATCCAGCACGATAATATTTATGGTGCTGGCCAGTGCCCCCATCGAGCCCAGTCGCCGTTTCTGTa                                                                                                                                >  1:3557068/1‑150 (MQ=255)
                        aTAACGGTGAACGCGAGCACACAGCCGACGATTTGCCCTAATGCGCTGTCGAGAAACTGACTGAAATGGAAAGTCATCGGGTTATCCAGCACGATAATATTTATGGTGCTGGCCAGTGCCCCCATAGAGCCCAGTCGCCGTTTCTGTACt                                                                                                                              <  2:3447612/150‑1 (MQ=255)
                                      gAGCACACAGCCGACGATTTGCCCTAATGCGCTGTCGAGAAACTGACTGAAATGGAAAGTCATCGGGTTATCCAGCACGATAATATTTATGGTGCTGGCCAGTGCCCCCATCGAGCCCAGTCGCCGTTTCTGTACTTCTATACCGAGgaa                                                                                                                >  2:3153206/1‑150 (MQ=255)
                                       aGCACACAGCCGACGATTTGCCCTAATGCGCTGTCGAGAAACTGACTGAAATGGAAAGTCATCGGGTTATCCAGCACGATAATATTTATGGTGCTGGCCAGTGCCCCCATCGAGCCCAGTCGCCGTTTCTGTACTTCTATACCGAGgaag                                                                                                               >  2:3206599/1‑150 (MQ=255)
                                        gCACACAGCCGACGATTTGCCCTAATGCGCTGTCGAGAAACTGACTGAAATGGAAAGTCATCGGGTTATCCAGCACGATAATATTTATGGTGCTGGCCAGTGCCCCCATCGAGCCCAGTCGCCGTTTCTGTACTTCTATACCGAGgaaga                                                                                                              >  2:1389254/1‑150 (MQ=255)
                                          acacaGCCGACGATTTGCCCTAATGCGCTGTCGAGAAACTGACTGAAATGGAAAGTCATCGGGTTATCCAGCACGATAATATTTATGGTGCTGGCCAGTGCCCCCATCGAGCCCAGTCGCCGTTTCTGTACTTCTATACCGAGGAAGAAt                                                                                                            <  2:3175502/150‑1 (MQ=255)
                                          acacaGCCGACGATTTGCCCTAATGCGCTGTCGAGAAACTGACTGAAATGGAAAGTCATCGGGTTATCCAGCACGATAATATTTATGGTGCTGGCCAGTGCCCCCATCGAGCCCAGTCGCCGTTTCTGTACTTCTATACCGAGGAAGAAt                                                                                                            <  2:938745/150‑1 (MQ=255)
                                              aGCCGACGATTTGCCCTAATGCGCTGTCGAGAAACTGACTGAAATGGAAAGTCATCGGGTTATCCAGCACGATAATATTTATGGTGCTGGCCAGTGCCCCCATAGAGCCCAGTCGCCGTTTCTGTACTTCTATACCGAGGAAGAATCCCa                                                                                                        >  2:982235/1‑150 (MQ=255)
                                                ccGACGATTTGCCCTAATGCGCTGTCGAGAAACTGACTGAAATGGAAAGTCATCGGGTTATCCAGCACGATAATATTTATGGTGCTGGCCAGTGCCCCCATCGAGCCCAGTCGCCGTTTCTGTACTTCTATACCGAGGAAGAATCCCAGc                                                                                                      >  1:4290881/1‑150 (MQ=255)
                                                  gacgaTTTGCCCTAATGCGCTGTCGAGAAACTGACTGAAATGGAAAGTCATCGGGTTATCCAGCACGATAATATTTATGGTGCTGGCCAGTGCCCCCATCGAGCCCAGTCGCCGTTTCTGTACTTCTATACCGAGGAAGAATCCCAGCAc                                                                                                    >  1:4849370/1‑150 (MQ=255)
                                                   acgaTTTGCCCTAATGCGCTGTCGAGAAACTGACTGAAATGGAAAGTCATCGGGTTATCCAGCACGATAATATTTATGGTGCTGGCCAGTGCCCCCATCGAGCCCAGTCGCCGTTTCTGTACTTCTATACCGAGGAAGAATCCCAGCACt                                                                                                   <  2:4172419/150‑1 (MQ=255)
                                                     gaTTTGCCCTAATGCGCTGTCGAGAAACTGACTGAAATGGAAAGTCATCGGGTTATCCAGCACGATAATATTTATGGTGCTGGCCAGTGCCCCCATCGAGCCCAGTCGCCGTTTCTGTACTTCTATACCGAGGAAGAATCCCAGCACTGc                                                                                                 >  2:4994773/1‑150 (MQ=255)
                                                      aTTTGCCCTAATGCGCTGTCGAGAAACTGACTGAAATGGAAAGTCATCGGGTTATCCAGCACGATAATATTTATGGTGCTGGCCAGTGCCCCCATAGAGCCCAGTCGCCGTTTCTGTACTTCTATACCGAGGAAGAATCCCAGCACTGcc                                                                                                <  1:1968877/150‑1 (MQ=255)
                                                        ttGCCCTAATGCGCTGTCGAGAAACTGACTGAAATGGAAAGTCATCGGGTTATCCAGCACGATAATATTTATGGTGCTGGCCAGTGCCCCCATCGAGCCCAGTCGCCGTTTCTGTACTTCTATACCGAGGAAGAATCCCAGCACTGCCAg                                                                                              <  1:4984945/150‑1 (MQ=255)
                                                          gCCCTAATGCGCTGTCGAGAAACTGACTGAAATGGAAAGTCATCGGGTTATCCAGCACGATAATATTTATGGTGCTGGCCAGTGCCCCCATAGAGCCCAGTCGCCGTTTCTGTACTTCTATACCGAGGAAGAATCCCAGCACTGCCAGGc                                                                                            >  2:74530/1‑150 (MQ=255)
                                                             cTAATGCGCTGTCGAGAAACTGACTGAACTGGAAAGTCATCGGGTTATCCAGCACGATAATATTTATGGTGCTGGCCAGTGCCCCCATCGAGCCCAGTCGCCGTTTCTGTACTTCTATACCGAGGAAGAATCCCAGCACTGCCAGGCTaa                                                                                         <  2:2261548/150‑1 (MQ=255)
                                                             cTAATGCGCTGTCGAGAAACTGACTGAAATGGAAAGTCATCGGGTTATCCAGCACGATAATATTTATGGTGCTGGCCAGTGCCCCCATCGAGCCCAGTCGCCGTTTCTGTACTTCTATACCGAGGAAGAATCCCAGCACTGCCAGGCTaa                                                                                         <  1:3895900/150‑1 (MQ=255)
                                                              tAATGCGCTGTCGAGAAACTGACTGAAATGGAAAGTCATCGGGTTATCCAGCACGATAATATTTATGGTGCTGGCCAGTGCCCCCATAGAGCCCAGTCGCCGTTTCTGTACTTCTATACCGAGGAAGAATCCCAGCACTGCCAGGCTAAt                                                                                        <  2:1698691/150‑1 (MQ=255)
                                                               aaTGCGCTGTCGAGAAACTGACTGAAATGGAAAGTCATCGGGTTATCCAGCACGATAATATTTATGGTGCTGGCCAGTGCCCCCATAGAGCCCAGTCGCCGTTTCTGTACTTCTATACCGAGGAAGAATCCCAGCACTGCCAGGCTAATg                                                                                       >  2:2324980/1‑150 (MQ=255)
                                                                aTGCGCTGTCGAGAAACTGACTGAAATGGAAAGTCATCGGGTTATCCAGCACGATAATATTTATGGTGCTGGCCAGTGCCCCCATCGAGCCCAGTCGCCGTTTCTGTACTTCTATACCGAGGAAGAATCCCAGCACTGCCAGGCTAATGc                                                                                      <  2:2430009/150‑1 (MQ=255)
                                                                 tGCGCTGTCGAGAAACTGACTGAAATGGAAAGTCATCGGGTTATCCAGCACGATAATATTTATGGTGCTGGCCAGTGCCCCCATCGAGCCCAGTCGCCGTTTCTGTACTTCTATACCGAGGAAGAATCCCAGCACTGCCAGGCTAATGca                                                                                     <  1:908371/150‑1 (MQ=255)
                                                                   cgcTGTCGAGAAACTGACTGAAATGGAAAGTCATCGGGTTATCCAGCACGATAATATTTATGGTGCTGGCCAGTGCCCCCATCGAGCCCAGTCGCCGTTTCTGTACTTCTATACCGAGGAAGAATCCCAGCACTGCCAGGCTAATGcaca                                                                                   >  2:3311982/1‑150 (MQ=255)
                                                                   cgcTGTCGAGAAACTGACTGAAATGGAAAGTCATCGGGTTATCCAGCACGATAATATTTATGGTGCTGGCCAGTGCCCCCATCGAGCCCAGTCGCCGTTTCTGTACTTCTATACCGAGGAAGAATCCCAGCACTGCCAGGCTAATGcaca                                                                                   <  1:1376099/150‑1 (MQ=255)
                                                                     cTGTCGAGAAACTGACTGAAATGGAAAGTCATCGGGTTATCCAGCACGATAATATTTATGGTGCTGGCCAGTGCCCCCATAGAGCCCAGTCGCCGTTTCTGTACTTCTATACCGAGGAAGAATCCCAGCACTGCCAGGCTAATGCAcagc                                                                                 >  1:2855207/1‑150 (MQ=255)
                                                                       gTCGAGAAACTGACTGAAATGGAAAGTCATCGGGTTATCCAGCACGATAATATTTATGGTGCTGGCCAGTGCCCCCATCGAGCCCAGTCGCCGTTTCTGTACTTCTATACCGAGGAAGAATCCCAGCACTGCCAGGCTAATGCAcagcag                                                                               >  1:2308853/1‑150 (MQ=255)
                                                                       gTCGAGAAACTGACTGAAATGGAAAGTCATCGGGTTATCCAGCACGATAATATTTATGGTGCTGGCCAGTGCCCCCATAGAGCCCAGTCGCCGTTTCTGTACTTCTATACCGAGGAAGAATCCCAGCACTGCCAAGCTAATGCAcagcag                                                                               >  2:1836326/1‑150 (MQ=255)
                                                                         cGAGAAACTGACTGAAATGGAAAGTCATCGGGTTATCCAGCACGATAATATTTATGGTGCTGGCCAGTGCCCCCATCGAGCCCAGTCGCCGTTTCTGTACTTCTATACCGAGGAAGAATCCCAGCACTGCCAGGCTAATGCACAgcagca                                                                             >  2:2957505/1‑150 (MQ=255)
                                                                         cGAGAAACTGACTGAAATGGAAAGTCATCGGGTTATCCAGCACGATAATATTTATGGTGCTGGCCAGTGCCCCCATCGAGCCCAGTCGCCGTTTCTGTACTTCTATACCGAGGAAGAATCCCAGCACTGCCAGGCTAATGCACAgcagca                                                                             >  1:2687391/1‑150 (MQ=255)
                                                                          gagaAACTGACTGAAATGGAAAGTCATCGGGTTATCCAGCACGATAATATTTATGGTGCTGGCCAGTGCCCCCATAGAGCCCAGTCGCCGTTTCTGTACTTCTATACCGAGGAAGAATCCCAGCACTGCCAGGCTAATGCACAGCAGCaa                                                                            >  1:3285628/1‑150 (MQ=255)
                                                                          gagaAACTGACTGAAATGGAAAGTCATCGGGTTATCCAGCACGATAATATTTATGGTGCTGGCCAGTGCCCCCATAGAGCCCAGTCGCCGTTTCTGTACTTCTATACCGAGGAAGAATCCCAGCACTGCCAGGCTAATGCACAGCAGCaa                                                                            <  1:1163592/150‑1 (MQ=255)
                                                                          gagaAACTGACTGAAATGGAAAGTCATCGGGTTATCCAGCACGATAATATTTATGGTGCTGGCCAGTGCCCCCATAGAGCCCAGTCGCCGTTTCTGTACTTCTATACCGAGGAAGAATCCCAGCACTGCCAGGCTAATGCACAGCAGCaa                                                                            >  2:2890221/1‑150 (MQ=255)
                                                                            gaAACTGACTGAAATGGAAAGTCATCGGGTTATCCAGCACGATAATATTTATGGTGCTGGCCAGTGCCCCCATAGAGCCCAGTCGCCGTTTCTGTACTTCTATACCGAGGAAGAATACCAGCACTGACAGGCTAATGCACAGCAGCAACa                                                                          >  1:2385333/1‑150 (MQ=255)
                                                                                 tgactgaAATGGAAAGTCATCGGGTTATCCAGCACGATAATATTTATGGTGCTGGCCAGTGCCCCCATCGAGCCCAGTCGCCGTTTCTGTACTTCTATACCGAGGAAGAATCCCAGCACTGCCAGGCTAATGCACAGCAGCAACATGctc                                                                     >  2:4260776/1‑150 (MQ=255)
                                                                                 tgactgaAATGGAAAGTCATCGGGTTATCCAGCACGATAATATTTATGGTGCTGGCCAGTGCCCCCATCGAGCCCAGTCGCCGTTTCTGTACTTCTATACCGAGGAAGAATCCCAGCACTGCCAGGCTAATGCACAGCAGCAACATGctc                                                                     >  2:4970846/1‑150 (MQ=255)
                                                                                   actgaAATGGAAAGTCATCGGGTTATCCAGCACGATAATATTTATGGTGCTGGCCAGTGCCCCCATCGAGCCCAGTCGCCGTTTCTGTACTTCTATACCGAGGAAGAATCCCAGCACTGCCAGGCTAATGCACAGCAGCAACATGCTCTg                                                                   >  2:2924189/1‑150 (MQ=255)
                                                                                    ctgaAATGGAAAGTCATCGGGTTATCCAGCACGATAATATTTATGGTGCTGGCCAGTGCCCCCATCGAGCCCAGTCGCCGTTTCTGTACTTCTATACCGAGGAAGAATCCCAGCACTGCCAGGCTAATGCACAGCAGCAACATGCTATGt                                                                  <  1:3210954/150‑1 (MQ=255)
                                                                                      gaAATGGAAAGTCATCGGGTTATCCAGCACGATAATCTTTATGGTGCTGGCCAGTGCCCCCATCGAGCCCAGTCGCCGTTTCTGTACTTCTATCCCGAGGAAGAATCCCAGCACTGCCAGGCTAATGCACAGCAGCAACATGCTCTGTTg                                                                >  2:2637860/1‑150 (MQ=255)
                                                                                      gaAATGGAAAGTCATCGGGTTATCCAGCACGATAATATTTATGGTGCTGGCCAGTGCCCCCATAGAGCCCAGTCGCCGTTTCTGTACTTCTATACCGAGGAAGAATCCCAGCACTGCCAGGCTAATGCACAGCAGCAACATGCTCTGTTg                                                                <  1:3417941/150‑1 (MQ=255)
                                                                                       aaaTGGAAAGTCATCGGGTTATCCAGCACGATAATATTTATGGTGCTGGCCAGTGCCCCCATCGAGCCCAGTCGCCGTTTCTGTACTTCTATACCGAGGAAGAATCCCAGCACTGCCAGGCTAATGCACAGCAGCAACATGCTCTGTTgg                                                               <  1:736051/150‑1 (MQ=255)
                                                                                       aaaTGGAAAGTCATCGGGTTATCCAGCACGATAATATTTATGGTGCTGGCCAGTGCCCCCATAGAGCCCAGTCGCCGTTTCTGTACTTCTATACCGAGGAAGAATCCCAGCACTGCCAGGCTAATGCACAGCAGCAACATGCTCTGTTgg                                                               >  2:254676/1‑150 (MQ=255)
                                                                                         aTGGAAAGTCATCGGGTTATCCAGCACGATAATATTTATGGTGCTGGCCAGTGCCCCCATCGAGCCCAGTCGCCGTTTCTGTACTTCTATACCGAGGAAGAATCCCAGCACTGCCAGGCTAATGCACAGCAGCAACATGCTCTGTTGGGt                                                             <  2:1021967/150‑1 (MQ=255)
                                                                                          tGGAAAGTCATCGGGTTATCCAGCACGATAATATTTATGGTGCTGGCCAGTGCCCCCATAGAGCCCAGTCGCCGTTTCTGTACTTCTATACCGAGGAAGAATCCCAGCACTGCCAGGCTAATGCACAGCAGCAACATGCTCTGTTGGGTa                                                            <  2:4698084/150‑1 (MQ=255)
                                                                                           ggAAAGTCATCGGGTTATCCAGCACGATAATATTTATGGTGCTGGCCAGTGCCCCCATCGAGCCCAGTCGCCGTTTCTGTACTTCTATACCGAGGAAGAATCCCAGCACTGCCAGGCTAATGCACAGCAGCAACATGCTCTGTTGGGTAt                                                           <  2:531963/150‑1 (MQ=255)
                                                                                           ggAAAGTCATCGGGTTATCCAGCACGATAATATTTATGGTGCTGGCCAGTGCCCCCATAGAGCCCAGTCGCCGTTTCTGTACTTCTATACCGAGGAAGAATCCCAGCACTGCCAGGCTAATGCACAGCAGCAACATGCTCTGTTGGGTAt                                                           >  1:348245/1‑150 (MQ=255)
                                                                                            gAAAGTCATCGGGTTATCCAGCACGATAATATTTATGGTGCTGGCCAGTGCCCCCATCGAGCCCAGTCGCCGTTTCTGTACTTCTATACCGAGGAAGAATCCCAGCACTGCCAGGCTAATGCACAGCAGCAACATGCTCTGTTGGGTAtt                                                          >  1:3331985/1‑150 (MQ=255)
                                                                                              aaGTCATCGGGTTATCCAGCACGATAATATTTATGGTGCTGGCCAGTGCCCCCATAGAGCCCAGTCGCCGTTTCTGTACTTCTATACCGAGGAAGAATCCCAGCACTGCCAGGCTAATGCACAGCAGCAACATGCTCTGTTGGGTATTAg                                                        <  1:4817536/150‑1 (MQ=255)
                                                                                               aGTCATCGGGTTATCCAGCACGATAATCGTGATGGTGCTGGCCCGTGCCCCCCTAGCGCCCCGTCGCCGTTTCTGTACTTCGCGCCCGAGGACGACTCCCCGCACTGCCAGGCTAATGCACAGCAGCCACAGGCTCTGTTGGGTATTAgg                                                       <  2:1574664/150‑1 (MQ=255)
                                                                                               aGTCATCGGGTTATCCAGCACGATAATATTTATGGTGCTGGCCAGTGCCCCCATCGAGCCCAGTCGCCGTTTCTGTACTTCTATACCGAGGAAGAATCCCAGCACTGCCAGGCTAATGCACAGCAGCAACATGCTCTGTTGGGTATTAgg                                                       <  1:39868/150‑1 (MQ=255)
                                                                                                gTCATCGGGTTATCCAGCACGATAATATTTATGGTGCTGGCCAGTGCCCCCATCGAGCCCAGTCGCCGTTTCTGTACTTCTATACCGAGGAAGAATCCCAGCACTGCCAGGCTAATGCACAGCAGCAACATGCTCTGTTGGGTATTAggg                                                      >  1:1991563/1‑150 (MQ=255)
                                                                                                gTCATCGGGTTATCCAGCACGATAATATTTATGGTGCTGGCCAGTGCCCCCATCGAGCCCAGTCGCCGTTTCTGTACTTCTATACCGAGGAAGAATCCCAGCACTGCCAGGCTAATGCACAGCAGCAACATGCTCTGTTGGGTATTAggg                                                      <  1:2083967/150‑1 (MQ=255)
                                                                                                gTCATCGGGTTATCCAGCACGATAATATTTATGGTGCTGGCCAGTGCCCCCATAGAGCCCAGTCGCCGTTTCTGTACTTCTATACCGAGGAAGAATCCCAGCACTGCCAGGCTAATGCACAGCAGCAACATGCTCTGTTGGGTATTAtgg                                                      <  2:512757/150‑1 (MQ=255)
                                                                                                 tCATCGGGTTATCCAGCACGATAATATTTATGGTGCTGGCCAGTGCCCCCATAGAGCCCAGTCGCCGTTTCTGTACTTCTATACCGAGGACGAATCCCAGCACTGCCAGGCTAATGCACAGCAGCAACATGCTCTGTTGGGTATTAGGGa                                                     >  1:2767250/1‑150 (MQ=255)
                                                                                                  cATCGGGTTATCCAGCACGATAATATTTAGGGTGCTGGCCAGTGCCCCCATCGAGCCCAGTCGCCGTTTCTGTACTTCTATACCGAGGAAGAATCCCAGCACTGCCAGGCTAATGCACAGCAGCAACATGCTCTGTTGGGTATTAGGGAt                                                    <  2:3773001/150‑1 (MQ=255)
                                                                                                  cATAGGGTTATCCAGCACGATAATATTTATGGTGCTGGCCAGTGCCCCCATAGAGCCCAGTCGCCGTTTCTGTACTTCTATACCGAGGAAGAATCCCAGCACTGCCAGGCTAATGCACAGCAGCAACATGCTCTGTTGGGTATTAGGGAt                                                    <  2:4855214/150‑1 (MQ=255)
                                                                                                     cGGGTTATCCAGCACGATAATATTTATGGTGCTGGCCAGTGCCCCCATAGAGCCCAGTCGCCGTTTCTGTACTTCTATACCGAGGAAGAATCCCAGCACTGCCAGGCTAATGCACAGCAGCAACATGCTCTGTTGGGTATTAGGGATAAt                                                 <  1:1027188/150‑1 (MQ=255)
                                                                                                       ggTTATCCAGCACGATAATATTTATGGTGCTGGCCAGTGCCCCCATCGAGCCCAGTCGCCGTTTCTGTACTTCTATACCGAGGAAGAATCCCAGCACTGCCAGGCTAATGCACAGCAGCAACATGCTCTGTTGGGTATTAGGGATAATCa                                               >  2:1165975/1‑150 (MQ=255)
                                                                                                        gTTATCCAGCACGATAATATTTATGGTGCTGGCCAGTGCCCCCATCGAGCCCAGTCGCCGTTTCTGTACTTCTATACCGAGGAAGAATCCCAGCACTGCCAGGCTAATGCACAGCAGCAACATGCTCTGTTGGGTATTAGGGATAATCAc                                              <  2:2119236/150‑1 (MQ=255)
                                                                                                            tCCAGCACGATAATATTTATGGTGCTGGCCAGTGCCCCCATAGAGCCCAGTCGCCGTTTCTGTACTTCTATACCGAGGAAGAATCCCAGCACTGCCAGGCTAATGCACAGCAGCAACATGCTCTGTTGGGTATTAGGGATAATCACCaaa                                          <  1:3691056/150‑1 (MQ=255)
                                                                                                              cAGCACGATAATATTTATGGTGCTGGCCAGTGCCCCCATCGAGCCCAGTCGCCGTTTCTGTACTTCTATACCGAGGAAGAATCCCAGCACTGCCAGGCTAATGCACAGCAGCAACATGCTCTGTTGGGTATTAGGGATAATCACCaaaaa                                        >  1:4593766/1‑150 (MQ=255)
                                                                                                              cAGCACGATAATATTTATGGTGCTGGCCAGTGCCCCCATCGAGCCCAGTCGCCGTTTCTGTACTTCTATACCGAGGAAGAATCCCAGCACTGCCAGGCTAATGCACAGCAGCAACATGCTCTGTTGGGTATTAGGGATAATCACCaaaaa                                        >  1:352712/1‑150 (MQ=255)
                                                                                                                 cACGATAATAATTATGGTGCTGGCCAGTGCCCCCATCGAGCCCAGTCGCCGTTTCTGTACTTCTATACCGAGGAAGAATCCCAGCACTGCCAGGCTAATGCACAGCAGCAACATGCTCTGTTGGGTATTAGGGATAATCACCAAAAAGTa                                     <  2:2232799/150‑1 (MQ=255)
                                                                                                                    gNTAATATTTATGGTGCTGGCCAGTGCCCCCATCGAGCCCAGTCGCCGTTTCAGTACTTCTATACCGAGGAAGAATCCCAGCACTGCCAGGCTAATGCACAGCAGCAACATGCTCTGTTGGGTATTAGGGATAATCACCAAAAAGTagag                                  >  1:3966702/1‑150 (MQ=255)
                                                                                                                    gATAATATTTATGGTGCTGGCCAGTGCCCCCATCGAGCCCAGTCGCCGTTTCTGTACTTCTATACCGAGGAAGAATCCCAGCACTGCCAGGCTAATGCACAGCAGCAACATGCTCTGTTGGGTATTAGGGATAATCACCAAAAAGTagag                                  >  2:3638107/1‑150 (MQ=255)
                                                                                                                    gATAATATTTATGGTGCTGGCCAGTGCCCCCATAGAGCCCAGTCGCCGTTTCTGTACTTCTATACCGAGGAAGAATCCCAGCACTGCCAGGCTAATGCACAGCAGCAACATGCTCTGTTGGGTATTAGGGATAATCACCAAAAAGTagag                                  >  2:2678109/1‑150 (MQ=255)
                                                                                                                     ataataTTTATGGTGCTGGCCAGTGCCCCCATCGAGCCCAGTCGCCGTTTCTGTACTTCTATACCGAGGAAGAATCCCAGCACTGCCAGGCTAATGCACAGCAGCAACATGCTCTGTTGGGTATTAGGGATAATCACCAAAAAGTAGagc                                 >  2:358000/1‑150 (MQ=255)
                                                                                                                      taataTTTATGGTGCTGGCCAGTGCCCCCATCGAGCCCAGTCGCCGTTTCTGTACTTCTATACCGAGGAAGAATCCCAGCACTGCCAGGCTAATGCACAGCAGCAACATGCTCTGTTGGGTATTAGGGATAATCACCAAAAAGTAGagca                                <  2:4043412/150‑1 (MQ=255)
                                                                                                                             tATGGTGCTGGCCAGTGCCCCCATCGAGCCCAGTCGCCGTTTCTGTACTTCTATACCGAGGAAGAATCCCAGCACTGCCAGGCTAATGCACAGCAGCAACATGCTCTGTTGGGTATTAGGGATAATCACCAAAAAGTAGAGCAGCCCTaa                         <  2:4066813/150‑1 (MQ=255)
                                                                                                                              aTGGTGCTGGCCAGTGCCCCCATCGAGCCCAGTCGCCGTTTCTGTACTTCTATACCGAGGAAGAATCCCAGCACTGCCAGGCTAATGCACAGCAGCAACATGCTCTGTTGGGTATTAGGGATAATCACCAAAAAGTAGAGCAGCCCTAAc                        <  2:4778058/150‑1 (MQ=255)
                                                                                                                              aTGGTGCTGGCCAGTGCCCCCATCGAGCCCAGTCGCCGTTTCTGTACTTCTATACCGAGGAAGAATCCCAGCACTGCCAGGCTAATGCACAGCAGCAACATGCTCTGTTGGGTATTAGGGATAATCACCAAAAAGTAGAGCAGCCCTAAc                        <  1:220091/150‑1 (MQ=255)
                                                                                                                               tGGTGCTGGCCAGTGCCCCCATAGAGCCCAGTCGCCGTTTCTGTACTTCTATACCGAGGAAGAATCCCAGCACTGCCAGGCTAATGCACAGCAGCAACATGCTCTGTTGGGTATTAGGGATAATCACCAAAAAGTAGAGCAGCCCTAACg                       <  1:2216770/150‑1 (MQ=255)
                                                                                                                                ggTGCTGGCCAGTGCCCCCATCGAGCCCAGTCGCCGTTTCTGTACTTCTATACCGAGGAAGAATCCCAGCACTGCCAGGCTAATGCACAGCAGCAACATGCTCTGTTGGGTATTAGGGATAATCACCAAAAAGTAGAGCAGCCCTAACgg                      >  1:1061630/1‑150 (MQ=255)
                                                                                                                                 gTGCTGGCCAGTGCCCCCATCGAGCCCAGTCGCCGTTTCTGTACTTCTATACCGAGGAAGAATCCCAGCACTGCCAGGCTAATGCACAGCAGCAACATGCTCTGTTGGGTATTAGGGATAATCACCAAAAAGTAGAGCAGCCCTAACGGc                     <  1:1637983/150‑1 (MQ=255)
                                                                                                                                  tGCTGGCCAGTGCCCCCATAGAGCCCAGTCGCCGTTTCTGTACTTCTATACCGAGGAAGAATCCCAGCACTGCCAGGCTAATGCACAGCAGCAACATGCTCTGTTGGGTATTAGGGATAATCACCAAAAAGTAGAGCAGCCCTAACGGCa                    >  1:1890735/1‑150 (MQ=255)
                                                                                                                                  tGCTGGCCAGTGCCCCCATAGAGCCCAGTCGCCGTTTCTGTACTTCTATACCGAGGAAGAATCCCAGCACTGCCAGGCTAATGCACAGCAGCAACATGCTCTGTTGGGTATTAGGGATAATCACCAAAAAGTAGAGCAGCCCTAACGGCa                    >  1:2619164/1‑150 (MQ=255)
                                                                                                                                      ggCCAGTGCCCCCATCGAGCCCAGTCGCCGTTTCTGTACTTCTATACCGAGGAAGAATCCCAGCACTGCCAGGCTAATGCACAGCAGCAACATGCTCTGTTGGGTATTAGGGATAATCACCAAAAAGTAGAGCAGCCCTAACGGCAgcgc                >  2:2483063/1‑150 (MQ=255)
                                                                                                                                       gCCAGTGCCCCCATCGAGCCCAGTCGCCGTTTCTGTACTTCTATACCGAGGAAGAATCCCAGCACTGCCAGGCTAATGCACAGCAGCAACATGCTCTGTTGGGTATTAGGGATAATCACCAAAAAGTAGAGCAGCCCTAACGGCAGcgcg               <  2:1300841/150‑1 (MQ=255)
                                                                                                                                       gCCAGTGCCCCCATCGAGCCCAGTCGCCGTTTCTGTACTTCTATACCGAGGAAGAATCCCAGCACTGCCAGGCTAATGCACAGCAGCAACATGCTCTGTTGGGTATTAGGGATAATCACCAAAAAGTAGAGCAGCCCTAACGGCAGcgcg               >  1:4246275/1‑150 (MQ=255)
                                                                                                                                         cAGTGCCCCCATCGAGCCCAGTCGCCGTTTCTGTACTTCTATACCGAGGAAGAATCCCAGCACTGCCAGGCTAATGCACAGCAGCAACATGCTCTGTTGGGTATTAGGGATAATCACCAAAAAGTAGAGCAGCCCTAACGGCAGCGCGGc             >  1:302893/1‑150 (MQ=255)
                                                                                                                                             gCCCCCATAGAGCCCAGTCGCCGTTTCTGTACTTCTATACCGAGGAAGAATCCCAGCACTGCCAGGCTAATGCACAGCAGCAACATGCTCTGTTGGGTATTAGGGATAATCACCAAAAAGTAGAGCAGCCCTAACGGCAGCGCGGCCAGc         >  2:4461686/1‑150 (MQ=255)
                                                                                                                                              cccccATCGAGCCCAGTCGCCGTTTCTGTACTTCTATACCGAGGAAGAATCCCAGCACTGCCAGGCTAATGCACAGCAGCAACATGCTCTGTTGGGTATTAGGGATAATCACCAAAAAGTAGAGCAGCCCTAACGGCAGCGCGGCCAGCg        >  2:5044050/1‑150 (MQ=255)
                                                                                                                                               cccaATAGAGCCCAGTCGCAGTTTATGTACTTCTATACCGAGGAAGAATCCCAGCACTGCCAGGCTAATGCACAGCAGCAACATGCTCTGTTGGGTATTAGGGATAATCACCAAAAAGTAGAGCAGCCCTAACGGCAGCGCGGCCAGCGt       <  2:2998110/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
CGCGATTTATCCCGCACCAGCAAAATAACGGTGAACGCGAGCACACAGCCGACGATTTGCCCTAATGCGCTGTCGAGAAACTGACTGAAATGGAAAGTCATCGGGTTATCCAGCACGATAATATTTATGGTGCTGGCCAGTGCCCCCATCGAGCCCAGTCGCCGTTTCTGTACTTCTATACCGAGGAAGAATCCCAGCACTGCCAGGCTAATGCACAGCAGCAACATGCTCTGTTGGGTATTAGGGATAATCACCAAAAAGTAGAGCAGCCCTAACGGCAGCGCGGCCAGCGTCCCGT  >  CP009273/3380025‑3380322

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