Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 169,092 T→G 8.3% L382V (TTG→GTG)  fhuB → fused iron‑hydroxamate transporter subunits of ABC superfamily: membrane components

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273169,0920TG 8.3% 769.9 / 18.9 257L382V (TTG→GTG) fhuBfused iron‑hydroxamate transporter subunits of ABC superfamily: membrane components
Reads supporting (aligned to +/- strand):  ref base T (75/155);  new base G (25/2);  total (100/157)
Fisher's exact test for biased strand distribution p-value = 1.18e-09
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 264 total aligned reads displayed.

GGATATGAAGGTCAACGATCGTGTCGCGGCTGAACGCCAACATGTGCTGGCGTTTGCCCTCGCGGGCGGCGTGCTGCTGTTGATGGCTGTGGTGGTGGCGCTGTCGTTTGGTCGTGATGCGCACGGCTGGACGTGGGCGAGCGGGGCGTTGCTCGAGGATTTAATGCCCTGGCGCTGGCCGCGAATTATGGCGGCGCTGTTTGCGGGCGTCATGCTGGCGGTGGCGGGCTGTATTATTCAGCGACTGACCGGAAACCCGATGGCAAGCCCGGAAGTGCTGGGGATTAGCTCCGGCGC  >  CP009273/168944‑169240
                                                                                                                                                    |                                                                                                                                                    
ggATATGAAGGTCAACGATCGTGTCGCGGCTGAACGCCAACATGTGCTGGCGTTTGCCCACGCGGGCGGCGTGCTGCTGTTGATGGCTGTGGTGGTGGCGCTGTCGTTTTGTAGTGTTTCGCACCGCGGGCCCTGGGCGAGCGGGGCGt                                                                                                                                                      >  1:4521890/1‑149 (MQ=255)
 gATATGAAGGTCAACGATCGTGTCGCGGCTGAACGCCAACATGTGCTGGCGTTTGCCCTCGCGGGCGGCGTGCTGCTGTTGATGGCTGTGGTGGTGGCGCTGTCGTTTGGTCGTGATGCGCACGGCTGGACGTGGGCGAGCGGGGCGtt                                                                                                                                                     <  2:1906830/149‑1 (MQ=255)
 gATATGAAGGTCAACGATCGTGTCGCGGCTGAACGCCAACATGTGCTGGCGTTTGCCCTCGCGGGCGGCGTGCTGCTGTTGATGGCTGTGGTGGTGGCGCTGTCGTTTGGTCGTGATGCGCACGGCTGGACGTGGGCGAGCGGGGCGtt                                                                                                                                                     <  2:4066299/149‑1 (MQ=255)
  atatGAAGGTCAACGATCGTGTCGCGGCTGAACGCCAACATGTGCTGGCGTTTGCCCTCGCGGGCGGCGTGCTGCTGTTGATGGCTGTGGTGGTGGCGCTGTCGTTTGGTCGTGATGCGCACCGCTGGACGTGGGCGAGCGGGGCCTTg                                                                                                                                                    >  2:4453115/1‑149 (MQ=255)
     tGAAGGTCAACGATCGTGTCGCGGCTGAACGCCAACATGTGCTGGCGTTTGCCCTCGCGGGCGGCGTGCTGCTGTTGATGGCTGTGGTGGTGGCGCTGTCGTTTGGTCGTTATGCGCACGGCTGGAAGTGGGCGAGCGGGGCGTTGCTc                                                                                                                                                 >  2:67783/1‑149 (MQ=255)
       aaGGTCAACGATCGTGTCGCGGCTGAACGCCAACATGTGCTGGCGTTTGCCCTCGCGGGCGGCGTGCTGCTGTTGATGGCTGTGGTGGTGGCGCTGTCGTTTGGTCGTGATGCGCACGGCTGGACGTGGGCGAGCGGGGCGTTGCTCGa                                                                                                                                               <  2:3699129/149‑1 (MQ=255)
       aaGGTCAACGATCGTGTCGCGGCTGAACGCCAACATGTGCTGGCGTTTGCCCACGCGGGCGGCGTGCTGCTGTTGATGGCGGTGGTGGTGGCGCTGTCGTTTGGTCGTGATGCGCACCGCTGGTCGTGGGCGAGCGGGGCGTTGCGCGa                                                                                                                                               >  1:745831/1‑149 (MQ=255)
        atgTCCACGGTCGGGTCGGGGCAGTACGGCAACCTGGGGTGGGGGTTGGCCACGGGGGGGGGGGGCCGGTGTTGTTGGGCGGTGTGGTGGCGCTGTCGTTTGGTCGTGATGCGCACGGCTGGACGTGGGCGAGCGGGGCGTTGCTCGAg                                                                                                                                              <  1:2253426/147‑1 (MQ=255)
           tCAACGATCGTGTCGCGGCTGAACGCCAACATGTGCTGGCGTTTGCCCTCGCGGGCGGCGTGCTGCTGTTGATGGCTGTGGTGGTGGCGCTGTCGTTTGGTCGTGATGCGCACGGCTGGACGTGGGCGAGCGGGGCGTTGCTCGAGGAt                                                                                                                                           <  2:1433326/149‑1 (MQ=255)
           tCAACGATCGTGTCGCGGCTGAACGCCAACATGTGCTGGCGTTTGCCCTCGCGGGCGGCGTGCTGCTGTTGATGGCTGTCGTGGTGGCGCTGTCGTTTGGTCGTGGTGCGCCCGGCGGGGCGTGGGCGCGCGGGGCGGTGCTCGAGGaa                                                                                                                                           >  2:450281/1‑148 (MQ=255)
            cAACGATCGTGTCGCGGCTGAACGCCACCAGGTGCTGGCGTTTGCCCTCGCGGGCGGCGTGCTGCTGTTGATGGCTGTGGTGGTGGCGCTGTCGTTTGGTCGTGATGCGCACGGCTGGACGTGGGCGAGCGGGGCGTTGCTCGAGGAtt                                                                                                                                          <  2:2945105/149‑1 (MQ=255)
                 aTCGTGTCGCGGCTGAACGCCAACATGTGCTGGCGTTTGCCCTCGCGGGCGGCGTGCTGCTGTTGATGGCTGTGGTGGTGGCGCTGTCGTTTGGGCGTGATGCGCCCGGCTGGACGTGGGCGAGCGGGGCGTGGCTCGCGGATTTTATg                                                                                                                                     >  2:2160540/1‑149 (MQ=255)
                 aTCGTGTCGCGGCTGAACGCCAACATGTGCTGGCGTTTGCCCTCGCGGGCGGCGTGCTGCTGTTGATGGCTGGGGTGGTGGCGCTGTCGTTTGGTCGTGATGCGCACGGCTGGACGTGGGCGAGCGGGGCGTTGCTCGAGGATTTAATg                                                                                                                                     >  2:662192/1‑149 (MQ=255)
                  tCGTGTCGCGGCTGAACGCCAACATGTGCTGGCGTTTGCCCTCGCGGGCGGCGTGCTGCTGTTGATGGCTGTGGTGGTGGCGCTGTCGTTTGGTCGTGATGCGCACGGCTGGACGTGGGCGAGCGGGGCGTTGCTCGAGGATTTAATGc                                                                                                                                    <  2:1347759/149‑1 (MQ=255)
                  tCGTGTCGCGGCTGAACGCCAACATGTGCTGGCGTTTGCCCTCGCGGGCGGCGTGCTGCTGTTGATGGCTGTGGTGGTGGCGCTGTCGTTTGGTCGTGATGCGCACGGCTGGACGTGGGCGAGCGGGGCGTTGCTCGAGGATTTAATGc                                                                                                                                    <  2:1579075/149‑1 (MQ=255)
                    gtgtCGCGGCTGAACGCCAACATGTGCTGGCGTTTGCCCTCGCGGGCGGCGTGCTGCTGTTGATGGCTGTGGTGGTGGCGCTGTCGTTTGGTGGTGATTCGCCCGGCGGGACGGTGGCGAGCGGGGCGGTGGTAGAGGATTTaaatgcc                                                                                                                                  >  2:3083873/1‑144 (MQ=255)
                      gtCGCGGCTGAACGCCAACATGTGCTGGCGTTTGCCCTCGCGGGCGGCGTGCTGCTGTTGATGGCTGTGGTGGTGGCGCTGTCGTTTGGTCGTGATGCGCACGGCTGGACGTGGGCGAGCGGGGCGTTGCTCGAGGATTTAATGCCCTg                                                                                                                                <  2:1202450/149‑1 (MQ=255)
                      gtCGCGGCTGAACGCCAACATGTGCTGGCGTTTGCCCTCGCGGGCGGCGTGCTGCTGTTGATGGCTGTGGTGGTGGCGCTGTCGTTTGGTCGTGATGCGCACGGCTGGACGTGGGCGAGCGGGGCGTTGCTCGAGGATTTAATGCCCTg                                                                                                                                <  2:205361/149‑1 (MQ=255)
                       tCGCGGCTGAACGCCAACATGTGCTGGCGTTTGCCCTCGCGGGCGGCGTGCTGCTGTTGATGGCTGTGGTGGTGGCGCTGTCGTTTGGTCGTGATGCGCACGGCTGGACGTGGGCGAGCGGGGCGTTGCTCGAGGATTTAATGCCCTgg                                                                                                                               <  1:3295850/149‑1 (MQ=255)
                           ggCTGAACGCCAACATGTGCTGGCGTTTGCCCTCGCGGGCGGCGTGCTGCTGTTGATGGCTGTGGTGGTGGCGCTGTCGTTTGGTCGTGATGCGCACGGCTGGACGTGGGCGAGCGGGGCGTTGCTCGAGGATTTAATGCCCTGGCGCt                                                                                                                           <  2:502911/149‑1 (MQ=255)
                            gCTGAACGCCAACATGTGCTGGCGTTTGCCCTCGCGGGCGGCGTGCTGCTGTTGATGGCGGTGGTGGTGGCGCTGTCGTTTTGTCGTGATGCGCACGGCTGGACGTGGGCGCGCGGGGCGTTGCGCGAGGATTTAATGCCCTGGCGCTg                                                                                                                          >  2:3980561/1‑149 (MQ=255)
                              tGAACGCCAACATGTGCTGGCGTTTGCCCTCGCGGGCGGCGTGCTGCTGTTGATGGCTGTGGTGGTGGCGCTGTCGTTTGGTCGTGATGCGCACGGCTGGACGTGGGCGAGCGGGGCGTTGCTCGAGGATTTAATGCCCTGGCGCTGGc                                                                                                                        <  1:1441424/149‑1 (MQ=255)
                              tGAACGCCAACATGTGCTGGCGTTTGCCCTCGCGGGCGGCGTGCTGCTGTTGATGGCGGTGGTGGTGGCGCTGTCGTTTGGTCGTGATGCGCACGGCTGGACGTGGGCGAGCGGGGCGTTGCTCGAGGATTTAATGCCCTGGCGCTGGc                                                                                                                        >  2:2246348/1‑149 (MQ=255)
                                   gCCAACATGTGCTGGCGTTTGCCCTCGCGGGCGGCGTGCTGCTGTTGATGGCTGTGGTGGTGGCGCTGTCGTTTGGTCGTGATGCGCACGGCTGGACGTGGGCGAGCGGGGCGTTGCTCGAGGATTTAATGCCCTGGCGCTGGCCGCGa                                                                                                                   <  2:2560979/149‑1 (MQ=255)
                                   gCCAACATGTGCTGGCGTTTGCCCTCGCGGGCGGCGTGCTGCTGTTGATGGCTGTGGTGGTGGCGCTGTCGTTTGGTCGTGATGCGCACGGCTGGACGTGGGCGAGCGGGGCGTTGCTCGAGGATTTAATGCCCTGGCGCTGGCCGCGa                                                                                                                   <  2:271121/149‑1 (MQ=255)
                                     cAACATGTGCTGGCGTTTGCCCTCGCGGGCGGCGTGCTGCTGTTGATGGCTGTGGTGGTGGCGCTGTCGTTTGGTCGTGATGCGCACGGCTGGACGTGGGCGAGCGGGGCGTTGCTCGAGGATTTAATGCCCTGGCGCTGGCCGCGAAt                                                                                                                 <  2:3623941/149‑1 (MQ=255)
                                       aCATGTGCTGGCGTTTGCCCTCGCGGGCGGCGTGCTGCTGTTGATGGCTGTGGTGGTGGCGCTGTCGTTTGGTCGTGATGCGCACGGCTGGACGTGGGCGAGCGGGGCGTTGCTCGAGGATTTAATGCCCTGGCGCTGGCCGCGAATTa                                                                                                               <  2:1514734/149‑1 (MQ=255)
                                       aCATGTGCTGGCGTTTGCCCTCGCGGGCGGCGTGCTGCTGTTGATGGCTGTGGTGGTGGCGCTGTCGTTTGGTCGTGATGCGCACGGCTGGACGTGGGCGAGCGGGGCGTTGCTCGAGGATTTAATGCCCTGGCGCTGGCCGCGAATTa                                                                                                               <  2:2637968/149‑1 (MQ=255)
                                        cATGTGCTGGCGTTTTCCCAAGCGGGCGGCGTGCTGCTGTTGATGGCTGTGGTGGTGGCGCTGTAGTTTTGTCGTGATTCCCACCGCAGGACGTGGGCGAGCCGGGCGGTGCTAGAGGAATTAATTCCCAGGCGCTTGCCGCCAAatat                                                                                                              >  2:394060/1‑149 (MQ=255)
                                           gtgCTGGCGTTTTCCCACGCGGGCGGCGTGCTGCTGTTGATGGCTGTGGTGGTGGCGCTGTCGTTTTGTCGTGATTCGCACCGCTGGAAGGGGGCGCGCGGGGCGGTGCTCGAGGAGTTAATTCCCAGGCGCTGGCCCCGAAATATtgc                                                                                                           >  2:2839225/1‑149 (MQ=255)
                                                ggCGTTTGCCCTCGCGGGCGGCGTGCTGCTGTTGATGGCTGTGGTGGTGGCGCTGTCGTTTTGTCGTGAAGCGCACGGCTGGACGTGGGCGAGCGGGGCGTTGCTCGAGGCTTTAATGCCCTGGCGCTGGCCGCGAATTATGGCGgcgc                                                                                                      >  1:1807881/1‑149 (MQ=255)
                                                ggCGTTTGCCCTCGCGGGCGGCGTGCTGCTGTTGATGGCTGTGGTGGTGGCGCTGTCGTTTGGTCGTGATGCGCACGGCTGGACGTGGGCGAGCGGGGCGTTGCTCGAGGATTTAATGCCCTGGCGCTGGCCGCGAATTATGGCGgcgc                                                                                                      <  1:399429/149‑1 (MQ=255)
                                                 gCGTTTGCCCTCGCGGGCGGCGTGCTGCTGTTGATGGCTGTGGTGGTGGCGCTGTCGTTTGGTCGTGATGCGCACGGCTGGACGTGGGCGAGCGGGGCGTTGCTCGAGGATTTAATGCCCTGGCGCTGGCCGCGAATTATGGCGGCGCt                                                                                                     <  2:3386885/149‑1 (MQ=255)
                                                 gCGTTTGCCCTCGCGGGCGGCGTGCTGCTGTTGATGGCTGTGGTGGTGGCGCTGTCGTTTGGTCGTGATGCGCACGGCTGGACGTGGGCGAGCGGGGCGGTGCTCGAGGATTTAATGCCCTGGCGCTGGCCGCGACTTATGGCGGCGCt                                                                                                     >  2:1275916/1‑149 (MQ=255)
                                                  cGTTTGCGCTCGCGGGCGGCGTGCTGCTGTTGTTGGCTGTGGTGGTGGCGCTGTCGTTTGGTCGTGATGCGCACGGCTGGACGTGGGCGAGCGGGGCGTTGCTCGAGGATTTAATGCCCTGGCGCTGGCCGCGAATTATGGCGGCGCTg                                                                                                    <  1:4019663/149‑1 (MQ=255)
                                                    tttGCCCTCGCGGGCGGCGTGCTGCTGTTGATGGCTGTGGTGGTGGCGCTGTCGTTTGGTCGTGTTGCGCACGGCTGGACGGGGGCGAGCGGGGCGTTGCTCGAGGATTTAATGCCCTGGCGCTGGCCGCGAATTATGGCGGCGCTGtt                                                                                                  >  2:1993133/1‑149 (MQ=255)
                                                    tttGCCCTCGCGGGCGGCGTGCTGCTGTTGATGGCTGTGGTGGTGGCGCTGTCGTTTGGTCGTGATGCGCACGGCTGGACGTGGGCGAGCGGGGCGTTGCTCGAGGATTTAATGCCCTGGCGCTGGCCGCGAATTATGGCGGCGCTGtt                                                                                                  >  2:611336/1‑149 (MQ=255)
                                                     ttGCCCTCGCGGGCGGCGTGCTGCTGTTGATGGCTGTGGTGGTGGCGCTGTCGTTTGGTCGTGATGCGCACGGCTGGACGTGGGCGAGCGGGGCGTTGCTCGAGGATTTAATGCCCTGGCGCTGGCCGCGAATTATGGCGGCGCTGttt                                                                                                 <  2:2977644/149‑1 (MQ=255)
                                                      tGCCCTCGCGGGCGGCGTGCTGCTGTTGATGGCTGTGGTGGTGGCGCTGTCGTTTGGTCGTGATGCGCACGGCTGGACGTGGGCGAGCGGGGCGTTGCTCGAGGATTTAATGCCCTGGCGCTGGCCGCGAATTATGGCGGGGCTGTTTg                                                                                                >  1:4536296/1‑149 (MQ=255)
                                                      tGCCCTCGCGGGCGGCGTGCTGCTGTTGATGGCTGTGGTGGTGGCGCTGTCGTTTGGTCGTGATGCGCACGGCTGGACGTGGGCGAGCGGGGCGTTGCTCGAGGATTTAATGCCCTGGCGCTGGCCGCGAATTATGGCGGCGCTGTTTg                                                                                                <  2:1559804/149‑1 (MQ=255)
                                                      tGCCCTCGCGGGCGGCGTGCTGCTGTTGATGGCTGTGGTGGTGGCGCTGTCGTTTGGTCGGGATGCGCACGGCGGGCCGTGGGCGCGCGGGGCGTTGGGCGAGGATTTTATGCCCAGGGGCTGGCCGCGAATTATGGCGGCGCTGTTTg                                                                                                >  1:2719516/1‑149 (MQ=255)
                                                        cccTCGCGGGCGGCGTGCTGCTGTTGATGGCTGTGGTGGTGGCGCTGTCGTTTGGTCGTGATGCGCACGGCTGGACGTGGGCGAGCGGGGCGTTGCTCGAGGATTTAATGCCCTGGCGCTGGCCGCGAATTATGGCGGCGCTGTTTGCg                                                                                              <  2:2291305/149‑1 (MQ=255)
                                                         ccTCGCGGGCGGCGTGCTGCTGTTGATGGCTGTGGTGGTGGCGCTGTCGTTTGGTCGTGATGCGCACGGCGGGACGTGGGCGAGCGGGGCGTTGCTCGAGGATTTAATGCCCTGGCGCTGGCCGCGAATTATGGCGGCGCTGTTTGCgg                                                                                             >  1:980849/1‑149 (MQ=255)
                                                          cTCGCGGGCGGCGTGCTGCTGTTGATGGCTGTGGTGGTGGCGCTGTCGTTTTGGCGGGATTCGCACGGCTGGACGTGGGCGTGCGGGGCGGTGCTCGAGGATTTAATGCCCTGGCGCTGGCCGCGCATTATGGGGGCGCTGTTTGCggg                                                                                            >  1:893963/1‑149 (MQ=255)
                                                          cTCGCGGGCGGCGTGCTGCTGTTGATGGCTGTGGTGGTGGCGCTGTCGTTTGGTCGTGATGCGCACGGCTGGACGTGGGCGAGCGGGGCGTTGCTCGAGGATTTAATTCCCTGGCGCTGGCCGCGAATTATTGCGGCGCTGTTTGCggg                                                                                            >  1:2167472/1‑149 (MQ=255)
                                                          cTCGCGGGCGGCGTGCTGCTGTTGATGGCTGTGGTGGTGGCGCTGTCGTTTGGTCGTGATGCGCACGGCTGGACGTGGGCGAGCGGGGCGTTGCTCGAGGATTTAATGCCCTGGCGCTGGCCGCGAATTATGGCGGCGCTGTTTGCggg                                                                                            <  1:2077910/149‑1 (MQ=255)
                                                           tCGCGGGCGGCGTGCTGCTGTTGATGGCTGTGGTGGTGGCGCTGTCGTTTGGTCGTGATGCGCACGGCTGGACGTGGGCGAGCGGGGCGTTGCTCGAGGATTTAATGCCCTGGCGCTGGCCGCGAATTATGGCGGCGCTGTTTGCGGGc                                                                                           <  2:1098332/149‑1 (MQ=255)
                                                             gcgggcggCGTGCTGCTGTTGATGGCTGTGGTGGTGGCGCTGTCGTTTGGTCGTGATGCGCACGGCTGGACGTGGGCGAGCGGGGCGTTGCTCGAGGATTTAATGCCCTGGCGCTGGCCGCGAATTATGGCGGCGCTGTTTGCGGGCGt                                                                                         <  1:2781870/149‑1 (MQ=255)
                                                              cgggcggCGTGCTGCTGTTGATGGCTGTGGTGGTGGCGCTGTCGTTTGGTCGTGATGCGCACGGCTGGACGTGGGCGAGCGGGGCGTTGCTCGAGGATTTAATGCCCTGGCGCTGGCCGCGAATTATGGCGGCGCTGTTTGCGGGCGTc                                                                                        <  1:1334824/149‑1 (MQ=255)
                                                              cgggcggCGTGCTGCTGTTGATGGCTGTGGTGGTGGCGCTGTCGTTTGGTCGTGATGCGCACGGCTGGACGTGGGCGAGCGGGGCGTTGCTCGAGGATTTAATGCCCTGGCGCTGGCCGCGAATTATGGCGGCGCTGTTTGCGGGCGTc                                                                                        <  1:1665672/149‑1 (MQ=255)
                                                               gggcggCGTGCTGCTGTTGATGGCTGTGGTGGTGGCGCTGTCGTTTGGTCGTGATGCGCACGGCTGGACGTGGGCGAGCGGGGCGTTGCTCGAGGATTTACTGCCCTGGCGCTGGCCGCGAATTATGGCGGCGCTGTTTGCGGGCGTCa                                                                                       <  1:4539687/149‑1 (MQ=255)
                                                               gggcggCGTGCTGCTGTTGATGGCTGTGGTGGTGGCGCTGTCGTTTGGTCGTGATGCGCACGGCTGGACGTGGGCGAGCGGGGCGTTGCTCGAGGATTTAATGCCCTGGCGCTGGCCGCGAATTATTGCGGCGCTGGTTGCGGGCGTCa                                                                                       >  1:4472005/1‑149 (MQ=255)
                                                                ggcggcGTGCTGCTGTTGATGGCTGTGGTGGTGGCGCTGTCGTTTGGTCGTGTTGCGCACGGCTGGACGTGGGCGAGCGGGGCGGTGCTCGAGGATTTAATGCCCGGGCGCGGGCCGCGAAATATGGCGGCGCTGGTTGCGGGCGTCaa                                                                                      >  1:3162350/1‑148 (MQ=255)
                                                                ggcggcGTGCTGCTGTTGATGGCTGTGGTGGTGGCGCTGTCGTTTGGTCGTGATGCGCACGGCTGGACGTGGGCGAGCGGGGCGTTGCTCGAGGATTTAATGCCCTGGCGCTGGCCGCGAATTATGGCGGCGCTGTTTGCGGGCGTCAt                                                                                      <  2:3895628/149‑1 (MQ=255)
                                                                 gcggcgTGCTGCTGTTGATGGCTGTGGTGGTGGCGCTGTCGTTTTGTCGTGATGCGCACGGCTGGACGTGGGCGATCGGGGCGTTGCTCGAGGATTTAATGCCCTTGCGCAGGCCGCGAAATAAGGCGGCGCAGTTTGCGGGCCTCaag                                                                                     >  1:305644/1‑147 (MQ=255)
                                                                  cggcgTGCTGCTGTTGATGGCTGTGGTGGTGGCGCTGTCGTTTGGTCGTGATGCGCACGGCTGGACGTGGGCGAGCGGGGCGTTGCTCGAGGATTTAATGCCCTGGCGCTGGCCGCGAATTATGGCGGCGCTGTTTGCGGGCGTCATGc                                                                                    <  2:1926426/149‑1 (MQ=255)
                                                                  cggcgTGCTGCTGTTGATGGCTGTGGTGGTGGCGCTGTCGTTTGGTCGTGATGCGCACGGCTGGACGTGGGCGAGCGGGGCGTTGCTCGAGGATTTAATGCCCTGGCGCTGGCCGCGAATTATGGCGGCGCTGTTTGCGGGCGTCATGc                                                                                    <  1:819879/149‑1 (MQ=255)
                                                                  cggcgTGCTGCTGTTGATGGCTGTGGTGGTGGCGCTGTCGTTTGGTCGTGATGCGCACGGCTGGACGTGGGCGAGCGGGGCGTTGCTCGAGGATTTAATGCCCTGGCGCTGGCCGCGAATTATGGCGGCGCTGTTTGCGGGCGTCATGc                                                                                    <  2:4318169/149‑1 (MQ=255)
                                                                   ggcgTGCTGCTGTTGATGGCTGTGGTGGTGGCGCTGTCGTTTGGTCGTGATGCGCACGGCTGGACGTGGGCGAGCGGGGCGTTGCTCGAGGATTTAATGCCCTGGCGCTGGCCGCGAATTATGGCGGCGCTGTTTGCGGGCGTCATGCt                                                                                   <  1:1765493/149‑1 (MQ=255)
                                                                   ggcgTGCTGCTGTTGATGGCTGTGGTGGTGGCGCTGTCGTTTGGTCGTGATGCGCACGGCTGGACGTGGGCGAGCGGGGCGTTGCTCGAGGATTTAATGCCCTGGCGCTGGCCGCGAATTATGGCGGCGCTGTTTGCGGGCGTCATGCt                                                                                   <  2:3312646/149‑1 (MQ=255)
                                                                    gcgTGCTGCTGTTGATGGCTGTGGTGGTGGCGCTGTCGTTTGGTCGTGATGCGCACGGCTGGACGTGGGCGAGCGGGGCGTTGCTCGAGGATTTAATGCCCTGGCGCTGGCCGCGAATTATGGCGGCGCTGTTTGCGGGCGTCATGCtg                                                                                  <  1:2347353/149‑1 (MQ=255)
                                                                       tgctgcTGTTGATGGCTGTGGTGGTGGCGCTGTCGTTTGGTCGTGATGCGCACGGCTGGACGTGGGCGAGCGGGGCGTTGCTCGAGGATTTAATGCCCTGGCGCTGGCCGCGAATTATGGCGGCGCTGTTTGCGGGCGTCATGCtggcg                                                                               <  1:1538982/149‑1 (MQ=255)
                                                                        gctgctGTTGATGGCTGTGGTGGTGGCGCTGTCGTTTGGTCGTGATGCGCACGGCTGGACGTGGGCGAGCGGGGCGTTGCTCGAGGATTTAATGCCCTGGCGCTGGCCGCGAATTATGGCGGCGCTGTTTGCGGGCGTCATGCtggcgg                                                                              <  1:831265/149‑1 (MQ=255)
                                                                          tgctgTTGATGGCTGTGGTGGTGGCGCTGTCGTTTGGTCGTGATGCGCACGGCTGGACGTGGGCGAGCGGGGCGTTGCTCGAGGATTTAATTCCCAGGCGCTGGCCGCGAATTATTGCGGCGCGGTTTTCGGGCGTCAATCtggcggtg                                                                            >  1:3641984/1‑149 (MQ=255)
                                                                          tgctgTTGATGGCTGTGGTGGTGGCGCTGTCGTTTGGTCGTGATGCGCACGGCTGGACGTGGGCGAGCGGGGCGTTGCTCGAGGATTTAATGCCCTGGCGCTGGCCGCGAATTATGGCGGCGCTGTTTGCGGGCGTCATGCtggcggtg                                                                            <  2:3685466/149‑1 (MQ=255)
                                                                           gctgTTGATGGCTGTGGTGGTGGCGCTGTCGTTTGGTCGTGATGCGCACGGCTGGACGTGGGCGAGCGGGGCGTTGCTCGAGGATTTAATGCCCTGGCGCTGGCCGCGAATTATGGCGGCGCTGTTTGCGGGCGTCATGCtggcggtgg                                                                           <  2:2872155/149‑1 (MQ=255)
                                                                              gTTGATGGCTGTGGTGGTGGCGCTGTCGTTTGGTCGTGATGCGCACGGCTGGACGTGGGCGAGCGGGGCGTTGCTCGAGGATTTAATGCCCTGGCGCTGGCCGCGAATTATGGCGGCGCTGTTTGCGGGCGTCATGCtggcggtggcgg                                                                        <  1:738891/149‑1 (MQ=255)
                                                                               ttGATGGCTGTGGTGGTGGCGCTGTCGTTTGGTCGTGATGCGCACGGCTGGACGTGGGCGAGCGGGGCGTTGCTCGAGGATTTAATGCCCTGGCGCTGGCCGCGAATTATGGCGGCGCTGTTTGCGGGCGTCATGCTGGCGGTGGCggg                                                                       <  1:4039216/149‑1 (MQ=255)
                                                                                   tGGCTGTGGTGGTGGCGCTGTCGTTTGGTGGTGATTCGCCAGGCTGGACGTGGGCGAGCGGGGCGGTTGTCGAGGATTTAATGACCACGGGCGGGCCGCGAAATATGGCGGCCCTGTTTTAGGGCGTCAATCGGGCGGTGGGGGGCGGt                                                                   >  1:4231550/1‑149 (MQ=255)
                                                                                    ggCTGTGGTGGTGGCGCTGTCGTTTGGTCGTGATGCGCACGGCTGGACGTGGGCGAGCGGGGCGTTGCTCGAGGATTTAATGCCCTGGCGCTGGCCGCGAATTATGGCGGCGCTGTTTGCGGGCGTCATGCTGGCGGTGGCGGGCTGta                                                                  <  1:3258775/149‑1 (MQ=255)
                                                                                    ggCTGTGGTGGTGGCGCTGTCGTTTGGTCGTGATGCGCACGGCTGGACGTGGGCGAGCGGGGCGTTGCTCGAGGATTTAATGCCCTGGCGCTGGCCGCGAATTATGGCGGCGCTGTTTGCGGGCGTCATGCTGGCGGTGGCGGGCTGta                                                                  <  2:3163950/149‑1 (MQ=255)
                                                                                    ggCTGTGGTGGTGGCGCTGTCGTTTGGTCGTGATGCGCACGGCTGGACGTGGGCGAGCGGGGCGTTGCTCGAGGAATTAATGCCCAGGCGCTGGCCGCGAATTATGGCGGCGCTGTTTGCGGGCGTCATGCTGGCGGTGGCGGGCTGta                                                                  >  2:4398277/1‑149 (MQ=255)
                                                                                      ccgtggtggtgGCGCTGTCGTTTGGTCGTGATGCGCACGGCTGGACGTGGGCGAGCGGGGCGTTGCTCGAGGATTTAATGCCCTGGCGCTGGCCGCGAATTATGGCGGCGCTGTTTGCGGGCGTCATGCTGGCGGTGGCGGGCTGtatt                                                                <  2:3534594/147‑1 (MQ=255)
                                                                                      cTGTGGTGGTGGCGCTGTCGTTTGGTCGTGATGCGCACGGCTGGACGTGGGCGAGCGGGGCGTTGCTCGAGGATTTAATGCCCTGGCGCTGGCCGCGAATTATGGCGGCGCTGTTTGCGGGCGTCATGCTGGCGGTGGCGGGCTGtatt                                                                <  2:251382/149‑1 (MQ=255)
                                                                                      cTGTGGTGGTGGCGCTGTCGTTTGGTCGTGATGCGCACGGCTGGACGTGGGCGAGCGGGGCGTTGCTCGAGGATTTAATGCCCTGGCGCTGGCCGCGAATTATGGCGGCGCTGTTTGCGGGCGTCATGCTGGCGGTGGCGGGCTGtatt                                                                <  1:4462972/149‑1 (MQ=255)
                                                                                       tgtgGTGGTGGCGCTGTCGTTTGGTCGTGATGCGCACGGCTGGACGTGGGCGAGCGGGGCGTTGCTCGAGGATTTAATGCCCTGGCGCTGGCCGCGAATTATGGCGGCGCTGTTTGCGGGCGTCATGCTGGCGGTGGCGGGCTGtatta                                                               <  2:1651466/149‑1 (MQ=255)
                                                                                        gtggtggtgGCGCTGTCGTTTGGTCGTGATGCGCACGGCTGGACGTGGGCGAGCGGGGCGTTGCTCGAGGATTTAATGCCCTGGCGCTGGCCGCGAATTATGGCGGCGCTGTTTGCGGGCGTCATGCTGGCGGTGGCGGGCTGtattat                                                              <  2:2031151/149‑1 (MQ=255)
                                                                                        gtggtggtgGCGCTGTCGTTTGGTCGTGATGCGCACGGCTGGACGTGGGCGAGCGGGGCGTTGCTCGAGGATTTAATGCCCTGGCGCTGGCCGCGAATTATGGCGGCGCTGTTTGCGGGCGTCATGCTGGCGGTGGCGGGCTGtattat                                                              <  1:1859758/149‑1 (MQ=255)
                                                                                            tggtggCGCTGTCGTTTGGTCGTGATGCGCACGGCTGGACGTGGGCGAGCGGGGCGTTGCTCGAGGATTTAATGCCCTGGCGCTGGCCGCGAATTATGGCGGCGCTGTTTGCGGGCGTCATGCTGGCGGTGGCGGGCTGTATTATTCAg                                                          <  1:3079697/149‑1 (MQ=255)
                                                                                            tggtggCGCTGTCGTTTGGTCGTGATGCGCACGGCTGGACGTGGGCGAGCGGGGCGTTGCTCGAGGATTTAATGCCCTGGCGCTGGCCGCGAATTATGGCGGCGCTGTTTGCGGGCGTCATGCTGGCGGTGGCGGGCTGTATTATTCAg                                                          <  1:2639015/149‑1 (MQ=255)
                                                                                             ggtggCGCTGTCGTTTGGTCGTGATGCGCACGGCTGGACGTGGGCGAGCGGGGCGTTGCTCGAGGATTTAATTCCCGGGCGCTGGCCGCGAATTATGGCGGCGCTGTTTGCGGGCGTGATTCTGGCGGTGGCGGGCTGTATTATTCAGc                                                         >  1:3219575/1‑149 (MQ=255)
                                                                                             ggtggCGCTGTCGTTTGGTCGTGATGCGCACGGCTGGACGTGGGCGAGCGGGGCGTTGCTCGAGGATTTAATGCCCTGGCGCTGGCCGCGAATTATGGCGGCGCTGTTTGCGGGCGTCATGCTGGCGGTGGCGGGCTGTATTATTAAGc                                                         >  1:1421884/1‑149 (MQ=255)
                                                                                               tggCGCTGTCGTTTGGTCGTGATGCGCACGGCTGGACGTGGGCGAGCGGGGCGTTGCTCGAGGATTTAATGCCCTGGCGCTGGCCGCGAATTATGGCGGCGCTGTTTGCGGGCGTCATGCTGGCGGTGGCGGGCTGTATTATTCAGCGa                                                       <  1:1333974/149‑1 (MQ=255)
                                                                                                ggCGCTGTCGTTTGGTCGTGATGCGCACGGCTGGACGTGGGCGAGCGGGGCGTTGCTCGAGGATTTAATGCCCTGGCGCTGGCCGCGAATTATGGCGGCGCTGTTTGCGGGCGTCATGCTGGCGGTGGCGGGCTGTATTATTCAGCGAc                                                      <  1:4299578/149‑1 (MQ=255)
                                                                                                ggCGCTGTCGTTTGGGCGTGATGCGCACGGCTGGACGTGGGCGAGCGGGGCGGTGCTCGAGGATTTAATGCCCTGGCGCTGGCCGCGAATTATGGCGGCGCTGTTTGCGGGCGTCATGCGGGCGGTGGCGGGCTGTATTATTCAGCGAc                                                      >  2:3617954/1‑149 (MQ=255)
                                                                                                 gcgcTGTCGTTTGGTCGTGATGCGCACGGCTGGACGTGGGCGAGCGGGGCGTTGCTCGAGGATTTAATGCCCTGGCGCTGGCCGCGAATTATGGCGGCGCTGTTTGCGGGCGTCATGCTGGCGGTGGCGGGCTGTATTATTCAGCGACt                                                     >  2:932763/1‑149 (MQ=255)
                                                                                                  cgcTGTCGTTTGGTCGTGATGCGCACGGCTGGACGTGGGCGAGCGGGGCGTTGCTCGAGGATTTAATGCCCTGGCGCTGGCCGCGAATTATGGCGGCGCTGTTTGCGGGCGTCATGCTGGCGGTGGCGGGCTGTATTATTCAGCGACTg                                                    <  2:3036642/149‑1 (MQ=255)
                                                                                                     tGTCGTTTGGTCGTGATGCGCACGGCTGGACGTGGGCGAGCGGGGCGTTGCTCGAGGATTTAATGCCCTGGCGCTGGCCGCGAATTATGGCGGCGCTGTTTGCGGGCGTCATGCTGGCGGTGGCGGGCTGTATTATTCAGCGACTGAcc                                                 >  1:657687/1‑149 (MQ=255)
                                                                                                          tttGGTCGTGATGCGCACGGCTGGACGTGGGCGAGCGGGGCGTTGCTCGAGGATTTAATGCCCTGGCGCTGGCCGCGAATTATGGCGGCGCTGTTTGCGGGCGTCATGCTGGCGGTGGCGGGCTGTATTATTCAGCGACTGACCGGaaa                                            <  1:1007134/149‑1 (MQ=255)
                                                                                                                             gCTGGACGTGGGCGAGCGGGGCGTTGCTCGAGGATTTAATGCCCTGGCGCTGGCCGCGAATTATGGCGGCGCTGTTTGCGGGCGTCATGCTGGCGGTGGCGGGCTGTATTATTCAGCGACTGACCGGAAACCCGATGGCAAGCCCGGaa                         <  2:1837142/149‑1 (MQ=255)
                                                                                                                                  tagTGGGCGTGCGGGGGGGTGGTAGTGGTATTAATGGCCTGGCGCTGGCCGCGAATTATGGCGGCGCTGTTTGCGGGCGTCATTGTGGCGGTGGCGGGCTGTATTATTCAGCGACTGACCGGAAACCCGATGGCAAGCCCGGAAGTGCt                    <  2:3305087/147‑1 (MQ=255)
                                                                                                                                       ggCGAGCGGGGCGTTGCTCGAGGATTTAATGCCCTGGCGCTGGCCGCGAATTATGGCGGCGCTGTTTGCGGGCGTCATGCTGGCGGTGGCGGGCTGTATTATTCAGCGACTGACCGGAAACCCGATGGCAAGCCCGGAAGTGCTGGGGa               <  1:4303627/149‑1 (MQ=255)
                                                                                                                                       ggCGAGCGGGGCGTTGCTCGAGGATTTAATGCCCTGGCGCTGGCCGCGAATTATGGCGGCGCTGTTTGCGGGCGTCATGCTGGCGGTGGCGGGCTGTATTATTCAGCGACTGACCGGAAACCCGATGGCAAGCCCGGAAGTGCTGGGGa               <  1:1119385/149‑1 (MQ=255)
                                                                                                                                            gCGGGGCGTTGCTCGAGGATTTAATGCCCTGGCGCTGGCCGCGAATTATGGCGGCGCTGTTTGCGGGCGTCATGCTGGCGGTGGCGGGCTGTATTATTCAGCGACTGACCGGAAACCCGATGGCAAGCCCGGAAGTGCTGGGGATTAGc          >  2:3730585/1‑149 (MQ=255)
                                                                                                                                            gCGGGGCGTTGCTCGAGGATTTAATGCCCTGGCGCTGGCCGCGAATTATGGCGGCGCTGTTTGCGGGCGTCATGCTGGCGGTGGCGGGCTGTATTATTCAGCGACTGACCGGAAACCCGATGGCAAGCCCGGAAGTGCTGGGGATTAGc          <  2:1020941/149‑1 (MQ=255)
                                                                                                                                            gCGGGGCGTTGCTCGAGGATTTAATGCCCTGGCGCTGGCCGCGAATTATGGCGGCGCTGTTTGCGGGCGTCATGCTGGCGGTGGCGGGCTGTATTATTCAGCGACTGACCGGAAACCCGATGGCAAGCCCGGAAGTGCTGGGGATTAGc          <  1:4406295/149‑1 (MQ=255)
                                                                                                                                               gggCGTTGCTCGAGGATTTAATGCCCTGGCGCTGGCCGCGAATTATGGCGGCGCTGTTTGCGGGCGTCATGCTGGCGGTGGCGGGCTGTATTATTCAGCGACTGACCGGAAACCCGATGGCAAGCCCGGAAGTGCTGGGTATTAGCTcc       >  1:3386315/1‑149 (MQ=255)
                                                                                                                                               gggCGTTGCTCGAGGATTTAATGCCCTGGCGCTGGCCGCGAATTATGGCGGCGCTGTTTGCGGGCGTCATGCTGGCGGTGGCGGGCTGTATTATTCAGCGACTGACCGGAAACCCGATGGCAAGCCCGGAAGTGCTGGGGATTAGCTcc       <  2:3813137/149‑1 (MQ=255)
                                                                                                                                                ggCGTTGCTCGAGGATTTAATGCCCTGGCGCTGGCCGCGAATTATGGCGGCGCTGTTTGCGGGCGTCATGCTGGCGGTGGCGGGCTGTATTATTCAGCGACTGACCGGAAACCCGATGGCAAGCCCGGAAGTGCTGGGGATTAGCTCcg      <  2:4262037/149‑1 (MQ=255)
                                                                                                                                                   gTTGCTCGAGGATTTAATGCCCTGGCGCTGGCCGCGAATTATGGCGGCGCTGTTTGCGGGCGTCATGCTGGCGGTGGCGGGCTGTATTATTCAGCGACTGACCGGAAACCCGATGGCAAGCCCGGAAGTGCTGGGGATTAGCTCcggcg   >  2:2166873/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
GGATATGAAGGTCAACGATCGTGTCGCGGCTGAACGCCAACATGTGCTGGCGTTTGCCCTCGCGGGCGGCGTGCTGCTGTTGATGGCTGTGGTGGTGGCGCTGTCGTTTGGTCGTGATGCGCACGGCTGGACGTGGGCGAGCGGGGCGTTGCTCGAGGATTTAATGCCCTGGCGCTGGCCGCGAATTATGGCGGCGCTGTTTGCGGGCGTCATGCTGGCGGTGGCGGGCTGTATTATTCAGCGACTGACCGGAAACCCGATGGCAAGCCCGGAAGTGCTGGGGATTAGCTCCGGCGC  >  CP009273/168944‑169240

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: