Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 968,415 G→T 94.7% S148* (TCG→TAG)  ycbC ← DUF218 superfamily protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273968,4150GT94.7% 869.7 / 31.3 259S148* (TCG→TAG) ycbCDUF218 superfamily protein
Reads supporting (aligned to +/- strand):  ref base G (7/7);  new base T (113/132);  total (120/139)
Fisher's exact test for biased strand distribution p-value = 7.90e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 259 total aligned reads displayed.

AAAAGATCATCGCCCGCGGCAGATGTGATGCGGAGGTCACCAGCAGGAAAGGGGCGTCACCAATCGCCTGTTTCACTGCTGCAGCTTCTTCTTCGGTATCTTTTGGCAAATCCAGGGTGATAATTTGCTCGCGCGGCACGCCCAGCGATTGCGCAACTCTGGCACCTACTTCCGCTGTACTCACCGTATTGGTTTTTGCTACGCCTCCCGTGAAGATCAGTTTTGATCCCGGATTTTCGCGCCATAAGCGAATACCTTCATTCAGGCGTGGCAGACTGTTATTGATTAAATT  >  CP009273/968269‑968560
                                                                                                                                                  |                                                                                                                                                 
aaaaGATCATCGCCCGCGGCAGATGTGATGCGGAGGTCACCAGCAGGAAAGGGGCGTCACCAATCGCCTGTTTCACTGCTGCAGCTTCTTCTTCGGTATCTTTTGGCAAATCCAGGGTGATAATTTGCTCGCGCGGCACGCCCAGCTAtt                                                                                                                                                >  1:2518360/1‑150 (MQ=255)
 aaaGATCATCGCCCGCGGCAGATGTGATGCGGAGGTCACCAGCAGGAAAGGGGCGTCACCAATCGCCTGTTTCACTGCTGCAGCTTCTTCTTCGGTATCTTTTGGCAAATCCAGGGTGATAATTTGCTCGCGCGGCACGCCCAGCGATTg                                                                                                                                               <  2:198585/150‑1 (MQ=255)
   aGATCATCGCCCGCGGCAGATGTGATGCGGAGGTCACCAGCAGGAAAGGGGCGTCACCAATCGCCTGTTTCACTGCTGCAGCTTCTTCTTCGGTATCTTTTGGCAAATCCAGGGTGATAATTTGCTCGCGCGGCACGCCCAGCGATTgcg                                                                                                                                             >  1:1632288/1‑150 (MQ=255)
       catcGCCCGCGGCAGATGTGATGCGGAGGTCACCAGCAGGAAAGGGGCGTCACCAATCGCCTGTTTCACTGCTGCAGCTTCTTCTTCGGTATCTTTTGGCAAATCCAGGGTGATAATTTGCTCGCGCGGCACGCCCAGCTATTGCGCAAc                                                                                                                                         >  1:1811058/1‑150 (MQ=255)
        atcGCCCGCGGCAGATGTGATGCGGAGGTCACCAGCAGGAAAGGGGCGTCACCAATCGCCTGTTTCACTGCTGCAGCTTCTTCTTCGGTATCTTTTGGCAAATCCAGGGTGATAATTTGCTCGCGCGGCACGCCCAGCTCTTGCGCAAct                                                                                                                                        >  1:377672/1‑150 (MQ=255)
        atcGCCCGCGGCAGATGTGATGCGGAGGTCACCAGCAGGAAAGGGGCGTCACCAATCGCCTGTTTCACTGCTGCAGCTTCTTCTTCGGTATCTTTTGGCAAATCCAGGGTGATAATTTGCTCGCGCGGCACGCCCAGCTATTGCGCAAct                                                                                                                                        >  2:3206861/1‑150 (MQ=255)
        atcGCCCGCGGCAGATGTGATGCGGAGGTCACCAGCAGGAAAGGGGCGTCACCAATCGCCTGTTTCACTGCTGCAGCTTCTTCTTCGGTATCTTTTGGCAAATCCAGGGTGATAATTTGCTCGCGCGGCACGCCCAGCTATTGCGCAAct                                                                                                                                        >  1:2758204/1‑150 (MQ=255)
         tcGCCCGCGGCAGATGTGATGCGGAGGTCACCAGCAGGAAAGGGGCGTCACCAATCGCCTGTTTCACTGCTGCAGCTTCTTCTTCGGTATCTTTTGGCAAATCCAGGGTGATAATTTGCTCGCGCGGCACGCCCAGCTATTGCGCAActc                                                                                                                                       <  1:960024/150‑1 (MQ=255)
         tcGCCCGCGGCAGATGTGATGCGGAGGTCACCAGCAGGAAAGGGGCGTCACCAATCGCCTGTTTCACTGCTGCAGCTTCTTCTTCGGTATCTTTTGGCAAATCCAGGGTGATAATTTGCTCGCGCGGCACGCCCAGCTATTGCGCAActc                                                                                                                                       <  1:4216209/150‑1 (MQ=255)
          cGCCCGCGGCAGATGTGATGCGGAGGTCACCAGCAGGAAAGGGGCGTCACCAATCGCCTGTTTCACTGCTGCAGCTTCTTCTTCGGTATCTTTTGGCAAATCCAGGGTGATAATTTGCTCGCGCGGCACGCCCAGCTATTGCGCAActct                                                                                                                                      >  2:3273880/1‑150 (MQ=255)
           gCCCGCGGCAGATGTGATGCGGAGGTCACCAGCAGGAAAGGGGCGTCACCAATCGCCTGTTTCACTGCTGCAGCTTCTTCTTCGGTATCTTTTGGCAAATCCAGGGTGATAATTTGCTCGCGCGGCACGCCCAGCTATTGCGCAACTCTg                                                                                                                                     >  2:5115925/1‑150 (MQ=255)
              cgcgGCAGATGTGATGCGGAGGTCACCAGCAGGAAAGGGGCGTCACCAATCGCCTGTTTCACTGCTGCAGCTTCTTCTTCGGTATCTTTTGGCAAATCCAGGGTGATAATTTGCTCGCGCGGCACGCCCAGCTATTGCGCAACTCTGGCa                                                                                                                                  <  2:1703440/150‑1 (MQ=255)
                 ggCAGATGTGATGCGGAGGTCACCAGCAGGAAAGGGGCGTCACCAATCGCCTGTTTCACTGCTGCAGCTTCTTCTTCGGTATCTTTTGGCAAATCCAGGGTGATAATTTGCTCGCGCGGCACGCCCAGCTATTGCGCAACTCTGGCACCt                                                                                                                               <  1:2439518/150‑1 (MQ=255)
                  gCAGATGTGATGCGGAGGTCACCAGCAGGAAAGGGGCGTCACCAATCGCCTGTTTCACTGCTGCAGCTTCTTCTTCGGTATCTTTTGGCAAATCCAGGGTGATAATTTGCTCGCGCGGCACGCCCAGCTATTGCGCAACTCTGGCACCTa                                                                                                                              <  2:4303637/150‑1 (MQ=255)
                   cAGATGTGATGCGGAGGTCACCAGCAGGAAAGGGGCGTCACCAATCGCCTGTTTCACTGCTGCAGCTTCTTCTTCGGTATCTTTTGGCAAATCCAGGGTGATAATTTGCTCGCGCGGCACGCCCAGCTATTGCGCAACTCTGGCACCTAc                                                                                                                             >  1:424947/1‑150 (MQ=255)
                   cAGATGTGATGCGGAGGTCACCAGCAGGAAAGGGGCGTCACCAATCGCCTGTTTCACTGCTGCAGCTTCTTCTTCGGTATCTTTTGGCAAATCCAGGGTGATAATTTGCTCGCGCGGCACGCCCAGCGATTGCGCAACTCTGGCACCTAc                                                                                                                             >  1:2057529/1‑150 (MQ=255)
                      aTGTGATGCGGAGGTCACCAGCAGGAAAGGGGCGTCACCAATCGCCTGTTTCACTGCTGCAGCTTCTTCTTCGGTATCTTTTGGCAAATCCAGGGTGATAATTTGCTCGCGCGGCACGCCCAGCTATTGCGCAACTCTGGCACCTACTTc                                                                                                                          >  2:3640443/1‑150 (MQ=255)
                         tgATGCGGAGGTCACCAGCAGGAAAGGGGCGTCACCAATCGCCTGTTTCACTGCTGCAGCTTCTTCTTCGGTATCTTTTGGCAAATCCAGGGTGATAATTTGCTCGCGCGGCACGCCCAGCTATTGCGCAACTCTGGCACCTACTTCCGc                                                                                                                       <  1:2098768/150‑1 (MQ=255)
                         tgATGCGGAGGTCACCAGCAGGAAAGGGGCGTCACCAATCGCCTGTTTCACTGCTGCAGCTTCTTCTTCGGTATCTTTTGGCAAATCCAGGGTGATAATTTGCTCGCGCGGCACGCCCAGCTATTGCGCAACTCTGGCACCTACTTCCGc                                                                                                                       >  2:797347/1‑150 (MQ=255)
                           aTGCGGAGGTCACCAGCAGGAAAGGGGCGTCACCAATCGCCTGTTTCACTGCTGCAGCTTCTTCTTCGGTATCTTTTGGCAAATCCAGGGTGATAATTTGCTCGCGCGGCACGCCCAGCTATTGCGCAACTCTGGCACCTACTTCCGCTg                                                                                                                     >  2:4027441/1‑150 (MQ=255)
                            tGCGGAGGTCACCAGCAGGAAAGGGGCGTCACCAATCGCCTGTTTCACTGCTGCAGCTTCTTCTTCGGTATCTTTTGGCAAATCCAGGGTGATAATTTGCTCGCGCGGCACGCCCAGCTATTGCGCAACTCTGGCACCTACTTCCGCTGt                                                                                                                    <  2:834727/150‑1 (MQ=255)
                            tGCGGAGGTCACCAGCAGGAAAGGGGCGTCACCAATCGCCTGTTTCACTGCTGCAGCTTCTTCTTCGGTATCTTTTGGCAAATCCAGGGTGATAATTTGCTCGCGCGGCACGCCCAGCTATTGCGCAACTCTGGCACCTACTTCCGCTGt                                                                                                                    <  2:2443318/150‑1 (MQ=255)
                             gCGGAGGTCACCAGCAGGAAAGGGGCGTCACCACTCGCCTGTTTCACTGCTGCAGCTTCTTCTTCGGTATCTTTTGGCAAATCCAGGGTGATAATTTGCTCGCGCGGCACGCCCAGCTATTGCGCAACTCTGGCACCTACTTCCGCTGTa                                                                                                                   <  2:1862825/150‑1 (MQ=255)
                             gCGGAGGTCACCAGCAGGAAAGGGGCGTCACCAATCGCCTGTTTCACTGCTGCAGCTTCTTCTTCGGTATCTTTTGGCAAATCCAGGGTGATAATTTGCTCGCGCGGCACGCCCAGCTATTGCGCAACTCTGGCACCTACTTCCGCTGTa                                                                                                                   <  2:1862953/150‑1 (MQ=255)
                             gCGGAGGTCACCAGCAGGAAAGGGGCGTCACCAATCGCCTGTTTCACTGCTGCAGCTTCTTCTTCGGTATCTTTTGGCAAATCCAGGGTGATAATTTGCTCGCGCGGCACGCCCAGCTATTGCGCAACTCTGGCACCTACTTCCGCTGTa                                                                                                                   >  1:3471989/1‑150 (MQ=255)
                               ggAGGTCACCAGCAGGAAAGGGGCGTCACCAATCGCCTGTTTCACTGCTGCAGCTTCTTCTTCGGTATCTTTTGGCAAATCCAGGGTGATAATTTGCTCGCGCGGCACGCCCAGCTATTGCGCAACTCTGGCACCTACTTCCGCTGTACt                                                                                                                 >  1:2614626/1‑150 (MQ=255)
                                  ggTCACCAGCAGGAAAGGGGCGTCACCAATCGCCTGTTTCACTGCTGCAGCTTCTTCTTCGGTATCTTTTGGCAAATCCAGGGTGATAATTTGCTCGCGCGGCACGCCCTGCTATTGCGCAACTCTGGCACCTACTTCCGCTGTACTCAc                                                                                                              <  2:4696436/150‑1 (MQ=255)
                                        cagcagGAAAGGGGCGTCACCAATCGCCTGTTTCACTGCTGCAGCTTCTTCTTCGGTATCTTTTGGCAAATCCAGGGTGATAATTTGCTCGCGCGGCACGCCCAGCTATTGCGCAACTCTGGCACCTACTTCCGCTGTACTCACCGTAtt                                                                                                        <  2:934570/150‑1 (MQ=255)
                                         agcagGAAAGGGGCGTCACCAATCGCCTGTTTCACTGCTGCAGCTTCTTCTTCGGTATCTTTTGGCAAATCCAGGGTGATAATTTGCTCGCGCGGCACGCCCAGCTATTGCGCAACTCTGGCACCTACTTCCGCTGTACTCACCGTATTg                                                                                                       >  2:3466133/1‑150 (MQ=255)
                                           cagGAAAGGGGCGTCACCAATCGCCTGTTTCACTGCTGCAGCTTCTTCTTCGGTATCTTTTGGCAAATCCAGGGTGATAATTTGCTCGCGCGGCACGCCCAGCTATTGCGCAACTCTGGCACCTACTTCCGCTGTACTCACCGTATTGGt                                                                                                     >  1:411557/1‑150 (MQ=255)
                                           cagGAAAGGGGCGTCACCAATCGCCTGTTTCACTGCTGCAGCTTCTTCTTCGGTATCTTTTGGCAAATCCAGGGTGATAATTTGCTCGCGCGGCACGCCCAGCTATTGCGCAACTCTGGCACCTACTTCCGCTGTACTCACCGTAGTGGt                                                                                                     >  2:1449879/1‑150 (MQ=255)
                                            agGAAAGGGGCGTCACCAATCGCCTGTTTCACTGCTGCAGCTTCTTCTTCGGTATCTTTTGGCAAATCCAGGGTGATAATTTGCTCGCGCGGCACGCCCAGCTATTGCGCAACTCTGGCACCTACTTCCGCTGTACTCACCGTATTGGtt                                                                                                    >  1:2152673/1‑150 (MQ=255)
                                            agGAAAGGGGCGTCACCAATCGCCTGTTTCACTGCTGCAGCTTCTTCTTCGGTATCTTTTGGCAAATCCAGGGTGATAATTTGCTCGCGCGGCACGCCCAGCTATTGCGCAACTCTGGCACCTACTTCCGCTGTACTCACCGTATTGGtt                                                                                                    >  1:4045481/1‑150 (MQ=255)
                                            agGAAAGGGGCGTCACCAATCGCCTGTTTCACTGCTGCAGCTTCTTCTTCGGTATCTTTTGGCAAATCCAGGGTGATAATTTGCTCGCGCGGCACGCCCAGCTATTGCGCAACTCTGGCACCTACTTCCGCTGTACTCACCGTATTGGtt                                                                                                    >  1:1925478/1‑150 (MQ=255)
                                              gAAAGGGGCGTCACCAATCGCCTGTTTCACTGCTGCAGCTTCTTCTTCGGTATCTTTTTGCAAATCCAGGGTGATAATTTGCTCGCGCGGCACGCCCAGCTATTGCGCAACTCTGGCACCTACTTCCGCTGTACTCACCGTATTGGtttt                                                                                                  >  2:1398716/1‑150 (MQ=255)
                                              gAAAGGGGCGTCACCAATCGCCTGTTTCACTGCTGCAGCTTCTTCTTCGGTATCTTTTGGCAAATCCAGGGTGATAATTTGCTCGCGCGGCACGCCCAGCTCTTGCGCAACTCTGGCACCTACTTCCGCTGTACTCACCGTATTGGtttt                                                                                                  >  2:4778798/1‑150 (MQ=255)
                                               aaaGGGGCGTCACCAATCGCCTGTTTCACTGCTGCAGCTTCTTCTTCGGTATCTTTTGGCAAATCCAGGGTGATAATTTGCTCGCGCGGCACGCCCAGCTATTGCGCAACTCTGGCACCTACTTCCGCTGTACTCACCGTATTGGttttt                                                                                                 <  2:1011469/150‑1 (MQ=255)
                                                aaGGGGCGTCACCAATCGCCTGTTTCACTGCTGCAGCTTCTTCTTCGGTATCTTTTGGCAAATCCAGGGTGATAATTTGCTCGCGCGGCACGCCCAGCTATTGCGCAACTCTGGCACCTACTTCCGCTGTACTCACCGTATTGGTTTTTg                                                                                                >  2:599653/1‑150 (MQ=255)
                                                    ggCGTCACCAATCGCCTGTTTCACTGCTGCAGCTTCTTCTTCGGTATCTTTTGGCAAATCCAGGGTGATAATTTGCTCGCGCGGCACGCCCAGCTATTGCGCAACTCTGGCACCTACTTCCGCTGTACTCACCGTATTGGTTTTTGCTAc                                                                                            <  1:5159608/150‑1 (MQ=255)
                                                    ggCGTCACCAATCGCCTGTTTCACTGCTGCAGCTTCTTCTTCGGTATCTTTTGGCAAATCCAGGGTGATAATTTGCTCGCGCGGCACGCCCAGCGATTGCGCAACTCTGGCACCTACTTCCGCTGTACTCACCGTATTGGTTTTTGCTAc                                                                                            >  1:494684/1‑150 (MQ=255)
                                                     gCGTCACCAATCGCCTGTTTCACTGCTGCAGCTTCTTCTTCGGTATCTTTTGGCAAATCCAGGGTGATAATTTGCTCGCGCGGCACGCCCAGCTATTGCGCAACTCTGGCACCTACTTCCGCTGTACTCACCGTATTGGTTTTTGCTACg                                                                                           <  1:2769338/150‑1 (MQ=255)
                                                       gTCACCAATCGCCTGTTTCACTGCTGCAGCTTCTTCTTCGGTATCTTTTGGCAAATCCAGGGTGATAATTTGCTCGCGCGGCACGCCCAGCTATTGCGCAACTCTGGCACCTACTTCCGCTGTACTCACCGTATTGGTTTTTGCTACGcc                                                                                         >  1:1036102/1‑150 (MQ=255)
                                                       gTCACCAATCGCCTGTTTCACTGCTGCAGCTTCTTCTTCGGTATCTTTTGGCAAATCCAGGGTGATAATTTGCTCGCGCGGCACGCCCAGCTATTGCGCAACTCTGGCACCTACTTCCGCTGTACTCACCGTATTGGTTTTTGCTACGcc                                                                                         >  1:2871223/1‑150 (MQ=255)
                                                              aTCGCCTGTTTCACTGCTGCAGCTTCTTCTTCGGTATCTTTTGGCAAATCCAGGGTGATAATTTGCTCGCGCGGCACGCCCAGCTATTGCGCAACTCTGGCACCTACTTCCGCTGTACTCACCGTATTGGTTTTTGCTACGCCTCCCGTg                                                                                  <  2:2193368/150‑1 (MQ=255)
                                                                cGCCTGTTTCACTGCTGCAGCTTCTTCTTCGGTATCTTTTGGCAAATCCAGGGTGATAATTTGCTCGCGCGGCACGCCCAGCTATTGCGCAACTCTGGCACCTGCTTCCGCTGTACTCACCGTATTGGTTTTTGCTACGCCTCCCGTGaa                                                                                <  1:2628707/150‑1 (MQ=255)
                                                                cGCCTGTTTCACTGCTGCAGCTTCTTCTTCGGTATCTTTTGGCAAATCCAGGGTGATAATTTGCTCGCGCGGCACGCCCAGCGATTGCGCAACTCTGGCACCTACTTCCGCTGTACTCACCGTATTGGTTTTTGCTACGCCTCCCGTGaa                                                                                <  1:1681940/150‑1 (MQ=255)
                                                                    tGTTTCACTGCTGCAGCTTCTTCTTCGGTATCTTTTGGCAAATCCAGGGTGATAATTTGCTCGCGCGGCACGCCCAGCTATTGCGCAACTCTGGCACCTACTTCCGCTGTACTCACCGTATTGGTTTTTGCTACGCCTCCCGTGAAGATc                                                                            <  2:3068473/150‑1 (MQ=255)
                                                                    tGTTTCACTGCTGCAGCTTCTTCTTCGGTATCTTTTGGCAAATCCAGGGTGATAATTTGCTCGCGCGGCACGCCCAGCTATTGCGCAACTCTGGCACCTACTTCCGCTGTACTCACCGTATTGGTTTTTGCTACGCCTCCCGTGAAGATc                                                                            <  2:5037105/150‑1 (MQ=255)
                                                                        tCACTGCTGCAGCTTCTTCTTCGGTATCTTTTGGCAAATCCAGGGTGATAATTTGCTCGCGCGGCACGCCCAGCTATTGCGCAACTCTGGCACCTACTTCCGCTGTACTCACCGTATTGGTTTTTGCTACGCCTCCCGTGAAGATCAGtt                                                                        <  1:436014/150‑1 (MQ=255)
                                                                          aCTGCTGCAGCTTCTTCTTCGGTATCTTTTGGCAAATCCAGGGTGATAATTTGCTCGCGCGGCACGCCCAGCTATTGCGCAACTCTGGCACCTACTTCCGCTGTACTCACCGTATTGGTTTTTGCTACGCCTCCCGTGAAGATCAGtttt                                                                      <  2:2511973/150‑1 (MQ=255)
                                                                               tgcAGCTTCTTCTTCGGTATCTTTTGGCAAATCCAGGGTGATAATTTGCTCGCGCGGCACGCCCAGCTATTGCGCAACTCTGGCACCTACTTCCGCTGTACTCACCGTATTGGTTTTTGCTACGCCTCCCGTGAAGATCAGTTTTGATcc                                                                 <  1:3846944/150‑1 (MQ=255)
                                                                                   gCTTCTTCTTCGGTATCTTTTGGCAAATCCAGGGTGATAATTTGCTCGCGCGGCACGCCCAGCTATTGCGCAACTCTGGCACCTACTTCCGCTGTACTCACCGTATTGGTTTTTGCTACGCCTCCCGTGAAGATCAGTTTTGATCCCGGa                                                             >  2:4631832/1‑150 (MQ=255)
                                                                                    cttcttcttCGGTATCTTTTGGCAAATCCAGGGTGATAATTTGCTCGCGCGGCACGCCCAGCTATTGCGCAACTCTGGCACCTACTTCCGCTGTACTCACCGTATTGGTTTTTGCTACGCCTCCCGTGAAGATCAGTTTTGATCCCGGAt                                                            <  1:1305905/150‑1 (MQ=255)
                                                                                      tcttcttcGGTATCTTTTGGCAAATCCAGGGTGATAATTTGCTCGCGCGGCACGCCCAGCTATTGCGCAACTCTGGCACCTACTTCCGCTGTACTCACCGTATTGGTTTTTGCTACGCCTCCCGTGAAGATCAGTTTTGATCCCGGAttt                                                          >  2:3148365/1‑150 (MQ=255)
                                                                                      tcttcttcGGTATCTTTTGGCAAATCCAGGGTGATAATTTGCTCGCGCGGCACGCCCAGCTATTGCGCAACTCGGGCACCTACTTCCGCTGTACTCACCGTATTGGTTTTTGCTACGCCTCCCGTGAAGATCAGTTTTGATCCCGGAttt                                                          <  2:2448010/150‑1 (MQ=255)
                                                                                       cttcttcGGTATCTTTTGGCAAATCCAGGGTGATAATTTGCTCGCGCGGCACGCCCAGCTATTGCGCAACTCTGGCACCTACTTCCGCTGTACTCACCGTATTGGTTTTTGCTACGCCTCCCGTGAAGATCAGTTTTGATCCCGGAtttt                                                         >  1:3648882/1‑150 (MQ=255)
                                                                                       cttcttcGGTATCTTTTGGCAAATCCAGGGTGATAATTTGCTCGCGCGGCACGCCCAGCTATTGCGCAACTCTGGCACCTACTTCCGCTGTACTCACCGTATTGGTTTTTGCTACGCCTCCCGTGAAGATCAGTTTTGATCCCGGAtttt                                                         >  2:3912623/1‑150 (MQ=255)
                                                                                       cttcttcGGTATCTTTTGGCAAATCCAGGGTGATAATTTGCTCGCGCGGCACGCCCAGCTATTGCGCAACTCTGGCACCTACTTCCGCTGTACTCACCGTATTGGTTTTTGCTACGCCTCCCGTGAAGATCAGTTTTGATCCCGGAtttt                                                         >  1:155629/1‑150 (MQ=255)
                                                                                        ttcttcGGTATCTTTTGGCAAATCCAGGGTGATAATTTGCTCGCGCGGCACGCCCAGCTATTGCGCAACTCTGGCACCTACTTCCGCTGTACTCACCGTATTGGTTTTTGCTACGCCTCCCGTGAAGATCAGTTTTGATCCCGGATTTTc                                                        >  1:1973919/1‑150 (MQ=255)
                                                                                         tcttcGGTATCTTTTGGCAAATCCAGGGTGATAATTTGCTCGCGCGGCACGCCCAGCTCTTGCGCAACTCTGGCACCTACTTCCGCTGTACTCACCGTATTGGTTTTTGCTACGCCTCCCGTGAAGATCAGTTTTGATCCCGGATTTTcg                                                       >  1:4450526/1‑150 (MQ=255)
                                                                                           ttcGGTATCTTTTGGCAAATCCAGGGTGATAATTTGCTCGCGCGGCACGCCCAGCTATTGCGCAACTCTGGCACCTACTTCCGCTGTACTCACCGTATTGGTTTTTGCTACGCCTCCCGTGAAGATCAGTTTTGATCCCGGATTTTcgcg                                                     <  1:3383555/150‑1 (MQ=255)
                                                                                           ttcGGTATCTTTTGGCAAATCCAGGGTGATAATTTGCTCGCGCGGCACGCCCAGCTATTGCGCAACTCTGGCACCTACTTCCGCTGTACTCACCGTATTGGTTTTTGCTACGCCTCCCGTGAAGATCAGTTTTGATCCCGGATTTTcgcg                                                     >  2:37417/1‑150 (MQ=255)
                                                                                              ggTATCTTTTGGCAAATCCAGGGTGATAATTTGCTCGCGCGGCACGCCCAGCTATTGCGCAACTCTGGCACCTACTTCCGCTGTACTCACCGTATTGGTTTTTGCTACGCCTCCCGTGAAGATCAGTTTTGATCCCGGATTTTCGCGCCa                                                  <  1:4724808/150‑1 (MQ=255)
                                                                                                tATCTTTTGGCAAATCCAGGGTGATAATTTGCTCGCGCGGCACGCCCAGCTATTGCGCAACTCTGGCACCTACTTCCGCTGTACTCACCGTATTGGTTTTTGCTACGCCTCCCGTGAAGATCAGTTTTGATCCCGGATTTTCGCGCCATa                                                >  2:298220/1‑150 (MQ=255)
                                                                                                tATCTTTTGGCAAATCCAGGGTGATAATTTGCTCGCGCGGCACGCCCAGCTATTGCGCAACTCTGGCACCTACTTCCGCTGTACTCACCGTATTGGTTTTTGCTACGCCTCCCGTGAAGATCAGTTTTGATCCCGGATTTTCGCGCCATa                                                <  2:4922848/150‑1 (MQ=255)
                                                                                                tATCTTTTGGCAAATCCAGGGTGATAATTTGCTCGCGCGGCACGCCCAGCGATTGCGCAACTCTGGCACCTACTTCCGCTGTACTCACCGTATTGGTTTTTGCTACGCCTCCCGTGAAGATCAGTTTTGATCCCGGATTTTCGCGCCATa                                                >  2:360919/1‑150 (MQ=255)
                                                                                                 aTCTTTTGGCAAATCCAGGGTGATAATTTGCTCGCGCGGCACGCCCAGCTATTGCGCAACTCTGGCACCTACTTCCGCTGTACTCACCGTATTGGTTTTTGCTACGCCTCCCGTGAAGATCAGTTTTGATCCCGGATTTTCGCGCCATaa                                               >  1:2252855/1‑150 (MQ=255)
                                                                                                 aTCTTTTGGCAAATCCAGGGTGATAATTTGCTCGCGCGGCACGCCCAGCTATTGCGCAACTCTGGCACCTACTTCCGCTGTACTCACCGTATTGGTTTTTGCTACGCCTCCCGTGAAGATCAGTTTTGATCCCGGATTTTCGCGCCATaa                                               <  2:448320/150‑1 (MQ=255)
                                                                                                  tCTTTTGGCAAATCCAGGGTGATAATTTGCTCGCGCGGCACGCCCAGCTATTGCGCAACTCTGGCACCTACTTCCGCTGTACTCACCGTATTGGTTTTTGCTACGCCTCCCGTGAAGATCAGTTTTGATCCCGGATTTTCGCGCCATAAg                                              <  2:4599499/150‑1 (MQ=255)
                                                                                                     tttGGCAAATCCAGGGTGATAATTTGCTCGCGCGGCACGCCCAGCTATTGCGCAACTCTGGCACCTACTTCCGCTGTACTCACCGTATTGGTTTTTGCTACGCCTCCCGTGAAGATCAGTTTTGATCCCGGATTTTCGCGCCATAAGCGa                                           >  1:136400/1‑150 (MQ=255)
                                                                                                       tGGCAAATCCAGGGTGATAATTTGCTCGCGCGGCACGCCCAGCTATTGCGCAACTCTGGCACCTACTTCCGCTGTACTCACCGTATTGGTTTTTGCTACGCCTCCCGTGAAGATCAGTTTTGATCCCGGATTTTCGCGCCATAAGCGAAt                                         >  1:119751/1‑150 (MQ=255)
                                                                                                       tGGCAAATCCAGGGTGATAATTTGCTCGCGCGGCACGCCCAGCTATTGCGCAACTCTGGCACCTACTTCCGCTGTACTCACCGTATTGGTTTTTGCTACGCCTCCCGTGAAGATCAGTTTTGATCCCGGATTTTCGCGCCATAAGCGAAt                                         >  1:1348116/1‑150 (MQ=255)
                                                                                                       tGGCAAATCCAGGGTGATAATTTGCTCGCGCGGCACGCCCAGCTATTGCGCAACTCTGGCACCTACTTCCGCTGTACTCACCGTATTGGTTTTTGCTACGCCTCCCGTGAAGATCAGTTTTGATCCCGGATTTTCGCGCCATAAGCGAAt                                         <  2:2879681/150‑1 (MQ=255)
                                                                                                       tGGCAAATCCAGGGTGATAATTTGCTCGCGCGGCACGCCCAGCTATTGCGCAACTCTGGCACCTACTTCCGCTGTACTCACCGTATTGGTTTTTGCTACGCCTCCCGTGAAGATCAGTTTTGATCCCGGATTTTCGCGCCATAAGCGAAt                                         <  1:1466955/150‑1 (MQ=255)
                                                                                                        ggCAAATCCAGGGTGATAATTTGCTCGCGCGGCACGCCCAGCTATTGCGCAACTCTGGCACCTACTTCCGCTGTACTCACCGTATTGGTTTTTGCTACGCCTCCCGTGAAGATCAGTTTTGATCCCGGATTTTCGCGCCATAAGCGAATa                                        <  2:4989308/150‑1 (MQ=255)
                                                                                                        ggCAAATCCAGGGTGATAATTTGCTCGCGCGGCACGCCCAGCTATTGCGCAACTCTGGCACCTACTTCCGCTGTACTCACCGTATTGGTTTTTGCTACGCCTCCCGTGAAGATCAGTTTTGATCCCGGATTTTCGCGCCATAAGCGAATa                                        >  2:1295129/1‑150 (MQ=255)
                                                                                                          cAAATCCAGGGTGATAATTTGCTCGCGCGGCACGCCCAGCTATTGCGCAACTCTGGCACCTACTTCCGCTGTACTCACCGTATTGGTTTTTGCTACGCCTCCCGTGAAGATCAGTTTTGATCCCGGATTTTCGCGCCATAAGCGAATAcc                                      <  1:3989172/150‑1 (MQ=255)
                                                                                                          cAAATCCAGGGTGATAATTTGCTCGCGCGGCACGCCCAGCTATTGCGCAACTCTGGCACCTACTTCCGCTGTACTCACCGTATTGGTTTTTGCTACGCCTCCCGTGAAGATCAGTTTTGATCCCGGATTTTCGCGCCATAAGCGAATAcc                                      <  2:4045481/150‑1 (MQ=255)
                                                                                                          cAAATCCAGGGTGATAATTTGCTCGCGCGGCACGCCCAGCTATTGCGCAACTCTGGCACCTACTTCCGCTGTACTCACCGTATTGGTTTTTGCTACGCCTCCCGTGAAGATCAGTTTTGATCCCGGATTTTCGCGCCATAAGCGAATAcc                                      <  2:538080/150‑1 (MQ=255)
                                                                                                          cAAATCCAGGGTGATAATTTGCTCGCGCGGCACGCCCAGCTATTGCGCAACTCTGGCACCTACTTCCGCTGTACTCACCGTATTGGTTTTTGCTACGCCTCCCGTGAAGATCAGTTTTGATCCCGGATTTTCGCGCCATAAGCGAATAcc                                      <  1:599653/150‑1 (MQ=255)
                                                                                                          cAAATCCAGGGTGATAATTTGCTCGCGCGGCACGCCCAGCTATTGCGCAACTCTGGCACCTACTTCCGCTGTACTCACCGTATTGGTTTTTGCTACGCCTCCCGTGAAGATCAGTTTTGATCCCGGATTTTCGCGCCATAAGCGAATAcc                                      <  1:1502882/150‑1 (MQ=255)
                                                                                                           aaaTCCAGGGTGATAATTTGCTCGCGCGGCACGCCCAGCTATTGCGCAACTCTGGCACCTACTTCCGCTGTACTCACCGTATTGGTTTTTGCTACGCCTCCCGTGAAGATCAGTTTTGATCCCGGATTTTCGCGCCATAAGCGAATACCt                                     >  1:2062875/1‑150 (MQ=255)
                                                                                                           aaaTCCAGGGTGATAATTTGCTCGCGCGGCACGCCCAGCTATTGCGCAACTCTGGCACCTACTTCCGCTGTACTCACCGTATTGGTTTTTGCTACGCCTCCCGTGAAGATCAGTTTTGATCCCGGATTTTCGCGCCATAAGCGAATACCt                                     >  2:2107490/1‑150 (MQ=255)
                                                                                                               ccAGGGTGATAATTTGCTCGCGCGGCACGCCCAGCTATTGCGCAACTCTGGCACCTACTTCCGCTGTACTCACCGTATTGGTTTTTGCTACGCCTCCCGTGAAGATCAGTTTTGATCCCGGATTTTCGCGCCATAAGCGAATACCttcat                                 <  1:3176803/150‑1 (MQ=255)
                                                                                                                 aGGGTGATAATTTGCTCGCGCGGCACGCCCAGCTATTGCGCAACTCTGGCACCTACTTCCGCTGTACTCACCGTATTGGTTTTTGCTACGCCTCCCGTGAAGATCAGTTTTGATCCCGGATTTTCGCGCCATAAGCGAATACCttcattc                               <  1:4922195/150‑1 (MQ=255)
                                                                                                                         aaTTTGCTCGCGCGGCACGCCCAGCTATTGCGCAACTCTGGCACCTACTTCCGCTGTACTCACCGTATTGGTTTTTGCTACGCCTCCCGTGAAGATCAGTTTTGATCCCGGATTTTCGCGCCATAAGCGAATACCTTCATTCAGGCGTgg                       <  2:3020673/150‑1 (MQ=255)
                                                                                                                          aTTTGCTCGCGCGGCACGCCCAGCTATTGCGCAACTCTGGCACCTACTTCCGCTGTACTCACCGTATTGGTTTTTGCTACGCCTCCCGTGAAGATCAGTTTTGATCCCGGATTTTCGCGCCATAAGCGAATACCTTCATTCAGGCGTGGc                      <  2:3144068/150‑1 (MQ=255)
                                                                                                                           tttGCTCGCGCGGCACGCCCAGCTATTGCGCAACTCTGGCACCTACTTCCGCTGTACTCACCGTATTGGTTTTTGCTACGCCTCCCGTGAAGATCAGTTTTGATCCCGGATTTTCGCGCCATAAGCGAATACCTTCATTCAGGCGTGGCa                     <  1:1474996/150‑1 (MQ=255)
                                                                                                                                 cgcgcgGCACGCCCAGCTATTGCGCAACTCTGGCACCTACTTCCGCTGTACTCACCGTATTGGTTTTTGCTACGCCTCCCGTGAAGATCAGTTTTGATCCCGGATTTTCGCGCCATAAGCGAATACCTTCATTCAGGCGTGGCAGACTGt               >  1:5187835/1‑150 (MQ=255)
                                                                                                                                    gcgGCACGCCCAGCTATTGCGCAACTCTGGCACCTACTTCCGCTGTACTCACCGTATTGGTTTTTGCTACGCCTCCCGTGAAGATCAGTTTTGATCCCGGATTTTCGCGCCATAAGCGAATACCTTCATTCAGGCGTGGCAGACTGTTAt            <  1:4604720/150‑1 (MQ=255)
                                                                                                                                    gcgGCACGCCCAGCGATTGCGCAACTCTGGCACCTACTTCCGCTGTACTCACCGTATTGGTTTTTGCTACGCCTCCCGTGAAGATCAGTTTTGATCCCGGATTTTCGCGCCATAAGCGAATACCTTCATTCAGGCGTGGCAGACTGTTAt            <  1:58830/150‑1 (MQ=255)
                                                                                                                                     cgGCACGCCCAGCTATTGCGCAACTCTGGCACCTACTTCCGCTGTACTCACCGTATTGGTTTTTGCTACGCCTCCCGTGAAGATCAGTTTTGATCCCGGATTTTCGCGCCATAAGCGAATACCTTCATTCAGGCGTGGCAGACTGTTAtt           <  1:5061699/150‑1 (MQ=255)
                                                                                                                                     cgGCACGCCCAGCTATTGCGCAACTCTGGCACCTACTTCCGCTGTACTCACCGTATTGGTTTTTGCTACGCCTCCCGTGAAGATCAGTTTTGATCCCGGATTTTCGCGCCATAAGCGAATACCTTCATTCAGGCGTGGCAGACTGTTAtt           >  2:4021299/1‑150 (MQ=255)
                                                                                                                                     cgGCACGCCCAGCTATTGCGCAACTCTGGCACCTACTTCCGCTGTACTCACCGTATTGGTTTTTGCTACGCCTCCCGTGAAGATCAGTTTTGATCCCGGATTTTCGCGCCATAAGCGAATACCTTCATTCAGGCGTGGCAGACTGTTAtt           <  1:716378/150‑1 (MQ=255)
                                                                                                                                        cACGCCCAGCTATTGCGCAACTCTGGCACCTACTTCCGCTGTACTCACCGTATTGGTTTTTGCTACGCCTCCCGTGAAGATCAGTTTTGATCCCGGATTTTCGCGCCATAAGCGAATACCTTCATTCAGGCGTGGCAGACTGTTATTGAt        >  1:561960/1‑150 (MQ=255)
                                                                                                                                          cGCCCAGCTATTGCGCAACTCTGGCACCTACTTCCGCTGTACTCACCGTATTGGTTTTTGCTACGCCTCCCGTGAAGATCAGTTTTGATCCCGGATTTTCGCGCCATAAGCGAATACCTTCATTCAGGCGTGGCAGACTGTTATTGATTa      <  1:3273880/150‑1 (MQ=255)
                                                                                                                                           gCCCAGCTATTGCGCAACTCTGGCACCTACTTCCGCTGTACTCACCGTATTGGTTTTTGCTACGCCTCCCGTGAAGATCAGTTTTGATCCCGGATTTTCGCGCCATAAGCGAATACCTTCATTCAGGCGTGGCAGACTGTTATTGATTaa     >  2:1634160/1‑150 (MQ=255)
                                                                                                                                             ccAGCTATTGCGCAACTCTGGCACCTACTTCCGCTGTACTCACCGTATTGGTTTTTGCTACGCCTCCCGTGAAGATCAGTTTTGATCCCGGATTTTCGCGCCATAAGCGAATACCTTCATTCAGGCGTGGCAGACTGTTATTGATTAAAt   <  1:2961794/150‑1 (MQ=255)
                                                                                                                                              cAGCTATTGCGCAACTCTGGCACCTACTTCCGCTGTACTCACCGTATTGGTTTTTGCTACGCCTCCCGTGAAGATCAGTTTTGATCCCGGATTTTCGCGCCATAAGCGAATACCTTCATTCAGGCGTGGCAGACTGTTATTGATTAAAtt  >  2:1046868/1‑150 (MQ=255)
                                                                                                                                              cAGCTATTGCGCAACTCTGGCACCTACTTCCGCTGTACTCACCGTATTGGTTTTTGCTACGCCTCCCGTGAAGATCAGTTTTGATCCCGGATTTTCGCGCCATAAGCGAATACCTTCATTCAGGCGTGGCAGACTGTTATTGATTAAAtt  >  2:2903606/1‑150 (MQ=255)
                                                                                                                                                  |                                                                                                                                                 
AAAAGATCATCGCCCGCGGCAGATGTGATGCGGAGGTCACCAGCAGGAAAGGGGCGTCACCAATCGCCTGTTTCACTGCTGCAGCTTCTTCTTCGGTATCTTTTGGCAAATCCAGGGTGATAATTTGCTCGCGCGGCACGCCCAGCGATTGCGCAACTCTGGCACCTACTTCCGCTGTACTCACCGTATTGGTTTTTGCTACGCCTCCCGTGAAGATCAGTTTTGATCCCGGATTTTCGCGCCATAAGCGAATACCTTCATTCAGGCGTGGCAGACTGTTATTGATTAAATT  >  CP009273/968269‑968560

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