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

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273557,5920GT96.7% 661.0 / 12.1 214D814Y (GAC→TAC) sfmDputative outer membrane export usher protein
Reads supporting (aligned to +/- strand):  ref base G (1/6);  new base T (106/100);  total (107/107)
Fisher's exact test for biased strand distribution p-value = 6.52e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.92e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

Only 100 of 217 total aligned reads displayed.

GACGCAAGGCGCGTTGGTTCGTGCCAATTTTGATACCCGCATAGGCGTGCGGGCGCTCATTACCGTTACCCAGGGCGGAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATGACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACTGGTTCAGTGGGGAGACGGCGCGAACTCACGCTGCATTGCGCACTATGTATTGCCGAAGCAAAGCTTACAGCAAGCCGTCACTGTTATTTCGGCAGTT  >  CP009273/557447‑557735
                                                                                                                                                  |                                                                                                                                               
gACGCAAGGCGCGTTGGTTCGTGCCAATTTTGATACCCGCATAGGCGTGCGGGCGCTCATTACCGTTACCCAGGGCGGAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATTAc                                                                                                                                               <  2:2977663/148‑1 (MQ=255)
  cGCAAGGCGCGTTGGTTCGTGCCAATTTTGATACCCGCATAGGCGTGCGGGCGCTCATTACCGTTACCCAGGGCGGAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATTACgg                                                                                                                                             >  1:418100/1‑148 (MQ=255)
     aaGGCGCGTTGGTTCGTGCCAATTTTGATACCCGCATAGGCGTGCGGGCGCTCATTACCGTTACCCAGGGCGGAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATTACGGGCa                                                                                                                                          >  1:3456141/1‑148 (MQ=255)
      aGGCGCGTTGGTTCGTGCCAATTTTGATACCCGCATAGGCGTGCGGGCGCTCATTACCGTTACCCAGGGCGGAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATTACGGGCaa                                                                                                                                         <  1:2571336/148‑1 (MQ=255)
        gcgcGTTGGTTCGTGCCAATTTTGATACCCGCATAGGCGTGCGGGCGCTCATTACCGTTACCCAGGGCGGAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATTACGGGCAAGt                                                                                                                                       <  1:2737455/148‑1 (MQ=255)
        gcgcGTTGGTTCGTGCCAATTTTGATACCCGCATAGGCGTGCGGGCGCTCATTACCGTTACCCAGGGCGGAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATTACGGGCAAGt                                                                                                                                       <  1:3205498/148‑1 (MQ=255)
         cgcgTTGGTTCGTGCCAATTTTGATACCCGCATAGGCGTGCGGGCGCTCATTACCGTTACCCAGGGCGGAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATTACGGGCAAGtt                                                                                                                                      <  1:2923339/148‑1 (MQ=255)
           cgTTGGTTCGTGCCAATTTTGATACCCGCATAGGCGTGCGGGCGCTCATTACCGTTACCCAGGGCGGAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATTACGGGCAAGttta                                                                                                                                    <  1:624201/148‑1 (MQ=255)
             ttGGTTCGTGCCAATTTTGATACCCGCATAGGCGTGCGGGCGCTCATTACCGTTACCCAGGGCGGAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATTACGGGCAAGtttatt                                                                                                                                  >  2:607785/1‑148 (MQ=255)
             ttGGTTCGTGCCAATTTTGATACCCGCATAGGCGTGCGGGCGCTCATTACCGTTACCCAGGGCGGAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATTACGGGAAATtttatt                                                                                                                                  >  2:3026634/1‑148 (MQ=255)
              tGGTTCGTGCCAATTTTGATACCCGCATAGGCGTGCGGGCGCTCATTACCGTTACCCAGGGCGGAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATTAAGGGCAAGtttattt                                                                                                                                 >  1:1276529/1‑148 (MQ=255)
              tGGTTCGTGCCAATTTTGATACCCGCATAGGCGTGCGGGCGCTCATTACCGTTACCCAGGGCGGAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGGGGGTGATTACGGGCAAGtttattt                                                                                                                                 >  1:1335285/1‑148 (MQ=255)
               ggTTCGTGCCAATTTTGATACCCGCATAGGCGTGCGGGCGCTCATTACCGTTACCCAGGGCGGAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATTACGGGCAAGtttattta                                                                                                                                <  2:3417865/148‑1 (MQ=255)
                gTTCGTGCCAATTTTGAGACCCGCAGAGGCGTGCGGGCGCTCATTACCGTTACCCAGGGCGGAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATTACGGGCAAGTTTATTTaa                                                                                                                               <  2:2070562/148‑1 (MQ=255)
                 ttCGTGCCAATTTTGATACCCGCATAGGCGTGCGGGCGCTCATTACCGTTACCCAGGGCGGAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATTACGGGCAAGTTTATTTAAg                                                                                                                              >  2:853789/1‑148 (MQ=255)
                 ttCGTGCCAATTTTGATACCCGCATAGGCGTGCGGGCGCTCATTACCGTTACCCAGGGCGGAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATTACGGGCAAGTTTATTTAAg                                                                                                                              >  2:174018/1‑148 (MQ=255)
                  tCGTGCCAATTTTGATACCCGCATAGGCGTGCGGGCGCTCATTACCGTTACCCAGGGCGGAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATTACGGGCACGTTTATTTAAgt                                                                                                                             >  1:3151133/1‑148 (MQ=255)
                  tCGTGCCAATTTTGATACCCGCATAGGCGTGCGGGCGCTCATTACCGTTACCCAGGGCGGAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATTACGGGCAAGTTTATTTAAgt                                                                                                                             >  1:1996493/1‑148 (MQ=255)
                  tCGTGCCAATTTTGATACCCGCATAGGCGTGCGGGCGCTCATTACCGTTACCCAGGGCGGAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATTACGGGCAAGTTTATTTAAgt                                                                                                                             >  1:201334/1‑148 (MQ=255)
                  tCGTGCCAATTTTGATACCCGCATAGGCGTGCGGGCGCTCATTACCGTTACCCAGGGCGGAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATTACGGGCAAGTTTATTTAAgt                                                                                                                             <  1:2740910/148‑1 (MQ=255)
                   cGTGCCAATTTTGATACCCGCATAGGCGTGCGGGCGCTCATTACCGTTACCCAGGGCGGAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATTACGGGCAAGTTTATTTAAgtg                                                                                                                            <  2:2884547/148‑1 (MQ=255)
                     tGCCAATTTTGATACCCGCATAGGCGTGCGGGCGCTCATTACCGTTACCCAGGGCGGAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATTACGGGCAAGTTTATTTAAgtggt                                                                                                                          <  2:3081645/148‑1 (MQ=255)
                          aTTTTGATACCCGCATAGGCGTGCGGGCGCTCATTACCGTTACCCAGGGCGGAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATGACGGGCAAGTTTATTTAAGTGGTGCGcc                                                                                                                     <  1:968284/148‑1 (MQ=255)
                           ttttGATACCCGCATAGGCGTGCGGGCGCTCATTACCGTTACCCAGGGCGGAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATTACGGGCAAGTTTATTTAAGTGGTGCGCCa                                                                                                                    <  2:2332879/148‑1 (MQ=255)
                              tGATACCCGCATAGGCGTGCGGGCGCTCATTACCGTTACCCAGGGCGGAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATTACGGGCAAGTTTATTTAAGTGGTGCGCCATTg                                                                                                                 >  2:423785/1‑148 (MQ=255)
                                aTACCCGCATAGGCGTGCGGGCGCTCATTACCGTTACCCAGGGCGGAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATTACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTc                                                                                                               >  2:2620417/1‑148 (MQ=255)
                                  aCCCGCATAGGCGTGCGGGCGCTCATTACCGTTACCCAGGGCGGAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATTACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTg                                                                                                             >  2:1213360/1‑148 (MQ=255)
                                       cATAGGCGTGCGGGCGCTCATTACCGTTACCCAGGGCGGAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCTGTATGGTGGGTGATTACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGaa                                                                                                        <  2:498497/148‑1 (MQ=255)
                                        aTAGGCGTGCGGGCGCTCATTACCGTTACCCAGGGCGGAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATTACGGGCAAGTTTATTTAAGTGGTGCGcc                                                                                                                     >  1:3557763/1‑134 (MQ=255)
                                        aTAGGCGTGCGGGCGCTCATTACCGTTACCCAGGGCGGAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATTACGGGCAAGTTTATTTAAGTGGTGCGcc                                                                                                                     <  2:3557758/134‑1 (MQ=255)
                                        aTAGGCGTGCGGGCGCTCATTACCGTTACCCAGGGCGGAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATGACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAAt                                                                                                       <  1:816698/148‑1 (MQ=255)
                                          aGGCGTGCGGGCGCTCATTACCGTTACCCAGGGCGGAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATTACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTa                                                                                                     <  1:1394198/148‑1 (MQ=255)
                                          aGGCGTGCGGGCGCTCATTACCGTTACCCAGGGCGGAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAAAAACCAGTATGGTGGGTGATTACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTa                                                                                                     >  1:1638549/1‑148 (MQ=255)
                                          aGGCGTGCGGGCGCTCATTACCGTTACCCAGGGCGGAAAACCGGTGCCGTTTGGATCACAGGTACGGGAAAAAAGTAACGGAATAACCCGGATGGTGGGTGATTAAGGGCAAGTTTATTTAAGTGGTGCGCCCGTGTCTGGGGAATTa                                                                                                     >  1:1523096/1‑148 (MQ=255)
                                          aGGCGTGCGGGCGCTCATTACCGTTACCCAGGGCGGAAAACCGGTGCCGGGTGGATCACTGGTACGGGAAAACAGTACCGGAAATAACCGTATGGTGGGTGATTACGGGCCAGTTTATTTAAGTGGTGCGCCAGTGTCTGGTGGATTa                                                                                                     >  1:3296319/1‑148 (MQ=255)
                                          aGGCGTGCGGGCGCTCAGTACCGTTACCCAGGGCGGAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATTACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTa                                                                                                     <  2:163193/148‑1 (MQ=255)
                                             cGTGCGGGCGCTCATTACCGTTACCCAGGGCGGAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATTACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACTg                                                                                                  <  2:632337/148‑1 (MQ=255)
                                                  gggCGCTCATTACCGTTACCCAGGGCGGAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATTACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACTGGTTCa                                                                                             <  1:2810712/148‑1 (MQ=255)
                                                   ggCGCTCATTACCGTTACCCAGGGCGGAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATTACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACTGGTTCAg                                                                                            >  1:454687/1‑148 (MQ=255)
                                                    gcgcTCATTACCGTTACCCAGGGCGGAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATTACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACTGGTTCAGt                                                                                           >  1:428359/1‑148 (MQ=255)
                                                     cgcTCATTACCGTTACCCAGGGCGGAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATTACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACTGGTTCAGTg                                                                                          <  1:119636/148‑1 (MQ=255)
                                                     cgcTCATTACCGTTACCCAGGGCGGAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATTACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACTGGTTCAGTg                                                                                          <  1:1714160/148‑1 (MQ=255)
                                                     cgcTCATTACCGTTACCCAGGGCGGAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATTACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACTGGTTCAGTg                                                                                          >  2:2221097/1‑148 (MQ=255)
                                                     cgcTCATTACCGTTACCCAGGGCGGAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATGACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACTGGTTCAGTg                                                                                          >  2:3458391/1‑148 (MQ=255)
                                                     cgcTCATTACCGGTACCCAGGGCGGAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATTACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACTGGTTCAGTg                                                                                          <  1:1919863/148‑1 (MQ=255)
                                                        tCATTACCGTTACCCAGGGCGGAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATTACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACTGGTTCAGTgggg                                                                                       >  1:720394/1‑148 (MQ=255)
                                                           ttACCGTTACCCAGGGCGGAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATTACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACTGGTTCAGTGGGgaga                                                                                    >  1:3459311/1‑148 (MQ=255)
                                                              ccGTTACCCAGGGCGGAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATTACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACTGGTTCAGTGGGGAGACgg                                                                                 >  2:1219005/1‑148 (MQ=255)
                                                              ccGTTACCCAGGGCGGAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATTACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACTGGTTCAGTGGGGAGACgg                                                                                 >  2:2277340/1‑148 (MQ=255)
                                                                gTTACCCAGGGCGGAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATTACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACTGGTTCAGTGGGGAGACGgcg                                                                               <  2:1319868/148‑1 (MQ=255)
                                                                gTTACCCAGGGCGGAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATTACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTATGGTGAATTACTGGTTCAGTGGGGAGACGgcg                                                                               <  1:845413/148‑1 (MQ=255)
                                                                          gCGGAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATTACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACTGGTTCAGTGGGGAGACGGCGCGAACTCACg                                                                     <  1:2785780/148‑1 (MQ=255)
                                                                          gCGGAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATTACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACTGGTTCAGTGGGGAGACGGCGCGAACTCACg                                                                     <  1:3414582/148‑1 (MQ=255)
                                                                             gAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATTACGGTCTAGTTTATTTAAGTGGTGCGCAAGTGTCTGGTGAATTACTGGTTCAGTGGGGAGACGGCGCGATCTCACGCTg                                                                  <  1:869977/148‑1 (MQ=255)
                                                                             gAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATTACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACTGGTTCAGTGGGGAGACGGCGCGAACTCACGCTg                                                                  >  1:1817689/1‑148 (MQ=255)
                                                                             gAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATTACGGGCAAGTTTATTTAAGTGGTGCGCCAGTGTCTGGTGAATTACTGGTTCAGTGGGGAGACGGCGCGAACTCACGCTg                                                                  <  2:2141708/148‑1 (MQ=255)
                                                                                aaCCGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATTACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACTGGTTCAGTGGGGAGACGGCGCGAACTCACGCTGCAt                                                               >  2:384768/1‑148 (MQ=255)
                                                                                    ggTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATTACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACTGGTTCAGTGGGGAGACGGCGCGAACTCACGCTGCATTgcg                                                           <  1:1791642/148‑1 (MQ=255)
                                                                                     gTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATTACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACTGGTTCAGTGGGGAGACGGCGCGAACTCACGCTGCATTgcgc                                                          <  2:558452/148‑1 (MQ=255)
                                                                                     gTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATTACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACTGGTTCAGTGGGGAGACGGCGCGAACTCACGCTGCATTgcgc                                                          <  2:739637/148‑1 (MQ=255)
                                                                                     gTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATGACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACTGGTTCAGTGGGGAGACGGCGCGAACTCACGCTGCATTgcgc                                                          <  1:2208176/148‑1 (MQ=255)
                                                                                      tGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATTACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACTGGTTCAGTGGGGAGACGGCGCGAACTCACGCTGCATTGCGCa                                                         <  2:542030/148‑1 (MQ=255)
                                                                                      tGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATTACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACTGGTTCAGTGGGGAGACGGCGCGAACTCACGCTGCATTGCGCa                                                         <  1:182509/148‑1 (MQ=255)
                                                                                          gTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATTACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACTGGTTCAGTGGGGAGACGGCGCGAACTCACGCTGCATTGCGCACTAt                                                     <  2:1989906/148‑1 (MQ=255)
                                                                                          gTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATTACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACTGGTTCAGTGGGGAGACGGCGCGAACTCACGCTGCATTGCGCACTAt                                                     <  2:716238/148‑1 (MQ=255)
                                                                                             tGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATTACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACTGGTTCAGTGGGGAGACGGCGCGAACTCACGCTGCATTGCGCACTATGTa                                                  <  1:2488614/148‑1 (MQ=255)
                                                                                               gATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATTACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACTGGTTCAGTGGGGAGACGGCGCGAACTCACGCTGCATTGCGCACTATGTAtt                                                >  1:1553264/1‑148 (MQ=255)
                                                                                               gATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATTACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACTGGTTCAGTGGGGAGACGGCGCGAACTCACGCTGCATTGCGCACTATGTAtt                                                >  2:3375449/1‑148 (MQ=255)
                                                                                                  cACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATTACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACTGGTTCAGTGGGGAGACGGCGCGAACTCACGCTGCATTGCGCACTATGTATTGcc                                             <  1:2145417/148‑1 (MQ=255)
                                                                                                  cACGGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATTACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACTGGTTCAGTGGGGAGACGGCGCGAACTCACGCTGCATTGCGCACTATGTATTGcc                                             <  2:3286984/148‑1 (MQ=255)
                                                                                                   aCTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATTACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACTGGTTCAGTGGGGAGACGGCGCGAACTCACGCTGCATTGCGCACTATGTATTGCCg                                            >  2:1414956/1‑148 (MQ=255)
                                                                                                      ggTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATTACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACTGGTTCAGTGGGGAGACGGCGCGAACTCACGCTGCATTGCGCACTATGTATTGCCGAAg                                         <  1:523770/148‑1 (MQ=255)
                                                                                                        tACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATTACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACTGGTTCAGTGGGGAGACGGCGCGAACTCACGCTGCATTGCGCACTATGTATTGCCGAAGCa                                       >  1:205059/1‑148 (MQ=255)
                                                                                                          cGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATTACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACTGGTTCAGTGGGGAGACGGCGCGAACTCACGCTGCATTGCGCACTATGTATTGCCGAAGCaaa                                     >  2:1552981/1‑148 (MQ=255)
                                                                                                           gggAAAACAGTACCGGAATAACCAGTATGGTGGGTGATACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACTGGTTCAGTGGGGAGACGGCGCGAACTCACGCTGCATTGCGCACTATGTATTGCCGAAGCAAAGc                                   <  2:319814/148‑1 (MQ=255)
                                                                                                            ggAAAACAGTACCGGAATAACCAGTATGGTGGGTGATTACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACTGGTTCAGTGGGGAGACGGCGCGAACTCACGCTGCATTGCGCACTATGTATTGCCGAAGCAAAGc                                   >  2:410329/1‑148 (MQ=255)
                                                                                                               aaaCAGTACCGGAATAACCAGTATGGTGGGTGATTACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACTGGTTCAGTGGGGAGACGGCGCGAACTCACGCTGCATTGCGCACTATGTATTGCCGAAGCAAAGCTTa                                >  2:3002980/1‑148 (MQ=255)
                                                                                                               aaaCAGTACCGGAATAACCAGTATGGTGGGTGATTACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACTGGTTCAGTGGGGAGACGGCGCGAACTCACGCTGCATTGCGCACTATGTATTGCCGAAGCAAAGCTTa                                <  2:201343/148‑1 (MQ=255)
                                                                                                                aaCAGTACCGGAATAACCAGTATGGTGGGTGATTACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACTGGTTCAGTGGGGAGACGGCGCGAACTCACGCTGCATTGCGCACTATGTATTGCCGAAGCAAAGCTTAc                               >  1:2175844/1‑148 (MQ=255)
                                                                                                                 aCAGTACCGGAATAACCAGTATGGTGGGTGATTACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACTGGTTCAGTGGGGAGACGGCGCGAACTCACGCTGCATTGCGCACTATGTATTGCCGAAGCAAAGCTTACa                              >  2:2949467/1‑148 (MQ=255)
                                                                                                                     tACCGGAATAACCAGTATGGTGGGTGATTACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACTGGTTCAGTGGGGAGACGGCGCGAACTCACGCTGCATTGCGCACTATGTATTGCCGAAGCAAAGCTTACAGCaa                          >  1:1676732/1‑148 (MQ=255)
                                                                                                                     tACCGGAATAACCAGTATGGTGGGTGATTACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACTGGTTCAGTGGGGAGACGGCGCGAACTCACGCTGCATTGCGCACTATGTATTGCCGAAGCAAAGCTTACAGCaa                          >  1:1079050/1‑148 (MQ=255)
                                                                                                                      aCCGGAATAACCAGTATGGTGGGTGATTACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACTGGTTCAGTGGGGAGACGGCGCGAACTCACGCTGCATTGCGCACTATGTATTGCCGAAGCAAAGCTTACAGCAAg                         <  2:1504837/148‑1 (MQ=255)
                                                                                                                       ccGGAATAACCAGTATGGTGGGTGATGACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACTGGTTCAGTGGGGAGACGGCGCGAACTCACGCTGCATTGCGCACTATGTATTGCCGAAGCAAAGCTTACAGCAAGc                        <  1:2769436/148‑1 (MQ=255)
                                                                                                                        cGGAATAACCAGTATGGTGGGTGATTACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACTGGTTCAGTGGGGAGACGGCGCGAACTCACGCTGCATTGCGCACTATGTATTGCCGAAGCAAAGCTTACAGCAAGcc                       >  1:2849723/1‑148 (MQ=255)
                                                                                                                            aaTAACCAGTATGGTGGGTGATTACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACTGGTTCAGTGGGGAGACGGCGCGAACTCACGCTGCATTGCGCACTATGTATTGCCGAAGCAAAGCTTACAGCAAGCCGTc                    >  2:837611/1‑148 (MQ=255)
                                                                                                                              tAACCAGTATGGTGGGTGATTACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACTGGTTCAGTGGGGAGACGGCGCGAACTCACGCTGCATTGCGCACTATGTATTGCCGAAGCAAAGCTTACAGCAAGCCGTCAc                  >  1:1417373/1‑148 (MQ=255)
                                                                                                                              tAACCAGTATGGTGGGTGATTACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACTGGTTCAGTGGGGAGACGGCGCGAACTCACGCTGCATTGCGCACTATGTATTGCCGAAGCAAAGCTTACAGCAAGCCGTCAc                  <  2:2646234/148‑1 (MQ=255)
                                                                                                                              tAACCAGTATGGTGGGTGATTACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACTGGTTCAGTGGGGAGACGGCGCGAACTCACGCTGCATTGCGCACTATGTATTGCCGAAGCAAAGCTTACAGCAAGCCGTCAc                  <  1:3525720/148‑1 (MQ=255)
                                                                                                                              tAACCAGTATGGTGGGTGATTACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACTGGTTCAGTGGGGAGACGGCGCGAACTCACGCTGCATTGCGCACTATGTATTGCCGAAGCAAAGCTTACAGCAAGCCGTCAc                  >  2:717602/1‑148 (MQ=255)
                                                                                                                                  cAGTATGGTGGGTGATTACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACTGGTTCAGTGGGGAGACGGCGCGAACTCACGCTGCATTGCGCACTATGTATTGCCGAAGCAAAGCTTACAGCAAGCCGTCACTGtt              >  2:2575235/1‑148 (MQ=255)
                                                                                                                                   aGTATGGTGGGTGATTACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACTGGTTCAGTGGGGAGACGGCGCGAACTCACGCTGCATTGCGCACTATGTATTGCCGAAGCAAAGCTTACAGCAAGCCGTCACTGTTa             >  1:2454737/1‑148 (MQ=255)
                                                                                                                                    gTATGGTGGGTGATTACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACTGGTTCAGTGGGGAGACGGCGCGAACTCACGCTGCATTGCGCACTATGTATTGCCGAAGCAAAGCTTACAGCAAGCCGTCACTGTTAt            <  2:1819398/148‑1 (MQ=255)
                                                                                                                                     tATGGTGGGTGATTACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACTGGTTCAGTGGGGAGACGGCGCGAACTCACGCTGCATTGCGCACTATGTATTGCCGAAGCAAAGCTTACAGCAAGCCGTCACTGTTAtt           <  1:3276129/148‑1 (MQ=255)
                                                                                                                                     tATGGTGGGTGATTACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACTGGTTCAGTGGGGAGACGGCGCGAACTCACGCTGCATTGCGCACTATGTATTGCCGAAGCAAAGCTTACAGCAAGCCGTCACTGTTAtt           >  2:276392/1‑148 (MQ=255)
                                                                                                                                            ggtgATTACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACTGGTTCAGTGGGGAGACGTCGCGAACTCACGCTGCATTGCGCACTATGTATTGCCGAAGCAAAGCTTACAGCAAGCCGTCACTGTTATTTCGGCAg    >  1:1758992/1‑148 (MQ=255)
                                                                                                                                             gtgATTACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACTGGTTCAGTGGGGAGACGGCGCGAACTCACGCTGCATTGCGCACTATGTATTGCCGAAGCAAAGCTTACAGCAAGCCGTCACTGTTATTTCGGCAGt   <  1:795805/148‑1 (MQ=255)
                                                                                                                                             gtgATTACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACTGGTTCAGTGGGGAGACGGCGCGAACTCACGCTGCATTGCGCACTATGTATTGCCGAAGCAAAGCTTACAGCAAGCCGTCACTGTTATTTCGGCAGt   <  1:2579752/148‑1 (MQ=255)
                                                                                                                                              tgattaCGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACTGGTTCAGTGGGGAGACGGCGCGAACTCACGCTGCATTGCGCACTATGTATTGCCGAAGCAAAGCTTACAGCAAGCCGTCACTGTTATTTCGGCAGtt  <  1:1585446/148‑1 (MQ=255)
                                                                                                                                              tgattaCGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACTGGTTCAGTGGGGAGACGGCGCGAACTCACGCTGCATTGCGCACTATGTATTGCCGAAGCAAAGCTTACAGCAAGCCGTCACTGTTATTTCGGCAGtt  >  2:1041461/1‑148 (MQ=255)
                                                                                                                                                  |                                                                                                                                               
GACGCAAGGCGCGTTGGTTCGTGCCAATTTTGATACCCGCATAGGCGTGCGGGCGCTCATTACCGTTACCCAGGGCGGAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATGACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACTGGTTCAGTGGGGAGACGGCGCGAACTCACGCTGCATTGCGCACTATGTATTGCCGAAGCAAAGCTTACAGCAAGCCGTCACTGTTATTTCGGCAGTT  >  CP009273/557447‑557735

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: