Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 608,575 T→A 25.7% L128Q (CTG→CAG)  fes → enterobactin/ferric enterobactin esterase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273608,5750TA25.7% 368.7 / 152.0 202L128Q (CTG→CAG) fesenterobactin/ferric enterobactin esterase
Reads supporting (aligned to +/- strand):  ref base T (71/79);  new base A (30/22);  total (101/101)
Fisher's exact test for biased strand distribution p-value = 2.60e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.98e-01

Only 100 of 202 total aligned reads displayed.

GACGACACAACTCAATGCCAACTGGCGCGGCAGCTACTGCTTTATTCCCACCGAACGCGATGACATTTTTTCTGTACCATCCCCCGATCGCCTCGAATTGCGCGAAGGCTGGCGAAAACTATTACCCCAGGCGATAGCCGATCCGCTGAACCTACAAAGCTGGAAAGGCGGGCGAGGGCACGCTGTTTCTGCACTCGAAATGCCGCAAGCGCCTCTGCAACCGGGATGGGATTGTCCGCAAGCGCCAGAAATACCTGCCAAAGAAATTATCTGGAAAAGTGAACGGTT  >  CP009273/608429‑608716
                                                                                                                                                  |                                                                                                                                             
 acgacACAACTCAATGCCAACTGGCGCGGCAGCTACTGCTTTATTCCCACCGAACGCGATGACATTTTTTCTGTACCATCCCCCGATCGCCTCGAATTGCGCGAAGGCTGGCGAAAACTATTACCCCAGGCGATAGCCGATCCGCAGa                                                                                                                                             >  1:2989988/1‑148 (MQ=255)
               tGCCAACTGGCGCGGCAGCTACTGCTTTATTCCCACCGAACGCGATGACATTTTTTCTGTACCATCCCCCGATCGCCTCGAATTGCGCGAAGGCTGGCGAAAACTATTACCCCAGGCGATAGCCGATCCGCTGAACCTACAAAGCTgg                                                                                                                               <  1:3568143/148‑1 (MQ=255)
                  cAACTGGCGCGGCAGCTACTGCTTTATTCCCACCGAACGCGATGACATTTTTTCTGTACCATCCCCCGATCGCCTCGAATTGCGCGAAGGCTGGCGAAAACTATTACCCCAGGCGATAGCCGATCCGCTGAACCTACAAAGCTGGaaa                                                                                                                            <  1:839862/148‑1 (MQ=255)
                   aaCTGGCGCGGCAGCTACTGCTTTATTCCCACCGAACGCGATGACATTTTTTCTGTACCATCCCCCGATCGCCTCGAATTGCGCGAAGGCTGGCGAAAACTATTACCCCAGGCGATAGCCGATCCGCTGAACCTACAAAGCTGGAAAg                                                                                                                           <  2:2540612/148‑1 (MQ=255)
                   aaCTGGCGCGGCAGCTACTGCTTTATTCCCACCGAACGCGATGACATTTTTTCTGTACCATCCCCCGATCGCCTCGAATTGCGCGAAGGCTGGCGAAAACTATTACCCCAGGCGATAGCCGATCCGCTGAACCTACAAAGCTGGAAAg                                                                                                                           >  2:751741/1‑148 (MQ=255)
                   aaCTGGCGCGGCAGCTACTGCTTTATTCCCACCGAACGCGATGACATTTTTTCTGTACCATCCCCCGATCGCCTCGAATTGCGCGAAGGCTGGCGAAAACTATTACCCCAGGCGATAGCCGATCCGCTGAACCTACAAAGCTGGAAAg                                                                                                                           >  2:2902215/1‑148 (MQ=255)
                    aCTGGCGCGGCAGCTACTGCTTTATTCCCACCGAACGCGATGACATTTTTTCTGTACCATCCCCCGATCGCCTCGAATTGCGCGAAGGCTGGCGAAAACTATTACCCCAGGCGATAGCCGATCCGCTGAACCTACAAAGCTGGAAAgg                                                                                                                          <  1:1849604/148‑1 (MQ=255)
                    aCTGGCGCGGCAGCTACTGCTTTATTCCCACCGAACGCGATGACATTTTTTCTGTACCATCCCCCGATCGCCTCGAATTGCGCGAAGGCTGGCGAAAACTATTACCCCAGGCGATAGCCGATCCGCTGAACCTACAAAGCTGGAAAgg                                                                                                                          <  1:290533/148‑1 (MQ=255)
                        gcgcGGCAGCTACTGCTTTATTCCCACCGAACGCGATGACATTTTTTCTGTACCATCCCCCGATCGCCTCGAATTGCGCGAAGGCTGGCGAAAACTATTACCCCAGGCGATAGCCGATCCGCTGAACCTACAAAGCTGGAAAggcggg                                                                                                                      <  2:3591854/148‑1 (MQ=255)
                          gcgGCAGCTACTGCTTTATTCCCACCGAACGCGATGACATTTTTTCTGTACCATCCCCCGATCGCCTCGAATTGCGCGAAGGCTGGCGAAAACTATTACCCCAGGCGATAGCCGATCCGCAGAACCTACAAAGCTGGAAAggcgggcg                                                                                                                    <  2:1144789/148‑1 (MQ=255)
                           cgGCAGCTACTGCTTTATTCCCACCGAACGCGATGACATTTTTTCTGTACCATCCCCCGATCGCCTCGAATTGCGCGAAGGCTGGCGAAAACTATTACCCCAGGCGATAGCCGATCCGCTGAACCTACAAAGCTGGAAAGGCGGGCGa                                                                                                                   <  2:325038/148‑1 (MQ=255)
                           cgGCAGCTACTGCTTTATTCCCACCGAACGCGATGACATTTTTTCTGTACCATCCCCCGATCGCCTCGAATTGCGCGAAGGCTGGCGAAAACTATTACCCCAGGCGATAGCCGATCCGCTGAACCTACAAAGCTGGAAAGGCGGGCGa                                                                                                                   <  1:3795874/148‑1 (MQ=255)
                            ggCAGCTACTGCTTTATTCCCACCGAACGCGATGACATTTTTTCTGTACCATCCCCCGATCGCCTCGAATTGCGCGAAGGCTGGCGAAAACTATTACCCCAGGCGATAGCCGATCCGCTGAACCTACAAAGCTGGAAAGGCGGGCGAg                                                                                                                  <  2:3716333/148‑1 (MQ=255)
                            ggCAGCTACTGCTTTATTCCCACCGAACGCGATGACATTTTTTCTGTACCATCCCCCGATCGCCTCGAATTGCGCGAAGGCTGGCGAAAACTATTACCCCAGGCGATAGCCGATCCGCTGAACCTACAAAGCTGGAAAGGCGGGCGAg                                                                                                                  <  2:3717105/148‑1 (MQ=255)
                            ggCAGCTACTGCTTTATTCCCACCGAACGCGATGACATTTTTTCTGTACCATCCCCCGATCGCCTCGAATTGCGCGAAGGCTGGCGAAAACTATTACCCCAGGCGATAGCCGATCCGCAGAACCTACAAAGCTGGAAAGGCGGGCGAg                                                                                                                  <  2:542163/148‑1 (MQ=255)
                             gCAGCTACTGCTTTATTCCCACCGAACGCGATGACATTTTTTCTGTACCATCCCCCGATCGCCTCGAATTGCGCGAAGGCTGGCGAAAACTATTACCCCAGGCGATAGCCGATCCGCTGAACCTACAAAGCTGGAAAGGCGGGCGAgg                                                                                                                 <  1:1093799/148‑1 (MQ=255)
                              cAGCTACTGCTTTATTCCCACCGAACGCGATGACATTTTTTCTGTACCATCCCCCGATCGCCTCGAATTGCGCGAAGGCTGGCGAAAACTATTACCCCAGGCGATAGCCGATCCGCAGAACCTACAAAGCTGGAAAGGCGGGCGAggg                                                                                                                >  2:1822223/1‑148 (MQ=255)
                              cAGCTACTGCTTTATTCCCACCGAACGCGATGACATTTTTTCTGTACCATCCCCCGATCGCCTCGAATTGCGCGAAGGCTGGCGAAAACTATTACCCCAGGCGATAGCCGATCCGCAGAACCTACAAAGCTGGAAAGGCGGGCGAggg                                                                                                                <  2:242686/148‑1 (MQ=255)
                                gCTACTGCTTTATTCCCACCGAACGCGATGACATTTTTTCTGTACCATCCCCCGATCGCCTCGAATTGCGCGAAGGCTGGCGAAAACTATTACCCCAGGCGATAGCCGATCCGCTGAACCTACAAAGCTGGAAAGGCGGGCGAGGGCa                                                                                                              <  1:1702275/148‑1 (MQ=255)
                                   aCTGCTTTATTCCCACCGAACGCGATGACATTTTTTCTGTACCATCCCCCGATCGCCTCGAATTGCGCGAAGGCTGGCGAAAACTATTACCCCAGGCGATAGCCGATCCGCTGAACCTACAAAGCTGGAAAGGCGGGCGAGGGCACGc                                                                                                           >  2:2085217/1‑148 (MQ=255)
                                    cTGCTTTATTCCCACCGAACGCGATGACATTTTTTCTGTACCATCCCCCGATCGCCTCGAATTGCGCGAAGGCTGGCGAAAACTATTACCCCAGGCGATAGCCGATCCGCTGAACCTACAAAGCTGGAAAGGCGGGCGAGGGCACGCt                                                                                                          >  2:2138630/1‑148 (MQ=255)
                                    cTGCTTTATTCCCACCGAACGCGATGACATTTTTTCTGTACCATCCCCCGATCGCCTCGAATTGCGCGAAGGCTGGCGAAAACTATTACCCCAGGCGATAGCCGATCCGCAGAACCTACAAAGCTGGAAAGGCGGGCGAGGGCACGCt                                                                                                          <  1:1018891/148‑1 (MQ=255)
                                      gCTTTATTCCCACCGAACGCGATGACATTTTTTCTGTACCATCCCCCGATCGCCTCGAATTGCGCGAAGGCTGGCGAAAACTATTACCCCAGGCGATAGCCGATCCGCTGAACCTACAAAGCTGGAAAGGCGGGCGAGGGCACGCTGt                                                                                                        >  2:495187/1‑148 (MQ=255)
                                         ttATTCCCACCGAACGCGATGACATTTTTTCTGTACCATCCCCCGATCGCCTCGAATTGCGCGAAGGCTGGCGAAAACTATTACCCCAGGCGATAGCCGATCCGCTGAACCTACAAAGCTGGAAAGGCGGGCGAGGGCACGCTGTTTc                                                                                                     <  1:2558489/148‑1 (MQ=255)
                                         ttATTCCCACCGAACGCGATGACATTTTTTCTGTACCATCCCCCGATCGCCTCGAATTGCGCGAAGGCTGGCGAAAACTATTACCCCAGGCGATAGCCGATCCGCTGAACCTACAAAGCTGGAAAGGCGGGCGAGGGCACGCTGTTTc                                                                                                     >  2:2536198/1‑148 (MQ=255)
                                             tCCCACCGAACGCGATGACATTTTTTCTGTACCATCCCCCGATCGCCTCGAATTGCGCGAAGGCTGGCGAAAACTATTACCCCAGGCGATAGCCGATCCGCTGAACCTACAAAGCTGGAAAGGCGGGCGAGGGCACGCTGTTTCTGCa                                                                                                 <  2:3622978/148‑1 (MQ=255)
                                                cACCGAACGCGATGACATTTTTTCTGTACCATCCCCCGATCGCCTCGAATTGCGCGAAGGCTGGCGAAAACTATTACCCCAGGCGATAGCCGATCCGCAGAACCTACAAAGCTGGAAAGGCGGGCGAGGGCACGCTGTTTCTGCACTc                                                                                              >  1:987949/1‑148 (MQ=255)
                                                    gAACGCGATGACATTTTTTCTGTACCATCCCCCGATCGCCTCGAATTGCGCGAAGGCTGGCGAAAACTATTACCCCAGGCGATAGCCGATCCGCTGAACCTACAAAGCTGGAAAGGCGGGCGAGGGCACGCTGTTTCTGCACTCGaaa                                                                                          >  2:3138626/1‑148 (MQ=255)
                                                     aaCGCGATGACATTTTTTCTGTACCATCCCCCGATCGCCTCGAATTGCGCGAAGGCTGGCGAAAACTATTACCCCAGGCGATAGCCGATCCGCAGAACCTACAAAGCTGGAAAGGCGGGCGAGGGCACGCTGTTTCTGCACTCGAAAt                                                                                         >  1:1227756/1‑148 (MQ=255)
                                                     aaCGCGATGACATTTTTTCTGTACCATCCCCCGATCGCCTCGAATTGCGCGAAGGCTGGCGAAAACTATTACCCCAGGCGATAGCCGATCCGCAGAACCTACAAAGCTGGAAAGGCGGGCGAGGGCACGCTGTTTCTGCACTCGAAAt                                                                                         >  1:3132416/1‑148 (MQ=255)
                                                       cgcgATGACATTTTTTCTGTACCATCCCCCGATCGCCTCGAATTGCGCGAAGGCTGGCGAAAACTATTACCCCAGGCGATAGCCGATCCGCTGAACCTACAAAGCTGGAAAGGCGGGCGAGGGCACGCTGTTTCTGCACTCGAAATGc                                                                                       >  1:1626794/1‑148 (MQ=255)
                                                       cgcgATGACATTTTTTCTGTACCATCCCCCGATCGCCTCGAATTGCGCGAAGGCTGGCGAAAACTATTACCCCAGGCGATAGCCGATCCGCTGAACCTACAAAGCTGGAAAGGCGGGCGAGGGCACGCTGTTTCTGCACTCGAAATGc                                                                                       >  1:1619393/1‑148 (MQ=255)
                                                         cgATGACATTTTTTCTGTACCATCCCCCGATCGCCTCGAATTGCGCGAAGGCTGGCGAAAACTATTACCCCAGGCGATAGCCGATCCGCTGAACCTACAAAGCTGGAAAGGCGGGCGAGGGCACGCTGTTTCTGCACTCGAAATGCCg                                                                                     >  1:2875383/1‑148 (MQ=255)
                                                           aTGACATTTTTTCTGTACCATCCCCCGATCGCCTCGAATTGCGCGAAGGCTGGCGAAAACTATTACCCCAGGCGATAGCCGATCCGCTGAACCTACAAAGCTGGAAAGGCGGGCGAGGGCACGCTGTTTCTGCACTCGAAATGCCGCa                                                                                   <  2:1205938/148‑1 (MQ=255)
                                                             gACATTTTTTCTGTACCATCCCCCGATCGCCTCGAATTGCGCGAAGGCTGGCGAAAACTATTACCCCAGGCGATAGCCGATCCGCAGAACCTACAAAGCTGGAAAGGCGGGCGAGGGCACGCTGTTTCTGCACTCGAAATGCCGCAAg                                                                                 >  2:2323076/1‑148 (MQ=255)
                                                              aCATTTTTTCTGTACCATCCCCCGATCGCCTCGAATTGCGCGAGGGCTGGCGAAAACTATTACCCCAGGCGATAGCCGATCCGCTGAACCTACAAAGCTGGAAAGGCGGGCGAGGGCACGCTGTTTCTGCACTCGAAATGCCGCAAgc                                                                                >  1:1845878/1‑148 (MQ=255)
                                                              aCATTTTTTCTGTACCATCCCCCGATCGCCTCGAATTGCGCGAAGGCTGGCGAAAACTATTACCCCAGGCGATAGCCGATCCGCTGAACCTACAAAGCTGGAAAGGCGGGCGAGGGCACGCTGTTTCTGCACTCGAAATGCCGCAAgc                                                                                >  1:1845693/1‑148 (MQ=255)
                                                              aCATTTTTTCTGTACCATCCCCCGATCGCCTCGAATTGCGCGAAGGCTGGCGAAAACTATTACCCCAGGCGATAGCCGATCCGCTGAACCTACAAAGCTGGAAAGGCGGGCGAGGGCACGCTGTTTCTGCACTCGAAATGCCGCAAgc                                                                                >  1:3735171/1‑148 (MQ=255)
                                                               cATTTTTTCTGTACCATCCCCCGATCGCCTCGAATTGCGCGAAGGCTGGCGAAAACTATTACCCCAGGCGATAGCCGATCCGCTGAACCTACAAAGCTGGAAAGGCGGGCGAGGGCACGCTGTTTCTGCACTCGAAATGCCGCAAgcg                                                                               >  2:2649837/1‑148 (MQ=255)
                                                               cATTTTTTCTGTACCATCCCCCGATCGCCTCGAATTGCGCGAAGGCTGGCGAAAACTATTACCCCAGGCGATAGCCGATCCGCAGAACCTACAAAGCTGGAAAGGCGGGCGAGGGCACGCTGTTTCTGCACTCGAAATGCCGCAAgcg                                                                               >  2:721240/1‑148 (MQ=255)
                                                                 ttttttCTGTACCATCCCCCGATCGCCTCGAATTGCGCGAAGGCTGGCGAAAACTATTACCCCAGGCGATAGCCGATCCGCTGAACCTACAAAGCTGGAAAGGCGGGCGAGGGCACGCTGTTTCTGCACTCGAAATGCCGCAAGCGcc                                                                             >  1:99022/1‑148 (MQ=255)
                                                                  tttttCTGTACCATCCCCCGATCGCCTCGAATTGCGCGAAGGCTGGCGAAAACTATTACCCCAGGCGATAGCCGATCCGCTGAACCTACAACGCTGGAAAGGCGGGCGAGGGCACGCTGTTTCTGCACTCGAAATGCCGCAAGCGCct                                                                            <  2:1263976/148‑1 (MQ=255)
                                                                  tttttCTGTACCATCCCCCGATCGCCTCGAATTGCGCGAAGGCTGGCGAAAACTATTACCCCAGGCGATAGCCGATCCGCAGAACCTACAAAGCTGGAAAGGCGGGCGAGGGCACGCTGTTTCTGCACTCGAAATGCCGCAAGCGCct                                                                            <  2:410272/148‑1 (MQ=255)
                                                                        tGTACCATCCCCCGATCGCCTCGAATTGCGCGAAGGCTGGCGAAAACTATTACCCCAGGCGATAGCCGATCCGCTGAACCTACAAAGCTGGAAAGGCGGGCGAGGGCACGCTGTTTCTGCACTCGAAATGCCGCAAGCGCCTCTGCaa                                                                      <  2:3552064/148‑1 (MQ=255)
                                                                        tGTACCATCCCCCGATCGCCTCGAATTGCGCGAAGGCTGGCGAAAACTATTACCCCAGGCGATAGCCGATCCGCTGAACCTACAAAGCTGGAAAGGCGGGCGAGGGCACGCTGTTTCTGCACTCGAAATGCCGCAAGCGCCTCTGCaa                                                                      <  2:3552032/148‑1 (MQ=255)
                                                                        tGTACCATCCCCCGATCGCCTCGAATTGCGCGAAGGCTGGCGAAAACTATTACCCCAGGCGATAGCCGATCCGCAGAACCTCCAAAGCTGGAAAGGCGGGCGAGGGCACGCTGTTTCTGCACTCGAAATGCCGCAAGCGCCTCTGCaa                                                                      <  2:1921561/148‑1 (MQ=255)
                                                                          tACCATCCCCCGATCGCCTCGAATTGCGCGAAGGCTGGCGAAAACTATTACCCCAGGCGATAGCCGATCCGCTGAACCTACAAAGCTGGAAAGGCGGGCGAGGGCACGCTTTTTCTGCACTCGAAATGCCGCAAGCGCCTCTGCAAcc                                                                    <  1:215137/148‑1 (MQ=255)
                                                                            ccATCCCCCGATCGCCTCGAATTGCGCGAAGGCTGGCGAAAACTATTACCCCAGGCGATAGCCGATCCGCAGAACCTACAAAGCTGGAAAGGCGGGCGAGGGCACGCTGTTTCTGCACTCGAAATGCCGCAAGCGCCTCTGCAACCgg                                                                  >  1:3261204/1‑148 (MQ=255)
                                                                              aTCCCCCGATCGCCTCGAATTGCGCGAAGGCTGGCGAAAACTATTACCCCAGGCGATAGCCGATCCGCAGAACCTACAAAGCTGGAAAGGCGGGCGAGGGCACGCTGTTTCTGCACTCGAAATGCCGCAAGCGCCTCTGCAACCggga                                                                >  2:2177620/1‑148 (MQ=255)
                                                                               tCCCCCGATCGCCTCGAATTGCGCGAAGGCTGGCGAAAACTATTACCCCAGGCGATAGCCGATCCGCTGAACCTACAAAGCTGGAAAGGCGGGCGAGGGCACGCTGTTTCTGCACTCGAAATGCCGCAAGCGCCTCTGCAACCgggat                                                               >  1:941716/1‑148 (MQ=255)
                                                                                cccccGATCGCCTCGAATTGCGCGAAGGCTGGCGAAAACTATTACCCCAGGCGATAGCCGATCCGCTGAACCTACAAAGCTGGAAAGGCGGGCGAGGGCACGCTGTTTCTGCACTCGAAATGCCGCAAGCGCCTCTGCAACCgggatg                                                              >  1:2285272/1‑148 (MQ=255)
                                                                                cccccGATCGCCTCGAATTGCGCGAAGGCTGGCGAAAACTATTACCCCAGGCGATAGCCGATCCGCAGAACCTACAAAGCTGGAAAGGCGGGCGAGGGCACGCTGTTTCTGCACTCGAAATGCCGCAAGCGCCTCTGCAACCgggatg                                                              >  2:465519/1‑148 (MQ=255)
                                                                                  cccGATCGCCTCGAATTGCGCGAAGGCTGGCGAAAACTATTACCCCAGGCGATAGCCGATCCGCTGAACCTACAAAGCTGGAAAGGCGGGCGAGGGCACGCTGTTTCTGCACTCGAAATGCCGCAAGCGCCTCTGCAACCgggatggg                                                            <  1:2487175/148‑1 (MQ=255)
                                                                                    cGATCGCCTCGAATTGCGCGAAGGCTGGCGAAAACTATTACCCCAGGCGATAGCCGATCCGCAGAACCTACAAAGCTGGAAAGGCGGGCGAGGGCACGCTGTTTCTGCACTCGAAATGCCGCAAGCGCCTCTGCAACCgggatgggat                                                          >  2:1189523/1‑148 (MQ=255)
                                                                                     gATCGCCTCGAATTGCGCGAAGGCTGGCGAAAACTATTACCCCAGGCGATAGCCGATCCGCTGAACCTACAAAGCTGGAAAGGCGGGCGAGGGCACGCTGTTTCTGCACTCGAAATGCCGCAAGCGCCTCTGCAACCGGGATGGGAtt                                                         >  1:1882923/1‑148 (MQ=255)
                                                                                      aTCGCCTCGAATTGCGCGAAGGCTGGCGAAAACTATTACCCCAGGCGATAGCCGATCCGCTGAACCTACAAAGCTGGAAAGGCGGGCGAGGGCACGCTGTTTCTGCACTCGAAATGCCGCAAGCGCCTCTGCAACCGGGATGGGATTg                                                        >  1:579899/1‑148 (MQ=255)
                                                                                       tCGCCTCGAATTGCGCGAAGGCTGGCGAAAACTATTACCCCAGGCGATAGCCGATCCGCTGAACCTACAAAGCTGGAAAGGCGGGCGAGGGCACGCTGTTTCTGCACTCGAAATGCCGCAAGCGCCTCTGCAACCGGGATGGGATTGt                                                       >  1:3210794/1‑148 (MQ=255)
                                                                                       tCGCCTCGAATTGCGCGAAGGCTGGCGAAAACTATTACCCCAGGCGATAGCCGATCCGCTGAACCTACAAAGCTGGAAAGGCGGGCGAGGGCACGCTGTTTCTGCACTCGAAATGCCGCAAGCGCCTCTGCAACCGGGATGGGATTGt                                                       >  2:833277/1‑148 (MQ=255)
                                                                                       tCGCCTCGAATTGCGCGAAGGCTGGCGAAAACTATTACCCCAGGCGATAGCCGATCCGCAGAACCTACAAAGCTGGAAAGGCGGGCGAGGGCACGCTGTTTCTGCACTCGAAATGCCGCAAGCGCCTCTGCAACCGGGATGGGATTGt                                                       >  1:2721369/1‑148 (MQ=255)
                                                                                        cGCCTCGAATTGCGCGAAGGCTGGCGAAAACTATTACCCCAGGCGATAGCCGATCCGCAGAACCTACAAAGCTGGAAAGGCGGGCGAGGGCACGCTGTTTCTGCACTCGAAATGCCGCAAGCGCCTCTGCAACCGGGATGGGATTGTc                                                      <  1:2232076/148‑1 (MQ=255)
                                                                                        cGCCTCGAATTGCGCGAAGGCTGGCGAAAACTATTACCCCAGGCGATAGCCGATCCGCAGAACCTACAAAGCTGGAAAGGCGGGCGAGGGCACGCTGTTTCTGCACTCGAAATGCCGCAAGCGCCTCTGCAACCGGGAGGGGATTGTc                                                      >  1:2861190/1‑148 (MQ=255)
                                                                                        cGCCTCGAATGGCGCGAAGGCTGGCTAAAACTATTCCCCCCGGCGATAGCCGATCCGCACAACCTACAAAGCTGGAATGGCGGGCGTGGGCACGCTGTTTCTGCACTCGAAATGCCGCACGCGCCTCTGCAACCGGGATGGGATTGTc                                                      <  1:50780/148‑1 (MQ=255)
                                                                                            tCGAATTGCGCGAAGGCTGGCGAAAACTATTACCCCAGGCGATAGCCGATCCGCTGAACCTACAAAGCTGGAAAGGCGGGCGAGGGCACGCTGTTTCTGCACTCGAAATGCCGCAAGCGCCTCTGCAACCGGGATGGGATTGTCCGCa                                                  <  2:1089081/148‑1 (MQ=255)
                                                                                             cGAATTGCGCGAAGGCTGGCGAAAACTATTACCCCAGGCGATAGCCGATCCGCTGAACCTACAAAGCTGGAAAGGCGGGCGAGGGCACGCTGTTTCTGCACTCGAAATGCCGCAAGCGCCTCTGCAACCGGGATGGGATTGTCCGCaa                                                 >  2:1666692/1‑148 (MQ=255)
                                                                                                aTTGCGCGAAGGCTGGCGAAAACTATTACCCCAGGCGATAGCCGATCCGCTGAACCTACAAAGCTGGAAAGGCGGTCGAGGGCACGCTGTTTCTGCACTCGAAATGCCGCAAGCGCCTCTGCAACCGGGATGGGATTGTCCGCAAgcg                                              <  2:1862209/148‑1 (MQ=255)
                                                                                                  tGCGCGAAGGCTGGCGAAAACTATTACCCCAGGCGATAGCCGATCCGCTGAACCTACAAAGCTGGAAAGGCGGGCGAGGGCACGCTGTTTCTGCACTCGAAATGCCGCAAGCGCCTCTGCAACCGGGATGGGATTGTCCGCAAGCGcc                                            >  2:1713255/1‑148 (MQ=255)
                                                                                                  tGCGCGAAGGCTGGCGAAAACTATTACCCCAGGCGATAGCCGATCCGCTGAACCTACAAAGCTGGAAAGGCGGGCGAGGGCACGCTGTTTCTGCACTCGAAATGCCGCAAGCGCCTCTGCAACCGGGAGGGGATTGTCCGCAAGCGcc                                            >  1:667980/1‑148 (MQ=255)
                                                                                                   gcgcGAAGGCTGGCGAAAACTATTACCCCAGGCGATAGCCGATCCGCTGAACCTACAAAGCTGGAAAGGCGGGCGAGGGCACGCTGTTTCTGCACTCGAAATGCCGCAAGCGCCTCTGCAACCGGGATGGGATTGTCCGCAAGCGCCa                                           <  2:272153/148‑1 (MQ=255)
                                                                                                          ggCTGGCGAAAACTATTACCCCAGGCGATAGCCGATCCGCTGAACCTACAAAGCTGGAAAGGCGGGCGAGGGCACGCTGTTTCTGCACTCGAAATGCCGCAAGCGCCTCTGCAACCGGGATGGGATTGTCCGCAAGCGCCAGAAATAc                                    <  2:3440253/148‑1 (MQ=255)
                                                                                                          ggCTGGCGAAAACTATTACCCCAGGCGATAGCCGATCCGCTGAACCTACAAAGCTGGAAAGGCGGGCGAGGGCACGCTGTTTCTGCACTCGAAATGCCGCAAGCGCCTCTGCAACCGGGATGGGATTGTCCGCAAGCGCCAGAAATAc                                    >  1:426711/1‑148 (MQ=255)
                                                                                                          ggCTGGCGAAAACTATTACCCCAGGCGATAGCCGATCCGCAGAACCTACAAAGCTGGAAAGGCGGGCGAGGGCACGCTGTTTCTGCACTCGAAATGCCGCAAGCGCCTCTGCAACCGGGATGGGATTGTCCGCAAGCGCCAGAAATAc                                    <  2:1227756/148‑1 (MQ=255)
                                                                                                           gCTGGCGAAAACTATTACCCCAGGCGATAGCCGATCCGCTGAACCTACAAAGCTGGAAAGGCGGGCGAGGGCACGATGTTTCTGCACTCGAAATGCCGCAAGCGCCTCTGCAACCGGGATGGGATTGTCCGCAAGCGCCAGAAATAcc                                   >  2:949326/1‑148 (MQ=255)
                                                                                                            cTGGCGAAAACTATTCCCCCCGGCGATAGCCGATCCGCTGAACCTACAAAGGTGGAAAGGCGGGCGAGGGCACGCTGTTTCTGCACTCGAAAGGCCGCAAGCGCCTCTGCAACCGGGATGGGATTGTCCGCAAGCGCCAGAAATACCt                                  <  1:348198/148‑1 (MQ=255)
                                                                                                             tGGCGAAAACTATTACCCCAGGCGATAGCCGATCCGATGAACCTACAAAGCTGGATAGGGGGGCGAGGGCACGCTTTTTCTGCACTCTAAATGCCGAAAGCGCCTCTGCAACCGGGATGGGATTGTCCGCAAGCGCCAGAAATACCTg                                 <  1:2536494/148‑1 (MQ=255)
                                                                                                               gCGAAAACTATTACCCCAGGCGATAGCCGATCCGCTGAACCTACAAAGCTGGAAAGGCGGGCGAGGGCACGCTGTTTCTGCACTCGAAATGCCGCAAGCGCCTCTGCAACCGGGATGGGATTGTCCGCAAGCGCCAGAAATACCTGcc                               >  1:116448/1‑148 (MQ=255)
                                                                                                                cGAAAACTATTACCCCAGGCGATAGCCGATCCGCTGAACCTACAAAGCTGGAAAGGCGGGCGAGGGCACGCTGTTTCTGCACTCGAAATGCCGCAAGCGCCTCTGCAACCGGGATGGGATTGTCCGCAAGCGCCAGAAATACCTGCCa                              >  2:1671903/1‑148 (MQ=255)
                                                                                                                  aaaaCTATTACCCCAGGCGATAGCCGATCCGCTGAACCTACAAAGCTGGAAAGGCGGGCGAGGGCACGCTGTTTCTGCACTCGAAATGCCGCAAGCGCCTCTGCAACCGGGATGGGATTGTCCGCAAGCGCCAGAAATACCTGCCaaa                            >  2:771932/1‑148 (MQ=255)
                                                                                                                   aaaCTATTACCCCAGGCGATAGCCGATCCGCAGAACCTACACAGCTGGAAAGGCGGGCGAGGGCACGCTGTTTCTGCACTCGAAATGCCGCAAGCGCCTCTGCAACCGGGATGGGATTGTCCGCAAGCGCCAGAAATACCTGCCAAAg                           <  2:2516447/148‑1 (MQ=255)
                                                                                                                    aaCTATTACCCCAGGCGATAGCCGATCCGCTGAACCTACAAAGCTGGAAAGGCGGGCGAGGGCACGCTGTTTCTGCACTCGAAATGCCGCAAGCGCCTCTGCAACCGGGATGGGATTGTCCGCAAGCGCCAGAAATACCTGCAAAAGa                          <  1:2206956/148‑1 (MQ=255)
                                                                                                                    aaCTATTACCCCAGGCGATAGCCGATCCGCAGAACCTACAAAGCTGGAAAGGCGGGCGAGGGCACGCTGTTTCTGCACTCGAAATGCCGCAAGCGCCTCTGCAACCGGGATGGGATTGTCCGCAAGCGCCAGAAATACCTGCCAAAGa                          >  1:88992/1‑148 (MQ=255)
                                                                                                                     aCTATTACCCCAGGCGATAGCCGATCCGCAGAACCTACAAAGCTGGAAAGGCGGGCGAGGGCACGCTGTTTCTGCACTCGAAATGCCGCAAGCGCCTCTGCAACCGGGATGGGATTGTCCGCAAGCGCCAGAAATACCTGCCAAAGaa                         <  2:3082296/148‑1 (MQ=255)
                                                                                                                     aCTATTACCCCAGGCGATAGCCGATCCGCAGAACCTACAAAGCTGGAAAGGCGGGCGAGGGCACGCTGTTTCTGCACTCGAAAGGCCGCAAGCGCCTCTGCAACCGGGATGGGATTGTCCGCAAGCGCCAGAAATACCTGCCAAAGaa                         <  2:3082363/148‑1 (MQ=255)
                                                                                                                       tATTCCCCCAGGCGATAGCCGATCCGCTGAACCTACAAAGCTGTAATGGCGGGCGAGGGCACGCTGTTTCTGCACTCGAAATGCCGCAAGCGCCTCTGCAACCGGGATGGGATTGTCCGCAAGCGCCAGAAATACCTGCCAAAGAAAt                       <  1:3355974/148‑1 (MQ=255)
                                                                                                                        aTTACCCCAGGCGATAGCCGATCCGCTGAACCTACAAAGCTGGAAAGGCGGGCGAGGGCACGCTGTTTCTGCACTCGAAATGCCGCAAGCGCCTCTGCAACCGGGATGGGATTGTCCGCAAGCGCCAGAAATACCTGCCAAAGAAAtt                      >  1:3335571/1‑148 (MQ=255)
                                                                                                                           aCCCCAGGCGATAGCCGATCCGCTGAACCTACAAAGCTGGAAAGGCGGGCGAGGGCACGCTGTTTCTGCACTCGAAATGCCGCAAGCGCCTCTGCAACCGGGATGGGATTGTCCGCAAGCGCCAGAAATACCTGCCAAAGAAATTATc                   <  2:2165763/148‑1 (MQ=255)
                                                                                                                               cAGGCGATAGCCGATCCGCTGAACCTACAAAGCTGGAAAGGCGGGCGAGGGCACGCTGTTTCTGCACTCGAAATGCCGCAAGCGCCTCTGCAACCGGGATGGGATTGTCCGCAAGCGCCAGAAATACCTGCCAAAGAAATTATCTGGa               <  2:2585312/148‑1 (MQ=255)
                                                                                                                               cAGGCGATAGCCGATCCGCTGAACCTACAAAGCTGGAAAGGCGGGCGAGGGCACGCTGTTTCTGCACTCGAAATGCCGCAAGCGCCTCTGCAACCGGGATGGGATTGTCCGCAAGCGCCAGAAATACCTGCCAAAGAAATTATCTGGa               <  2:2584490/148‑1 (MQ=255)
                                                                                                                                aGGCGATAGCCGATCCGCTGAACCTACAAAGCTGGAAAGGCGGGCGAGGGCACGCTGTTTCTGCACTCGAAATGCCGCAAGCGCCTCTGCAACCGGGATGGGATTGTCCGCAAGCGCCAGAAATACCTGCCAAAGAAATTATCTGGaa              <  1:176646/148‑1 (MQ=255)
                                                                                                                                 tgCGATTGCCGATCCGCTGAACCTACAAAGCTGGTATGGGTGGTGAGGGAACGTTGTTTCTGCACTCGATATGCCGCAAGCGCCTCTGCATCCGGGATGGGATTGTACGCAAGCGCCAGAAATACCTGCCAAAGAAATTATCTGGaaa             <  1:2486131/147‑1 (MQ=255)
                                                                                                                                 ggCGATAGCCGATCCGCTGAACCTACAAAGCTGGAAAGGCGGGCGAGGGCACGCTGTTTCTGCACTCGAAATGCCGCAAGCGCCTCTGCAACCGGGATGGGATTGTCCGCAAGCGCCAGAAATACCTGCCAAAGAAATTATCTGGaaa             <  2:29781/148‑1 (MQ=255)
                                                                                                                                  gCGATAGCCGATCCGCTGAACCTACAAAGCTGGAAAGGCGGGCGAGGGCACGCTGTTTCTGCACTCGAAATGCCGCAAGCGCCTCTGCAACCGGGATGGGATTGTCCGCAAGCGCCAGAAATACCTGCCAAAGAAATTATCTGGaaaa            >  1:3058659/1‑148 (MQ=255)
                                                                                                                                   cGATAGCCGATCCGCTGAACCTACAAAGCTGGAAAGGCGGGCGAGGGCACGCTGTTTCTGCACTCTAAATGCCGCAAGCGCCTCTGCAACCGGGGTGGGATTGTCCGCAAGCGCCAGAAATACCTGCCAAAGAAATTATCTGGAAAAg           <  1:2317609/148‑1 (MQ=255)
                                                                                                                                     aTAGCCGATCCGCTGAACCTACAATGCTGGAATGGCGGGCGAGGGCACGCTGTTTATGCACTCGAAATGCCGCAAGCGCCTCTGCAACCGGGATGGGATTGTCCGCAAGCGCCAGAAATACCTGCCAAAGAAATTATCTGGAAAAGTg         <  1:835458/148‑1 (MQ=255)
                                                                                                                                       tgCCGATCCGCAGAACCTACAAAGCTGGAATGGGGGGCGCGGGCACGCTTTTTCTGCACTCGAAAGCCCGCAAGCGCCTCTGCAACCGGGATGGGATTGTCCGCAAGCGCCAGAAATACCTGCCAAAGAAATTATCTGGAAAAGTGaa       <  1:465519/147‑1 (MQ=255)
                                                                                                                                       aGCCGATCCGCAGAACCTACAAAGCTGGAAAGGCGGGCGAGGGCACGCTGTTTCTGCACTCGAAATGCCGCAAGCGCCTCTGCAACCGGGATGGGATTGTCCGCAAGCGCCAGAAATACCTGCCAAAGAAATTATCTGGAAAAGTGaa       >  2:1561881/1‑148 (MQ=255)
                                                                                                                                          cGATCCGCTGAACCTACAAAGCTGGAAAGGCGGGCGAGGGCACGCTGTTTCTGCACTCGAAATGCCGCAAGCGCCTCTGCAACCGGGATGGGATTGTCCGCAAGCGCCAGAAATACCTGCCAAAGAAATTATCTGGAAAAGTGAACgg    <  1:490373/148‑1 (MQ=255)
                                                                                                                                           gATCCGCTGAACCTACAAAGCTGGAAAGGCGGGCGAGGGCACGCTGTTTCTGCACTCGAAATGCCGCAAGCGCCTCTGCAACCGGGATGGGATTGTCCGCAAGCGCCAGAAATACCTGCCAAAGAAATTATCTGGAAAAGTGAACGGt   <  2:1845698/148‑1 (MQ=255)
                                                                                                                                           gATCCGCAGAACCTACAAAGCTGGAAAGGCGGGCGAGGGCACGCTGTTTCTGCACTCGAAATGCCGCAAGCGCCTCTGCAACCGGGATGGGATTGTCCGCAAGCGCCAGAAATACCTGCCAAAGAAATTATCTGGAAAAGTGAACGGt   <  2:3392133/148‑1 (MQ=255)
                                                                                                                                            aTCCGCTGAACCTACAAAGCTGGAAAGGCGGGCGAGGGCACGCTGTTTCTGCACTCGAAATGCCGCAAGCGCCTCTGCAACCGGGATGGGATTGTCCGCAAGCGCCAGAAATACCTGCCAAAGAAATTATCTGGAAAAGTGAACGGtt  >  1:1539942/1‑148 (MQ=255)
                                                                                                                                            aTCCGCTGAACCTACAAAGCTGGAAAGGCGGGCGAGGGCACGCTGTTTCTGCACTCGAAATGCCGCAAGCGCCTCTGCAACCGGGATGGGATTGTCCGCAAGCGCCAGAAATACCTGCCAAAGAAATTATCTGGAAAAGTGAACGGtt  >  1:3518791/1‑148 (MQ=255)
                                                                                                                                                  |                                                                                                                                             
GACGACACAACTCAATGCCAACTGGCGCGGCAGCTACTGCTTTATTCCCACCGAACGCGATGACATTTTTTCTGTACCATCCCCCGATCGCCTCGAATTGCGCGAAGGCTGGCGAAAACTATTACCCCAGGCGATAGCCGATCCGCTGAACCTACAAAGCTGGAAAGGCGGGCGAGGGCACGCTGTTTCTGCACTCGAAATGCCGCAAGCGCCTCTGCAACCGGGATGGGATTGTCCGCAAGCGCCAGAAATACCTGCCAAAGAAATTATCTGGAAAAGTGAACGGTT  >  CP009273/608429‑608716

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