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,5920GT98.2% 689.7 / 3.3 219D814Y (GAC→TAC) sfmDputative outer membrane export usher protein
Reads supporting (aligned to +/- strand):  ref base G (2/2);  new base T (100/115);  total (102/117)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.57e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 221 total aligned reads displayed.

GACGCAAGGCGCGTTGGTTCGTGCCAATTTTGATACCCGCATAGGCGTGCGGGCGCTCATTACCGTTACCCAGGGCGGAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATGACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACTGGTTCAGTGGGGAGACGGCGCGAACTCACGCTGCATTGCGCACTATGTATTGCCGAAGCAAAGCTTACAGCAAGCCGTCACTGTTATTTCGGCAGTTT  >  CP009273/557447‑557736
                                                                                                                                                 |                                                                                                                                                
 aCGCAAGGCGCGTTGGTTCGTGCCAATTTTGATACCCGCATAGGCGTGCGGGCGCTCATTACCGTTACCCAGGGCGGAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATTACg                                                                                                                                               <  1:3483245/148‑1 (MQ=255)
     aaGGCGCGTTGGTTCGTGCCAATTTTGATACCCGCATAGGCGTGCGGGCGCTCATTACCGTTACCCAGGGCGGAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATTACGGGCa                                                                                                                                           >  1:2551757/1‑148 (MQ=255)
      aGGCGCGTTGGTTCGTGCCAATTTTGATACCCGCATAGGCGTGCGGGCGCTCATTACCGTTACCCAGGGCGGAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATTACGGGCaa                                                                                                                                          <  1:854252/148‑1 (MQ=255)
        gcgcGTTGGTTCGTGCCAATTTTGATACCCGCATAGGCGTGCGGGCGCTCATTACCGTTACCCAGGGCGGAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATTACGGGCAAGt                                                                                                                                        <  2:203861/148‑1 (MQ=255)
        gcgcGTTGGTTCGTGCCAATTTTGATACCCGCATAGGCGTGCGGGCGCTCATTACCGTTACCCAGGGCGGAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATTACGGGCAAGt                                                                                                                                        <  2:1771603/148‑1 (MQ=255)
        gcgcGTTGGTTCGTGCCAATTTTGATACCCGCATAGGCGTGCGGGCGCTCATTACCGTTACCCAGGGCGGAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATTACGGGCAAGt                                                                                                                                        <  2:711313/148‑1 (MQ=255)
        gcgcGTTGGTTCGTGCCAATTTTGATACCCGCATAGGCGTGCGGGCGCTCATTACCGTTACCCAGGGCGGAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATTACGGGCAAGt                                                                                                                                        <  1:41966/148‑1 (MQ=255)
          gcgTTGGTTCGTGCCAATTTTGATACCCGCATAGGCGTGCGGGCGCTCATTACCGTTACCCAGGGCGGAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATTACGGGCAAGttt                                                                                                                                      <  1:1409383/148‑1 (MQ=255)
          gcgTTGGTTCGTGCCAATTTTGATACCCGCATAGGCGTGCGGGCGCTCATTACCGTTACCCAGGGCGGAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATTACGGGCAAGttt                                                                                                                                      <  1:198224/148‑1 (MQ=255)
            gTTGGTTCGTGCCAATTTTGATACCCGCATAGGCGTGCGGGCGCTCATTACCGTTACCCAGGGCGGAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATTACGGGCAAGtttat                                                                                                                                    <  1:3167698/148‑1 (MQ=255)
             ttGGTTCGTGCCAATTTTGATACCCGCATAGGCGTGCGGGCGCTCATTACCGTTACCCAGGGCGGAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATTACGGGCAAGtttatt                                                                                                                                   >  1:2628997/1‑148 (MQ=255)
             ttGGTTCGTGCCAATTTTGATACCCGCATAGGCGTGCGGGCGCTCATTACCGTTACCCAGGGCGGAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATTACGGGCAAGtttatt                                                                                                                                   <  1:901682/148‑1 (MQ=255)
             ttGGTTCGTGCCAATTTTGATACCCGCATAGGCGTGCGGGCGCTCATTACCGTTACCCAGGGCGGAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATTACGGGCAAGtttatt                                                                                                                                   <  1:2762058/148‑1 (MQ=255)
               ggTTCGTGCCAATTTTGATACCCGCATAGGCGTGCGGGCGCTCATTACCGTTACCCAGGGCGGAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATTACGGGCAAGtttattta                                                                                                                                 <  2:558033/148‑1 (MQ=255)
               ggTTCGTGCCAATTTTGATACCCGCATAGGCGTGCGGGCGCTCATTACCGTTACCCAGGGCGGAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATTACGGGCAAGtttattta                                                                                                                                 <  1:1472821/148‑1 (MQ=255)
                gTTCGTGCCAATTTTGATACCCGCATAGGCGTGCGGGCGCTCATTACCGTTACCCAGGGCGGAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATTACGGGCAAGTTTATTTaa                                                                                                                                <  2:1956753/148‑1 (MQ=255)
                gTTCGTGCCAATTTTGATACCCGCATAGGCGTGCGGGCGCTCATTACCGTTACCCAGGGCGGAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATTACGGGCAAGTTTATTTaa                                                                                                                                >  2:2377187/1‑148 (MQ=255)
                gTTCGTGCCAATTTTGATACCCGCATAGGCGTGCGGGCGCTCATTACCGTTACCCAGGGCGGAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATTACGGGCAAGTTTATTTaa                                                                                                                                <  1:2493328/148‑1 (MQ=255)
                  tCGTGCCAATTTTGATACCCGCATAGGCGTGCGGGCGCTCATTACCGTTACCCAGGGCGGAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATTACGGGCAAGTTTATTTAAgt                                                                                                                              >  2:3433965/1‑148 (MQ=255)
                       ccAATTTTGATACCCGCATAGGCGTGCGGGCGCTCATTACCGTTACCCAGGGCGGAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGGTTACGGGCAAGTTTATTTTAgtggtgg                                                                                                                         >  1:3233890/1‑147 (MQ=255)
                          aTTTTGATACCCGCATAGGCGTGCGGGCGCTCATTACCGTTACCCAGGGCGGAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTCGGGTGGGTGATTACGGGCCAGTTTATTGAAGTGGTGCGcc                                                                                                                      >  2:2163407/1‑148 (MQ=255)
                          aTTTTGATACCCGCATAGGCGTGCGGGCGCTCATTACCGTTACCCAGGGCGGAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATTACGGGCAAGTTTATTTAAGTGGTGCGcc                                                                                                                      >  2:2163425/1‑148 (MQ=255)
                           ttttGATACCCGCATAGGCGTGCGGGCGCTCATTACCGTTACCCAGGGCGGAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATTACGGGCAAGTTTATTTAAGTGGTGCGCCa                                                                                                                     >  1:2828559/1‑148 (MQ=255)
                           ttttGATACCCGCATAGGCGTGCGGGCGCTCATTACCGTTACCCAGGGCGGAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATTACGGGCAAGTTTATTTAAGTGGTGCGCCa                                                                                                                     >  1:2827039/1‑148 (MQ=255)
                           ttttGATACCCGCATAGGCGTGCGGGCGCTCATTACCGTTACCCAGGGCGGAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATTACGGGCAAGTTTATTTAAGTGGTGCGCCa                                                                                                                     >  1:2827059/1‑148 (MQ=255)
                           ttttGATACCCGCATAGGCGTGCGGGCGCTCATTACCGTTACCCAGGGCGGAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATTACGGGCAAGTTTATTTAAGTGGTGCGCCa                                                                                                                     >  1:1830911/1‑148 (MQ=255)
                           ttttGATACCCGCATAGGCGTGCGGGCGCTCATTACCGTTACCCAGGGCGGAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATGACGGGCAAGTTTATTTAAGTGGTGCGCCa                                                                                                                     <  1:1955214/148‑1 (MQ=255)
                           ttttGAGACCCGCATAGGCGTGCGGGCGCTCATTACCGTTACCCAGGGCGGAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATTACGGGCAAGTTTATTTAAGTGGTGCGCCa                                                                                                                     <  2:2396181/148‑1 (MQ=255)
                             ttGATACCCGCATAGGCGTGCGGGCGCTCATTACCGTTACCCAGGGCGGAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATTACGGGCAAGTTTATTTAAGTGGTGCGCCAtt                                                                                                                   <  2:453951/148‑1 (MQ=255)
                               gATACCCGCATAGGCGTGCGGGCGCTCATTACCGTTACCCAGGGCGGAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATTACGGGCAAGTTTATTTAAGTGGTGCGCCATTGt                                                                                                                 <  1:486998/148‑1 (MQ=255)
                               gATACCCGCATAGGCGTGCGGGCGCTCATTACCGTTACCCAGGGCGGAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATTACGGGCAAGTTTATTTAAGTGGTGCGCCATTGt                                                                                                                 <  1:426476/148‑1 (MQ=255)
                                aTACCCGCATAGGCGTGCGGGCGCTCATTACCGTTACCCAGGGCGGAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATTACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTc                                                                                                                <  2:248652/148‑1 (MQ=255)
                                aTACCCGCATAGGCGTGCGGGCGCTCATTACCGTTACCCAGGGCGGAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATTACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTc                                                                                                                <  2:3489285/148‑1 (MQ=255)
                                aTACCCGCATAGGCGTGCGGGCGCTCATTACCGTTACCCAGGGCGGAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATTACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTc                                                                                                                >  2:2579790/1‑148 (MQ=255)
                                 tACCCGCATAGGCGTGCGGGCGCTCATTACCGTTACCCAGGGCGGAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATTACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCt                                                                                                               <  1:1802773/148‑1 (MQ=255)
                                   cccGCATAGGCGTGCGGGCGCTCATTACCGTTACCCAGGGCGGAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATTACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTgg                                                                                                             >  2:2281800/1‑148 (MQ=255)
                                    ccGCATAGGCGTGCGGGCGCTCATTACCGTTACCCAGGGCGGAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATTACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGt                                                                                                            >  1:182106/1‑148 (MQ=255)
                                    ccGCATAGGCGTGCGGGCGCTCATTAACGTTACCCAGGGCGGAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATTACTGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGt                                                                                                            >  2:3416112/1‑148 (MQ=255)
                                      gCATAGGCGTGCGGGCGCTCATTACCGTTACCCAGGTAGTAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTAACGGAATAACCAGTATGGTGGGTGATTACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGa                                                                                                          <  2:2759041/148‑1 (MQ=255)
                                        aTAGGCGTGCGGGCGCTCATTACCGTTACCCAGGGCGGAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATTACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAAt                                                                                                        <  2:40916/148‑1 (MQ=255)
                                         tAGGCGTGCGGGCGCTCATTACCGTTACCCAGGGCGGAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATTACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAAtt                                                                                                       >  2:2895942/1‑148 (MQ=255)
                                          aGGCGTGCGGGCGCTCATTACCGTTACCCAGGGCGGAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATTACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTa                                                                                                      >  2:3376434/1‑148 (MQ=255)
                                          aGGCGTGCGGGCGCTCATTACCGTTACCCAGGGCGGAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATTACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTa                                                                                                      >  2:910785/1‑148 (MQ=255)
                                          aGGCGTGCGGGCGCTCATTACCGTTACCCAGGGCGGAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTAACGGAATAACCAGTATGGTGGGTGATTACGGGCAAGTTTATTTAATTGGTGCGCCAGTGTCTGGTGAATTa                                                                                                      >  2:2710872/1‑148 (MQ=255)
                                           ggCGTGCGGGCGCTCATTACCGTTACCCAGGGCGGAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATTACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTAc                                                                                                     >  2:729319/1‑148 (MQ=255)
                                           ggCGTGCGGGCGCTCATTACCGTTACCCAGGGCGGAAAACCGGTGCCGTTTGGATCACAGGTAAGGGAAAACAGTACCGGAATAACCAGTATGGTTGGTGATTACGGGCAACTTTATTTAAATGGTGCGCCATTGTCTGGTGAAATaa                                                                                                     >  1:596144/1‑147 (MQ=255)
                                                gCGGGCGCTCATTACCGTTACCCAGGGCGGAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATTACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACTGGtt                                                                                                <  2:1167025/148‑1 (MQ=255)
                                                  gggCGCTCATTACCGTTACCCAGGGCGGAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATTACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTTAATTACTGGTTCa                                                                                              <  1:3383092/148‑1 (MQ=255)
                                                  gggCGCTCATTACCGTTACCCAGGGCGGAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATTACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACTGGTTCa                                                                                              <  1:1217989/148‑1 (MQ=255)
                                                   ggCGCTCATTACCGTTACCCAGGGCGGAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATTACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACTGGTTCAg                                                                                             <  2:2996494/148‑1 (MQ=255)
                                                         cATTACCGTTACCCAGGGCGGAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATTACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACTGGTTCAGTGGGga                                                                                       >  1:3272299/1‑148 (MQ=255)
                                                             aCCGTTACCCAGGGCGGAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATTACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACTGGTTCAGTGGGGAGACg                                                                                   <  2:1087833/148‑1 (MQ=255)
                                                                  tACCCAGGGCGGAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATTACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACTGGTTCAGTGGGGAGACGGcgcg                                                                              >  1:2767906/1‑148 (MQ=255)
                                                                      cAGGGCGGAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATTACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACTGGTTCAGTGGGGAGACGGCGCGAACt                                                                          <  2:846965/148‑1 (MQ=255)
                                                                       aGGGCGGAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAAGAACCAGTATGGTGGGTGATTACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACTTGTTCAGTGGGGAGACGGCGCGAACTc                                                                         <  2:3324791/148‑1 (MQ=255)
                                                                        gggCGGAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATTACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACTGGTTCAGTGGGGAGACGGCGCGAACTCa                                                                        >  1:1343716/1‑148 (MQ=255)
                                                                          gCGGAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATTACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACTGGTTCAGTGGGGAGACGGCGCGAACTCACg                                                                      <  2:182106/148‑1 (MQ=255)
                                                                             gAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATTACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACTGGTTCAGTGGGGAGACGGCGCGAACTCACGCTg                                                                   <  2:3074902/148‑1 (MQ=255)
                                                                              aaaaCCGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATTACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACTGGTTCAGTGGGGAGACGGCGCGAACTCACGCTGc                                                                  <  1:2646765/148‑1 (MQ=255)
                                                                                  ccGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATTACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACTGGTTCAGTGGGGAGACGGCGCGAACTCACGCTGCATTg                                                              >  2:1876391/1‑148 (MQ=255)
                                                                                  ccGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATTACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACTGGTTCAGTGGGGAGACGGCGCGAACTCACGCTGCATTg                                                              >  1:2359265/1‑148 (MQ=255)
                                                                                      tGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATTACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACTGGTTCAGTGGGGAGACGGCGCGAACTCACGCTGCATTGCGCa                                                          <  2:2881818/148‑1 (MQ=255)
                                                                                      tGCCGTTTGGATCACTGGTAAGGGAAAACAGTACCGGAAAAAACAGTATGGTGGGTGATTACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTAATGGTTCAGTGGGGAGACGGCGCGAACTCACGCTGCATTGCGCa                                                          >  1:2682310/1‑148 (MQ=255)
                                                                                       gCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATTACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACTGGTTCAGTGGGGAGACGGCGCGAACTCACGCTGCATTGCGCAc                                                         >  2:2597107/1‑148 (MQ=255)
                                                                                       gCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATTACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACTGGTTCAGTGGGGAGACGGCGCGAACTCACGCTGCATTGCGCAc                                                         >  2:659609/1‑148 (MQ=255)
                                                                                         cGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATTACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACTGGTTCAGTGGGGAGACGGCGCGAACTCACGCTGCATTGCGCACTa                                                       >  1:165390/1‑148 (MQ=255)
                                                                                          gTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATTACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACTGGTTCAGTGGGGAGACGGCGCGAACTCACGCTGCATTGCGCACTAt                                                      >  1:431956/1‑148 (MQ=255)
                                                                                            ttGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATTACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACTGGTTCAGTGGGGAGACGGCGCGAACTCACGCTGCATTGCGCACTATGt                                                    <  2:2774489/148‑1 (MQ=255)
                                                                                                aTCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATTACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACTGGTTCAGTGGGGAGACGGCGCGAACTCACGCTGCATTGCGCACTATGTATTg                                                <  2:3128679/148‑1 (MQ=255)
                                                                                                aTCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATTACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACTGGTTCAGTGGGGAGACGGCGCGAACTCACGCTGCATTGCGCACTATGTATTg                                                <  2:2827040/148‑1 (MQ=255)
                                                                                                aTCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATTACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACTGGTTCAGTGGGGAGACGGCGCGAACTCACGCTGCATTGCGCACTATGTATTg                                                <  2:2828559/148‑1 (MQ=255)
                                                                                                  cACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATTACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACTGGTTCAGTGGGGAGACGGCGCGAACTCACGCTGCATTGCGCACTATGTATTGcc                                              <  2:2373937/148‑1 (MQ=255)
                                                                                                   aCTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATTACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACTGGTTCAGTGGGGAGACGGCGCGAACTCACGCTGCATTGCGCAATATGTATTGCCg                                             <  1:2579353/148‑1 (MQ=255)
                                                                                                           gggAAAACAGTACCGGAATAACCAGTATGGTGGGTGATTACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACTGGTTCAGTGGGGAGACGGCGCGAACTCACGCTGCATTGCGCACTATGTATTGCCGAAGCAAAg                                     <  2:2628997/148‑1 (MQ=255)
                                                                                                            ggAAAACAGTACCGGAATAACCAGTATGGTGGGTGATTACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACTGGTTCAGTGGGGAGACGGCGCGAACTCACGCTGCATTGCGCACTATGTATTGCCGAAGCAAAGc                                    >  1:1716316/1‑148 (MQ=255)
                                                                                                            ggAAAACAGTACCGGAATAACCAGTATGGTGGGTGATTACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACTGGTTCAGTGGGGAGACGGCGCGAACTCACGCTGCATTGCGCACTATGTATTGCCGAAGCAAAGc                                    >  2:1027321/1‑148 (MQ=255)
                                                                                                             gAAAACAGTACCGGAATAACCAGTATGGTGGGTGATTACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACTGGTTCAGTGGGGAGACGGCGCGAACTCACGCTGCATTGCGCACTATGTATTGCCGAAGCAAAGCt                                   <  2:680555/148‑1 (MQ=255)
                                                                                                             gAAAACAGTACCGGAATAACCAGTATGGTGGGTGATTACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACTGGTTCAGTGGGGAGACGGCGCGAACTCACGCTGCATTGCGCACTATGTATTGCCGAAGCAAAGCt                                   <  2:289185/148‑1 (MQ=255)
                                                                                                              aaaaCAGTACCGGAATAACCAGTATGGTGGGTGATTACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACTGGTTCAGTGGGGAGACGGCGCGAACTCACGCTGCATTGCGCACTATGTATTGCCGAAGCAAAGCtt                                  >  2:3291862/1‑148 (MQ=255)
                                                                                                               aaaCAGTACCGGAATAACCAGTATGGTGGGTGATTACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACTGGTTCAGTGGGGAGACGGCGCGAACTCACGCTGCATTGCGCACTATGTATTGCCGAAGCAAAGCTTa                                 <  1:1893963/148‑1 (MQ=255)
                                                                                                                aaCAGTACCGGAATAACCAGTATGGTGGGTGATTACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACTGGTTCAGTGGGGAGACGGCGCGAACTCACGCTGCATTGCGCACTATGTATTGCCGAAGCAAAGCTTAc                                >  2:1899634/1‑148 (MQ=255)
                                                                                                                 aCAGTACCGGAATAACCAGTATGGTGGGTGATTACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACTGGTTCAGTGGGGAGACGGCGCGAACTCACGCTGCATTGCGCACTATGTATTGCCGAAGCAAAGCTTACa                               >  1:1001404/1‑148 (MQ=255)
                                                                                                                  cAGTACCGGAATAACCAGTATGGTGGGTGATTACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACTGGTTCAGTGGGGAGACGGCGCGAACTCACGCTGCATTGCGCACTATGTATTGCCGAAGCAAAGCTTACAg                              <  2:1178086/148‑1 (MQ=255)
                                                                                                                  cAGTACCGGAATAACCAGTATGGTGGGTGATTACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACTGGTTCAGTGGGGAGACGGCGCGAACTCACGCTGCATTGCGCACTATGTATTGCCGAAGCAAAGCTTACAg                              >  2:2740265/1‑148 (MQ=255)
                                                                                                                     tACCGGAATAACCAGTATGGTGGGTGATTACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACTGGTTCAGTGGGGAGACGGCGCGAACTCACGCTGCATTGCGCACTATGTATTGCCGAAGCAAAGCTTACAGCaa                           >  2:2544820/1‑148 (MQ=255)
                                                                                                                      aCCGGAATAACCAGTATGGTGGGTGATTACGGTCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACTGGTTCAGTGGGGAGACGGCGCGAACTCACGCTGCATTGCGCACTATGTATTGCCGAAGCAAAGCTTACAGCAAg                          <  1:1137721/148‑1 (MQ=255)
                                                                                                                      aCCGGAATAACCAGTATGGTGGGTGATTACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACTGGTTCAGTGGGGAGACGGCGCGAACTCACGCTGCATTGCGCACTATGTATTGCCGAAGCAAAGCTTACAGCAAg                          <  2:1236842/148‑1 (MQ=255)
                                                                                                                      aCCGGAATAACCAGTATGGTGGGTGATTACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACTGGTTCAGTGGGGAGACGGCGCGAACTCACGCTGCATTGCGCACTATGTATTGCCGAAGCAAAGCTTACAGCAAg                          <  2:545799/148‑1 (MQ=255)
                                                                                                                            aTAACCAGTATGGTGGGTGATTACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACTGGTTCAGTGGGGAGACGGCGCGAACTCACGCTGCATTGCGCACTATGTATTGCCGAAGCAAAGCTTACAGCAAGCCGTCa                    <  2:2551807/148‑1 (MQ=255)
                                                                                                                            aTAACCAGTATGGTGGGTGATTACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACTGGTTCAGTGGGGAGACGGCGCGAACTCACGCTGCATTGCGCACTATGTATTGCCGAAGCAAAGCTTACAGCAAGCCGTCa                    <  1:3115307/148‑1 (MQ=255)
                                                                                                                             tAACCAGTATGGTGGGTGATTACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACTGGTTCAGTGGGGAGACGGCGCGAACTCACGCTGCATTGCGCACTATGTATTGCCGAAGCAAAGCTTACAGCAAGCCGTCAc                   >  1:1532779/1‑148 (MQ=255)
                                                                                                                              aaCCAGTATGGTGGGTGATTACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACTGGTTCAGTGGGGAGACGGCGCGAACTCACGCTGCATTGCGCACTATGTATTGCCGAAGCAAAGCTTACAGCAAGCCGTCACt                  >  1:2616304/1‑148 (MQ=255)
                                                                                                                              aaCCAGTATGGTGGGTGATTACGGGCAAGTTTATTTAAGTGGTGCGCCATAGTCTGGTGAATTACTGGTTCAGTGGGGAGACGGCGCGAACTCACGCTGCATTGCGCACTATGTATTGCCGAAGCAAATCTTACAGCAAGCCGTCACt                  >  1:2987527/1‑148 (MQ=255)
                                                                                                                               aCCAGTATGGTGGGTGATTACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACTGGTTCAGTGGGGAGACGGCGCGAACTCACGCTGCATTGCGCACTATGTATTGCCGAAGCAAAGCTTACAGCAAGCCGTCACTg                 <  2:139161/148‑1 (MQ=255)
                                                                                                                                ccAGTATGGTGGGTGATTACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACTGGTTCAGTGGGGAGACGGCGCGAACTCACGCTGCATTGCGCACTATGTATTGCCGAAGCAAAGCTTACAGCAAGCCGTCACTGt                >  1:380791/1‑148 (MQ=255)
                                                                                                                                ccAGTATGGTGGGTGATTACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACTGGTTCAGTGGGGAGACGGCGCGAACTCACGCTGCATTGCGCACTATGTATTGCCGAAGCAAAGCTTACAGCAAGCCGTCACTGt                <  1:968576/148‑1 (MQ=255)
                                                                                                                                    tATGGTGGGTGATTACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACTGGTTCAGTGGGGAGACGGCGCGAACTCACGCTGCATTGCGCACTATGTATTGCCGAAGCAAAGCTTACAGCAAGCCGTCACTGTTAtt            >  1:3493146/1‑148 (MQ=255)
                                                                                                                                       ggtgggtgATTACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACTGGTTCAGTGGGGAGACGGCGCGAACTCACGCTGCATTGCGCACTATGTATTGCCGAAGCAAAGCTTACAGCAAGCCGTCACTGTTATTTCg         >  1:1513450/1‑148 (MQ=255)
                                                                                                                                             tgattaCGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACTGGTTCAGTGGGGAGACGGCGCGAACTCACGCTGCATTGCGCACTATGTATTGCCGAAGCAAAGCTTACAGCAAGCCGTCACTGTTATTTCGGCAGtt   <  1:2921622/148‑1 (MQ=255)
                                                                                                                                              gattaCGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACTGGTTCAGTGGGGAGACGGCGCGAACTCACGCTGCATTGCGCACTATGTATTGCCGAAGCAAAGCTTACAGCAAGCCGTCACTGTTATTTCGGCAGttt  >  2:209822/1‑148 (MQ=255)
                                                                                                                                                 |                                                                                                                                                
GACGCAAGGCGCGTTGGTTCGTGCCAATTTTGATACCCGCATAGGCGTGCGGGCGCTCATTACCGTTACCCAGGGCGGAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATGACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACTGGTTCAGTGGGGAGACGGCGCGAACTCACGCTGCATTGCGCACTATGTATTGCCGAAGCAAAGCTTACAGCAAGCCGTCACTGTTATTTCGGCAGTTT  >  CP009273/557447‑557736

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: