Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,582,677 G→T 9.0% E402* (GAA→TAA)  tsr → methyl‑accepting chemotaxis protein I, serine sensor receptor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,582,6770GT 9.0% 1010.7 / 47.8 273E402* (GAA→TAA) tsrmethyl‑accepting chemotaxis protein I, serine sensor receptor
Reads supporting (aligned to +/- strand):  ref base G (133/115);  new base T (15/10);  total (148/125)
Fisher's exact test for biased strand distribution p-value = 6.75e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 279 total aligned reads displayed.

GCGATATCTCCACCAGTTCGCAGAAAATCGCCGATATTATCAGCGTAATTGACGGCATTGCCTTCCAGACCAATATTCTGGCTTTGAACGCGGCGGTTGAGGCTGCGCGTGCGGGTGAGCAAGGGCGCGGTTTTGCGGTGGTCGCGGGAGAAGTGCGTAATCTGGCCCAGCGCAGCGCCCAGGCGGCTCGTGAAATTAAAAGCCTGATTGAAGACTCGGTGGGGAAAGTGGATGTTGGCTCTACGCTGGTCGAAAGCGCCGGGGAAACAATGGCGGAGATTGTCAGCGCCGTGACCCGC  >  CP009273/4582528‑4582826
                                                                                                                                                     |                                                                                                                                                     
 cgATATCTCCACCAGTTCGCAGAAAATCGCCGATATTATCAGCGTAATTGACGGCATTGCCTTCCAGACCAATATTCTGGCTTTGAACGCGGCGGTTGAGGCTGCGCGTGCGGGTGAGCAAGGGCGCGGTTTTGCGGTGGTCGCGGgaga                                                                                                                                                      >  2:4299390/1‑150 (MQ=255)
    tatCTCCACCAGTTCGCAGAAAATCGCCGATATTATCAGCGTAATTGACGGCATTGCCTTCCAGACCAATATTCTGGCTTTGAACGCGGCGGTTGAGGCTGCGCGTGCGGGTGAGCAAGGGCGCGGTTTTGCGGTGGTCGCGGGAGAAGt                                                                                                                                                   >  2:2477170/1‑150 (MQ=255)
     atCTCCACCAGTTCGCAGAAAATCGCCGATATTATCAGCGTAATTGACGGCATTGCCTTCCAGACCAATATTCTGGCTTTGAACGCGGCGGTTGAGGCTGCGCGTGCGGGTGAGCAAGGGCGCGGTTTTGCGGTGGTCGCGGGAGAAGTg                                                                                                                                                  >  1:276893/1‑150 (MQ=255)
      tctcCACCAGTTCGCAGAAAATCGCCGATATTATCAGCGTAATTGACGGCATTGCCTTCCAGACCAATATTCTGGCTTTGAACGCGGCGGTTGAGGCTGCGCGTGCGGGTGAGCAAGGGCGCGGTTTTGCGGTGGTCGCGGGAGAAGTGc                                                                                                                                                 >  1:4315818/1‑150 (MQ=255)
        tcCACCAGTTCGCAGAAAATCGCCGATATTATCAGCGTAATTGACGGCATTGCCTTCCAGACCAATATTCTGGCTTTGAACGCGGCGGTTGAGGCTGCGCGTGCGGGTGAGCAAGGGCGCGGTTTTGCGGTGGTCGCGGGAGAAGTGCGt                                                                                                                                               >  2:1145712/1‑150 (MQ=255)
        tcCACCAGTTCGCAGAAAATCGCCGATATTATCAGCGTAATTGACGGCATTGCCTTCCAGACCAATATTCTGGCTTTGAACGCGGCGGTTGAGGCTGCGCGTGCGGGTGAGCAAGGGCGCGGTTTTGCGGTGGTCGCGGGAGAAGTGCGt                                                                                                                                               >  2:2512242/1‑150 (MQ=255)
           accaGTTCGCAGAAAATCGCCGATATTATCAGCGTAATTGACGGCATTGCCTTCCAGACCAATATTCTGGCTTTGAACGCGGCGGTTGAGGCTGCGCGTGCGGGTGAGCAAGGGCGCGGTTTTGCGGTGGTCGCGGGAGAAGTGCGTAAt                                                                                                                                            <  1:392037/150‑1 (MQ=255)
               gTTCGCAGAAAATCGCCGATATTATCAGCGTAATTGACGGCATTGCCTTCCAGACCAATATTCTGGCTTTGAACGCGGCGGTTGAGGCTGCGCGTGCGGGTGAGCAAGGGCGCGGTTTTGCGGTGGTCGCGGGAGAAGTGCGTAATCTgg                                                                                                                                        <  2:1015628/150‑1 (MQ=255)
                 tCGCAGAAAATCGCCGATATTATCAGCGTAATTGACGGCATTGCCTTCCAGACCAATATTCTGGCTTTGAACGCGGCGGTTGAGGCTGCGCGTGCGGGTGAGCAAGGGCGCGGTTTTGCGGTGGTCGCGGGAGAAGTGCGTAATCTGGcc                                                                                                                                      <  1:1087718/150‑1 (MQ=255)
                 tCGCAGAAAATCGCCGATATTATCAGCGTAATTGACGGCATTGCCTTCCAGACCAATATTCTGGCTTTGAACGCGGCGGTTGAGGCTGCGCGTGCGGGTGAGCAAGGGCGCGGTTTTGCGGTGGTCGCGGGAGAAGTGCGTAATCTGGcc                                                                                                                                      >  2:5048555/1‑150 (MQ=255)
                  cGCAGAAAATCGCCGATATTATCAGCGTAATTGACGGCATTGCCTTCCAGACCAATATTCTGGCTTTGAACGCGGCGGTTGAGGCTGCGCGTGCGGGTGAGCAAGGGCGCGGTTTTGCGGTGGTCGCGGGAGAAGTGCGTAATCTGGccc                                                                                                                                     <  2:638251/150‑1 (MQ=255)
                  cGCAGAAAATCGCCGATATTATCAGCGTAATTGACGGCATTGCCTTCCAGACCAATATTCTGGCTTTGAACGCGGCGGTTGAGGCTGCGCGTGCGGGTGAGCAAGGGCGCGGTTTTGCGGTGGTCGCGGGAGAAGTGCGTAATCTGGccc                                                                                                                                     <  2:682876/150‑1 (MQ=255)
                   gCAGAAAATCGCCGATATTATCAGCGTAATTGACGGCATTGCCTTCCAGACCAATATTCTGGCTTTGAACGCGGCGGTTGAGGCTGCGCGTGCGGGTGAGCAAGGGCGCGGTTTTGCGGTGGTCGCGGGAGAAGTGCGTAATCTGGCCca                                                                                                                                    <  1:250843/150‑1 (MQ=255)
                     aGAAAATCGCCGATATTATCAGCGTAATTGACGGCATTGCCTTCCAGACCAATATTCTGGCTTTGAACGCGGCGGTTGAGGCTGCGCGTGCGGGTGAGCAAGGGCGCGGTTTTGCGGTGGTCGCGGGAGAAGTGCGTAATCTGGCCcagc                                                                                                                                  <  2:4499561/150‑1 (MQ=255)
                       aaaaTCGCCGATATTAACAGCGTAATTGACGGCATTGCCTTCCAGACCAATATTCTGGCTTTGAACGCGGCGGTTGAGGCTGCGCGTGCGGGTGAGCAAGGGCGCGGTTTTGCGGTGGTCGCGGGAGAAGTGCGTAATCTGGCCcagcgc                                                                                                                                <  1:660559/150‑1 (MQ=255)
                        aaaTCGCATATATTATAAGCGTAATTGACGGCATTGCCTTCCAGACCAATATTCTGTCTTTGAACGCGGCGGTTGAGGCAGCGAGTGCGGGTGAGAAAGGGCGCGGTTTTGCGGTGGTCGCGGGAGAAGTGAGTAATCTGGCCaagctca                                                                                                                               >  1:3948822/1‑150 (MQ=255)
                         aaTCGCCGATATTATCAGCGTAATTGACGGCATTGCCTTCCAGACCAATATTCTGGCTTTGAACGCGGCGGTTGAGGCTGCGCGTGCGGGTGAGCAAGGGCGCGGTTTTGCGGTGGTCGCGGGAGAAGTGCGTAATCTGGCCcagcgcag                                                                                                                              <  1:3498355/150‑1 (MQ=255)
                           tCGCCGATATTATCAGCGTAATTGACGGCATTGCCTTCCAGACCAATATTCTGGCTTTGAACGCGGCGGGTGAGGCTGCGCGTGCGGGTGAGCAAGGGCGCGGTTTTGCGGTGGTCGCGGGAGAAGTGCGTAATCTGGCCcagcgcagcg                                                                                                                            >  2:5072611/1‑150 (MQ=255)
                            cGCCGATATTATCAGCGTAATTGACGGCATTGCCTTCCAGACCAATATTCTGGCTTTGAACGCGGCGGTTGAGGCTGCGCGTGCGGGTGAGCAAGGGCGCGGTTTTGCGGTGGTCGCGGGATAAGTGCGTAATCTGGCCCagcgcagcgc                                                                                                                           >  1:155343/1‑150 (MQ=255)
                             gCCGATATTATCAGCGTAATTGACGGCATTGCCTTCCAGACCAATATTCTGGCTTTGAACGCGGCGGTTGAGGCTGCGCGTGCGGGTGAGCAAGGGCGCGGTTTTGCGGTGGTCGCGGGAGAAGTGCGTAATCTGGCCCAGCGCAGCGcc                                                                                                                          <  2:3394353/150‑1 (MQ=255)
                              ccGATAATATCAGCGTAATTGACGGCATTGCCTTCCAGACCAATATTCTGGCTTTGAACGCGGCGGTTGAGGCTGCGCGTGCGGGTGAGCAAGGGCGCGGTTTTGCGGTGGTCGCGGGAGAAGTGCGTAATCTGGCCCAGCGCAGCGccc                                                                                                                         >  2:76964/1‑150 (MQ=255)
                               cGATATTATCAGCGTAATTGACGGCATTGCCTTCCAGACCAATATTCTGGCTTTGAACGCGGCGGTTGAGGCTGCGCGTGCGGGTGAGCAAGGGCGCTGTTTTGCGGTGGTCGCGGGAGAAGTGAGTAATCTGGCCCAGCGCAGCGCCCa                                                                                                                        <  2:1852247/150‑1 (MQ=255)
                               cGATATTATCAGCGTAATTGACGGCATTGCCTTCCAGACCAATATTCTGGCTTTGAACGCGGCGGTTGAGGCTGCGCGTGCGGGTGAGCAAGGGCGCGGTTTTGCGGTGGTCGCGGGAGAAGTGCGTAATCTGGCCCAGCGCAGCGCCCa                                                                                                                        >  1:275181/1‑150 (MQ=255)
                                gATATTATCAGCGTAATTGACGGCATTGCCTTCCAGACCAATATTCTGGCTTTGAACGCGGCGGTTGAGGCTGCGCGTGCGGGTGAGCAAGGGCGCGGTTTTGCGGTGGTCGCGGGAGAAGTGCGTAATCTGGCCCAGCGCAGCGCCCAg                                                                                                                       >  1:2409577/1‑150 (MQ=255)
                                 atatTATCAGCGTAATTGACGGCATTGCCTTCCAGACCAATATTCTGGCTTTGAACGCGGCGGTTGAGGCTGCGCGTGCGGGTGAGCAAGGGCGCGGTTTTGCGGTGGTCGCGGGAGAAGTGCGTAATCTGGCCCAGCGCAGCGCCCAgg                                                                                                                      <  1:3613478/150‑1 (MQ=255)
                                  tattatCAGCGTAATTGACGGCATTGCCTTCCAGACCAATATTCTGGCTTTGAACGCGGCGGTTGAGGCTGCGCGTGCGGGTGAGCAAGGGCGCGGTTTTGCGGTGGTCGCGGGATAAGTGCGTAATCTGGCCCAGCGCAGCGCCCAggc                                                                                                                     <  1:4741964/150‑1 (MQ=255)
                                     tatCAGCGTAATTGACGGCATTGCCTTCCAGACCAATATTCTGGCTTTGAACGCGGCGGTTGAGGCTGCGCGTGCGGGTGAGCAAGGGCGCGGTTTTGCGGTGGTCGCGGGATAAGTGCGTAATCTGGCCCAGCGCAGCGCCCAggcggc                                                                                                                  <  2:976210/150‑1 (MQ=255)
                                     tatCAGCGTAATTGACGGCATTGCCTTCCAGACCAATATTCTGGCTTTGAACGCGGCGGTTGAGGCTGCGCGTGCGGGTGAGCAAGGGCGCGGTTTTGCGGTGGTCGCGGGAGAAGTGCGTAATCTGGCCCAGCGCAGCGCCCAggcggc                                                                                                                  >  2:710886/1‑150 (MQ=255)
                                      atCAGCGTAATTGACGGCATTGCCTTCCAGACCAATATTCTGGCTTTGAACGCGGCGGTTGAGGCTGCGCGTGCGGGTGAGCAAGGGCGCGGTTTTGCGGTGGTCGCGGGAGAAGTGCGTAATCTGGCCCAGCGCAGCGCCCAGGCGGCt                                                                                                                 >  2:2217462/1‑150 (MQ=255)
                                       tCAGCGTAATTGACGGCATTGCCTTCCAGACCAATATTCTGGCTTTGAACGCGGCGGTTGAGGCTGCGCGTGCGGGTGAGCAAGGGCGCGGTTTTGCGGTGGTCGCGGGAGAAGTGCGTACTCTGGCCCAGCGCAGCGCCCAGGCGGCTc                                                                                                                >  2:1419187/1‑150 (MQ=255)
                                       tCAGCGTAATTGACGGCATTGCCTTCCAGACCAATATTCTGGCTTTGAACGCGGCGGTTGAGGCTGCGCGTGCGGGTGAGCAAGGGCGCGGTTTTGCGGTGGTCGCGGGAGAAGTGCGTAATCTGGCCCAGCGCAGCGCCCAGGCGGCTc                                                                                                                >  2:4535475/1‑150 (MQ=255)
                                        cAGCGTAATTGACGGCATTGCCTTCCAGACCAATATTCTGGCTTTGAACGCGGCGGTTGAGGCTGCGCGTGCGGGTGAGCAAGGGCGCGGTTTTGCGGTGGTCGCGGGAGAAGTGCGTAATCTGGCCCAGCGCAGCGCCCAGGCGGCTCg                                                                                                               <  2:4246181/150‑1 (MQ=255)
                                              aaTTGACGGCATTGCCTTCCAGACCAATATTCTGGCTTTGAACGCGGCGGTTGAGGCTGCGCGTGCGGGTGAGCAAGGGCGCGGTTTTGCGGTGGTCGCGGGATAAGTGCGTAATCTGGCCCAGCGCAGCGCCCAGGCGGCTCGTGAAAt                                                                                                         >  1:4285887/1‑150 (MQ=255)
                                               aTTGACGGCATTGCCTTCCAGACCAATATTCTGGCTTTGAACGCGGCGGTTGAGGCTGCGCGTGCGGGTGAGCAAGGGCGCGGTTTTGCGGTGGTCGCGGGAGAAGTGCGTAATCTGGCCCAGCGCAGCGCCCAGGCGGCTCGTGAAAtt                                                                                                        >  2:3904365/1‑150 (MQ=255)
                                               aTTGACGGCATTGCCTTCCAGACCAATATTCTGGCTTTGAACGCGGCGGTTGAGGCTGCGCGTGCGGGTGAGCAAGGGCGCGGTTTTGCGGTGGTCGCGGGAGAAGTGCGTAATCTGGCCCAGCGCAGCGCCCAGGCGGCTCGTGAAAtt                                                                                                        >  1:3366804/1‑150 (MQ=255)
                                                  gACGGCATTGCCTTCCAGACCAATATTCTGGCTTTGAACGCGGCGGTTGAGGCTGCGCGTGCGGGTGAGCAAGGGCGCGGTTTTGCGGTGGTCGCGGGAGAAGTGCGTAATCTGGCCCAGCGCAGCGCCCAGGCGGCTCGTGAAATTaaa                                                                                                     <  2:1088512/150‑1 (MQ=255)
                                                   aCGGCATTGCCTTCCAGACCAATATTCTGGCTTTGAACGCGGCGGTTGAGGCTGCGCGTGCGGGTGAGCAAGGGCGCGGTTTTGCGGTGGTCGCGGGAGAAGTGCGTAATCTGGCCCAGCGCAGCGCCCAGGCGGCTCGTGAAATTaaaa                                                                                                    >  1:4023140/1‑150 (MQ=255)
                                                   aCGGCATTGCCTTCCAGACCAATATTCTGGCTTTGAACGCGGCGGTTGAGGCTGCGCGTGCGGGTGAGCAAGGGCGCGGTTTTGCGGTGGTCGCGGGAGAAGTGCGTAATCTGGCCCAGCGCAGCGCCCAGGCGGCTCGTGAAATTaaaa                                                                                                    >  2:1252135/1‑150 (MQ=255)
                                                   aCGGCATTGCCTTCCAGACCAATATTCTGGCTTTGAACGCGGCGGTTGAGGCTGCGCGTGCGGGTGAGCAAGGGCGCGGTTTTGCGGTGGTCGCGGGAGAAGTGCGTAATCTGGCCCAGCGCAGCGCCCAGGCGGCTCGTGAAATTaaaa                                                                                                    >  1:4731773/1‑150 (MQ=255)
                                                    cGGCATTGCCTTCCAGACCAATATTCTGGCTTTGAACGCGGCGGTTGAGGCTGCGCGTGCGGGTGAGCAAGGGCGCGGTTTTGCGGTGGTCGCGGGAGAAGTGCGTAATCTGGCCCAGCGCAGCGCCCAGGCGGCTCGTGAAATTAAAAg                                                                                                   >  2:5135242/1‑150 (MQ=255)
                                                       cATTGCCTTCCAGACCAATATTCTGGCTTTGAACGCGGCGGTTGAGGCTGCGCGTGCGGGTGAGCAAGGGCGCGGTTTTGCGGTGGTCGCGGGATAAGTGCGTAATCTGGCCCAGCGCAGCGCCCAGGCGGCTCGTGAAATTAAAAGCCt                                                                                                >  1:5051418/1‑150 (MQ=255)
                                                       cATTGCCTTCCAGACCAATATTCTGGCTTTGAACGCGGCGGTTGAGGCTGCGCGTGCGGGTGAGCAAGGGCGCGGTTTTGCGGTGGTCGCGGGAGAAGTGCGTAATCTGGCCCAGCGCAGCGCCCAGGCGGCTCGTGAAATTAAAAGCCt                                                                                                >  1:3097198/1‑150 (MQ=255)
                                                        aTTGCCTTCCAGACCAATATTCTGGCTTTGAACGCGGCGGTTGAGGCTGCGCGTGCGGGTGAGCAAGGGCGCGGTTTTGCGGTGGTCGCGGGAGAAGTGCGTAATCTGGCCCAGCGCAGCGCCCAGGCGGCTCGTGAAATTAAAAGCCTg                                                                                               <  1:1736711/150‑1 (MQ=255)
                                                          tGCCTTCCAGACCAATATTCTGGCTTTGAACGCGGCGGTTGAGGCTGCGCGTGCGGGTGAGCAAGGGCGCGGTTTTGCGGTGGTCGCGGGAGAAGTGCGTAATCTGGCCCAGCGCAGCGCCCAGGCGGCTCGTGAAATTAAAAGCCTGAt                                                                                             >  1:1339094/1‑150 (MQ=255)
                                                              ttCCAGACCAATATTCTGGCTTTGAACGCGGCGGTTGAGGCTGCGCGTGCGGGTGAGCAAGGGCGCGGTTTTGCGGTGGTCGCGGGAGAAGTGCGTAATCTGGCCCAGCGCAGCGCCCAGGCGGCTCGTGAAATTAAAAGCCTGATTGaa                                                                                         >  1:3188744/1‑150 (MQ=255)
                                                                ccAGACCAATATTCTGGCTTTGAACGCGGCGGTTGAGGCTGCGCGTGCGGGTGAGCAAGGGCGCGGTTTTGCGGTGGTCGCGGGAGAAGTGCGTAATCTGGCCCAGCGCAGCGCCCAGGCGGCTCGTGAAATTAAAAGCCTGATTGAAGa                                                                                       >  1:4453507/1‑150 (MQ=255)
                                                                 cAGACCAATATTCTGGCTTTGAACGCGGCGGTTGAGGCTGCGCGTGCGGGTGAGCAAGGGCGCGGTTTTGCGGTGGTCGCGGGAGAAGTGCGTAATCTGGCCCAGCGCAGCGCCCAGGCGGCTCGTGAAATTAAAAGCCTGATTGAAGAc                                                                                      >  2:3304047/1‑150 (MQ=255)
                                                                  aGACCAATATTCTGGCTTTGAACGCGGCGGTTGAGGCTGCGCGTGCGGGTGAGCAAGGGCGCGGTTTTGCGGTGGTCGCGGGAGAAGTGCGTAATCTGGCCCAGCGCAGCGCCCAGGCGGCTCGTGAAATTAAAAGCCTGATTGAAGACt                                                                                     >  2:465253/1‑150 (MQ=255)
                                                                    aCCAATATTCTGGCTTTGAACGCGGCGGTTGAGGCTGCGCGTGCGGGTGAGCAAGGGCGCGGTTTTGCGGTGGTCGCGGGAGAAGTGCGTAATCTGGCCCAGCGCAGCGCCCAGGCGGCTCGTGAAATTAAAAGCCTGATTGAAGACTCg                                                                                   >  1:3678175/1‑150 (MQ=255)
                                                                    aCCAATATTCTGGCTTTGAACGCGGCGGTTGAGGCTGCGCGTGCGGGTGAGCAAGGGCGCGGTTTTGCGGTGGTCGCGGGAGAAGTGCGTAATCTGGCCCAGCGCAGCGCCCAGGCGGCTCGTGAAATTAAAAGCCTGATTGAAGACTCg                                                                                   <  2:1612372/150‑1 (MQ=255)
                                                                     ccAATATTCTGGCTTTGAACGCGGCGGTTGAGGCTGCGCGTGCGGGTGAGCAAGGGCGCGGTTTTGCGGTGGTCGCGGGATAAGTGCGTAATCTGGCCCAGCGCAGCGCCCAGGCGGCTCGTGAAATTAAAAGCCTGATTGAAGACTCgg                                                                                  >  1:2777179/1‑150 (MQ=255)
                                                                     ccAATATTCTGGCTTTGAACGCGGCGGTTGAGGCTGCGCGTGCGGGTGAGCAAGGGCGCGGTTTTGCGGTGGTCGCGGGAGAAGTGCGTAATCTGGCCCAGCGCAGCGCCCAGGCGGCTCGTGAAATTAAAAGCCTGATTGACGACTCgg                                                                                  >  2:3097419/1‑150 (MQ=255)
                                                                        atatTCTGGCTTTGAACGCGGCGGTTGAGGCTGCGCGTGCGGGTGAGCAAGGGCGCGGTTTTGCGGTGGTCGCGGGAGAAGTGCGTAATCTGGCCCAGCGCAGCGCCCAGGCGGCTCGTGAAATTAAAAGCCTGATTGAAGACTCGGTgg                                                                               <  2:4912832/150‑1 (MQ=255)
                                                                        atatTCTGGCTTTGAACGCGGCGGTTGAGGCTGCGCGTGCGGGTGAGCAAGGGCGCGGTTTTGCGGTGGTCGCGGGAGAAGTGCGTAATCTGGCCCAGCGCAGCGCCCAGGCGGCTCGTGAAATTAAAAGCCTGATTGAAGACTCGGTgg                                                                               >  1:2592747/1‑150 (MQ=255)
                                                                         tatTCTGGCTTTGAACGCGGCGGTTGAGGCTGCGCGTGCGGGTGAGCAAGGGCGCGGTTTTGCGGTGGTCGCGGGAGAAGTGCGTAATCTGGCCCAGCGCAGCGCCCAGGCGGCTCGTGAAATTAAAAGCCTGATTGAAGACTCGGTggg                                                                              >  2:3890460/1‑150 (MQ=255)
                                                                         tatTCTGGCTTTGAACGCGGCGGTTGAGGCTGCGCGTGCGGGTGAGCAAGGGCGCGGTTTTGCGGTGGTCGCGGGAGAAGTGCGTAATCTGGCCCAGCGCAGCGCCCAGGCGGCTCGTGAAATTAAAAGCCTGATTGAAGACTCGGTggg                                                                              >  1:4748266/1‑150 (MQ=255)
                                                                          atTCTGGCTTTGAACGCGGCGGTTGAGGCTGCGCGTGCGGGTGAGCAAGGGCGCGGTTTTGCGGTGGTCGCGGGAGAAGTGCGTAATCTGGCCCAGCGCAGCGCCCAGGCGGCTCGTGAAATTAAAAGCCTGATTGAAGACTCGGTgggg                                                                             >  2:3785322/1‑150 (MQ=255)
                                                                           ttCTGGCTTTGAACGCGGCGGTTGAGGCTGCGCGTGCGGGTGAGCAAGGGCGCGGTTTTGCGGTGGTCGCGGGAGAAGTGCGTAATCTGGCCCAGCGCAGCGCCCAGGCGGCTCGTGAAATTAAAAGCCTGATTGAAGACTCGGTGGGGa                                                                            <  2:368310/150‑1 (MQ=255)
                                                                           ttCTGGCTTTGAACGCGGCGGTTGAGGCTGCGCGTGCGGGTGAGCAAGGGCGCGGTTTTGCGGTGGTCGCGGGAGAAGTGCGTAATCTGGCCCAGCGCAGCGCCCAGGCGGCTCGTGAAATTAAAAGCCTGATTGAAGACTCGGTGGGGa                                                                            <  2:1387004/150‑1 (MQ=255)
                                                                            tCTGGCTTTGAACGCGGCGGTTGAGGCTGCGCGTGCGGGTGAGCAAGGGCGCGGTTTTGCGGTGGTCGCGGGAGAAGTGCGTAATCTGGCCCAGCGCAGCGCCCAGGCGGCTCGTGAAATTAAAAGCCTGATTGAAGACTCGGTGGGGaa                                                                           <  1:177380/150‑1 (MQ=255)
                                                                            tCTGGCTTTGAACGCGGCGGTTGAGGCTGCGCGTGCGGGTGAGCAAGGGCGCGGTTTTGCGGTGGTCGCGGGAGAAGTGCGTAATCTGGCCCAGCGCAGCGCCCAGGCGGCTCGTGAAATTAAAAGCCTGATTGAAGACTCGGTGGGGaa                                                                           >  1:1567549/1‑150 (MQ=255)
                                                                            tCTGGCTTTGAACGCGGCGGTTGAGGCTGCGCGTGCGGGTGAGCAAGGGCGCGGTTTTGCGGTGGTCGCGGGAGAAGTGCGTAATCTGGCCCAGCGCAGCGCCCAGGCGGCTCGTGAAATTAAAAGCCTGATTGAAGACTCGGTGGGGaa                                                                           <  1:85658/150‑1 (MQ=255)
                                                                             cTGGCTTTGAACGCGGCGGTTGAGGCTGCGCGTGCGGGTGAGCAAGGGCGCGGTTTTGCGGTGGTCGCGGGAGAAGTGCGTAATCTGGCCCAGCGCAGCGCCCAGGCGGCTCGTGAAATTAAAAGCCTGATTGAAGACTCGGTGGGGaaa                                                                          >  2:3335492/1‑150 (MQ=255)
                                                                                gCTTTGAACGCGGCGGTTGAGGCTGCGCGTGCGGGTGAGCAAGGGCGCGGTTTTGCGGTGGTCGCGGGAGAAGTGCGTAATCTGGCCCAGCGCAGCGCCCAGGCGGCTCGTGAAATTAAAAGCCTGATTGAAGACTCGGTGGGGAAAGTg                                                                       <  1:492849/150‑1 (MQ=255)
                                                                                 cTTTGAACGCGGCGGTTGAGGCTGCGCGTGCGGGTGAGCAAGGGCGCGGTTTTGCGGTGGTCGCGGGAGAAGTGCGTAATCTGGCCCAGCGCAGCGCCCAGGCGGCTCGTGAAATTAAAAGCCTGATTGAAGACTCGGTGGGGAAAGTgg                                                                      >  2:2167970/1‑150 (MQ=255)
                                                                                   ttGAACGCGGCGGTTGAGGCTGCGCGTGCGGTTGAGCAAGGGCGCGGTTTTGCTGTGGTCGCGGGAGAAGTGCGTAATCTGGCCCAGCGCAGCGCCCAGGCGGCTCTTGAAATTTAAAGCCTGTTTGAAGACTCGGTGGGGATTGTGGAt                                                                    <  2:4709519/150‑1 (MQ=255)
                                                                                      aaCGCGGCGGTTGAGGCTGCGCGTGCGGGTGAGCAAGGGCGCGGTTTTGCGGTGGTCGCGGGAGAAGTGCGTAATCTGGCCCAGCGCAGCGCCCAGGCGGCTCGTGAAATTAAAAGCCTGATTGAAGACTCGGTGGGGAAAGTGGATGtt                                                                 <  1:3752838/150‑1 (MQ=255)
                                                                                       aCGCGGCGGTTGAGGCTGCGCGTGCGGGTGAGCAAGGGCGCGGTTTTGCGGTGGTCGCGGGAGAAGTGCGTAATCTGGCCCAGCGCAGCGCCCAGGCGGCTCGTGAAATTAAAAGCCTGATTGAAGACTCGGTGGGGAAAGTGGATGTTg                                                                <  2:3949751/150‑1 (MQ=255)
                                                                                           ggcggTTGAGGCTGCGCGTGCGGGTGAGCAAGGGCGCGGTTTTGCGGTGGTCGCGGGAGAAGTGCGTAATCTGGCCCAGCGCAGCGCCCAGGCGGCTCGTGAAATTAAAAGCCTGATTGAAGACTCGGTGGGGAAAGTGGATGTTGGctc                                                            >  1:1668794/1‑150 (MQ=255)
                                                                                             cggTTGAGGCTGCGCGTGCGGGTGAGCAAGGGCGCGGTTTTGCGGTGGTCGCGGGAGAAGTGCGTAATCTGGCCCAGCGCAGCGCCCAGGCGGCTCGTGAAATTAAAAGCCTGATTGAAGACTCGGTGGGGAAAGTGGATGTTGGCTCTa                                                          <  1:1398478/150‑1 (MQ=255)
                                                                                              ggTTGAGGCTGCGCGTGCGGGTGAGCAAGGGCGCGGTTTTGCGGTGGTCGCGGGAGAAGTCCGTAATCTGGCCCAGCGCAGCGCCCAGGCGGCTCGTGAAATTAAAAGCCTGATTGAAGACTCGGTGGGGAAAGTGGATGTTGGCTCTAc                                                         >  2:1449091/1‑150 (MQ=255)
                                                                                                   aGGCTGCGCGTGCGGGTGAGCAAGGGCGCGGTTTTGCGGTGGTCGCGGGAGAAGTGCGTAATCTGGCCCAGCGCAGCGCCCAGGCGGCTCGTGAAATTAAAAGCCTGATTGAAGACTCGGTGGGGAAAGTGGATGTTGGCTCTACGCTgg                                                    <  1:4276659/150‑1 (MQ=255)
                                                                                                    ggCTGCGCGTGCGGGTGAGCTAGGGCGCGGTTTTGCGGTGGTCGCGGGAGAAGTGCGTAATCTGGCCCAGCGCAGCGCCCAGGCGGCTCGTGAAATTAAAAGCCTGATTGAAGACTCGGTGGGGAAAGTGGATGTTGGCTCTACGCTGGt                                                   <  2:2855342/150‑1 (MQ=255)
                                                                                                    ggCTGCGCGTGCGGGTGAGCAAGGGCGCGGTTTTGCGGTGGTCGCGGGAGAAGTGCGTAATCTGGCCCAGCGCAGCGCCCAGGCGGCTCGTGAAATTAAAAGCCTGATTGAAGACTCGGTGGGGAAAGTGGATGTTGGCTCTACGCTGGt                                                   <  2:1377324/150‑1 (MQ=255)
                                                                                                          gcgTGCGGGTGAGCAAGGGCGCGGTTTTGCGGTGGTCGCGGGAGAAGTGCGTAATCTGGCCCAGCGCAGCGCCCAGGCGGCTCGTGAAATTAAAAGCCTGATTGAAGACTCGGTGGGGAAAGTGGATGTTGGCTCTACGCTGGTCGAAAg                                             >  2:2329841/1‑150 (MQ=255)
                                                                                                           cgTGCGGGTGAGCAAGGGCGCGGTTTTGCGGTGGTCGCGGGAGAAGTGCGTAATCTGGCCCAGCGCAGCGCCCAGGCGGCTCGTGAAATTAAAAGCCTGATTGAAGACTCGGTGGGGAAAGTGGATGTTGGCTCTACGCTGGTCGAAAgc                                            <  1:2229260/150‑1 (MQ=255)
                                                                                                           cgTGCGGGTGAGCAAGGGCGCGGTTTTGCGGTGGTCGCGGGAGAAGTGCGTAATCTGGCCCAGAGCAGCGCCCAGGCGGCTCGTGAAATTAAAAGCCTGATTGAAGACTCGGTGGGGAAAGTGGATGTTGGCTCTACGCTGGTCGAAAgc                                            <  1:4823133/150‑1 (MQ=255)
                                                                                                            gTGCGGGTGAGCATGGGGGGGGTTTTGCGGTGGTCGCGGGAGAAGTGCGTATTCTGGCGCTGCGCAGGTCCGAGGCGGCTCGTGATATTAATAGCCTGATTGACGACTCGGTGGGGAAAGTGGATGTTGGCTCTACGCTGGTCGTTAgcg                                           <  2:4832465/150‑1 (MQ=255)
                                                                                                            gTGCGGGTGAGCAAGGGCGCGGTTTTGCGGTGGTCGCGGGAGAAGTGCGTAATCTGGCCCAGCGCAGCGCCCAGGCGGCTCGTGAAATTAAAAGCCTGATTGAAGACTCGGTGGGGAAAGTGGATGTTGGCTCTACGCTGGTCGAAAgcg                                           <  2:4085416/150‑1 (MQ=255)
                                                                                                             tGCGGGTGAGCAAGGGCGCGGTTTTGCGGTGGTCGCGGGAGAAGTGCGTAATCTGGCCCAGCGCAGCGCCCAGGCGGCTCGTGAAATTAAAAGCCTGATTGAAGACTCGGTGGGGAAAGTGGATGTTGGCTCTACGCTGGTCGAAAgcgc                                          >  1:1402970/1‑150 (MQ=255)
                                                                                                              gCGGGTGAGCAAGGGCGCGGTTTTGCGGTGGTCGCGGGATAAGTGCGTAATCTGGCCCAGCGCAGCGCCCAGGCGGCTCGTGAAATTAAAAGCCTGATTGAAGACTCGGTGGGGAAAGTGGATGTTGGCTCTACGCTGGTCGAAAGCGcc                                         <  2:155343/150‑1 (MQ=255)
                                                                                                                  gTGAGCAAGGGCGCGGTTTTGCGGTGGTCGCGGGAGAAGTGCGTAATCTGGCCCAGCGCAGCGCCCAGGCGGCTCGTGAAATTAAAAGCCTGATTGAAGACTCGGTGGGGAAAGTGGATGTTGGCTCTACGCTGGTCGAAAGCGCCgggg                                     <  2:4305593/150‑1 (MQ=255)
                                                                                                                          gggCGCGGTTTTGCGGTGGTCGCGGGAGAAGTGCGTAATCTGGCCCAGCGCAGCGCCCAGGCGGCTCGTGAAATTAAAAGCCTGATTGAAGACTCGGTGGGGAAAGTGGATGTTGGCTCTACGCTGGTCGAAAGCGCCGGGGAAACAATg                             <  1:4297354/150‑1 (MQ=255)
                                                                                                                           ggCGCGGTTTTGCGGTGGTCGCGGGATAAGTGCGTAATCTGGCCCAGCGCAGCGCCCAGGCGGCTCGTGAAATTAAAAGCCTGATTGAAGACTCGGTGGGGAAAGTGGATGTTGGCTCTACGCTGGTCGAAAGCGCCGGGGAAACAATgg                            >  1:4136832/1‑150 (MQ=255)
                                                                                                                            gcgcGGTTTTGCGGTGGTCGCGGGAGAAGTGCGTAATCTGGCCCAGCGCAGCGCCCAGGCGGCTCGTGAAATTAAAAGCCTGATTGAAGACTCGGTGGGGAAAGTGGATGTTGGCTCTACGCTGGTCGAAAGCGCCGGGGAAACAATgac                           <  2:4733127/150‑3 (MQ=255)
                                                                                                                            gcgcGGTTTTGCGGTGGTCGCGGGAGAAGTGCGTAATCTGGCCCAGCGCAGCGCCCAGGCGGCTCGTGAAATTAAAAGCCTGATTGAAGACTCGGTGGGGAAAGTGGATGTTGGCTCTACGCTGGTCGAAAGCGCCGGGGAAACAATGGc                           <  2:130686/150‑1 (MQ=255)
                                                                                                                                   tttGCGGTGGTCGCGGGAGAAGTGCGTAATCTGGCCCAGCGCAGCGCCCAGGCGGCTCGTGAAATTAAAAGCCTGATTGAAGACTCGGTGGGGAAAGTGGATGTTGGCTCTACGCTGGTCGAAAGCGCCGGGGAAACAATGGCGGAGAtt                    >  1:3942207/1‑150 (MQ=255)
                                                                                                                                   tttGCGGTGGTCGCGGGAGAAGTGCGTAATCTGGCCCAGCGCAGCGCCCAGGCGGCTCGTGAAATTAAAAGCCTGATTGAAGACTCGGTGGGGAAAGTGGATGTTGGCTCTACGCTGGTCGAAAGCGCCGGGGAAACAATGGCGGAGAtt                    >  1:940566/1‑150 (MQ=255)
                                                                                                                                    ttGCGGTGGTCGCGGGAGAAGTGCGTAATCTGGCCCAGCGCAGCGCCCAGGCGGCTCGTGAAATTAAAAGCCTGATTGAAGACTCGGTGGGGAAAGTGGATGTTGGCTCTACGCTGGTCGAAAGCGCCGGGGAAACAATGGCGGAGATTg                   <  1:1145712/150‑1 (MQ=255)
                                                                                                                                     tGCGGTGGTCGCGGGAGAAGTGCGTAATCTGGCCCAGCGCAGCGCCCAGGCGGCTCGTGAAATTAAAAGCCTGATTGAAGACTCGGTGGGGAAAGTGGATGTTGGCTCTACGCTGGTCGAAAGCGCCGGGGAAACAATGGCGGAGATTGt                  <  2:3122240/150‑1 (MQ=255)
                                                                                                                                         gtggtCGCGGGAGAAGTGCGTAATCTGGCCCAGCGCAGCGCCCAGGCGGCTCGTGAAATTAAAAGCCTGATTGAAGACTCGGTGGGGAAAGTGGATGTTGGCTCTACGCTGGTCGAAAGCGCCGGGGAAACAATGGCGGAGATTGTCAgc              >  1:2093798/1‑150 (MQ=255)
                                                                                                                                           ggtCGCGGGAGAAGTGCGTAATCTGGCCCAGCGCAGCGCCCAGGCGGCTCGTGAAATTAAAAGCCTGATTGAAGACTCGGTGGGGAAAGTGGATGTTGGCTCTACGCTGGTCGAAAGCGCCGGGGAAACAATGGCGGAGATTGTCAgcgc            >  1:1495761/1‑150 (MQ=255)
                                                                                                                                             tCGCGGGAGAAGTGCGTAATCTGGCCCAGCGCAGCGCCCAGGCGGCTCGTGAAATTAAAAGCCTGATTGAAGACTCGGTGGGGAAAGTGGATGTTGGCTCTACGCTGGTCGAAAGCGCCGGGGAAACAATGGCGGAGATTGTCAGCGCCg          <  1:3324461/150‑1 (MQ=255)
                                                                                                                                               gcgGGATAAGTGCGTAATCTGGCCCAGCGCAGCGCCCAGGCGGCTCGTGAAATTAAAAGCCTGATTGAAGACTCGGTGGGGAAAGTGGATGTTGGCTCTACGCTGGTCGAAAGCGCCGGGGAAACAATGGCGGAGATTGTCAGCGCCGTg        >  2:2073220/1‑150 (MQ=255)
                                                                                                                                                cgGGAGAAGTGCGTAATCTGGCCCAGCGCAGCGCCCAGGCGGCTCGTGAAATTAAAAGCCTGATTGAAGACTCGGTGGGGAAAGTGGATGTTGGCTCTACGCTGGTCGAAAGCGCCGGGGAAACAATGGCGGAGATTGTCAGCGCCGTGa       >  2:551040/1‑150 (MQ=255)
                                                                                                                                                 gggAGAAGTGCGTAATCTGGCCCAGCGCAGCGCCCAGGCGGCTCGTGAAATTAAAAGCCTGATTGAAGACTCGGTGGGGAAAGTGGATGTTGGCTCTACGCTGGTCGAAAGCGCCGGGGAAACAATGGCGGAGATTGTCAGCGCCGTGAc      >  2:2716548/1‑150 (MQ=255)
                                                                                                                                                 gggAGAAGTGCGTAATCTGGCCCAGCGCAGCGCCCAGGCGGCTCGTGAAATTAAAAGCCTGATTGAAGACTCGGTGGGGAAAGTGGATGTTGGCTCTACGCTGGTCGAAAGCGCCGGGGAAACAATGGCGGAGATTGTCAGCGCCGTGAc      >  2:4968478/1‑150 (MQ=255)
                                                                                                                                                 gggAGAAGTGCGTAATCTGGCCCAGCGCAGCGCCCAGGCGGCTCGTGAAATTAAAAGCCTGATTGAAGACTCGGTGGGGAAAGTGGATGTTGGCTCTACGCTGGTCGAAAGCGCCGGGGAAACAATGGCGGAGATTGTCAGCGCCGTGAc      >  1:4619080/1‑150 (MQ=255)
                                                                                                                                                  ggAGAAGTGCGTAATCTGGCCCAGCGCAGCGCCCAGGCGGCTCGTGAAATTAAAAGCCTGATTGAAGACTCGGTGGGGAAAGTGGATGTTGGCTCTACGCTGGTCGAAAGCGCCGGGGAAACAATGGCGGAGATTGTCAGCGCCGTGAcc     >  2:5077919/1‑150 (MQ=255)
                                                                                                                                                    agaAGTGCGTAATCTGGCCCAGCGCAGCGCCCAGGCGGCTCGTGAAATTAAAAGCCTGATTGAAGACTCGGTGGGGAAAGTGGATGTTGGCTCTACGCTGGTCGAAAGCGCCGGGGAAACAATGGCGGAGATTGTCAGCGCCGTGACCcg   >  1:592610/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GCGATATCTCCACCAGTTCGCAGAAAATCGCCGATATTATCAGCGTAATTGACGGCATTGCCTTCCAGACCAATATTCTGGCTTTGAACGCGGCGGTTGAGGCTGCGCGTGCGGGTGAGCAAGGGCGCGGTTTTGCGGTGGTCGCGGGAGAAGTGCGTAATCTGGCCCAGCGCAGCGCCCAGGCGGCTCGTGAAATTAAAAGCCTGATTGAAGACTCGGTGGGGAAAGTGGATGTTGGCTCTACGCTGGTCGAAAGCGCCGGGGAAACAATGGCGGAGATTGTCAGCGCCGTGACCCGC  >  CP009273/4582528‑4582826

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: