Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,959,915 C→A 100% P178P (CCG→CCT flhB ← flagellin export apparatus, substrate specificity protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,959,9150CA98.9% 587.1 / ‑1.8 181P178P (CCG→CCTflhBflagellin export apparatus, substrate specificity protein
Reads supporting (aligned to +/- strand):  ref base C (2/0);  new base A (89/90);  total (91/90)
Fisher's exact test for biased strand distribution p-value = 4.97e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.92e-01
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

Only 100 of 181 total aligned reads displayed.

GAATATCCTGCCGTGACATACGCAGCTTTTTCAGGTGGCTGAATATTTGGAAAAAGACGTCAAATCCCACCATTGGAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATCGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAGAAAAAACCCCGTCACGCTGCCAACCAGGATGGTTTTCAAAATTGCTTTAAGCAACTCCGCGCCAGTCTGAGCCGAGAACATCCGTTTAATGCCCGG  >  CP009273/1959770‑1960058
                                                                                                                                                 |                                                                                                                                               
  atatCCTGCCGTGACATACGCAGCTTTTTCAGGTGGCTGAATATTTGGAAAAAGACGTCAAATCCCACCATTGGAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATAGgag                                                                                                                                             <  1:2499969/148‑1 (MQ=255)
  atatCCTGCCGTGACATACGCAGCTTTTTCAGGTGGCTGAATATTTGGAAAAAGACGTCAAATCCCACCATTGGAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATAGgag                                                                                                                                             <  2:678801/148‑1 (MQ=255)
   tatCCTGCCGTGACATACGCAGCTTTTTCAGGTGGCTGAATATTTGGAAAAAGACGTCAAATCCCACCATTGGAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATAGgaga                                                                                                                                            <  1:1003813/148‑1 (MQ=255)
   tatCCTGCCGTGACATACGCAGCTTTTTCAGGTGGCTGAATATTTGGAAAAAGACGTCAAATCCCACCATTGGAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATAGgaga                                                                                                                                            <  2:1010059/148‑1 (MQ=255)
    atCCTGCCGTGACATACGCAGCTTTTTCAGGTGGCTGAATATTTGGAAAAAGACGTCAAATCCCACCATTGGAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATAGGAGAc                                                                                                                                           >  1:323303/1‑148 (MQ=255)
    atCCTGCCGTGACATACGCAGCTTTTTCAGGTGGCTGAATATTTGGAAAAAGACGTCAAATCCCACCATTGGAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATAGGAGAc                                                                                                                                           >  2:3247118/1‑148 (MQ=255)
      ccTGCCGTGACATACGCAGCTTTTTCAGGTGGCTGAATATTTGGAAAAAGACGTCAAATCCCACCATTGGAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATAGGAGACTc                                                                                                                                         >  1:2916178/1‑148 (MQ=255)
        tGCCGTGACATACGCAGCTTTTTCAGGTGGCTGAATATTTGGAAAAAGACGTCAAATCCCACCATTGGAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATAGGAGACTCgg                                                                                                                                       <  1:964288/148‑1 (MQ=255)
        tGCCGTGACATACGCAGCTTTTTCAGGTGGCTGAATATTTGGAAAAAGACGTCAAATCCCACCATTGGAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATAGGAGACTCgg                                                                                                                                       >  1:218160/1‑148 (MQ=255)
               aCATACGCAGCTTTTTCAGGTGGCTGAATATTTGGAAAAAGACGTCAAATCCCACCATTGGAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATAGGAGACTCGGCCATCaa                                                                                                                                <  1:841365/148‑1 (MQ=255)
               aCATACGCAGCTTTTTCAGGTGGCTGAATATTTGGAAAAAGACGTCAAATCCCACCATTGGAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATAGGAGACTCGGCCATCaa                                                                                                                                >  1:1286691/1‑148 (MQ=255)
                 aTACGCAGCTTTTTCAGGTGGCTGAATATTTGGAAAAAGACGTCAAATCCCACCATTGGAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATAGGAGACTCGGCCATCAAgc                                                                                                                              >  2:1917997/1‑148 (MQ=255)
                     gCAGCTTTTTCAGGTGGCTGAATATTTGGAAAAAGACGTCAAATCCCACCATTGGAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATAGGAGACTCGGCCATCAAGCGcat                                                                                                                          <  2:432929/148‑1 (MQ=255)
                      cAGCTTTTTCAGGTGGCTGAATATTTGGAAAAAGACGTCAAATCCCACCATTGGAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATAGGAGACTCGGCCATCAAGCGcatc                                                                                                                         <  2:475850/148‑1 (MQ=255)
                       aGCTTTTTCAGGTGGCTGAATATTTGGAAAAAGACGTCAAATCCCACCATTGGAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATAGGAGACTCGGCCATCAAGCGcatca                                                                                                                        <  1:2780239/148‑1 (MQ=255)
                        gCTTTTTCAGGTGGCTGAATATTTGGAAAAAGACGTCAAATCCCACCATTGGAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATAGGAGACTCGGCCATCAAGCGcatcat                                                                                                                       <  2:2895532/148‑1 (MQ=255)
                        gCTTTTTCAGGTGGCTGAATATTTGGAAAAAGACGTCAAATCCCACCATTGGAATGACAACAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATAGGAGACTCGGCCATCAAGCGcatcat                                                                                                                       <  2:2976396/148‑1 (MQ=255)
                         cTTTTTCAGGTGGCTGAATATTTGGAAAAAGACGTCAAATCCCACCATTGGAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATAGGAGACTCGGCCATCAAGCGCatcatc                                                                                                                      <  2:2209054/148‑1 (MQ=255)
                           ttttCAGGTGGCTGAATATTTGGAAAAAGACGTCAAATCCCACCATTGGAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATAGGAGACTCGGCCATCAAGCGCATCATCTg                                                                                                                    <  1:2094777/148‑1 (MQ=255)
                           ttttCAGGTGGCTGAATATTTGGAAAAAGACGTCAAATCCCACCATTGGAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATAGGAGACTCGGCCATCAAGCGCATCATCTg                                                                                                                    >  2:586432/1‑148 (MQ=255)
                                 ggTGGCTGAATATTTGGAAAAAGACGTCAAATCCAACCATTGGAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATAGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCa                                                                                                              >  1:586799/1‑148 (MQ=255)
                                  gTGGCTGAATATTTGGAAAAAGACGTCAAATCCCACCATTGGAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATAGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAg                                                                                                             <  1:2950522/148‑1 (MQ=255)
                                  gTGGCTGAATATTTGGAAAAAGACGTCAAATCCCACCATTGGAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATAGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAg                                                                                                             >  1:1113659/1‑148 (MQ=255)
                                   tGGCTGAATATTTGGAAAAAGACGTCAAATCCCACCATTGGAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCAGCGCATTACCCATGGCGGTAATAGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGt                                                                                                            <  2:2538729/148‑1 (MQ=255)
                                    ggCTGAATATTTGGAAAAAGACGTCAAATCCCACCATTGGAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATAGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTg                                                                                                           >  1:976756/1‑148 (MQ=255)
                                    ggCTGAATATTTGGAAAAAGACGTCAAATCCCACCATTGGAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATAGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTg                                                                                                           >  1:1923325/1‑148 (MQ=255)
                                       tGAATATTTGGAAAAAGACGTCAAATCCCACCATTGGAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATAGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATg                                                                                                        >  1:2168859/1‑148 (MQ=255)
                                       tGAATATTTGGAAAAAGACGTCAAATCCCACCATTGGAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATAGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATg                                                                                                        >  1:3540560/1‑148 (MQ=255)
                                       tGAATATTTGGAAAAAGACGTCAAATCCCACCATTGGAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATAGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATg                                                                                                        <  2:1205051/148‑1 (MQ=255)
                                           tatTTGGAAAAAGACGTCAAATCCCACCATTGGAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATAGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCag                                                                                                    <  2:1862187/148‑1 (MQ=255)
                                           tatTTGGAAAAAGACGTCAAATCCCACCATTGGAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATAGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCag                                                                                                    >  2:2538608/1‑148 (MQ=255)
                                           tatTTGGAAAAAGACGTCAAATCCCACCATTGGAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATAGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCag                                                                                                    >  2:985901/1‑148 (MQ=255)
                                           tatTTGGAAAAAGACGTCAAATCCCACCATTGGAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATAGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCag                                                                                                    >  1:3366994/1‑148 (MQ=255)
                                            atTTGGAAAAAGACGTCAAATCCCACCATTGGAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATAGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCaga                                                                                                   >  2:2705440/1‑148 (MQ=255)
                                             tttGGAAAAAGACGTCAAATCCCACCATTGGAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATAGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCagag                                                                                                  <  2:3190375/148‑1 (MQ=255)
                                             tttGGAAAAAGACGTCAAATCCCACCATTGGAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATAGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCagag                                                                                                  >  1:1544334/1‑148 (MQ=255)
                                             tttGGAAAAAGACGTCAAATCCCACCATTGGAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATAGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCagag                                                                                                  <  2:1762757/148‑1 (MQ=255)
                                              ttGGAAAAAGACGTCAAATCCCACCATTGGAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATAGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAgaga                                                                                                 >  2:274301/1‑148 (MQ=255)
                                                 gAAAAAGACGTCAAATCCCACCATTGGAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATAGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAGaaaa                                                                                              >  1:3100236/1‑148 (MQ=255)
                                                   aaaaGACGTCAAATCCCACCATTGGAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATAGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAGaaaaaa                                                                                            <  1:3524764/148‑1 (MQ=255)
                                                   aaaaGACGTCAAATCCCACCATTGGAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATAGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAGaaaaaa                                                                                            >  2:522845/1‑148 (MQ=255)
                                                   aaaaGACGTCAAATCCCACCATTGGAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATAGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAGaaaaaa                                                                                            >  1:1750862/1‑148 (MQ=255)
                                                      aGACGTCAAATCCCACCATTGGAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATAGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAGAAAAAAccc                                                                                         <  1:2167582/148‑1 (MQ=255)
                                                       gACGTCAAATCCCACCATTGGAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATAGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAGAAAAAAcccc                                                                                        >  2:3303722/1‑148 (MQ=255)
                                                        aCGTCAAATCCCACCATTGGAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATAGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAGAAAAAACCCCg                                                                                       >  2:3415474/1‑148 (MQ=255)
                                                        aCGTCAAATCCCACCATTGGAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATAGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAGAAAAAACCCCg                                                                                       >  2:863166/1‑148 (MQ=255)
                                                               aTCCCACCATTGGAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATCGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAGAAAAAACCCCGTCACGCt                                                                                >  2:2785304/1‑148 (MQ=255)
                                                                tCCCACCATTGGAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATAGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAGAAAAAACCCCGTCACGCTg                                                                               <  2:1544334/148‑1 (MQ=255)
                                                                 cccACCATTGGAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCCTTACCCATGGCGGTAATAGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAGAAAAAACCCCGTCACGCTGc                                                                              <  1:1061276/148‑1 (MQ=255)
                                                                    accaTTGGAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATAGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAGAAAAAACCCCGTCACGCTGCCaa                                                                           <  2:402773/148‑1 (MQ=255)
                                                                     ccaTTGGAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATAGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAGAAAAAACCCCGTCACGCTGCCAAc                                                                          <  2:2432509/148‑1 (MQ=255)
                                                                     ccaTTGGAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATAGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAGAAAAAACCCCGTCACGCTGCCAAc                                                                          >  2:1944670/1‑148 (MQ=255)
                                                                      caTTGGAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATAGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAGAAAAAACCCCGTCACGCTGCCAAcc                                                                         >  1:2043866/1‑148 (MQ=255)
                                                                      caTTGGAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATAGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAGAAAAAACCCCGTCACGCTGCCAAcc                                                                         <  1:1743688/148‑1 (MQ=255)
                                                                       aTTGGAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATAGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAGAAAAAACCCCGTCACGCTGCCAACCa                                                                        <  1:1406285/148‑1 (MQ=255)
                                                                       aTTGGAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATAGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAGAAAAAACCCCGTCACGCTGCCAACCa                                                                        <  1:1406293/148‑1 (MQ=255)
                                                                        ttGGAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATAGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAGAAAAAACCCCGTCACGCTGCCAACCAg                                                                       >  2:2042921/1‑148 (MQ=255)
                                                                           gAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATAGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAGAAAAAACCCCGTCACGCTGCCAACCAGGAt                                                                    <  2:1899794/148‑1 (MQ=255)
                                                                           gAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATAGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAGAAAAAACCCCGTCACGCTGCCAACCAGGAt                                                                    <  2:726529/148‑1 (MQ=255)
                                                                           gAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATAGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAGAAAAAACCCCGTCACGCTGCCAACCAGGAt                                                                    <  2:726469/148‑1 (MQ=255)
                                                                             aTGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATAGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAGAAAAAACCCCGTCACGCTGCCAACCAGGATgg                                                                  >  1:2891367/1‑148 (MQ=255)
                                                                                acacCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATAGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAGAAAAAACCCCGTCACGCTGCCAACCAGGATGGttt                                                               <  1:2861786/148‑1 (MQ=255)
                                                                                 cacCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATAGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAGAAAAAACCCCGTCACGCTGCCAACCAGGATGGtttt                                                              >  2:1437042/1‑148 (MQ=255)
                                                                                       gCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATAGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAGAAAAAACCCCGTCACGCTGCCAACAAGGATGGTTTTCAAAAt                                                        <  1:3118699/148‑1 (MQ=255)
                                                                                         accaccAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATAGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAGAAAAAACCCCGTCACGCTGCCAACCAGGATGGTTTTCAAAATTg                                                      <  1:985901/148‑1 (MQ=255)
                                                                                         accaccAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATAGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAGAAAAAACCCCGTCACGCTGCCAACCAGGATGGTTTTCAAAATTg                                                      <  1:314528/148‑1 (MQ=255)
                                                                                         accaccAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATAGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAGAAAAAACCCCGTCACGCTGCCAACCAGGATGGTTTTCAAAATTg                                                      <  1:314729/148‑1 (MQ=255)
                                                                                         accaccAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATAGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAGAAAAAACCCCGTCACGCTGCCAACCAGGATGGTTTTCAAAATTg                                                      <  1:2538608/148‑1 (MQ=255)
                                                                                            accaGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATAGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAGAAAAAACCCCGTCACGCTGCCAACCAGGATGGTTTTCAAAATTGCtt                                                   <  2:2535669/148‑1 (MQ=255)
                                                                                            accaGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATAGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAGAAAAAACCCCGTCACGCTGCCAACCAGGATGGTTTTCAAAATTGCtt                                                   <  1:3303722/148‑1 (MQ=255)
                                                                                                 cagTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATAGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAGAAAAAACCCCGTCACGCTGCCAACCAGGATGGTTTTCAAAATTGCTTTAAGc                                              <  2:781219/148‑1 (MQ=255)
                                                                                                 cagTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATAGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAGAAAAAACCCCGTCACGCTGCCAACCAGGATGGTTTTCAAAATTGCTTTAAGc                                              <  2:2563307/148‑1 (MQ=255)
                                                                                                  agTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATAGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAGAAAAAACCCCGTCACGCTGCCAACCAGGATGGTTTTCAAAATTGCTTTAAGCa                                             <  2:1651418/148‑1 (MQ=255)
                                                                                                     gcgcATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATAGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAGAAAAAACCCCGTCACGCTGCCAACCAGGATGGTTTTCAAAATTGCTTTAAGCAACt                                          <  2:1644574/148‑1 (MQ=255)
                                                                                                       gcATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATAGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAGAAAAAACCCCGTCACGCTGCCAACCAGGATGGTTTTCAAAATTGCTTTAAGCAACTcc                                        <  1:1756215/148‑1 (MQ=255)
                                                                                                        cATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATAGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAGAAAAAACCCCGTCACGCTGCCAACCAGGATGGTTTTCAAAATTGCTTTAAGCAACTCcg                                       >  2:2379583/1‑148 (MQ=255)
                                                                                                          tAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATAGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAGAAAAAACCCCGTCACGCTGCCAACCAGGATGGTTTTCAAAATTGCTTTAAGCAACTCcgcg                                     >  2:1956338/1‑148 (MQ=255)
                                                                                                          tAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATAGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAGAAAAAACCCCGTCACGCTGCCAACCAGGATGGTTTTCAAAATTGCTTTAAGCAACTCcgcg                                     <  2:3544133/148‑1 (MQ=255)
                                                                                                             cccTACCAAATCCATCGCATTACCCATGGCGGTAATAGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAGAAAAAACCCCGTCACGCTGCCAACCAGGATGGTTTTCAAAATTGCTTTAAGCAACTCCGCGCCa                                  >  2:2380239/1‑148 (MQ=255)
                                                                                                             cccTACCAAATCCATCGCATTACCCATGGCGGTAATAGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAGAAAAAACCCCGTCACGCTGCCAACCAGGATGGTTTTCAAAATTGCTTTAAGCAACTCCGCGCCa                                  >  1:2667556/1‑148 (MQ=255)
                                                                                                               cTACCAAATCCATCGAATTACCCATGGCGGTAATAGGAGACTGGGCCATAAAGCGCATCATCTGCGGCCTGTGATGCCAGAGAAAAAACCCCGTCACGCTGCAAACCTGTTTGTTTTTAAAAATTGCTTTAAGCAACTCCGCGCCAGt                                <  1:10185/148‑1 (MQ=255)
                                                                                                                  ccAAATCCATCGCATTACCCATGGCGGTAATAGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAGAAAAAACCCCGTCACGCTGCCAACCAGGATGGTTTTCAAAATTGCTTTAAGCAACTCCGCGCCAGTCTg                             <  2:2848698/148‑1 (MQ=255)
                                                                                                                  ccAAATCCATCGCATTACCCATGGCGGTAATAGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAGAAAAAACCCCGTCACGCTGCCAACCAGGATGGTTTTCAAAATTGCTTTAAGCAACTCCGCGCCAGTCTg                             <  2:785308/148‑1 (MQ=255)
                                                                                                                   cAAATCCATCGCATTACCCATGGCGGTAATAGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAGAAAAAACCCCGTCACGCTGCCAACCAGGATGGTTTTCAAAATTGCTTTAAGCAACTCCGCGCCAGTCTGa                            >  1:2117113/1‑148 (MQ=255)
                                                                                                                   cAAATCCATCGCATTACCCATGGCGGTAATAGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAGAAAAAACCCCGTCACGCTGCCAACCAGGATGGTTTTCAAAATTGCTTTAAGCAACTCCGCGCCAGTCTGa                            >  1:2117133/1‑148 (MQ=255)
                                                                                                                   cAAATCCATCGCATTACCCATGGCGGTAATAGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAGAAAAAACCCCGTCACGCTGCCAACCAGGATGGTTTTCAAAATTGCTTTAAGCAACTCCGCGCCAGTCTGa                            >  1:239487/1‑148 (MQ=255)
                                                                                                                    aaaTCCATCGCATTACCCATGGCGGTAATAGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAGAAAAAACCCCGTCACGCTGCCAACCAGGATGGTTTTCAAAATTGCTTTAAGCAACTCCGCGCCAGTCTGAg                           >  2:1167252/1‑148 (MQ=255)
                                                                                                                    aaaTCCATCGCATTACCCATGGCGGTAATAGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAGAAAAAACCCCGTCACGCTGCCAACCAGGATGGTTTTCAAAATTGCTTTAAGCAACTCCGCGCCAGTCTGAg                           >  2:491031/1‑148 (MQ=255)
                                                                                                                    aaaTCCATCGCATTACCCATGGCGGTAATAGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAGAAAAAACCCCGTCACGCTGCCAACCAGGATGGTTTTCAAAATTGCTTTAAGCAACTCCGCGCCAGTCTGAg                           >  2:492594/1‑148 (MQ=255)
                                                                                                                     aaTCCATCGCATTACCCATGGCGGTAATAGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAGAAAAAACCCCGTCACGCTGCCAACCAGGATGGTTTTCAAAATTGCTTTAAGCAACTCCGCGCCAGTCTGAGc                          >  2:2673187/1‑148 (MQ=255)
                                                                                                                       tCCATCGCATTACCCATGGCGGTAATAGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAGAAAAAACCCCGTCACGCTGCCAACCAGGATGGTTTTCAAAATTGCTTTAAGCAACTCCGCGCCAGTCTGAGCCg                        >  2:605046/1‑148 (MQ=255)
                                                                                                                         cATCGCATTACCCATGGCGGTAATAGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAGAAAAAACCCCGTCACGCTGCCAACCAGGATGGTTTTCAAAATTGCTTTAAGCAACTCCGCGCCAGTCTGAGCCgag                      >  2:1087329/1‑148 (MQ=255)
                                                                                                                          aTCGCATTACCCATGGCGGTAATAGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAGAAAAAACCCCGTCACGCTGCCAACCAGGATGGTTTTCAAAATTGCTTTAAGCAACTCCGCGCCAGTCTGAGCCgaga                     <  2:2261162/148‑1 (MQ=255)
                                                                                                                                ttACCCATGGCGGTAATAGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAGAAAAAACCCCGTCACGCTGCCAACCAGGATGGTTTTCAAAATTGCTTTAAGCAACTCCGCGCCAGTCTGAGCCGAGAACATcc               >  1:730103/1‑148 (MQ=255)
                                                                                                                                     cATGGCGGTAATAGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAGAAAAAACCCCGTCACGCTGCCAACCAGGATGGTTTTCAAAATTGCTTTAAGCAACTCCGCGCCAGTCTGAGCCGAGAACATCCGTTTa          >  2:1743348/1‑148 (MQ=255)
                                                                                                                                     cATGGCGGTAATAGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAGAAAAAACCCCGTCACGCTGCCAACCAGGATGGTTTTCAAAATTGCTTTAAGCAACTCCGCGCCAGTCTGAGCCGAGAACATCCGTTTa          <  2:2942093/148‑1 (MQ=255)
                                                                                                                                     cATGGCGGTAATAGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAGAAAAAACCCCGTCACGCTGCCAACCAGGATGGTTTTCAAAATTGCTTTAAGCAACTCCGCGCCAGTCTGAGCCGAGAACATCCGTTTa          <  2:2418371/148‑1 (MQ=255)
                                                                                                                                     cATGGCGGTAATAGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAGAAAAAACCCCGTCACGCTGCCAACCAGGATGGTTTTCAAAATTGCTTTAAGCAACTCCGCGCCAGTCTGAGCCGAGAACATCCGTTTa          >  2:1743381/1‑148 (MQ=255)
                                                                                                                                           ggTAATAGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAGAAAAAACCCCGTCACGCTGCCAACCAGGATGGTTTTCAAAATTGCTTTAAGCAACTCCGCGCCAGTCTGAGCCGAGAACATCCGTTTAATGccc    <  2:3500088/148‑1 (MQ=255)
                                                                                                                                             tAATAGGAGACTCGGCCATCAAGCGCATCATCGGCGGCCAGTGATGCCAGAGAAAAAACCCCGTCACGCTGCCAACCCGGATGGTTTTCAAAATTGCTTTAAGCAACTCCGCGCCAGTCTGAGCCGAGAACATCCGTTTAATGCCCgg  <  1:1653499/148‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                               
GAATATCCTGCCGTGACATACGCAGCTTTTTCAGGTGGCTGAATATTTGGAAAAAGACGTCAAATCCCACCATTGGAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATCGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAGAAAAAACCCCGTCACGCTGCCAACCAGGATGGTTTTCAAAATTGCTTTAAGCAACTCCGCGCCAGTCTGAGCCGAGAACATCCGTTTAATGCCCGG  >  CP009273/1959770‑1960058

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: