Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,328,078 A→T 6.3% Y121N (TAC→AAC)  yfaT ← DUF1175 family protein, putative host defense protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,328,0780AT 6.3% 1312.6 / 62.2 363Y121N (TAC→AAC) yfaTDUF1175 family protein, putative host defense protein
Reads supporting (aligned to +/- strand):  ref base A (150/190);  new base T (14/9);  total (164/199)
Fisher's exact test for biased strand distribution p-value = 1.33e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.83e-01

Only 100 of 368 total aligned reads displayed.

CCGTTGTCAGTTTTCGTGGCGCTTCCGGTGTGGTAGATGACGTAACGCCCCATCCAGACCATTAAGTGCTGGGCATCGCCCTGATCGAAAAAAATCATATCGCCAGGCAGCGCCTGGTTTATGTCCTGGCCAATAAACTGGCTGTTGTACTGAATCAAATTAATCGCGGTCACATAGGGGCCGGTTTTCCCGTTCCCCTGATTCCAGTTTTGCGCCAGTTGACGCTGTTCAGGTGTTAGCGTCATCTCTGGCGGCAAATACTGGCTAGATAAACCGTTACTTTTTAGCCATTTACT  >  CP009273/2327930‑2328225
                                                                                                                                                    |                                                                                                                                                   
ccGTTGTCAGTATTCGTGGCGCTTCCGGAGTGGTAGATGACGTAACGCCCCATCCATACCATTAAGTGCTGGGCATCGCCCTGATCGAAAAAAATCATATCGACAGGCAGCGCCTGGTTTATGTCCTGGCCAATAAACTGGCTGTTGTa                                                                                                                                                     >  2:3445172/1‑149 (MQ=255)
     gTCAGTTTTCGTGGCGCTTCCGGTGTGGTAGATGACGTAACGCCCCATCCAGACCATTAAGTGCTGGGCATCGCCCTGATCGAAAAAAATCATATCGCCAGGCAGCGCCTGGTTTATGTCCTGGCCAATAAACTGGCTGTTGTACTGaa                                                                                                                                                <  1:4992100/149‑1 (MQ=255)
      tCAGTTTTCGTGGCGCTTCCGGTGTGGTAGATGACGTAACGCCCCATCCAGACCATTAAGTGCTGGGCATCGCCCTGATCGAAAAAAATCATATCGCCAGGCAGCGCCTGGTTTATGTCCTGGCCAATAAACTGGCTGTTGTACTGAAt                                                                                                                                               >  1:2567677/1‑149 (MQ=255)
       cAGTTTTCGTGGCGCTTCCGGTGTGGTAGATGACGTAACGCCCCATCCAGACCATTAAGTGCTGGGCATCGCCCTGATCGAAAAAAATCATATCGCCAGGCAGCGCCTGGTTTATGTCCTGGCCAATAAACTGGCTGTTGTACTGAATc                                                                                                                                              <  1:932537/149‑1 (MQ=255)
       cAGTTTTCGTGGCGCTTCCGGTGTGGTAGATGACGTAACGCCCCATCCAGACCATTAAGTGCTGGGCATCGCCCTGATCGAAAAAAATCATATCGCCAGGCAGCGCCTGGTTTATGTCCTGGCCAATAAACTGGCTGTTGTACTGAATc                                                                                                                                              <  1:2300696/149‑1 (MQ=255)
              cGTGGCGCTTCCGGTGTGGTAGATGACGTAACGCCCCATCCAGACCATTAAGTGCTGGGCATCGCCCTGATCGAAAAAAATCATATCGCCAGGCAGCGCCTGGTTTATGTCCTGGCCAATAAACTGGCTGTTGTACTGAATCAAATTaa                                                                                                                                       <  1:1755426/149‑1 (MQ=255)
               gTGGCGCTTCCGGTGTGGTAGATGACGTAACGCCCCATCCAGACCATTAAGTGCTGGGCATCGCCCTGATCGAAAAAAATCATATCGCCAGGCAGCGCCTGGTTTATGTCCTGGCCAATAAACTGGCTGTTGTACTGAATCAAATTAAt                                                                                                                                      <  1:4998748/149‑1 (MQ=255)
                tGGCGCTTCCGGTGTGGTAGATGACGTAACGCCCCATCCAGACCATTAAGTGCTGGGCATCGCCCTGATCGAAAAAAATCATATCGCCAGGCAGCGCCTGGTTTATGTCCTGGCCAATAAACTGGCTGTTGTACTGAATCAAATTAATc                                                                                                                                     <  2:294343/149‑1 (MQ=255)
                tGGCGCTTCCGGTGTGGTAGATGACGTAACGCCCCATCCAGACCATTAAGTGCTGGGCATCGCCCTGATCGAAAAAAATCATATCGCCAGGCAGCGCCTGGTTTATGTCCTGGCCAATAAACTGGCTGTTGTACTGAATCAAATTAATc                                                                                                                                     >  1:5053144/1‑149 (MQ=255)
                    gcTTCCGGTGTGGTAGATGACGTAACGCCCCATCCAGACCATTAAGTGCTGGGCATCGCCCTGATCGAAAAAAATCATATCGCCAGGCAGCGCCTGGTTTATGTCCTGGCCAATAAACTGGCTGTTGTACTGAATCAAATTAATCGCgg                                                                                                                                 <  2:493853/149‑1 (MQ=255)
                    gcTTCCGGTGTGGTAGATGACGTAACGCCCCATCCAGACCATTAAGTGCTGGGCATCGCCCTGATCGAAAAAAATCATATCGCCAGGCAGCGCCTGGTTTATGTCCTGGCCAATAAACTGGCTGTTGTACTGAATCAAATTAATCGCgg                                                                                                                                 <  1:2442100/149‑1 (MQ=255)
                     cTTCCTGTGTGGTAGATGACGTAACGACCCATCCAGACCATTAAGTGCTGGGCATCGCCCTGATCGAAAAAAATCATATCGCCAGGCAGCGCCTGGTTTATGTCCTGGCCAATAAACTGGCTGTTGTTCTGAATCAAATTAATCGCGGt                                                                                                                                <  2:2825447/149‑1 (MQ=255)
                      ttCCGGTGTGGTAGATGACGTAACGCCCCATCCAGACCATTAAGTGCTGGGCATCGCCCTGATCGAAAAAAATCATATCGCCAGGCAGCGCCTGGTTTATGTCCTGGCCAATAAACTGGCTGTTGTACTGAATCAAATTAATCGCGGTc                                                                                                                               >  1:2088710/1‑149 (MQ=255)
                        ccGGTGTGGTAGATGACGTAACGCCCCATCCAGACCATTAAGTGCTGGGCATCGCCCTGATCGAAAAAAATCATATCGCCAGGCAGCGCCTGGTTTATGTCCTGGCCAATAAACTGGCTGTTGTACTGAATCAAATTAATCGCGGTcac                                                                                                                             <  1:5141308/149‑1 (MQ=255)
                          ggTGTGGTAGATGACGTAACGCCCCATCCAGACCATTAAGTGCTGGGCATCGCCCTGATCGAAAAAAATCATATCGCCAGGCAGCGCCTGGTTTATGTCCTGGCCAATAAACTGGCTGTTGTACTGAATCAAATTAATCGCGGTCACAt                                                                                                                           >  2:3172337/1‑149 (MQ=255)
                           gtgtGGTAGATGACGTAACGCCCCATCCAGACCATTAAGTGCTGGGCATCGCCCTGATCGAAAAAAATCATATCGCCAGGCAGCGCCTGGTTTATGTCCTGGCCAATAAACTGGCTGTTGTACTGAATCAAATTAATCGCGGTCACATa                                                                                                                          <  1:3998784/149‑1 (MQ=255)
                                 tAGATGACGTAACGCCCCATCCAGACCATTAAGTGCTGGGCATCGCCCTGATCGAAAAAAATCATATCGCCAGGCAGCGCCTGGTTTATGTCCTGGCCAATAAACTGGCTGTTGTACTGAATCAAATTAATCGCGGTCACATAGGGGcc                                                                                                                    <  2:2983381/149‑1 (MQ=255)
                                 tAGATGACGTAACGCCCCATCCAGACCATTAAGTGCTGGGCATCGCCCTGATCGAAAAAAATCATATCGCCAGGCAGCGCCTGGTTTATGTCCTGGCCAATAAACTGGCTGTTGTACTGAATCAAATTAATCGCGGTCACATAGGGGcc                                                                                                                    >  2:2185858/1‑149 (MQ=255)
                                       aCGTAACGCCCCATCCAGACCATTAAGTGCTGGGCATCGCCCTGATCGAAAAAAATCATATCGCCAGGCAGCGCCTGGTTTATGTCCTGGCCAATAAACTGGCTGTTGTACTGAATCAAATTAATCGCGGTCACATAGGGGCCGGtttt                                                                                                              >  2:1872165/1‑149 (MQ=255)
                                              gCCCCATCCAGACCATTAAGTGCTGGGCATCGCCCTGATCGAAAAAAATCATATCGCCAGGCAGCGCCTGGTTTATGTCCTGGCCAATAAACTGGCTGTTGTACTGAATCAAATTAATCGCGGTCACATAGGGGCCGGTTTTCCCGTTc                                                                                                       >  1:103507/1‑149 (MQ=255)
                                                 ccATCCAGACCATTAAGTGCTGGGCATCGCCCTGATCGAAAAAAATCATATCGCCAGGCAGCGCCTGGTTTATGTCCTGGCCAATAAACTGGCTGTTGTACTGAATCAAATTAATCGCGGTCACATAGGGGCCGGTTTTCCCGTTcccc                                                                                                    >  2:3898348/1‑149 (MQ=255)
                                                  cATCCAGACCATTAAGTGCTGGGCATCGCCCTGATCGAAAAAAATCATATCGCCAGGCAGCGCCTGGTTTATGTCCTGGCCAATAAACTGGCTGTTGTACTGAATCAAATTAATCGCGGTCACATAGGGGCCGGTTTTCCCGTTCCCCt                                                                                                   <  1:892269/149‑1 (MQ=255)
                                                   aTCCAGACCATTAAGTGCTGGGCATCGCCCTGATCGAAAAAAATCATATCGCCAGGCAGCGCCTGGTTTATGTCCTGGCCAATAAACTGGCTGTTGTACTGAATCAAATTAATCGCGGTCACATAGGGGCCGGTTTTCCCGTTCCCCTg                                                                                                  >  1:1701870/1‑149 (MQ=255)
                                                   aTCCAGACCATTAAGTGCTGGGCATCGCCCTGATCGAAAAAAATCATATCGCCAGGCAGCGCCTGGTTTATGTCCTGGCCAATAAACTGGCTGTTGTACTGAATCAAATTAATCGCGGTCACATAGGGGCCGGTTTTCCCGTTCCCCTg                                                                                                  <  1:3066125/149‑1 (MQ=255)
                                                    tCCAGACCATTAAGTGCTGGGCATCGCCCTGATCGAAAAAAATCATATCGCCAGGCAGCGCCTGGTTTATGTCCTGGCCAATAAACTGGCTGTTGTACTGAATCAAATTAATCGCGGTCACATAGGGGCCGGTTTTCCCGTTCCCCTGa                                                                                                 >  2:4284572/1‑149 (MQ=255)
                                                       aGACCATTAAGTGCTGGGCATCGCCCTGATCGAAAAAAATCATATCGCCAGGCAGCGCCTGGTGTATGTCCTGGCCAATAAACTGGCTGTTGTACTGAATCAAATTAATCGCGGTCACATAGGGGCCGGTTTTCCCGTTCCCCTGAGTc                                                                                              >  2:3306302/1‑149 (MQ=255)
                                                        gACCATTAAGTGCTGGGCATCGCCCTGATCGAAAAAAATCATATCGCCAGGCAGCGCCTGGTTTATGTCCTGGCCAATAAACTGGCTGTTGTACTGAATCAAATTAATCGCGGTCACATAGGGGCCGGTTTTCCCGTTCCCCTGATTcc                                                                                             >  1:692099/1‑149 (MQ=255)
                                                        gACCATTAAGTGCTGGGCATCGCCCTGATCGAAAAAAATCATATCGCCAGGCAGCGCCTGGTTTATGTCCTGGCCAATAAACTGGCTGTTGTACTGAATCAAATTAATCGCGGTCACATAGGGGCCGGTTTTCCCGTTCCCCTGATTcc                                                                                             >  2:5157932/1‑149 (MQ=255)
                                                           cATTAAGTGCTGGGCATCGCCCTGATCGAAAAAAATCATATCGCCAGGCAGCGCCTGGTTTATGTCCTGGCCAATAAACTGGCTGTTGTACTGAATCAAATTAATCGCGGTCACATAGGGGCCGGTTTTCCCGTTCCCCTGATTCCAGt                                                                                          <  1:1314584/149‑1 (MQ=255)
                                                           cATTAAGTGCTGGGCATCGCCCTGATCGAAAAAAATCATATCGCCAGGCAGCGCCTGGTTTATGTCCTGGCCAATAAACTGGCTGTTGTACTGAATCAAATTAATCGCGGTCACATAGGGGCCGGTTTTCCCGTTCCCCTGATTCCAGt                                                                                          <  1:4624047/149‑1 (MQ=255)
                                                            aTTAAGTGCTGGGCATCGCCCTGATCGAAAAAAATCATATCGCCAGGCAGCGCCTGGTTTATGTCCTGGCCAATAAACTGGCTGTTGTACTGAATCAAATTAATCGCGGTCACATAGGGGCCGGTTTTCCCGTTCCCCTGATTCCAGtt                                                                                         >  2:3505578/1‑149 (MQ=255)
                                                            aTTAAGTGCTGGGCATCGCCCTGATCGAAAAAAATCATATCGCCAGGCAGCGCCTGGTTTATGTCCTGGCCAATAAACTGGCTGTTGTACTGAATCAAATTAATCGCGGTCACATAGGGGCCGGTTTTCCCGTTCCCCTGATTCCAGtt                                                                                         >  2:1408117/1‑149 (MQ=255)
                                                               aaGTGCTGGGCATCGCCCTGATCGAAAAAAATCATATCGCCAGGCAGCGCCTGGTTTATGTCCTGGCCAATAAACTGGCTGTTGTACTGAATCAAATTAATCGCGGTCACATAGGGGCCGGTTTTCCCGTTCCCCTGATTCCAGTTTTg                                                                                      <  2:3758897/149‑1 (MQ=255)
                                                                 gTGCTGGGCATCGCCCTGATCGAAAAAAATCATATCGCCAGGCAGCGCCTGGTTTATGTCCTGGCCAATAAACTGGCTGTTGTACTGAATCAAATTAATCGCGGTCACATAGGGGCCGGTTTTCCCGTTCCCCTGATTCCAGTTTTgcg                                                                                    >  2:4681440/1‑149 (MQ=255)
                                                                   gCTGGGCATCGCCCTGATCGAAAAAAATCATATCGCCAGGCAGCGCCTGGTTTATGTCCTGGCCAATAAACTGGCTGTTGTACTGAATCAAATTAATCGCGGTCACATAGGGGCCGGTTTTCCCGTTCCCCTGATTCCAGTTTTGCGcc                                                                                  <  1:2938077/149‑1 (MQ=255)
                                                                      ggtcATCGCCCTGATCGAAAAAAATCATATCGCCAGGCAGCGCCTGGTTTATGTCCTGGCCAATAAACTGGCTGTTGTACTGAATCAAATAAATCGCGGTCACATAGGGGCCGGTTTTCCCGTTCCCCTGATTCAAGTTTTGCGCCAGt                                                                               >  2:2595088/4‑149 (MQ=255)
                                                                      gggCATCGCCCTGATCGAAAAAAATCATATCGCCAGGCAGCGCCTGGTTTATGTCCTGGCCAATAAACTGGCTGTTGTACTGAATCAAATTAATCGCGGTCACATAGGGGCCGGTTTTCCCGTTCCCCTGATTCCAGTTTTGCGCCAGt                                                                               >  2:1320753/1‑149 (MQ=255)
                                                                        gCATCGCCCTGATCGAAAAAAATCATATCGCCAGGCAGCGCCTGGTTTATGTCCTGGCCAATAAACTGGCTGTTGTACTGAATCAAATTAATCGCGGTCACATAGGGGCCGGTTTTCCCGTTCCCCTGATTCCAGTTTTGCGCCAGTTg                                                                             >  1:574286/1‑149 (MQ=255)
                                                                         cATCGCCCTGATCGAAAAAAATCATATCGCCAGGCAGCGCCTGGTTTATGTCCTGGCCAATAAACTGGCTGTTGTACTGAATCAAATTAATCGCGGTCACATAGGGGCCGGTTTTCCCGTTCCCCTGATTCCAGTTTTGCGCCAGTTGa                                                                            <  1:4845451/149‑1 (MQ=255)
                                                                         cATCGCCCTGATCGAAAAAAATCATATCGCCAGGCAGCGCCTGGTTTATGTCCTGGCCAATAAACTGGCTGTTGTACTGAATCAAATTAATCGCGGTCACATAGGGGCCGGTTTTCCCGTTCCCCTGATTCCAGTTTTGCGCCAGTTGa                                                                            <  2:3109403/149‑1 (MQ=255)
                                                                            cGCCCTGATCGAAAAAAATCATATCGCCAGGCAGCGCCTGGTTTATGTCCTGGCCAATAAACTGGCTGTTGTACTGAATCAAATTAATCGCGGTCACATAGGGGCCGGTTTTCCCGTTCCCCTGATTCCAGTTTTGCGCCAGTTGACGc                                                                         >  2:3541295/1‑149 (MQ=255)
                                                                            cGCCCTGATCGAAAAAAATCATATCGCCAGGCAGCGCCTGGTTTATGTCCTGGCCAATAAACTGGCTGTTGTACTGAATCAAATTAATCGCGGTCACATAGGGGCCGGTTTTCCCGTTCCCCTGATTCCAGTTTTGCGCCAGTTGACGc                                                                         <  2:2253214/149‑1 (MQ=255)
                                                                              cccTGATCGAAAAAAATCATATCGCCAGGCAGCGCCTGGTTTATGTCCTGGCCAATAAACTGGCTGTTGTACTGAATCAAATTAATCGCGGTCACATAGGGGCCGGTTTTCCCGTTCCCCTGATTCCAGTTTTGCGCCAGTTGACGCTg                                                                       >  2:1569184/1‑149 (MQ=255)
                                                                              cccTGATCGAAAAAAATCATATCGCCAGGCAGCGCCTGGTTTATGTCCTGGCCAATAAACTGGCTGTTGTACTGAATCAAATTAATCGCGGTCACATAGGGGCCGGTTTTCCCGTTCCCCTGATTCCAGTTTTGCGCCAGTTGACGCTg                                                                       >  2:1811464/1‑149 (MQ=255)
                                                                              cccTGATCGAAAAAAATCATATCGCCAGGCAGCGCCTGGTTTATGTCCTGGCCAATAAACTGGCTGTTGTACTGAATCAAATTAATCGCGGTCACATAGGGGCCGGTTTTCCCGTTCCCCTGATTCCAGTTTTGCGCCAGTTGACGCTg                                                                       >  2:749549/1‑149 (MQ=255)
                                                                               ccTGATCGAAAAAAATCATATCGCCAGGCAGCGCCTGGTTTATGTCCTGGCCAATAAACTGGCTGTTGTACTGAATCAAATTAATCGCGGTCACATAGGGGCCGGTTTTCCCGTTCCCCTGATTCCAGTTTTGCGCCAGTTGACGCTGt                                                                      >  2:4371736/1‑149 (MQ=255)
                                                                                 tGATCGAAAAAAATCATATCGCCAGGCAGCGCCTGGTTTATGTCCTGGCCAATAAACTGGCTGTTGTACTGAATCAAATTAATCGCGGTCACATAGGGGCCGGTTTTCCCGTTCCCCTGATTCCAGTTTTGCGCCAGTTGACGCTGTTc                                                                    >  1:537315/1‑149 (MQ=255)
                                                                                  gATCGAAAAAAATCATATCGCCAGGCAGCGCCTGGTTTATGTCCTGGCCAATAAACTGGCTGTTGTTCTGAATCAAATTAATCGCGGTCACATAGGGGCCGGTTTTCCCGTTCCCCTGATTCCAGTTTTGCGCCAGTTGACGCTGTTCa                                                                   >  2:2984818/1‑149 (MQ=255)
                                                                                   aTCGAAAAAAATCATATCGCCAGGCAGCGCCTGGTTTATGTCCTGGCCAATAAACTGGCTGTTGTACTGAATCAAATTAATCGCGGTCACATAGGGGCCGGTTTTCCCGTTCCCCTGATTCCAGTTTTGCGCCAGTTGACGCTGTTCAg                                                                  >  1:1654362/1‑149 (MQ=255)
                                                                                      gAAAAAAATCATATCGCCAGGCAGCGCCTGGTTTATGTCCTGGCCAATAAACTGGCTGTTGTACTGAATCAAATTAATCGCGGTCACATAGGGGCCGGTTTTCCCGTTCCCCTGATTCCAGTTTTGCGCCAGTTGACGCTGTTCAGgtg                                                               >  1:1450317/1‑149 (MQ=255)
                                                                                       aaaaaaaTCATATCGCCAGGCAGCGCCTGGTTTATGTCCTGGCCAAGAAACTGGCTGTTGTACTGAATCAAATTAATCGCGGTCACATAGGGGCCGGTTTTCCCGTTCCCCTGATTCCAGTTTTGCGCCAGTTGACGCTGTTCAGgtgt                                                              <  2:1470576/149‑1 (MQ=255)
                                                                                           aaaTCATATCGCCAGGCAGCGCCTGGTTTATGTCCTGGCCAATAAACTGGCTGTTGTTCTGAATCAAATTAATCGCGGTCACATAGGGGCCGGTTTTCCCGTTCCCCTGATTCCAGTTTTGCGCCAGTTGACGCTGTTCAGGTGTTAGc                                                          >  1:61659/1‑149 (MQ=255)
                                                                                                atatCGCCAGGCAGCGCCTGGTTTATGTCCTGGCCAATAAACTGGCTGTTGTACTGAATCAAATTAATCGCGGTCACATAGGGGCCGGTTTTCCCGTTCCCCTGATTCCAGTTTTGCGCCAGTTGACGCTGTTCAGGTGTTAGCGTCAt                                                     <  1:4302663/149‑1 (MQ=255)
                                                                                                atatCGCCAGGCAGCGCCTGGTTTATGTCCTGGCCAATAAACTGGCTGTTGTACTGAATCAAATTAATCGCGGTCACATAGGGGCCGGTTTTCCCGTTCCCCTGATTCCAGTTTTGCGCCAGTTGACGCTGTTCAGGTGTTAGCGTCAt                                                     >  1:3215084/1‑149 (MQ=255)
                                                                                                  atagCCAGGCAGCGCCTGGTTTATGTCCTGGCCAATAAACTGGCTGTTGTACTGAATCAAATTAATCGCGGTCACATAGGGGCCGGTTTTCCCGTTCCCCTGATTCCAGTTTTGCGCCAGTTGACGCTGTTCAGGTGTTAGCGTCAtct                                                   <  2:587254/146‑1 (MQ=255)
                                                                                                   tCGCCAGGCAGCGCCTGGTTTATGTCCTGGCCAATAAACTGGCTGTTGTACTGAATCAAATTAATCGCGGTCACATAGGGGCCGGTTTTCCCGTTCCCCTGATTCCAGTTTTGCGCCAGTTGACGCTGTTCAGGTGTTAGCGTCAtctc                                                  <  2:4413617/149‑1 (MQ=255)
                                                                                                   tCGCCAGGCAGCGCCTGGTTTATGTCCTGGCCAAGAAACTGGCTGTTGTACTGAATCAAATTAATCGCGGTCACATAGGGGCCGGTTTTCCCGTTCCCCTGATTCCAGTTTTGCGCCAGTTGACGCTGTTCAGGTGTTAGCGTCAtctc                                                  <  2:33790/149‑1 (MQ=255)
                                                                                                     gCCAGGCAGCGCCTGGTTTATGTCCTGGCCAATAAACTGGCTGTTGTACTGAATCAAATTAATCGCGGTCACATAGGGGCCGGTTTTCCCGTTCCCCTGATTCCAGTTTTGCGCCAGTTGACGCTGTTCAGGTGTTAGCGTCATCTCTg                                                >  1:3612683/1‑149 (MQ=255)
                                                                                                       cAGGCAGCGCCTGGTTTATGTCCTGGCCAATAAACTGGCTGTTGTACTGAATCAAATTAATCGCGGTCACATAGGGGCCGGTTTTCCCGTTCCCCTGATTCCAGTTTTGCGCCAGTTGACGCTGTTCAGGTGTTAGCGTCATCTCTggc                                              <  1:4479131/149‑1 (MQ=255)
                                                                                                       cAGGCAGCGCCTGGTTTATGTCCTGGCCAATAAACTGGCTGTTGTACTGAATCAAATTAATCGCGGTCACATAGGGGCCGGTTTTCCCGTTCCCCTGATTCCAGTTTTGCGCCAGTTGACGCTGTTCAGGTGTTAGCGTCATCTCTggc                                              <  1:43805/149‑1 (MQ=255)
                                                                                                          gCAGCGCCTGGTTTATGTCCTGGCCAATACACTGGCTGTTGTACTGAATCAAATTAATCGCGGTCACATAGGGGCCGGTTTTCCCGTTCCCCTGATTCCAGTTTTGCGCCAGTTGACGCTGTTCAGGTGTTAGCGTCATCTCTggcggc                                           <  2:21588/149‑1 (MQ=255)
                                                                                                           cAGCGCCTGGTTTATGTCCTGGCCAATAAACTGGCTGTTGTACTGAATCAAATTAATCGCGGTCACATAGGGGCCGGTTTTCCCGTTCCCCTGATTCCAGTTTTGCGCCAGTTGACGCTGTTCAGGTGTTAGCGTCATCTCTGGCGGCa                                          >  2:1084268/1‑149 (MQ=255)
                                                                                                            aGCGCCTGGTTTATGTCCTGGCCAATAAACTGGCTGTTGTACTGAATCAAATTAATCGCGGTCACATAGGGGCCGGTTTTCCCGTTCCCCTGATTCCAGTTTTGCGCCAGTTGACGCTGTTCAGGTGTTAGCGTCATCTCTGGCGGCaa                                         >  2:4238496/1‑149 (MQ=255)
                                                                                                            aGCGCCTGGTTTATGTCCTGGCCAATAAACTGGCTGTTGTACTGAATCAAATTAATCGCGGTCACATAGGGGCCGGTTTTCCCGTTCCCCTGATTCCAGTTTTGCGCCAGTTGACGCTGTTCAGGTGTTAGCGTCATCTCTGGCGGCaa                                         <  1:4238496/149‑1 (MQ=255)
                                                                                                             gcgcCTGGTTTATGTCCTGGCCAATAAACTGGCTGTTGTACTGAATCAAATTAATCGCGGTCACATAGGGGCCGGTTTTCCCGTTCCCCTGATTCCAGTTTTGCGCCAGTTGACGCTGTTCAGGTGTTAGCGTCATCTCTGGCGGCaaa                                        <  2:4509644/149‑1 (MQ=255)
                                                                                                             gcgcCTGGTTTATGTCCTGGCCAATAAACTGGCTGTTGTACTGAATCAAATTAATCGCGGTCACATAGGGGCCGGTTTTCCCGTTCCCCTGATTCCAGTTTTGCGCCAGTTGACGCTGTTCAGGTGTTAGCGTCATCTCTGGCGGCaaa                                        >  1:3337446/1‑149 (MQ=255)
                                                                                                                ccTGGTTTATGTCCTGGCCAATAAACTGGCTGTTGTACTGAATCAAATTAATCGCGGTCACATAGGGGCCGGTTTTCCCGTTCCCCTGATTCCAGTTTTGCGCCAGTTGACGCTGTTCAGGTGTTAGCGTCATCTCTGGCGGCAAATAc                                     <  1:724752/149‑1 (MQ=255)
                                                                                                                     tttATGTCCTGGCCAATAAACTGGCTGTTGTACTGAATCAAATTAATCGCGGTCACATAGGGGCCGGTTTTCCCGTTCCCCTGATTCCAGTTTTGCGCCAGTTGACGCTGTTCAGGTGTTAGCGTCATCTCTGGCGGCAAATACTGGCt                                <  2:626225/149‑1 (MQ=255)
                                                                                                                     tttATGTCCTGGCCAATAAACTGGCTGTTGTACTGAATCAAATTAATCGCGGTCACATAGGGGCCGGTTTTCCCGTTCCCCTGATTCCAGTTTTGCGCCAGTTGACGCTGTTCAGGTGTTAGCGTCATCTCTGGCGGCAAATACTGGCt                                <  1:1959215/149‑1 (MQ=255)
                                                                                                                     tttATGTCCTGGCCAATAAACTGGCTGTTGTACTGAATCAAATTAATCGCGGTCACATAGGGGCCGGTTTTCCCGTTCCCCTGATTCCAGTTTTGCGCCAGTTGACGCTGTTCAGGTGTTAGCGTCATCTCTGGCGGCAAATACTGGCt                                >  2:495185/1‑149 (MQ=255)
                                                                                                                     tttATGTCCTGGCCAATAAACTGGCTGTTGTACTGAATCAAATTAATCGCGGTCACATAGGGGCCGGTTTTCCCGTTCCCCTGATTCCAGTTTTGCGCCAGTTGACGCTGTTCAGGTGTTAGCGTCATCTCTGGCGGCAAATACTGGCt                                <  1:3455487/149‑1 (MQ=255)
                                                                                                                       tATGTCCTGGCCAATAAACTGGCTGTTGTACTGAATCAAATTAATCGCGGTCACATAGGGGCCGGTTTTCCCGTTCCCCTGATTCCAGTTTTGCGCCAGTTGACGCTGTTCAGGTGTTAGCGTCATCTCTGGCGGCAAATACTGGCTAg                              >  1:3265695/1‑149 (MQ=255)
                                                                                                                        aTGTCCTGGCCAATAAACTGGCTGTTGTACTGAATCAAATTAATCGCGGTCACATAGGGGCCGGTTTTCCCGTTCCCCTGATTCCAGTTTTGCGCCAGTTGACGCTGTTCAGGTGTTAGCGTCATCTCTGGCGGCAAATACTGGCTAGa                             <  1:2679689/149‑1 (MQ=255)
                                                                                                                         tGTCCTGGCCAATAAACTGGCTGTTGTACTGAATCAAATTAATCGCGGTCACATAGGGGCCGGTTTTCCCGTTCCCCTGATTCCAGTTTTGCGCCAGTTGACGCTGTTCAGGTGTTAGCGTCATCTCTGGCGGCAAATACTGGCTAGAt                            >  1:3856618/1‑149 (MQ=255)
                                                                                                                          gTCCTGGCCAATAAACTGGCTGTTGTACTGAATCAAATTAATCGCGGTCACATAGGGGCCGGTTTTCCCGTTCCCCTGATTCCAGTTTTGCGCCAGTTGACGCTGTTCAGGTGTTAGCGTCATCTCTGGCGGCAAATACTGGCTAGATa                           <  1:1261956/149‑1 (MQ=255)
                                                                                                                           tCCTGGCCAATAAACTGGCTGTTGTACTGAATCAAATTAATCGCGGTCACATAGGGGCCGGTTTTCCCGTTCCCCTGATTCCAGTTTTGCGCCAGTTGACGCTGTTCAGGTGTTAGCGTCATCTCTGGCGGCAAATACTGGCTAGATaa                          >  2:1683069/1‑149 (MQ=255)
                                                                                                                             cTGGCCAATAAACTGGCTGTTGTACTGAATCAAATTAATCGCGGTCACATAGGGGCCGGTTTTCCCGTTCCCCTGATTCCAGTTATGCGCCAGTTGACGCTGTTCAGGTGTTAGCGTCATCTCTGGCGGCAAATACTGTCTAGATAAAc                        >  1:1185129/1‑149 (MQ=255)
                                                                                                                              tGGCCAATAAACTGGCTGTTGTACTGAATCAAATTAATCGCGGTCACATAGGGGCCGGTTTTCCCGTTCCCCTGATTCCAGTTTTGCGCCAGTTGACGCTGTTCAGGTGTTAGCGTCATCTCTGGCGGCAAATACTGGCTAGATAAAcc                       <  2:3250343/149‑1 (MQ=255)
                                                                                                                                  cAATAAACTGGCTGTTGTACTGAATCAAATTAATCGCGGTCACATAGGGGCCGGTTTTCCCGTTCCCCTGATTCCAGTTTTGCGCCAGTTGACGCTGTTCAGGTGTTAGCGTCATCTCTGGCGGCAAATACTGGCTAGATAAACCGTTa                   >  2:464449/1‑149 (MQ=255)
                                                                                                                                  cAATAAACTGGCTGTTGTACTGAATCAAATTAATCGCGGTCACATAGGGGCCGGTTTTCCCGTTCCCCTGATTCCAGTTTTGCGCCAGTTGACGCTGTTCAGGTGTTAGCGTCATCTCTGGCGGCAAATACTGGCTAGATAAACCGTTa                   >  1:3787484/1‑149 (MQ=255)
                                                                                                                                   aaTAAACTGGCTGTTGTACTGAATCAAATTAATCGCGGTCACATAGGGGCCGGTTTTCCCGTTCCCCTGATTCCAGTTTTGCGCCAGTTGACGCTGTTCAGGTGTTAGCGTCATCTCTGGCGGCAAATACTGGCTAGATAAACCGTTAc                  <  2:3203088/149‑1 (MQ=255)
                                                                                                                                   aaTAAACTGGCTGTTGTACTGAATCAAATTAATCGCGGTCACATAGGGGCCGGTTTTCCCGTTCCCCTGATTCCAGTTTTGCGCCAGTTGACGCTGTTCAGGTGTTAGCGTCATCTCTGGCGGCAAATACTGGCTAGATAAACCGTTAc                  <  1:3093793/149‑1 (MQ=255)
                                                                                                                                    aTAAACTGGCTGTTGTACTGAATCAAATTAATCGCGGTCACATAGGGGCCGGTTTTCCCGTTCCCCTGATTCCAGTTTTGCGCCAGTTGACGCTGTTCAGGTGTTAGCGTCATCTCTGGCGGCAAATACTGGCTAGATAAACCGTTACt                 <  2:4980332/149‑1 (MQ=255)
                                                                                                                                     tAAACTGGCTGTTGTTCTGAATCAAATTAATCGCGGTCACATAGGGGCCGGTTTTCCCGTTCCCCTGATTCCAGTTTTGCGCCAGTTGACGCTGTTCAGGTGTTAGCGTCATCTCTGGCGGCAAATACTGGCTAGATAAACCGTTACtt                >  1:3993370/1‑149 (MQ=255)
                                                                                                                                     tAAACTGGCTGTTGTACTGAATCAAATTAATCGCGGTCACATAGGGGCCGGTTTTCCCGTTCCCCTGATTCCAGTTTTGCGCCAGTTGACGCTGTTCAGGTGTTAGCGTCATCTCTGGCGGCAAATACTGGCTAGATAAACCGTTACtt                >  1:2238950/1‑149 (MQ=255)
                                                                                                                                      aaaCTGGCTGTTGTACTGAATCAAATTAATCGCGGTCACATAGGGGCCGGTTTTCCCGTTCCCCTGATTCCAGTTTTGCGCCAGTTGACGCTGTTCAGGTGTTAGCGTCATCTCTGGCGGCAAATACTGGCTAGATAAACCGTTACttt               <  2:4178122/149‑1 (MQ=255)
                                                                                                                                       aaCTGGCTGTTGTACTGAATCAAATTAATCGCGGTCACATAGGGGCCGGTTTTCCCGTTCCCCTGATTCCAGTTTTGCGCCAGTTGACGCTGTTCAGGTGTTAGCGTCATCTCTGGCGGCAAATACTGGCTAGATAAACCGTTACtttt              <  1:3247070/149‑1 (MQ=255)
                                                                                                                                           ggCTGTTGTTCTGAATCAAATTAATCGCGGTCACATAGGGGCCGGTTTTCCCGTTCCCCTGATTCCAGTTTTGCGCCAGTTGACGCTGTTCAGGTGTTAGCGTCATCTCTGGCGGCAAATACTGGCTAGATAAACCGTTACTTTTTAGc          <  2:2296289/149‑1 (MQ=255)
                                                                                                                                           ggCTGTTGTACTGAATCAAATTAATCGCGGTCACATAGGGGCCGGTTTTCCCGTTCCCCTGATTCCAGTTTTGCGCCAGTTGACGCTGTTCAGGTGTTAGCGTCATCTCTGGCGGCAAATACTGGCTAGATAAACCGTTACTTTTTAGc          <  1:1872165/149‑1 (MQ=255)
                                                                                                                                           ggCTGTTGTACTGAATCAAATTAATCGCGGTCACATAGGGGCCGGTTTTCCCGTTCCCCTGATTCCAGTTTTGCGCCAGTTGACGCTGTTCAGGTGTTAGCGTCATCTCTGGCGGCAAATACTGGCTAGATAAACCGTTACTTTTTAGc          <  1:812701/149‑1 (MQ=255)
                                                                                                                                           ggCTGTTGTACTGAATCAAATTAATCGCGGTCACATAGGGGCCGGTTTTCCCGTTCCCCTGATTCCAGTTTTGCGCCAGTTGACGCTGTTCAGGTGTTAGCGTCATCTCTGGCGGCAAATACTGGCTAGATAAACCGTTACTTTTTAGc          >  2:234781/1‑149 (MQ=255)
                                                                                                                                            gCTGTTGTACTGAATCAAATTAATCGCGGTCACATAGGGGCCGGTTTTCCCGTTCCCCTGATTCCAGTTTTGCGCCAGTTGACGCTGTTCAGGTGTTAGCGTCATCTCTGGCGGCAAATACTGGCTAGATAAACCGTTACTTTTTAGcc         <  2:835678/149‑1 (MQ=255)
                                                                                                                                            gCTGTTGTACTGAATCAAATTAATCGCGGTCACATAGGGGCCGGTTTTCCCGTTCCCCTGATTCCAGTTTTGCGCCAGTTGACGCTGTTCAGGTGTTAGCGTCATCTCTGGCGGCAAATACTGGCTAGATAAACCGTTACTTTTTAGcc         <  2:3833830/149‑1 (MQ=255)
                                                                                                                                            gCTGTTGTACTGAATCAAATTAATCGCGGTCACATAGGGGCCGGTTTTCCCGTTCCCCTGATTCCAGTTTTGCGCCAGTTGACGCTGTTCAGGTGTTAGCGTCATCTCTGGCGGCAAATACTGGCTAGATAAACCGTTACTTTTTAGcc         <  2:2654720/149‑1 (MQ=255)
                                                                                                                                            gCTGGTGTACTGAATCAAATTAATCGCGGTCACATAGGGGCCGGTTTTCCCGTTCCCCTGATTCCAGTTTTGCGCCAGTTGACGCTGTTCAGGTGTTAGCGTCATCTCTGGCGGCAAATACTGGCTAGATAAACCGTTACTTTTTAGcc         <  2:4050884/149‑1 (MQ=255)
                                                                                                                                              tgttgtACTGAATCAAATTAATCGCGGTCACATAGGGGCCGGTTTTCCCGTTCCCCTGATTCCAGTTTTGCGCCAGTTGACGCTGTTCAGGTGTTAGCGTCATCTCTGGCGGCAAATACTGGCTAGATAAACCGTTACTTTTTAGCCAt       >  1:1002194/1‑149 (MQ=255)
                                                                                                                                              tgttgtACTGAATCAAATTAATCGCGGTCACATAGGGGCCGGTTTTCCCGTTCCCCTGATTCCAGTTTTGCGCCAGTTGACGCTGTTCAGGTGTTAGCGTCATCTCTGGCGGCAAATACTGGCTAGATAAACCGTTACTTTTTAGCCAt       <  2:2016973/149‑1 (MQ=255)
                                                                                                                                                ttgtACTGAATCAAATTAATCGCGGTCACATAGGGGCCGGTTTTCCCGTTCCCCTGATTCCAGTTTTGCGCCAGTTGACGCTGTTCAGGTGTTAGCGTCATCTCTGGCGGCAAATACTGGCTAGATAAACCGTTACTTTTTAGCCAttt     <  1:3727739/149‑1 (MQ=255)
                                                                                                                                                ttgtACTGAATCAAATTAATCGCGGTCACATAGGGGCCGGTTTTCCCGTTCCCCTGATTCCAGTTTTGCGCCAGTTGACGCTGTTCAGGTGTTAGCGTCATCTCTGGCGGCAAATACTGGCTAGATAAACCGTTACTTTTTAGCCAttt     <  1:4053216/149‑1 (MQ=255)
                                                                                                                                                   tACTGAATCAAATTAATCGCGGTCACATAGGGGCCGGTTTTCCCGTTCCCCTGATTCCAGTTTTGCGCCAGTTGACGCTGTTCAGGTGTTAGCGTCATCTCTGGCGGCAAATACTGGCTAGATAAACCGTTACTTTTTAGCCATTTACt  <  2:3129906/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
CCGTTGTCAGTTTTCGTGGCGCTTCCGGTGTGGTAGATGACGTAACGCCCCATCCAGACCATTAAGTGCTGGGCATCGCCCTGATCGAAAAAAATCATATCGCCAGGCAGCGCCTGGTTTATGTCCTGGCCAATAAACTGGCTGTTGTACTGAATCAAATTAATCGCGGTCACATAGGGGCCGGTTTTCCCGTTCCCCTGATTCCAGTTTTGCGCCAGTTGACGCTGTTCAGGTGTTAGCGTCATCTCTGGCGGCAAATACTGGCTAGATAAACCGTTACTTTTTAGCCATTTACT  >  CP009273/2327930‑2328225

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: