Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 557,592 G→T 100% D814Y (GAC→TAC)  sfmD → putative outer membrane export usher protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273557,5920GT98.4% 887.1 / 3.6 267D814Y (GAC→TAC) sfmDputative outer membrane export usher protein
Reads supporting (aligned to +/- strand):  ref base G (1/4);  new base T (126/136);  total (127/140)
Fisher's exact test for biased strand distribution p-value = 3.73e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
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 276 total aligned reads displayed.

CGACGCAAGGCGCGTTGGTTCGTGCCAATTTTGATACCCGCATAGGCGTGCGGGCGCTCATTACCGTTACCCAGGGCGGAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATGACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACTGGTTCAGTGGGGAGACGGCGCGAACTCACGCTGCATTGCGCACTATGTATTGCCGAAGCAAAGCTTACAGCAAGCCGTCACTGTTATTTCGGCAGTTTG  >  CP009273/557446‑557737
                                                                                                                                                  |                                                                                                                                                 
  aCGCAAGGCGCGTTGGTTCGTGCCAATTTTGATACCCGCATAGGCGTGCGGGCGCTCATTACCGTTACCCAGGGCGGAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATTACgg                                                                                                                                               <  1:4082153/149‑1 (MQ=255)
     cAAGGCGCGTTGGTTCGTGCCAATTTTGATACCCGCATAGGCGTGCGGGCGCTCATTACCGTTACCCAGGGCGGAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATTACGGGCa                                                                                                                                            >  1:839272/1‑149 (MQ=255)
      aaGGCTCCTTGGTTCGAGCCAATTTTGATACCCGCATAGGCTTGCGGGCGCTCATTACCGTTAACCATTGAGGAAAACCTGTGCCGAATGTATCAAAGGTAAGGGAAAACAGTACAGGAATAACCAGTAGGGTGGGTGATTAAGGGCaa                                                                                                                                           >  2:1291059/1‑149 (MQ=255)
      aaGGCGCGTTGGTTCGTGCCAATTTTGATACCCGCATAGGCGTGCGGGCGCTCATTACCGTTACCCAGGGCGGAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATTACGGGCaa                                                                                                                                           <  2:414250/149‑1 (MQ=255)
       aGGCGCGTTGGTTCGTGCCAATTTTGATACCCGCATAGGCGTGCGGGCGCTCATTACCGTTACCCAGGGCGGAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATTACGGGCAAg                                                                                                                                          >  2:1743616/1‑149 (MQ=255)
       aGGCGCGTTGGTTCGTGCCAATTTTGATACCCGCATAGGCGTGCGGGCGCTCATTACCGTTACCCAGGGCGGAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATTACGGGCAAg                                                                                                                                          <  1:2753623/149‑1 (MQ=255)
       aGGCGCGTTGGTTCGTGCCAATTTTGATACCCGCATAGGCGTGCGGGCGCTCATTACCGTTACCCAGGGCGGAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATTACGGGCAAg                                                                                                                                          <  1:4624367/149‑1 (MQ=255)
          cgcgTTGGTTCGTGCCAATTTTGATACCCGCATAGGCGTGCGGGCGCTCATTACCGTTACCCAGGGCGGAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATTACGGGCAAGttt                                                                                                                                       <  1:2323262/149‑1 (MQ=255)
           gcgTTGGTTCGTGCCAATTTTGATACCCGCATAGGCGTGCGGGCGCTCATTACCGTTACCCAGGGCGGAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATTACGGGCAAGttta                                                                                                                                      <  2:65207/149‑1 (MQ=255)
             gTTGGTTCGTGCCAATTTTGATACCCGCATAGGCGTGCGGGCGCTCATTACCGTTACCCAGGGCGGAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATTACGGGCAAGtttatt                                                                                                                                    <  1:415099/149‑1 (MQ=255)
              ttGGTTCGTGCCAATTTTGATACCCGCATAGGCGTGCGGGCGCTCATTACCGTTACCCAGGGCGGAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATTACGGGCAAGtttattt                                                                                                                                   >  2:3029780/1‑149 (MQ=255)
              ttGGTTCGTGCCAATTTTGATACCCGCATAGGCGTGCGGGCGCTCATTACCGTTACCCAGGGCGGAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATTACGGGCAAGtttattt                                                                                                                                   <  1:1311420/149‑1 (MQ=255)
                ggTTCGTGCCAATTTTGATACCCGCATAGGCGTGCGGGCGCTCATTACCGTTACCCAGGGCGGAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATTACGGGCAAGTTTATTTaa                                                                                                                                 >  2:2836828/1‑149 (MQ=255)
                ggTTCGTGCCAATTTTGATACCCGCATAGGCGTGCGGGCGCTCATTACCGTTACCCAGGGCGGAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATTACGGGCAAGTTTATTTaa                                                                                                                                 >  1:736971/1‑149 (MQ=255)
                ggTTCGTGCCAATTTTGATACCCGCATAGGCGTGCGGGCGCTCATTACCGTTACCCAGGGCGGAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGAGTACGGGCAAGTTTATTTaa                                                                                                                                 >  2:1809470/1‑149 (MQ=255)
                   tCGTGCCAATTTTGATACCCGCATAGGCGTGCGGGCGCTCATTACCGTTACCCAGGGCGGAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATTACGGGCAAGTTTATTTAAgtg                                                                                                                              >  1:1570345/1‑149 (MQ=255)
                       gCCAATTTTGATACCCGCATAGGCGTGCGGGCGCTCATTACCGTTACCCAGGGCGGAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATTACGGGCAAGTTTATTTAAGTGGTgc                                                                                                                          >  2:788051/1‑149 (MQ=255)
                       gCCAATTTTGATACCCGCATAGGCGTGCGGGCGCTCATTACCGTTACCCAGGGCGGAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATTACGGGCAAGTTTATTTAAGTGGTgc                                                                                                                          <  1:2837295/149‑1 (MQ=255)
                       gCCAATTTTGATACCCGCATAGGCGTGCGGGCGCTCATTACCGTTACCCAGGGCGGAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATTACGGGCAAGTTTATTTAAGTGGTgc                                                                                                                          >  2:1499912/1‑149 (MQ=255)
                         cAATTTTGATACCCGCATAGGCGTGCGGGCGCTCATTACCGTTACCCAGGGCGGAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATTACGGGCAAGTTTATTTAAGTGGTgcgc                                                                                                                        >  1:2520352/1‑149 (MQ=255)
                          aaTTTTGATACCCGCATAGGCGTGCGGGCGCTCATTACCGTTACCCAGGGCGGAAAACCGGTGCCGTTTGGATCACTGGTCCGGGAAAACAGTACCGGAATAACCAGTCTGGTGGGTGATTACGGGCAAGTTTATTTAAGTGGTGCGcc                                                                                                                       >  2:3798221/1‑149 (MQ=255)
                          aaTTTTGATACCCGCATAGGCGTGCGGGCGCTCATTACCGTTACCCAGGGCGGAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATTACGGGCAAGTTTATTTAAGTGGTGCGcc                                                                                                                       >  1:544124/1‑149 (MQ=255)
                          aaTTTTGATACCCGCATAGGCGTGCGGGCGCTCATTACCGTTACCCAGGGCGGAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATTACGGGCAAGTTTATTTAAGTGGTGCGcc                                                                                                                       <  2:3772768/149‑1 (MQ=255)
                            ttttGATACCCGCATAGGCGTGCGGGCGCTCATTACCGTTACCCAGGGCGGAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATTACGGGCAAGTTTATTTAAGTGGTGCGCCAt                                                                                                                     <  1:4170642/149‑1 (MQ=255)
                             tttGATACCCGCATAGGCGTGCGGGCGCTCATTACCGTTACCCAGGGCGGAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATTACGGGCAAGTTTATTTAAGTGGTGCGCCAtt                                                                                                                    <  2:541750/149‑1 (MQ=255)
                               tGATACCCGCATAGGCGTGCGGGCGCTCATTACCGTTACCCAGGGCGGAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATTACGGGCAAGTTTATTTAAGTGGTGCGCCATTGt                                                                                                                  >  1:4474472/1‑149 (MQ=255)
                               tGATACCCGCACAGGCGTGCGGGCGCTCATTACCGTTACCCAGGGCGGAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATTACGGGCAAGTTTATTTAAGTGGTGCGCCATTGt                                                                                                                  >  2:2238245/1‑149 (MQ=255)
                                    cccGCATAGGCGTGCGGGCGCTCATTACCGTTACCCAGGGCGGAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATTACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGt                                                                                                             >  1:2205303/1‑149 (MQ=255)
                                          tAGGCGTGCGGGCGCTCATTACCGTTACCCAGGGCGGAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATTACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTa                                                                                                       >  1:3722430/1‑149 (MQ=255)
                                           aGGCGTGCGGGCGCTCATTACCGTTACCCAGGGCGGAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATTACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTAc                                                                                                      >  1:1857778/1‑149 (MQ=255)
                                           aGGCGTGCGGGCGCTCATTACCGTTACCCAGGGCGGAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATTACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTAc                                                                                                      >  1:1187291/1‑149 (MQ=255)
                                             gCGTGCGGGCGCTCATTACCGTTACCCAGGGCGGAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATTACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACTg                                                                                                    <  1:700108/149‑1 (MQ=255)
                                                 gCGGGCGCTCATTACCGTTACCCAGGGCGGAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATTACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACTGGTTc                                                                                                >  2:1316367/1‑149 (MQ=255)
                                                 gCGGGCGCTCATTACCGTTACCCAGGGCGGAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATTACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACTGGTTc                                                                                                <  2:3157505/149‑1 (MQ=255)
                                                  cGGGCGCTCATTACCGTTACCCAGGGCGGAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATTACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACTGGTTCa                                                                                               <  2:3727337/149‑1 (MQ=255)
                                                   gggCGCTCATTACCGTTACCCAGGGCGGAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATTACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACTGGTTCAg                                                                                              <  2:2765033/149‑1 (MQ=255)
                                                    ggCGCTCATTACCGTTACCCAGGGCGGAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATTACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACTGGTTCAGt                                                                                             >  2:3099228/1‑149 (MQ=255)
                                                    ggCGCTCATTACCGTTACCCAGGGCGGAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATTACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACTGGTTCAGt                                                                                             >  2:707034/1‑149 (MQ=255)
                                                     gcgcTCATTACCGTTACCCAGGGCGGAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATTACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACTGGTTCAGTg                                                                                            <  2:2173657/149‑1 (MQ=255)
                                                        cTCATTACCGTTACCCAGGGCGGAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATTACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACTGGTTCAGTgggg                                                                                         >  2:3523599/1‑149 (MQ=255)
                                                          cATTACCGTTACCCAGGGCGGAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATTACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACTGGTTCAGTGGGgag                                                                                       >  1:4740825/1‑149 (MQ=255)
                                                             tACCGTTACCCAGGGCGGAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATTACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACTGGTTCAGTGGGGAGACg                                                                                    >  2:2183114/1‑149 (MQ=255)
                                                               ccGTTACCCAGGGCGGAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATTACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACTGGTTCAGTGGTGAGACGgc                                                                                  >  2:2549842/1‑149 (MQ=255)
                                                               ccGTTACCCAGGGCGGAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATTACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACTGGTTCAGTGGGGAGACGgc                                                                                  <  1:3965548/149‑1 (MQ=255)
                                                                cGTTACCCAGGGCGGAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATTACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACTGGTTCAGTGGGGAGACGgcg                                                                                 >  1:337773/1‑149 (MQ=255)
                                                                 gTTACCCAGGGCGGAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATTACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACTGGTTCAGTGGGGAGACGgcgc                                                                                >  2:2015549/1‑149 (MQ=255)
                                                                 gTTACCCAGGGCGGAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATTACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACTGGTTCAGTGGGGAGACGgcgc                                                                                >  2:1458129/1‑149 (MQ=255)
                                                                    aCCCAGGGCGGAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATTACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACTGGTTCAGTGGGGAGACGGCGCGaa                                                                             <  1:3700841/149‑1 (MQ=255)
                                                                         gggCGGAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATTACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACTGGTTCAGTGGGGAGACGGCGCGAACTCAc                                                                        >  2:4336714/1‑149 (MQ=255)
                                                                           gCGGAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATTACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACTGGTTCAGTGGGGAGACGGCGCGAACTCACGc                                                                      >  1:802502/1‑149 (MQ=255)
                                                                              gAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATTACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACTGGTTCAGTGGGGAGACGGCGCGAACTCACGCTGc                                                                   >  1:2654273/1‑149 (MQ=255)
                                                                              gAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATTACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACTGGTTCAGTGGGGAGACGGCGCGAACTCACGCTGc                                                                   <  1:1148450/149‑1 (MQ=255)
                                                                                   ccGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATTACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACTGGTTCAGTGGGGAGACGGCGCGAACTCACGCTGCATTgc                                                              <  1:30634/149‑1 (MQ=255)
                                                                                        gCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATTACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACTGGTTCAGTGGGGAGACGGCGCGAACTCACGCTGCATTGCGCACt                                                         >  2:3891225/1‑149 (MQ=255)
                                                                                        gCCGTGTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATTACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACTGGTTCAGTGGGGAGACGGCGCGAACTCACGCTGCATTGCGCACt                                                         <  1:3307714/149‑1 (MQ=255)
                                                                                          cGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATTACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACTGGTTCAGTGGGGAGACGGCGCGAACTCACGCTGCATTGCGCACTAt                                                       <  2:459896/149‑1 (MQ=255)
                                                                                           gTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATTACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACTGGTTCAGTGGGGAGACGGCGCGAACTCACGCTGCATTGCGCACTATg                                                      >  2:3226863/1‑149 (MQ=255)
                                                                                             ttGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATTACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACTGGTTCAGTGGGGAGACGGCGCGAACTCACGCTGCATTGCGCACTATGTa                                                    <  1:1419784/149‑1 (MQ=255)
                                                                                                gATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATTACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACTGGTTCAGTGGGGAGACGGCGCGAACTCACGCTGCATTGCGCACTATGTATTg                                                 >  1:995834/1‑149 (MQ=255)
                                                                                                gATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATTACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACTGGTTCAGTGGGGAGACGGCGCGAACTCACGCTGCATTGCGCACTATGTATTg                                                 >  1:793992/1‑149 (MQ=255)
                                                                                                gATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATTACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACTGGTTCAGTGGGGAGACGGCGCGAACTCACGCTGCATTGCGCACTATGTATTg                                                 >  1:3549147/1‑149 (MQ=255)
                                                                                                 aTCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATTACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACTGGTTCAGTGGGGAGACGGCGCGAACTCACGCTGCATTGCGCACTATGTATTGc                                                <  2:544124/149‑1 (MQ=255)
                                                                                                  tCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATTACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACTGGTTCAGTGGGGAGACGGCGCGAACTCACGCTGCATTGCGCACTATGTATTGcc                                               >  1:150919/1‑149 (MQ=255)
                                                                                                   cACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATTACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACTGGTTCAGTGGGGAGACGGCGCGAACTCACGCTGCATTGCGCACTATGTATTGCCg                                              >  1:1258637/1‑149 (MQ=255)
                                                                                                   cACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATTACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACTGGTTCAGTGGGGAGACGGCGCGAACTCACGCTGCATTGCGCACTATGTATTGCCg                                              >  1:500149/1‑149 (MQ=255)
                                                                                                     cTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATTACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACTGGTTCAGTGGGGAGACGGCGCGAACTCACGCTGCATTGCGCACTATGTATTGCCGaa                                            <  2:1622542/149‑1 (MQ=255)
                                                                                                     cTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATTACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACTGGTTCAGTGGGGAGACGGCGCGAACTCACGCTGCATTGCGCACTATGTATTGCCGaa                                            <  2:4481954/149‑1 (MQ=255)
                                                                                                     cTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATTACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACTGGTTCAGTGGGGAGACGGCGCGAACTCACGCTGCATTGCGCACTATGTATTGCCGaa                                            >  1:1086828/1‑149 (MQ=255)
                                                                                                     cTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATTACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACTGGTTCAGTGGGGAGACGGCGCGAACTCACGCTGCATTGCGCACTATGTATTGCCGaa                                            <  1:2482083/149‑1 (MQ=255)
                                                                                                     cTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATTACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACTGGTTCAGTGGGGAGACGGCGCGAACTCACGCTGCATTGCGCACTATGTATTGCCGaa                                            >  1:3475295/1‑149 (MQ=255)
                                                                                                      tGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATTACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACTGGTTCAGTGGGGAGACGGCGCGAACTCACGCTGCATTGCGCACTATGTATTGCCGAAg                                           >  1:180351/1‑149 (MQ=255)
                                                                                                       tgTACGGGAAAACAGTACCGCAATAACCAGTATGGTGGGTGATGACGGTCAAGTTTTTTTAAGTGGTGCGCCATGGTCTGGTGATTTAGTGGTTCAGTGGGGAGACGTCGCTAACTCACGCTGCATTGCGAACTATGTATTGCCGAAGc                                          <  2:3343633/148‑1 (MQ=255)
                                                                                                       ggTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATTACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACTGGTTCAGTGGGGAGACGGCGCGAACTCACGCTGCATTGCGCACTATGTATTGCCGAAGc                                          <  2:955515/149‑1 (MQ=255)
                                                                                                        gTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATTACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACTGGTTCAGTGGGGAGACGGCGCGAACTCACGCTGCATTGCGCACTATGTATTGCCGAAGCa                                         >  2:4176475/1‑149 (MQ=255)
                                                                                                           cGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATTACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACTGGTTCAGTGGGGAGACGGCGCGAACTCACGCTGCATTGCGCACTATGTATTGCCGAAGCAAAg                                      >  2:1950508/1‑149 (MQ=255)
                                                                                                            gtgAAAACAGTACCGGAATAACCAGTATGGTGGGTGATTACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACTGGTTCAGTGGGGAGACGGCGCGAACTCACGCTGCATTGCGCACTATGTATTGCCGAAGCAAAGc                                     <  2:1673820/147‑1 (MQ=255)
                                                                                                                  aCAGTACCGGAATAACCAGTATGGTGGGTGATTACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACTGGTTCAGTGGGGAGACGGCGCGAACTCACGCTGCATTGCGCACTATGTATTGCCGAAGCAAAGCTTACAg                               <  1:4089407/149‑1 (MQ=255)
                                                                                                                  aCAGTACCGGAATAACCAGTATGGTGGGTGATTACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACTGGTTCAGTGGGGAGACGGCGCGAACTCACGCTGCATTGCGCACTATGTATTGCCGAAGCAAAGCTTACAg                               >  2:4244534/1‑149 (MQ=255)
                                                                                                                   cAGTACCGGAATAACCAGTATGGTGGGTGATTACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACTGGTTCAGTGGGGAGACGGCGCGAACTCACGCTGCATTGCGCACTATGTATTGCCGAAGCAAAGCTTACAGc                              <  1:2039304/149‑1 (MQ=255)
                                                                                                                    aGTACCGGAATAACCAGTATGGTGGGTGATTACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACTGGTTCAGTGGGGAGACGGCGCGAACTCACGCTGCATTGCGCACTATGTATTGCCGAAGCAAAGCTTACAGCa                             >  2:2717926/1‑149 (MQ=255)
                                                                                                                     gTACCGGAATAACCAGTATGGTGGGTGATTACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACTGGTTCAGTGGGGAGACGGCGCGAACTCACGCTGCATTGCGCACTATGTATTGCCGAAGCAAAGCTTACAGCaa                            >  1:2385987/1‑149 (MQ=255)
                                                                                                                      tACCGGAATAACCAGTATGGTGGGTGATTACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACTGGTTCAGTGGGGAGACGGCGCGAACTCACGCTGCATTGCGCACTATGTATTGCCGAAGCAAAGCTTACAGCAAg                           >  2:4539735/1‑149 (MQ=255)
                                                                                                                       aCCGGAATAACCAGTATGGTGGGTGATTACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACTGGTTCAGTGGGGAGACGGCGCGAACTCACGCTGCATTGCGCACTATGTATTGCCGAAGCAAAGCTTACAGCAAGc                          <  1:2852919/149‑1 (MQ=255)
                                                                                                                       aCCGGAATAACCAGTATGGTGGGTGATTACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACTGGTTCAGTGGGGAGACGGCGCGAACTCACGCTGCATTGCGCACTATGTATTGCCGAAGCAAAGCTTACAGCAAGc                          <  2:3603500/149‑1 (MQ=255)
                                                                                                                          ggAATAACCAGTATGGTGGGTGATTACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACTGGTTCAGTGGGGAGACGGCGCGAACTCACGCTGCATTGCGCACTATGTATTGCCGAAGCAAAGCTTACAGCAAGCCGt                       <  2:2357566/149‑1 (MQ=255)
                                                                                                                            aaTAACCAGTATGGTGGGTGATTACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACTGGTTCAGTGGGGAGACGGCGCGAACTCACGCTGCATTGCGCACTATGTATTGCCGAAGCAAAGCTTACAGCAAGCCGTCa                     >  1:4785338/1‑149 (MQ=255)
                                                                                                                            aaTAACCAGTATGGTGGGTGATTACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACTGGTTCAGTGGGGAGACGGCGCGAACTCACGCTGCATTGCGCACTATGTATTGCCGAAGCAAAGCTTACAGCAAGCCGTCa                     >  2:1556855/1‑149 (MQ=255)
                                                                                                                            aaTAACCAGTATGGTGGGTGATTACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACTGGTTCAGTGGGGAGACGGCGCGAACTCACGCTGCATTGCGCACTATGTATTGCCGAAGCAAAGCTTACAGCAAGCCGTCa                     >  1:4013257/1‑149 (MQ=255)
                                                                                                                            aaTAACCAGTATGGTGGGTGATTACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACTGGTTCAGTGGGGAGACGGCGCGAACTCACGCTGCATTGCGCACTATGTATTGCCGAAGCAAAGCTTACAGCAAGCCGTCa                     >  2:1918240/1‑149 (MQ=255)
                                                                                                                              tAACCAGTATGGTGGGTGATTACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACTGGTTCAGTGGGGAGACGGCGCGAACTCACGCTGCATTGCGCACTATGTATTGCCGAAGCAAAGCTTACAGCAAGCCGTCACt                   <  1:3443968/149‑1 (MQ=255)
                                                                                                                                aCCAGTATGGTGGGTGATTACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACTGGTTCAGTGGGGAGACGGCGCGAACTCACGCTGCATTGCGCACTATGTATTGCCGAAGCAAAGCTTACAGCAAGCCGTCACTGt                 >  2:2578144/1‑149 (MQ=255)
                                                                                                                                 ccAGTATGGTGGGTGATTACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACTGGTTCAGTGGGGAGACGGCGCGAACTCACGCTGCATTGCGCACTATGTATTGCCGAAGCAAAGCTTACAGCAAGCCGTCACTGtt                >  2:1072381/1‑149 (MQ=255)
                                                                                                                                    gTATGGTGGGTGATTACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACTGGTTCAGTGGGGAGACGGCGCGAACTCACGCTGCATTGCGCACTATGTATTGCCGAAGCAAAGCTTACAGCAAGCCGTCACTGTTAtt             <  2:1297629/149‑1 (MQ=255)
                                                                                                                                        ggtgggtgATTACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACTGGTTCAGTGGGGAGACGGCGCGAACTCACGCTGCATTGCGCACTATGTATTGCCGAAGCAAAGCTTACAGCAAGCCGTCACTGTTATTTCgg         >  1:2038542/1‑149 (MQ=255)
                                                                                                                                         gtgggtgATTACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACTGGTTCAGTGGGGAGACGGCGCGAACTCACGCTGCATTGCGCACTATGTATTGCCGAAGCAAAGCTTACAGCAAGCCGTCACTGTTATTTCGGc        <  1:4318036/149‑1 (MQ=255)
                                                                                                                                         gtgggtgATTACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACTGGTTCAGTGGGGAGACGGCGCGAACTCACGCTGCATTGCGCACTATGTATTGCCGAAGCAAAGCTTACAGCAAGCCGTCACTGTTATTTCGGc        <  1:4317984/149‑1 (MQ=255)
                                                                                                                                           gggtgATTACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACTGGTTCAGTGGGGAGACGGCGCGAACTCACGCTGCATTGCGCACTATGTATTGCCGAAGCAAAGCTTACAGCAAGCCGTCACTGTTATTTCGGCAg      <  2:2391627/149‑1 (MQ=255)
                                                                                                                                              tgattaCGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACTGGTTCAGTGGGGAGACGGCGCGAACTCACGCTGCATTGCGCACTATGTATTGCCGAAGCAAAGCTTACAGCAAGCCGTCACTGTTATTTCGGCAGttt   >  1:1353330/1‑149 (MQ=255)
                                                                                                                                              tgattaCGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACTGGTTCAGTGGGGAGACGGCGCGAACTCACGCTGCATTGCGCACTATGTATTGCCGAAGCAAAGCTTACAGCAAGCCGTCACTGTTATTTCGGCAGttt   >  1:3742838/1‑149 (MQ=255)
                                                                                                                                               gattaCGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACTGGTTCAGTGGGGAGACGGCGCGAACTCACGCTGCATTGCGCACTATGTATTGCCGAAGCAAAGCTTACAGCAAGCCGTCACTGTTATTTCGGCAGTTTg  >  2:4574535/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                 
CGACGCAAGGCGCGTTGGTTCGTGCCAATTTTGATACCCGCATAGGCGTGCGGGCGCTCATTACCGTTACCCAGGGCGGAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATGACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACTGGTTCAGTGGGGAGACGGCGCGAACTCACGCTGCATTGCGCACTATGTATTGCCGAAGCAAAGCTTACAGCAAGCCGTCACTGTTATTTCGGCAGTTTG  >  CP009273/557446‑557737

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