Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,133,311 T→A 18.0% L162Q (CTG→CAG)  flgJ → flagellar rod assembly protein and murein hydrolase; flagellum‑specific muramidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,133,3110TA18.0% 433.4 / 81.7 172L162Q (CTG→CAG) flgJflagellar rod assembly protein and murein hydrolase; flagellum‑specific muramidase
Reads supporting (aligned to +/- strand):  ref base T (69/72);  new base A (17/14);  total (86/86)
Fisher's exact test for biased strand distribution p-value = 6.92e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 173 total aligned reads displayed.

GATGAAATTCCCGCTCGAAACTGTGGTGCGTTATCAAAATCAGGCGCTTTCGCAGCTGGTGCAAAAGGCCGTGCCACGTAACTACGATGATTCGCTGCCGGGTGACAGTAAAGCATTCCTCGCGCAACTCTCGCTGCCCGCCCAACTGGCAAGCCAGCAAAGCGGTGTGCCACATCATTTGATCCTCGCTCAGGCGGCACTGGAATCTGGTTGGGGGCAACGGCAAATCCGCCGCGAAAACGGCGAGCCGAGCTATAACCTGTTTGGTGTCAAAGCCTCTGGCAACTGGAAAGG  >  CP009273/1133165‑1133458
                                                                                                                                                  |                                                                                                                                                   
gATGAAATTCCCGCTCGAAACTGTGGTGCGTTATCAAAATCAGGCGCTTTCGCAGCTGGTGCAAAAGGCCGTGCCACGTAACTACGATGATTCGCTGCCGGGTGACAGTAAAGCATTCCTCGCGCAACTCTCGCTGCCCGCCCAACTg                                                                                                                                                    >  1:2543422/1‑148 (MQ=255)
gATGAAATTCCCGCTCGAAACTGTGGTGCGTTATCAAAATCAGGCGCTTTCGCAGCTGGTGCAAAAGGCCGTGCCACGTAACTACGATGATTCGCTGCCGGGTGACAGTAAAGCATTCCTCGCGCAACTCTCGCTGCCCGCCCAACTg                                                                                                                                                    >  1:2367123/1‑148 (MQ=255)
 aTGAAATTCCCGCTCGAAACTGTGGTGCGTTATCAAAATCAGGCGCTTTCGCAGCTGGTGCAAAAGGCCGTGCCACGTAACTACGATGATTCGCTGCCGGGTGACAGTAAAGCATTCCTCGCGCAACTCTCGCTGCCCGCCCAACTgg                                                                                                                                                   <  1:2956925/148‑1 (MQ=255)
      aTTCCCGCTCGAAACTGTGGTGCGTTATCAAAATCAGGCGCTTTCGCAGCTGGTGCAAAAGGCCGTGCCACGTAACTACGATGATTCGCTGCCGGGTGACAGTAAAGCATTCCTCGCGCAACTCTCGCTGCCCGCCCAACAGGCAAGc                                                                                                                                              >  1:1064196/1‑148 (MQ=255)
      aTTCCCGCTCGAAACTGTGGTGCGTTATCAAAATCAGGCGCTTTCGCAGCTGGTGCAAAAGGCCGTGCCACGTAACTACGATGATTCGCTGCCGGGTGACAGTAAAGCATTCCTCGCGCAACTCTCGCTGCCCGCCCAACAGGCAAGc                                                                                                                                              >  1:1064188/1‑148 (MQ=255)
          ccGCTCGAAACTGTGGTGCGTTATCAAAATCAGGCGCTTTCGCAGCTGGTGCAAAAGGCCGTGCCACGTAACTACGATGATTCGCTGCCGGGTGACAGTAAAGCATTCCTCGCGCAACTCTCGCTGCCCGCCCAACTGGCAAGCCAGc                                                                                                                                          >  2:2510261/1‑148 (MQ=255)
          ccGCTCGAAACTGTGGTGCGTTATCAAAATCAGGCGCTTTCGCAGCTGGTGCAAAAGGCCGTGCCACGTAACTACGATGATTCGCTGCCGGGTGACAGTAAAGCATTCCTCGCGCAACTCTCGCTGCCCGCCCAACTGGCAAGCCAGc                                                                                                                                          >  1:1620717/1‑148 (MQ=255)
             cTCGAAACTGTGGTGCGTTATCAAAATCAGGCGCTTTCGCAGCTGGTGCAAAAGGCCGTGCCACGTAACTACGATGATTCGCTGCCGGGTGACAGTAAAGCATTCCTCGCGCAACTCTCGCTGCCCGCCCAACTGGCACGCCAGCaca                                                                                                                                       >  1:2136478/1‑146 (MQ=255)
                  aaCTGTGGTGCGTTATCAAAATCAGGCGCTTTCGCAGCTGGTGCAAAAGGCCGTGCCACGTAACTACGATGATTCGCTGCCGGGTGACAGTAAAGCATTCCTCGCGCAACTCTCGCTGCCCGCCCAACTGGCAAGCCAGCAAAGCGgt                                                                                                                                  >  1:2855755/1‑148 (MQ=255)
                    cTGTGGTGCGTTATCAAAATCAGGCGCTTTCGCAGCTGGTGCAAAAGGCCGTGCCACGTAACTACGATGATTCGCTGCCGGGTGACAGTAAAGCATTCCTCGCGCAACTCTCGCTGCCCGCCCAACTGGCAAGCCAGCAAAGCGgtgt                                                                                                                                >  2:180011/1‑148 (MQ=255)
                    cTGTGGTGCGTTATCAAAATCAGGCGCTTTCGCAGCTGGTGCAAAAGGCCGTGCCACGTAACTACGATGATTCGCTGCCGGGTGACAGTAAAGCATTCCTCGCGCAACTCTCGCTGCCCGCCCAACTGGCAAGCCAGCAAAGCGgtgt                                                                                                                                >  1:388449/1‑148 (MQ=255)
                      gtggtgCGTTATCAAAATCAGGCGCTTTCGCAGCTGGTGCAAAAGGCCGTGCCACGTAACTACGATGATTCGCTGCCGGGTGACAGTAAAGCATTCCTCGCGCAACTCTCGCTGCCCGCCCAACTGGCAAGCCAGCAAAGCGGTGTGc                                                                                                                              >  1:606535/1‑148 (MQ=255)
                      gtggtgCGTTATCAAAATCAGGCGCTTTCGCAGCTGGTGCAAAAGGCCGTGCCACGTAACTACGATGATTCGCTGCCGGGTGACAGTAAAGCATTCCTCGCGCAACTCTCGCTGCCCGCCCAACTGGCAAGCCAGCAAAGCGGTGTGc                                                                                                                              >  1:3494704/1‑148 (MQ=255)
                           gCGTTATCAAAATCAGGCGCTTTATCAGATGGTGCAAAAGGCCGTGCCACGTAACTACGATGATTCGCTGCCGGGTGACAGTAAAGCATTCCTCGCGCAACTCTCGCTGCCCGCCCAACTGGCAAGCCAGCAAAGCGGTGTGCCAcat                                                                                                                         >  1:545563/1‑148 (MQ=255)
                            cGTTATCAAAATCAGGCGCTTTCGCAGCTGGTGCAAAAGGCCGTGCCACGTAACTACGATGATTCGCTGCCGGGTGACAGTAAAGCATTCCTCGCGCAACTCTCGCTGCCCGCCCAACTGGCAAGCCAGCAAAGCGGTGTGCCAcatc                                                                                                                        >  1:1694284/1‑148 (MQ=255)
                             gTTATCAAAATCAGGCGCTTTCGCAGCTGGTGCAAAAGGCCGTGCCACGTAACTACGATGATTCGCTGCCGGGTGACAGTAAAGCATTCCTCGCGCAACTCTCGCTGCCCGCCCAACTGGCAAGCCAGCAAAGCGGTGTGCCAcatca                                                                                                                       <  1:1833915/148‑1 (MQ=255)
                              ttATCAAAATCAGGCGCTTTCGCAGCTGGTGCAAAAGGCCGTGCCACGTAACTACGATGATTCGCTGCCGGGTGACAGTAAAGCATTCCTCGCGCAACTCTCGCTGCCCGCCCAACTGGCAAGCCAGCAAAGCGGTGTGCCAcatcat                                                                                                                      <  2:2394445/148‑1 (MQ=255)
                              ttATCAAAATCAGGCGCTTTCGCAGCTGGTGCAAAAGGCCGTGCCACGTAACTACGATGATTCGCTGCCGGGTGACAGTAAAGCATTCCTCGCGCAACTCTCGCTGCCCGCCCAACAGGCAAGCCAGCAAAGCGGTGTGCCAcatcat                                                                                                                      <  1:2305798/148‑1 (MQ=255)
                                 tCAAAATCAGGCGCTTTCGCAGCTGGTGCAAAAGGCCGTGCCACGTAACTACGATGATTCGCTGCCGGGTGACAGTAAAGCATTCCTCGCGCAACTCTCGCTGCCCGCCCAACTGGCAAGCCAGCAAAGCGGTGTGCCACATCATTTg                                                                                                                   >  1:662878/1‑148 (MQ=255)
                                 tCAAAATCAGGCGCTTTCGCAGCTGGTGCAAAAGGCCGTGCCACGTAACTACGATGATTCGCTGCCGGGTGACAGTAAAGCATTCCTCGCGCAACTCTCGCTGCCCGCCCAACTGGCAAGCCAGCAAAGCGGTGTGCCACATCATTTg                                                                                                                   >  2:3351301/1‑148 (MQ=255)
                                  cAAAATCAGGCTCTTTCGCAGCTTGTTCAAAAGGCCGTTCCACGTAACTCAGATGATTCGCTGCCGGGTGACAGTAAAGCATTACTAGCGCAACTCTCGCTGCCCGACCAACTGGCAAGCCCGCAAAGCGGTGTGCCACATCATTTGa                                                                                                                  <  1:1975525/148‑1 (MQ=255)
                                   aaaaTCAGGCGCTTTCGCAGCTGGTGCAAAAGGCCGTGCCACGTAACTACGATGATTCGCTGCCGGGTGACAGTAAAGCATTCCTCGCGCAACTCTCGCTGCCCGCCCAACTGGCAAGCCAGCAAAGCGGTGTGCCACATCATTTGAt                                                                                                                 <  1:718671/148‑1 (MQ=255)
                                    aaaTCAGGCGCTTTCGCAGCTGGTGCAAAAGGCCGTGCCACGTAACTACGATGATTCGCTGCCGGGTGACAGTAAAGCATTCCTCGCGCAACTCTCGCTGCCCGCCCAACTGGCAAGCCAGCAAAGCGGTGTGCCACATCATTTGATc                                                                                                                <  1:3756203/148‑1 (MQ=255)
                                    aaaTCAGGCGCTTTCGCAGCTGGTGCAAAAGGCCGTGCCACGTAACTACGATGATTCGCTGCCGGGTGACAGTAAAGCATTCCTCGCGCAACTCTCGCTGCCCGCCCAACTGGCAAGCCAGCAAAGCGGTGTGCCACATCATTTGATc                                                                                                                >  2:205035/1‑148 (MQ=255)
                                    aaaTCAGGCGCTTTCGCAGCTGGTGCAAAAGGCCGTGCCACGTAACTACGATGATTCGCTGCCGGGTGACAGTAAAGCATTCCTCGCGCAACTCTCGCTGCCCGCCCAACTGGCAAGCCAGCAAAGCGGTGTGCCACATCATTTGATc                                                                                                                >  2:205068/1‑148 (MQ=255)
                                    aaaTCAGGCGCTTTCGCAGCTGGTGCAAAAGGCCGTGCCACGTAACTACGATGATTCGCTGCCGGGTGACAGTAAAGCATTCCTCGCGCAACTCTCGCTGCCCGCCCAACTGGCAAGCCAGCAAAGCGGTGTGCCACATCATTTGATc                                                                                                                <  1:398569/148‑1 (MQ=255)
                                     aaTCAGGCGCTTTCGCAGCTGGTGCAAAAGGCCGTGCCACGTAACTACGATGATTCGCTGCCGGGTGACAGTAAAGCATTCCTCGCGCAACTCTCGCTGCCCGCCCAACTGGCAAGCCAGCAAAGCGGTGTGCCACATCATTTGATcc                                                                                                               >  1:2083890/1‑148 (MQ=255)
                                     aaTCAGGCGCTTTCGCAGCTGGTGCAAAAGGCCGTGCCACGTAACTACGATGATTCGCTGCCGGGTGACAGTAAAGCATTCCTCGCGCAACTCTCGCTGCCCGCCCAACTGGCAAGCCAGCAAAGCGGTGTGCCACATCATTTGATcc                                                                                                               <  2:1907332/148‑1 (MQ=255)
                                      aTCAGGCGCTTTCGCAGCTGGTGCAAAAGGCCGTGCCACGTAACTACGATGATTCGCTGCCGGGTGACAGTAAAGCATTCCTCGCGCAACTCTCGCTGCCCGCCCAACTGGCAAGCCAGCAAAGCGGTGTGCCACATCATTTGATCCt                                                                                                              <  1:138380/148‑1 (MQ=255)
                                           gcgcTTTCGCAGCTGGTGCAAAAGGCCGTGCCACGTAACTACGATGATTCGCTGCCGGGTGACAGTAAAGCATTCCTCGCGCAACTCTCGCTGCCCGCCCAACTGGCAAGCCAGCAAAGCGGTGTGCCACATCATTTGATCCTCGCTc                                                                                                         >  1:1162074/1‑148 (MQ=255)
                                            cgcTTTCGCAGCTGGTGCAAAAGGCCGTGCCACGTAACTACGATGATTCGCTGCCGGGTGACAGTAAAGCATTCCTCGCGCAACTCTCGCTGCCCGCCCAACTGGCAAGCCAGCAAAGCGGTGTGCCACATCATTTGATCCTCGCTCa                                                                                                        <  1:1520020/148‑1 (MQ=255)
                                              cTTTCGCAGCTGGTGCAAAAGGCCGTGCCACGTAACTACGATGATTCGCTGCCGGGTGACAGTAAAGCATTCCTCGCGCAACTCTCGCTGCCCGCCCAACTGGCAAGCCAGCAAAGCGGTGTGCCACATCATTTGATCCTCGCTCAgg                                                                                                      <  1:297525/148‑1 (MQ=255)
                                              cTTTCGCAGCTGGTGCAAAAGGCCGTGCCACGTAACTACGATGATTCGCTGCCGGGTGACAGTAAAGCATTCCTCGCGCAACTCTCGCTGCCCGCCCAACTGGCAAGCCAGCAAAGCGGTGTGCCACATCATTTGATCCTCGCTCAgg                                                                                                      <  1:298007/148‑1 (MQ=255)
                                              cTTTCGCAGCTGGTGCAAAAGGCCGTGCCACGTAACTACGATGATTCGCTGCCGGGTGACAGTAAAGCATTCCTCGCGCAACTCTCGCTGCCCGCCCAACTGGCAAGCCAGCAAAGCGGTGTGCCACATCATTTGATCCTCGCTCAgg                                                                                                      >  1:697825/1‑148 (MQ=255)
                                                ttCGCAGCTGGTGCAAAAGGCCGTGCCACGTAACTACGATGATTCGCTGCCGGGTGACAGTAAAGCATTCCTCGCGCAACTCTCGCTGCCCGCCCAACTGGCAAGCCAGCAAAGCGGTGTGCCACATCATTTGATCCTCGCTCAggcg                                                                                                    >  2:726673/1‑148 (MQ=255)
                                                ttCGCAGCTGGTGCAAAAGGCCGTGCCACGTAACTACGATGATTCGCTGCCGGGTGACAGTAAAGCATTCCTCGCGCAACTCTCGCTGCCCGCCCAACTGGCAAGCCAGCAAAGCGGTGTGCCACATCATTTGATCCTCGCTCAggcg                                                                                                    >  2:2643240/1‑148 (MQ=255)
                                                   gCAGCTGGTGCAAAAGGCCGTGCCACGTAACTACGATGATTCGCTGCCGGGTGACAGTAAAGCATTCCTCGCGCAACTCTCGCTGCCCGCCCAACAGGCAAGCCAGCAAAGCGGTGTGCCACATCATTTGATCCTCGCTCAGGCGGCa                                                                                                 <  1:2202381/148‑1 (MQ=255)
                                                      gCTGGTGCAAAAGGCCGTGCCACGTAACTACGATGATTCGCTGCCGGGTGACAGTAAAGCATTCCTCGCGCAACTCTCGCTGCCCGCCCAACTGGCAAGCCAGCAAAGCGGTGTGCCACATCATTTGATCCTCGCTCAGGCGGCACTg                                                                                              >  2:2701968/1‑148 (MQ=255)
                                                        tGGTGCAAAAGGCCGTGCCACGTAACTACGATGATTCGCTGCCGGGTGACAGTAAAGCATTCCTCGCGCAACTCTCGCTGCCCGCCCAACTGGCAAGCCAGCAAAGCGGTGTGCCACATCATTTGATCCTCGCTCAGGCGGCACTGGa                                                                                            <  1:2599336/148‑1 (MQ=255)
                                                          gTGCAAAAGGCCGTGCCACGTAACTACGATGATTCGCTGCCGGGTGACAGTAAAGCATTCCTCGCGCAACTCTCGCTGCCCGCCCAACAGGCAAGCCAGCAAAGCGGTGTGCCACATCATTTGATCCTCGCTCAGGCGGCACTGGAAt                                                                                          >  1:3202591/1‑148 (MQ=255)
                                                          gTGCAAAAGGCCGTGCCACGTAACTACGATGATTCGCTGCCGGGTGACAGTAAAGCATTCCTCGCGCAACTCTCGCTGCCCGCCCAACAGGCAAGCCAGCAAAGCGGTGTGCCACATCATTTGATCCTCGCTCAGGCGGCACTGGAAt                                                                                          >  1:2669157/1‑148 (MQ=255)
                                                          gTGCAAAAGGCCGTGCCACGTAACTACGATGATTCGCTGCCGGGTGACAGTAAAGCATTCCTCGCGCAACTCTCGCTGCCCGCCCAACAGGCAAGCCAGCAAAGCGGTGTGCCACATCATTTGATCCTCGCTCAGGCGGCACTGGAAt                                                                                          >  1:2669136/1‑148 (MQ=255)
                                                             cAAAAGGCCGTGCCACGTAACTACGATGATTCGCTGCCGGGTGACAGTAAAGCATTCCTCGCGCAACTCTCGCTGCCCGCCCAACAGGCAAGCCAGCAAAGCGGTGTGCCACATCATTTGATCCTCGCTCAGGCGGCACTGGAATCTg                                                                                       <  2:2946355/148‑1 (MQ=255)
                                                              aaaaGGCCGTGCCACGTAACTACGATGATTCGCTGCCGGGTGACAGTAAAGCATTCCTCGCGCAACTCTCGCTGCCCGCCCAACAGGCAAGCCAGCAAAGCGGTGTGCCACATCATTTGATCCTCGCTCAGGCGGCACTGGAATCTgg                                                                                      >  2:3471289/1‑148 (MQ=255)
                                                               aaaGGCCGTGCCACGTAACTACGATGATTCGCTGCCGGGTGACAGTAAAGCATTCCTCGCGCAACTCTCGCTGCCCGCCCAACTGGCAAGCCAGCAAAGCGGTGTGCCACATCATTTGATCCTCGCTCAGGCGGCACTGGAATCTGGt                                                                                     >  1:3547290/1‑148 (MQ=255)
                                                                   gCCGTGCCACGTAACTACGATGATTCGCTGCCGGGTGACAGTAAAGCATTCCTCGCGCAACTCTCGCTGCCCGCCCAACTGGCAAGCCAGCAAAGCGGTGTGCCACATCATTTGATCCTCGCTCAGGCGGCACTGGAATCTGGTTggg                                                                                 <  2:632349/148‑1 (MQ=255)
                                                                      gTGCCACGTAACTACGATGATTCGCTGCCGGGTGACAGTAAAGCAGTCCTCGCGCAACTCTCGCTGCCCGCCCAACAGGCAAGCCAGCAAAGCGGTGTGCCACATCATTTGATCCTCGCTCAGGCGGCACTGGAATCTGGTTGGGGGc                                                                              <  2:1789809/148‑1 (MQ=255)
                                                                       tGCCACGTAACTACGATGATTCGCTGCCGGGTGACAGTAAAGCATTCCTCGCGCAACTCTCGCTGCCCGCCCAACTGGCAAGCCAGCAAAGCGGTGTGCCACATCATTTGATCCTCGCTCAGGCGGCACTGGAATCTGGTTGGGGGCa                                                                             <  1:261269/148‑1 (MQ=255)
                                                                       tGCCACGTAACTACGATGATTCGCTGCCGGGTGACAGTAAAGCATTCCTCGCGCAACTCTCGCTGCCCGCCCAACTGGCAAGCCAGCAAAGCGGTGTGCCACATCATTTGATCCTCGCTCAGGCGGCACTGGAATCTGGTTGGGGGCa                                                                             <  1:260709/148‑1 (MQ=255)
                                                                       tGCCACGTAACTACGATGATTCGCTGCCGGGTGACAGTAAAGCATTCCTCGCGCAACTCTCGCTGCCCGCCCAACTGGCAAGCCAGCAAAGCGGTGTGCCACATCATTTGATCCTCGCTCAGGCGGCACTGGAATCTGGTTGGGGGCa                                                                             >  2:3697147/1‑148 (MQ=255)
                                                                           aCGTAACTACGATGATTCGCTGCCGGGTGACAGTAAAGCATTCCTCGCGCAACTCTCGCTGCCCGCCCAACTGGCAAGCCAGCAAAGCGGTGTGCCACATCATTTGATCCTCGCTCAGGCGGCACTGGAATCTGGTTGGGGGCAACgg                                                                         >  1:2275763/1‑148 (MQ=255)
                                                                           aCGTAACTACGATGATTCGCTGCCGGGTGACAGTAAAGCATTCCTCGCGCAACTCTCGCTGCCCGCCCAACAGGCAAGCCAGCAAAGCGGTGTGCCACATCATTTGATCCTCGCTCAGGCGGCACTGGAATCTGGTTGGGGGCAACgg                                                                         >  2:324989/1‑148 (MQ=255)
                                                                            cGTAACTACGATGATTCGCTGCCGGGTGACAGTAAAGCATTCCTCGCGCAACTCTCGCTGCCCGCCCAACAGGCAAGCCAGCAAAGCGGTGTGCCACATCATTTGATCCTCGCTCAGGCGGCACTGGAATCTGGTTGGGGGCAACGGc                                                                        >  1:532985/1‑148 (MQ=255)
                                                                               aaCTACGATGATTCGCTGCCGGGTGACAGTAAAGCATTCCTCGCGCAACTCTCGCTGCCCGCCCAACTGGCAAGCCAGCAAAGCGGTGTGCCACATCATTTGATCCTCGCTCAGGCGGCACTGGAATCTGGTTGGGGGCAACGGCaaa                                                                     >  2:1680354/1‑148 (MQ=255)
                                                                               aaCTACGATGATTCGCTGCCGGGTGACAGTAAAGCATTCCTCGCGCAACTCTCGCTGCCCGCCCAACTGGCAAGCCAGCAAAGCGGTGTGCCACATCATTTGATCCTCGCTCAGGCGGCACTGGAATCTGGTTGGGGGCAACGGCaaa                                                                     >  1:1991566/1‑148 (MQ=255)
                                                                                aCTACGATGATTCGCTGCCGGGTGACAGTAAAGCATTCCTCGCGCAACTCTCGCTGCCCGCCCAACTGGCAAGCCAGCAAAGCGGTGTGCCACATCATTTGATCCTCGCTCAGGCGGCACTGGAATCTGGTTGGGGGCAACGGCAAAt                                                                    >  1:944974/1‑148 (MQ=255)
                                                                                 cTACGATGATTCGCTGCCGGGTGACAGTAAAGCATTCCTCGCGCAACTCTCGCTGCCCGCCCAACTGGCAAGCCAGCAAAGCGGTGTGCCACATCATTTGATCCTCGCTCAGGCGGCACTGGAATCTGGTTGGGGGCAACGGCAAATc                                                                   <  1:814497/148‑1 (MQ=255)
                                                                                 cTACGATGATTCGCTGCCGGGTGACAGTAAAGCATTCCTCGCGCAACTCTCGCTGCCCGCCCAACTGGCAAGCCAGCAAAGCGGTGTGCCACATCATTTGATCCTCGCTCAGGCGGCACTGGAATCTGGTTGGGGGCAACGGCAAATc                                                                   <  1:892809/148‑1 (MQ=255)
                                                                                       tgatTCGCTGCCGGGTGACAGTAAAGCATTCCTCGCGCAACTCTCGCTGCCCGCCCAACTGGCAAGCCAGCAAAGCGGTGTGCCACATCATTTGATCCTCGCTCAGGCGGCACTGGAATCTGGTTGGGGGCAACGGCAAATCcgccgc                                                             >  2:3531034/1‑148 (MQ=255)
                                                                                        gatTCGCTGCCGGGTGACAGTAAAGCATTCCTCGCGCAACTCTCGCTGCCCGCCCAACTGGCAAGCCAGCAAAGCGGTGTGCCACATCATTTGATCCTCGCTCAGGCGGCACTGGAATCTGGTTGGGGGCAACGGCAAATCCGCcgcg                                                            >  2:2834677/1‑148 (MQ=255)
                                                                                        gatTCGCTGCCGGGTGACAGTAAAGCATTCCTCGCGCAACTCTCGCTGCCCGCCCAACTGGCAAGCCAGCAAAGCGGTGTGCCACATCATTTGATCCTCGCTCAGGCGGCACTGGAATCTGGTTGGGGGCAACGGCAAATCCGCcgcg                                                            >  2:2834804/1‑148 (MQ=255)
                                                                                         atTCGCTGCCGGGTGACAGTAAAGCATTCCTCGCGCAACTCTCGCTGCCCGCCCAACAGGCAAGCCAGCAAAGCGGTGTGCCACATCATTTGATCCTCGCTCAGGCGGCACTGGAATCTGGTTGGGGGCAACGGCAAATCCGCCGCGa                                                           >  2:1177771/1‑148 (MQ=255)
                                                                                          ttCGCTGCCGGGTGACAGTAAAGCATTCCTCGCGCAACTCTCGCTGCCCGCCCAACTGGCAAGCCAGCAAAGCGGTGTGCCACATCATTTGATCCTCGCTCAGGCGGCACTGGAATCTGGTTGGGGGCAACGGCAAATCCGCCGCGaa                                                          >  2:1884099/1‑148 (MQ=255)
                                                                                           tCGCTGCCGGGTGACAGTAAAGCATTCCTCGCGCAACTCTCGCTGCCCGCCCAACTGGCAAGCCAGCAAAGCGGTGTGCCACATCATTTGATCCTCGCTCAGGCGGCACTGGAATCTGGTTGGGGGCAACGGCAAATCCGCCGCGaaa                                                         <  2:1560606/148‑1 (MQ=255)
                                                                                                      tGACAGTAAAGCATTCCTCGCGCAACTCTCGCGGCCCGCCCAACTGGCAAGCCAGCAAAGCGGTGTGCCACATCATTTGATCCTCGCTCAGGCGGCACTGGAATCTGGTTGGGGGCAACGGCAAATCCGCCGCGAAAACGGcgagccg                                              <  2:865001/148‑1 (MQ=255)
                                                                                                        aCAGTAAAGCATTCCTCGCGCAACTCTCGCTGCCCGCCCAACTGGCAAGCCAGCAAAGCGGTGTGCCACATCATTTGATCCTCGCTCAGGCGGCACTGGAATCTGGTTGGGGGCAACGGCAAATCCGCCGCGAAAACGGcgagccgag                                            <  1:3481324/148‑1 (MQ=255)
                                                                                                         cAGTAAAGCATTCCTCGCGCAACTCTCGCTGCCCGCCCAACAGGCAAGCCAGCAAAGCGGTGTGCCACATCATTTGATCCTCGCTCAGGCGGCACTGGAATCTGGTTGGGGGCAACGGCAAATCCGCCGCGAAAACGGcgagccgagc                                           >  1:3453117/1‑148 (MQ=255)
                                                                                                          aGTAAAGCATTCCTCGCGCAACTCTCGCTGCCCGCCCAACTGGCAAGCCAGCAAAGCGGTGTGCCACATCATTTGATCCTCGCTCAGGCGGCACTGGAATCTGGTTGGGGGCAACGGCAAATCCGCCGCGAAAACGGCGAGCCGAGCt                                          <  2:1408556/148‑1 (MQ=255)
                                                                                                           gTAAAGCATTCCTCGCGCAACTCTCTCTGCCCGCCCAACTGGCAAGCCAGCAAAGCGGTGTGCCACATCATTTGATCCTCGCTCAGTCTGCACTGGATTCTGGTTGGGGGGAACGGCAAATCCGCCGCGAAAACGGCGAGCCGATCta                                         <  2:545563/148‑1 (MQ=255)
                                                                                                           gTAAAGCATTCCTCGCGCAACTCTCGCTGCCCGCCCAACTGGCAAGCCAGCAAAGCGGTGTGCCACATCATTTGATCCTCGCTCAGGCGGCACTGGAATCTGGTTGGGGGCAACGGCAAATCCGCCGCGAAAACGGCGAGCCGAGCta                                         <  2:432652/148‑1 (MQ=255)
                                                                                                            tAAAGCATTCCTCGCGCAACTCTCGCTGCCCGCCCAACTGGCAAGCCAGCAAAGCGGTGTGCCACATCATTTGATCCTCGCTCAGGCGGCACTGGAATCTGGTTGGGGGCAACGGCAAATCCGCCGCGAAAACGGCGAGCCGAGCtat                                        >  1:1792358/1‑148 (MQ=255)
                                                                                                               aGCATTCCTCGCGCAACTCTCGCTGCCCGCCCAACTGGCAAGCCAGCAAAGCGGTGTGCCACATCATTTGATCCTCGCTCAGGCGGCACTGGAATCTGGTTGGGGGCAACGGCAAATCCGCCGCGAAAACGGCGAGCCGAGCTATAAc                                     >  2:1396269/1‑148 (MQ=255)
                                                                                                               aGCATTCCTCGCGCAACTCTCGCTGCCCGCCCAACTGGCAAGCCAGCAAAGCGGTGTGCCACATCATTTGATCCTCGCTCAGGCGGCACTGGAATCTGGTTGGGGGCAACGGCAAATCCGCCGCGAAAACGGCGAGCCGAGCTATAAc                                     >  2:1901551/1‑148 (MQ=255)
                                                                                                                gCATTCCTCGCGCAACTCTCGCTGCCCGCCCAACTGGCAAGCCAGCAAAGCGGTGTGCCACATCATTTGATCCTCGCTCAGGCGGCACTGGAATCTGGTTGGGGGCAACGGCAAATCCGCCGCGAAAACGGCGAGCCGAGCTATAAcc                                    >  1:3608743/1‑148 (MQ=255)
                                                                                                                gCATTCCTCGCGCAACTCTCGCTGCCCGCCCAACAGGCAAGCCAGCAAAGCGGTGTGCCACATCATTTGATCCTCGCTCAGGCGGCACTGGAATCTGGTTGGGGGCAACGGCAAATCCGCCGCGAAAACGGCGAGCCGAGCTATAAcc                                    <  1:46838/148‑1 (MQ=255)
                                                                                                                  aTTCCTCGCGCAACTCTCGCTGCCCGCCCAACTGGCAAGCCAGCAAAGCGGTGTGCCACATCATTTGATCCTCGCTCAGGCGGCACTGGAATCTGGTTGGGGGCAACGGCAAATCCGCCGCGAAAACGGCGAGCCGAGCTATAACCTg                                  <  2:2862766/148‑1 (MQ=255)
                                                                                                                   ttCCTCGCGCAACTCTCGCTGCCCGCCCAACTGGCAAGCCAGCAAAGCGGTGTGCCACATCATTTGATCCTCGCTCAGGCGGCACTGGAATCTGGTTGGGGGCAACGGCAAATCCGCCGCGAAAACGGCGAGCCGAGCTATAACCTGt                                 <  2:621526/148‑1 (MQ=255)
                                                                                                                   ttCCTCGCGCAACTCTCGCTGCCCGCCCAACTGGCAAGCCAGCAAAGCGGTGTGCCACATCATTTGATCCTCGCTCAGGCGGCACTGGAATCTGGTTGGGGGCAACGGCAAATCCGCCGCGAAAACGGCGAGCCGAGCTATAACCTGt                                 >  1:647376/1‑148 (MQ=255)
                                                                                                                    tCCTCGCGCAACTCTCGCTGCCCGCCCAACTGGCAAGCCAGCAAAGCGGTGTGCCACATCATTTGATCCTCGCTCAGGCGGCACTGGAATCTGGTTGGGGGCAACGGCAAATCCGCCGCGAAAACGGCGAGCCGAGCTATAACCTGtt                                <  2:2126061/148‑1 (MQ=255)
                                                                                                                    tCCTCGCGCAACTCTCGCTGCCCGCCCAACTGGCAAGCCAGCAAAGCGGTGTGCCACATCATTTGATCCTCGCTCAGGCGGCACTGGAATCTGGTTGGGGGCAACGGCAAATCCGCCGCGAAAACGGCGAGCCGAGCTATAACCTGtt                                <  2:2862169/148‑1 (MQ=255)
                                                                                                                    tCCTCGCGCAACTCTCGCTGCCCGCCCAACTGGCAAGCCAGCAAAGCGGTGTGCCACATCATTTGATCCTCGCTCAGGCGGCACTGGAATCTGGTTGGGGGCAACGGCAAATCCGCCGCGAAAACGGCGAGCCGAGCTATAACCTGtt                                <  2:197574/148‑1 (MQ=255)
                                                                                                                    tCCTCGCGCAACTCTCGCTGCCCGCCCAACAGGCAAGCCAGCAAAGCGGTGTGCCACATCATTTGATCCTCGCTCAGGCGGCACTGGAATCTGGTTGGGGGCAACGGCAAATCCGCCGCGAAAACGGCGAGCCGAGCTATAACCTGtt                                <  1:537221/148‑1 (MQ=255)
                                                                                                                    tCCTCGCGCAACTCTCGCGGCCCGCCCAACAGGCAAGCCAGAAATGCGTGTTGCAACATCATTTTTTCCTCGCACTGGCGGAACTGGAATCTGGTTGGGGGAAACGCCAAATCCGCCGACAAAACGGCGAGCCGAGCTATAACCTGtt                                <  1:324989/148‑1 (MQ=255)
                                                                                                                      cTCGCGCAACTCTCGCTGCCCGCCCAACTGGCAAGCCAGCAAAGCGGTGTGCCACATCATTTGATCCTCGCTCAGGCGGCACTGGAATCTGGTTGGGGGCAACGGCAAATCCGCCGCGAAAACGGCGAGCCGAGCTATAACCTGTTTg                              <  1:607045/148‑1 (MQ=255)
                                                                                                                          cgcAACTCTCGCTGCCCGCCCAACAGGCAAGCCAGCAAAGCGGTGTGCCACATCATTTGATCCTCGCTCTGGCGGCACTGGAATCTGGTTGGGGGCAACGGCAAATCCGCCGCGAAAACGGCGAGCCGAGCTATAACCTGTTTGgtgt                          <  1:2328363/148‑1 (MQ=255)
                                                                                                                            cAACTCTCGCTGCCCGCCCAACTGGCAAGCCAGCAAAGCGGTGTGCCACATCATTTGATCCTCGCTCAGGCGGCACTGGAATCTGGTTGGGGGCAACGGCAAATCCGCCGCGAAAACGGCGAGCCGAGCTATAACCTGTTTGGTGTCa                        >  1:1295267/1‑148 (MQ=255)
                                                                                                                            cAACTCTCGCTGCCCGCCCAACTGGCAAGCCAGCAAAGCGGTGTGCCACATCATTTGATCCTCGCTCAGGCGGCACTGGAATCTGGTTGGGGGCAACGGCAAATCCGCCGCGAAAACGGCGAGCCGAGCTATAACCTGTTTGGTGTCa                        >  1:1295248/1‑148 (MQ=255)
                                                                                                                                 ctcGCTGCCCGCCCAACTGGCAAGCCAGCAAAGCGGTGTGCCACATCATTTGATCCTCGCTCAGGCGGCACTGGAATCTGGTTGGGGGCAACGGCAAATCCGCCGCGAAAACGGCGAGCCGAGCTATAACCTGTTTGGTGTCAAAGcc                   >  2:2628578/1‑148 (MQ=255)
                                                                                                                                 ctcGCTGCCCGCCCAACTGGCAAGCCAGCAAAGCGGTGTGCCACATCATTTGATCCTCGCTCAGGCGGCACTGGAATCTGGTTGGGGGCAACGGCAAATCCGCCGCGAAAACGGCGAGCCGAGCTATAACCTGTTTGGTGTCAAAGcc                   >  2:1001278/1‑148 (MQ=255)
                                                                                                                                  tcGCTGCCCGCCCAACTGGCAAGCCAGCAAAGCGGTGTGCCACATCATTTGATCCTCGCTCAGGCGGCACTGGAATCTGGTTGGGGGCAACGGCAAATCCGCCGCGAAAACGGCGAGCCGAGCTATAACCTGTTTGGTGTCAAAGCct                  <  1:2729426/148‑1 (MQ=255)
                                                                                                                                  tcGCTGCCCGCCCAACTGGCAAGCCAGCAAAGCGGTGTGCCACATCATTTGATCCTCGCTCAGGCGGCACTGGAATCTGGTTGGGGGCAACGGCAAATCCGCCGCGAAAACGGCGAGCCGAGCTATAACCTGTTTGGTGTCAAAGCct                  <  1:2729116/148‑1 (MQ=255)
                                                                                                                                  tcGCTGCCCGCCCAACTGGCAAGCCAGCAAAGCGGTGTGCCACATCATTTGATCCTCGCTCAGGCGGCACTGGAATCTGGTTGGGGGCAACGGCAAATCCGCCGCGAAAACGGCGAGCCGAGCTATAACCTGTTTGGTGTCAAAGCct                  <  2:3608743/148‑1 (MQ=255)
                                                                                                                                   cGCTGCCCGCCCAACTGGCAAGCCAGCAAAGCGGTGTGCCACATCATTTGATCCTCGCTCAGGCGGCACTGGAATCTGGTTGGGGGCAACGGCAAATCCGCCGCGAAAACGGCGAGCCGAGCTATAACCTGTTTGGTGTCAAAGCctc                 <  1:205035/148‑1 (MQ=255)
                                                                                                                                   cGCTGCCCGCCCAACTGGCAAGCCAGCAAAGCGGTGTGCCACATCATTTGATCCTCGCTCAGGCGGCACTGGAATCTGGTTGGGGGCAACGGCAAATCCGCCGCGAAAACGGCGAGCCGAGCTATAACCTGTTTGGTGTCAAAGCctc                 <  1:205068/148‑1 (MQ=255)
                                                                                                                                         ccgcccAACTGGCAAGCCAGCAAAGCGGTGTGCCACATCATTTGATCCTCGCTCAGGCGGCACTGGAATCTGGTTGGGGGCAACGGCAAATCCGCCGCGAAAACGGCGAGCCGAGCTATAACCTGTTTGGTGTCAAAGCCTCTGGCaa           <  1:2178238/148‑1 (MQ=255)
                                                                                                                                         ccgcccAACTGGCAAGCCAGCAAAGCGGTGTGCCACATCATTTGATCCTCGCTCAGGCGGCACTGGAATCTGGTTGGGGGCAACGGCAAATCCGCCGCGAAAACGGCGAGCCGAGCTATAACCTGTTTGGTGTCAAAGCCTCTGGCaa           <  2:1694476/148‑1 (MQ=255)
                                                                                                                                         ccgcccAACTGGCAAGCCAGCAAAGCGGTGTGCCACATCATTTGATCCTCGCTCAGGCGGCACTGGAATCTGGTTGGGGGCAACGGCAAATCCGCCGCGAAAACGGCGAGCCGAGCTATAACCTGTTTGGTGTCAAAGCCTCTGGCaa           <  2:1694284/148‑1 (MQ=255)
                                                                                                                                              cAACAGGCAAGCCAGCAAAGCGGTGTGCCACATCATTTGATCCTCGCTCAGGCGGCACTGGAATCTGGTTGGGGGCAACGGCAAATCCGCCGCGAAAACGGCGAGCCGAGCTATAACCTGTGTGGTGTCAAAGCCTCTGGCAACTGGa      >  2:2351226/1‑148 (MQ=255)
                                                                                                                                               aaCTGGCAAGCCAGCAAAGCGGTGTGCCACATCATTTGATCCTCGCTCAGGCGGCACTGGAATCTGGTTGGGGGCAACGGCAAATCCGCCGCGAAAACGGCGAGCCGAGCTATAACCTGTTTGGTGTCAAAGCCTCTGGCAACTGGaa     <  1:180011/148‑1 (MQ=255)
                                                                                                                                                  tGGCAAGCCAGCAAAGCGGTGTGCCACATCATTTGATCCTCGCTCAGGCGGCACTGGAATCTGGTTGGGGGCAACGGCAAATCCGCCGCGAAAACGGCGAGCCGAGCTATAACCTGTTTGGTGTCAAAGCCTCTGGCAACTGGAAAgg  <  1:2487641/148‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
GATGAAATTCCCGCTCGAAACTGTGGTGCGTTATCAAAATCAGGCGCTTTCGCAGCTGGTGCAAAAGGCCGTGCCACGTAACTACGATGATTCGCTGCCGGGTGACAGTAAAGCATTCCTCGCGCAACTCTCGCTGCCCGCCCAACTGGCAAGCCAGCAAAGCGGTGTGCCACATCATTTGATCCTCGCTCAGGCGGCACTGGAATCTGGTTGGGGGCAACGGCAAATCCGCCGCGAAAACGGCGAGCCGAGCTATAACCTGTTTGGTGTCAAAGCCTCTGGCAACTGGAAAGG  >  CP009273/1133165‑1133458

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: