Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 557,592 G→T 88.0% 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,5920GT88.0% 460.0 / 55.4 188D814Y (GAC→TAC) sfmDputative outer membrane export usher protein
Reads supporting (aligned to +/- strand):  ref base G (12/11);  new base T (75/90);  total (87/101)
Fisher's exact test for biased strand distribution p-value = 6.57e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 188 total aligned reads displayed.

GACGCAAGGCGCGTTGGTTCGTGCCAATTTTGATACCCGCATAGGCGTGCGGGCGCTCATTACCGTTACCCAGGGCGGAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATGACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACTGGTTCAGTGGGGAGACGGCGCGAACTCACGCTGCATTGCGCACTATGTATTGCCGAAGCAAAGCTTACAGCAAGCCGTCACTGTTATTTCGGCAGT  >  CP009273/557447‑557734
                                                                                                                                                 |                                                                                                                                              
gACGCAAGGCGCGTTGGTTCGTGCCAATTTTGATACCCGCATAGGCGTGCGGGCGCTCATTACCGTTACCCAGGGCGGAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATTAc                                                                                                                                              <  1:1027887/148‑1 (MQ=255)
gACGCAAGGCGCGTTGGTTCGTGCCAATTTTGATACCCGCATAGGCGTGCGGGCGCTCATTACCGTTACCCAGGGCGGAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATTAc                                                                                                                                              >  1:183559/1‑148 (MQ=255)
 aCGCAAGGCGCGTTGGTTCGTGCCAATTTTGATACCCGCATAGGCGTGCGGGCGCTCATTACCGTTACCCAGGGCGGAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATTACg                                                                                                                                             <  2:681610/148‑1 (MQ=255)
   gCAAGGCGCGTTGGTTCGTGCCAATTTTGATACCCGCATAGGCGTGCGGGCGCTCATTACCGTTACCCAGGGCGGAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATTACggg                                                                                                                                           <  2:2141866/148‑1 (MQ=255)
    cAAGGCGCGTTGGTTCGTGCCAATTTTGATACCCGCATAGGCGTGCGGGCGCTCATTACCGTTACCCAGGGCGGAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATGACGGGc                                                                                                                                          >  2:688285/1‑148 (MQ=255)
     aaGGCGCGTTGGTTCGTGCCAATTTTGATACCCGCATAGGCGTGCGGGCGCTCATTACCGTTACCCAGGGCGGAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATTACGGGCa                                                                                                                                         <  1:863786/148‑1 (MQ=255)
       ggCGCGTTGGTTCGTGCCAATTTTGATACCCGCATAGGCGTGCGGGCGCTCATTACCGTTACCCAGGGCGGAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATTACGGGCAAg                                                                                                                                       >  2:949280/1‑148 (MQ=255)
       ggCGCGTTGGTTCGTGCCAATTTTGATACCCGCATAGGCGTGCGGGCGCTCATTACCGTTACCCAGGGCGGAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATTACGGGCAAg                                                                                                                                       <  1:3257102/148‑1 (MQ=255)
        gcgcGTTGGTTCGTGCCAATTTTGATACCCGCATAGGCGTGCGGGCGCTCATTACCGTTACCCAGGGCGGAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATTACGGGCAAGt                                                                                                                                      <  2:3236189/148‑1 (MQ=255)
        gcgcGTTGGTTCGTGCCAATTTTGATACCCGCATAGGCGTGCGGGCGCTCATTACCGTTACCCAGGGCGGAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATTACGGGCAAGt                                                                                                                                      <  1:1984945/148‑1 (MQ=255)
        gcgcGTTGGTTCGTGCCAATTTTGATACCCGCATAGGCGTGCGGGCGCTCATTACCGTTACCCAGGGCGGAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATTACGGGCAAGt                                                                                                                                      >  1:2407288/1‑148 (MQ=255)
           cgTTGGTTCGTGCCAATTTTGATACCCGCATAGGCGTGCGGGCGCTCATTACCGTTACCCAGGGCGGAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATTACGGGCAAGttta                                                                                                                                   >  1:2810418/1‑148 (MQ=255)
           cgTTGGTTCGTGCCAATTTTGATACCCGCATAGGCGTGCGGGCGCTCATTACCGTTACCCAGGGCGGAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATGACGGGCAAGttta                                                                                                                                   >  2:3082255/1‑148 (MQ=255)
            gTTGGTTCGTGCCAATTTTGATACCCGCATAGGCGTGCGGGCGCTCATTACCGTTACCCAGGGCGGAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATTACGGGCAAGtttat                                                                                                                                  <  1:2556301/148‑1 (MQ=255)
             ttGGTTCGTGCCAATTTTGATACCCGCATAGGCGTGCGGGCGCTCATTACCGTTACCCAGGGCGGAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATTACGGGCAAGtttatt                                                                                                                                 >  1:1217855/1‑148 (MQ=255)
             ttGGTTCGTGCCAATTTTGATACCCGCATAGGCGTGCGGGCGCTCATTACCGTTACCCAGGGCGGAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATTACGGGCAAGtttatt                                                                                                                                 >  2:1117256/1‑148 (MQ=255)
             ttGGTTCGTGCCAATTTTGATACCCGCATAGGCGTGCGGGCGCTCATTACCGTTACCCAGGGCGGAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATGACGGGCAAGtttatt                                                                                                                                 >  1:2176623/1‑148 (MQ=255)
               ggTTCGTGCCAATTTTGATACCCGCATAGGCGTGCGGGCGCTCATTACCGTTACCCAGGGCGGAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATGACGGGCAAGtttattta                                                                                                                               <  1:55177/148‑1 (MQ=255)
                 ttCGTGCCAATTTTGATACCCGCATAGGCGTGCGGGCGCTCATTACCGTTACCCAGGGCGGAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATTACGGGCAAGTTTATTTAAg                                                                                                                             <  2:2203632/148‑1 (MQ=255)
                  tCGTGCCAATTTTGATACCCGCATAGGCGTGCGGGCGCTCATTACCGTTACCCAGGGCGGAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATTACGGGCAAGTTTATTTAAgt                                                                                                                            <  2:2529758/148‑1 (MQ=255)
                  tCGTGCCAATTTTGATACCCGCATAGGCGTGCGGGCGCTCATTACCGTTACCCAGGGCGGAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATTACGGGCAAGTTTATTTAAgt                                                                                                                            >  1:2694160/1‑148 (MQ=255)
                      gCCAATTTTGATACCCGCATAGGCGTGCGGGCGCTCATTACCGTTACCCAGGGCGGAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATTACGGGCAAGTTTATTTAAgtggtg                                                                                                                        >  2:2391827/1‑148 (MQ=255)
                      gCCAATTTTGATACCCGCATAGGCGTGCGGGCGCTCATTACCGTTACCCAGGGCGGAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATTACGGGCAAGTTTATTTAAgtggtg                                                                                                                        >  2:732979/1‑148 (MQ=255)
                           ttttGATACCCGCATAGGCGTGCGGGCGCTCATTACCGTTACCCAGGGCGGAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATTACGGGCAAGTTTATTTAAGTGGTGCGCCa                                                                                                                   <  2:3120045/148‑1 (MQ=255)
                           ttttGATACCCGCATAGGCGTGCGGGCGCTCATTACCGTTACCCAGGGCGGAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATTACGGGCAAGTTTATTTAAGTGGTGCGCCa                                                                                                                   <  1:1210027/148‑1 (MQ=255)
                           ttttGATACCCGCATAGGCGTGCGGGCGCTCATTACCGTTACCCAGGGCGGAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATGACGGGCAAGTTTATTTAAGTGGTGCGCCa                                                                                                                   <  1:672017/148‑1 (MQ=255)
                            tttGATACCCGCATAGGCGTGCGGGCGCTCATTACCGTTACCCAGGGCGGAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATTACGGGCAAGTTTATTTAAGTGGTGCGCCAt                                                                                                                  <  2:2594363/148‑1 (MQ=255)
                             ttGATACCCGCATAGGCGTGCGGGCGCTCATTACCGTTACCCAGGGCGGAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATTACGGGCAAGTTTATTTAAGTGGTGCGCCAtt                                                                                                                 <  2:3035878/148‑1 (MQ=255)
                             ttGATACCCGCATAGGCGTGCGGGCGCTCATTACCGTTACCCAGGGCGGAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATTACGGGCAAGTTTATTTAAGTGGTGCGCCAtt                                                                                                                 <  1:313166/148‑1 (MQ=255)
                              tGATACCCGCATAGGCGTGCGGGCGCTCATTACCGTTACCCAGGGCGGAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATTACGGGCAAGTTTATTTAAGTGGTGCGCCATTg                                                                                                                <  2:2141421/148‑1 (MQ=255)
                              tGATACCCGCATAGGCGTGCGGGCGCTCATTACCGTTACCCAGGGCGGAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATGACGGGCAAGTTTATTTAAGTGGTGCGCCATTg                                                                                                                <  1:966761/148‑1 (MQ=255)
                               gATACCCGCATAGGCGTGCGGGCGCTCATTACCGTTACCCAGGGCGGAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATTACGGGCAAGTTTATTTAAGTGGTGCGCCATTGt                                                                                                               <  2:2746888/148‑1 (MQ=255)
                                aTACCCGCATAGGCGTGCGGGCGCTCATTACCGTTACCCAGGGCGGAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATGACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTc                                                                                                              >  1:1374599/1‑148 (MQ=255)
                                 tACCCGCATAGGCGTGCGGGCGCTCATTACCGTTACCCAGGGCGGAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATTACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCt                                                                                                             <  1:27180/148‑1 (MQ=255)
                                 tACCCGCATAGGCGTGCGGGCGCTCATTACCGTTACCCAGGGCGGAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATTACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCt                                                                                                             <  2:781769/148‑1 (MQ=255)
                                  aCCCGCATAGGCGTGCGGGCGCTCATTACCGTTACCCAGGGCGGAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGGATGGTGGGTGATTACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTg                                                                                                            >  2:1419437/1‑148 (MQ=255)
                                   cccGCATAGGCGTGCGGGCGCTCATTACCGTTACCCAGGGCGGAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATTACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTgg                                                                                                           <  2:606688/148‑1 (MQ=255)
                                    ccGCATAGGCGTGCGGGCGCTCATTACCGTTACCCAGGGCGGAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATTACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGt                                                                                                          <  1:2821575/148‑1 (MQ=255)
                                    ccGCATAGGCGTGCGGGCGCTCATTACCGTTACCCAGGGCGGAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATTACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGt                                                                                                          <  1:1700566/148‑1 (MQ=255)
                                     cGCATAGGCGTGCGGGCGCTCATTACCGTTACCCAGGGCGGAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGCATAACCCGTATGGTGGGTGATTACGGGGAAGTTTATTTAATTGGTGCGCCCTTGTCTGGTg                                                                                                         >  1:303155/1‑148 (MQ=255)
                                     cGCATAGGCGTGCGGGCGCTCATTACCGTTACCCAGGGCGGAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATTACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTg                                                                                                         >  1:188394/1‑148 (MQ=255)
                                     cGCATAGGCGTGCGGGCGCTCATTACCGTTAACCAGGGCGGAAAACCGGTGCCGTTTTGATAACTGGGAAGGGAAAAAAGTTCCGGAAAAACCAGTATGGTGGGTGATTAACGGCAAGGTTATTTAAGTTGTGCGCCATTGTgtggtg                                                                                                         >  1:304137/1‑148 (MQ=255)
                                        aTAGGCGTGCGGGCGCTCATTACCGTTACCCAGGGCGGAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATTACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAAt                                                                                                      >  1:2292864/1‑148 (MQ=255)
                                         tAGGCGTGCGGGCGCTCATTACCGTTACCCAGGGCGGAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATTACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAAtt                                                                                                     >  2:375757/1‑148 (MQ=255)
                                         tAGGCGTGCGGGCGCTCATTACCGTTACCCAGGGCGGAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATGACGGGCAAGTTTATTTAAGTGGTGCGCCCTTGTCAGGTGAAtt                                                                                                     >  1:291559/1‑148 (MQ=255)
                                         tAGGCGTGCGGGCGCTCATTACCGTTACCCAGGGCGGAAAACAGGTGCCGTTTGGCTCACTGGTACGGTAAATCAGTTCCGTAATAACCATTTTTGTTGGTGATTACGGGCAAGTTTTTTTAAGTGGTGCGCCATTTTCTGGTgtttt                                                                                                     <  2:1680057/148‑5 (MQ=255)
                                          aGGCGTTCGTGCGCTCATTACCGTTACCCAGAGCGGAAAACCCGTTCCGTTTGGATCACTGGTACGGGAAAACATTACAGCAATAACCAGTAATGAGGGAGATGACGGTCAAGATTATTAAATTGGAGCGAAATTGTCTGGTGAATTa                                                                                                    >  2:1327489/1‑148 (MQ=255)
                                          aGGCGTGCGGGCGCTCATTACCGTTACCCAGGGCGGAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATTACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTa                                                                                                    <  1:2256888/148‑1 (MQ=255)
                                          aGGCGTGCGGGCGCTCATTACCGTTACCCAGGGCGGAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATGACGGGCAAGGTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTa                                                                                                    >  2:2823492/1‑148 (MQ=255)
                                           ggCGTGCGGGCGCTCATTACCGTTACCCAGGGCGGAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATTACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTAc                                                                                                   <  1:1246545/148‑1 (MQ=255)
                                             cGTGCGGGCGCTCATTACCGTTACCCAGGGCGGAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATTACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACTg                                                                                                 <  2:2694160/148‑1 (MQ=255)
                                              gTGCGGGCGCTCATTACCGTTACCCAGGGCGGAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATTACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACTgg                                                                                                <  2:1097917/148‑1 (MQ=255)
                                               tGCGGGCGCTCATTACCGTTACCCAGGGCGGAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATTACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACTGGt                                                                                               <  2:841169/148‑1 (MQ=255)
                                                gCGGGCGCTCATTACCGTTACCCAGGGCGGAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATTACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACTGGtt                                                                                              <  2:1138234/148‑1 (MQ=255)
                                                   ggCGCTCATTACCGTTACCCAGGGCGGAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATGACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACTGGTTCAg                                                                                           >  1:3083100/1‑148 (MQ=255)
                                                          aTTACCGTTACCCAGGGCGGAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATTACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACTGGTTCAGTGGGgag                                                                                    <  2:2854600/148‑1 (MQ=255)
                                                          aTTACCGTTACCCAGGGCGGAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATGACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACTGGTTCAGTGGGgag                                                                                    <  2:1374599/148‑1 (MQ=255)
                                                            tACCGTTACCCCGGGCGGAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATTACGGGCAAGTTTAGTTAAGTGGTGCGCCATTGTCTGGTGAATTACTGGTTCAGTGGGGAGAc                                                                                  >  2:242243/1‑148 (MQ=255)
                                                            tACCGTTACCCAGGGCGGAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATTACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACTGGTTCAGTGGGGAGAc                                                                                  >  2:629686/1‑148 (MQ=255)
                                                            tACCGTTACCCAGGGCGGAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATTACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACTGGTTCAGTGGGGAGAc                                                                                  >  2:1782647/1‑148 (MQ=255)
                                                               cGTTACCCAGGGCGGAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATTACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACTGGTTCAGTGGGGAGACGgc                                                                               >  2:1952143/1‑148 (MQ=255)
                                                                     ccAGGGCGGAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATTACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACTGGTTCAGTGGGGAGACGGCGCGAAc                                                                         <  2:3127057/148‑1 (MQ=255)
                                                                     ccAGGGCGGAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATTACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACTGGTTCAGTGGGGAGACGGCGCGAAc                                                                         >  1:767121/1‑148 (MQ=255)
                                                                      cAGGGCGGAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATTACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACTGGTTCAGTGGGGAGACGGCGCGAACt                                                                        >  1:1725667/1‑148 (MQ=255)
                                                                      cAGGGCGGAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATTACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACTGGTTCAGTGGGGAGACGGCGCGAACt                                                                        >  2:2666057/1‑148 (MQ=255)
                                                                         ggCGGAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATTACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACTGGTTCAGTGGGGAGACGGCGCGAACTCAc                                                                     >  1:2394690/1‑148 (MQ=255)
                                                                          gCGGAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATTACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACTGGTTCAGTGGGGAGACGGCGCGAACTCACg                                                                    <  2:508466/148‑1 (MQ=255)
                                                                              aaaaCCGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATTACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACTGGTTCAGTGGGGAGACGGCGCGAACTCACGCTGc                                                                <  2:1393181/148‑1 (MQ=255)
                                                                                  ccGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATTACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACTGGTTCAGTGGGGAGACGGCGCGAACTCACGCTGCATTg                                                            >  2:1467575/1‑148 (MQ=255)
                                                                                         cGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATTACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACTGGTTCAGTGGGGAGACGGCGCGAACTCACGCTGCATTGCGCACTa                                                     <  1:2210611/148‑1 (MQ=255)
                                                                                            ttGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATTACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACTGGTTCAGTGGGGAGACGGCGCGAACTCACGCTGCATTGCGCACTATGt                                                  >  1:2225610/1‑148 (MQ=255)
                                                                                             tGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATTACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACTGGTTCAGTGGGGAGACGGCGCGAACTCACGCTGCATTGCGCACTATGTa                                                 <  2:2111583/148‑1 (MQ=255)
                                                                                              ggATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATTACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACTGGTTCAGTGGGGAGACGGCGCGAACTCACGCTGCATTGCGCACTATGTAt                                                <  1:949280/148‑1 (MQ=255)
                                                                                                aTCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATTACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACTGGTTCAGTGGGGAGACGGCGCGAACTCACGCTGCATTGCGCACTATGTATTg                                              <  2:2269905/148‑1 (MQ=255)
                                                                                                aTCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATTACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACTGGTTCAGTGGGGAGACGGCGCGAACTCACGCTGCATTGCGCACTATGTATTg                                              <  2:609221/148‑1 (MQ=255)
                                                                                                 tCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATTACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACTGGTTCAGTGGGGAGACGGCGCGAACTCACGCTGCATTGCGCACTATGTATTGc                                             >  1:1006668/1‑148 (MQ=255)
                                                                                                 tCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATTACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACTGGTTCAGTGGGGAGACGGCGCGAACTCACGCTGCATTGCGCACTATGTATTGc                                             >  1:192576/1‑148 (MQ=255)
                                                                                                         aCGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATTACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACTGGTTCAGTGGGGAGACGGCGCGAACTCACGCTGCATTGCGCACTATGTATTGCCGAAGCaa                                     <  1:360459/148‑1 (MQ=255)
                                                                                                          cGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATTACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACTGGTTCAGTGGGGAGACGGCGCGAACTCACGCTGCATTGCGCACTATGTATTGCCGAAGCaaa                                    <  2:2326294/148‑1 (MQ=255)
                                                                                                          cGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATTACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACTGGTTCAGTGGGGAGACGGCGCGAACTCACGCTGCATTGCGCACTATGTATTGCCGAAGCaaa                                    >  1:1319361/1‑148 (MQ=255)
                                                                                                           gggAAAACAGTACCGGAATAACCAGTATGGTGGGTGATTACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACTGGTTCAGTGGGGAGACGGCGCGAACTCACGCTGCATTGCGCACTATGTATTGCCGAAGCAAAg                                   <  2:1046297/148‑1 (MQ=255)
                                                                                                           gggAAAACAGTACCGGAATAACCAGTATGGTGGGTGATTACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACTGGTTCAGTGGGGAGACGGCGCGAACTCACGCTGCATTGCGCACTATGTATTGCCGAAGCAAAg                                   <  2:1919953/148‑1 (MQ=255)
                                                                                                             gAAAACAGTACCGGAATAACCAGTATGGTGGGTGATTACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACTGGTTCAGTGGGGAGACGGCGCGAACTCACGCTGCATTGCGCACTATGTATTGCCGAAGCAAAGCt                                 >  1:434414/1‑148 (MQ=255)
                                                                                                                aaCAGTACCGGAATAACCAGTATGGTGGGTGATTACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACTGGTTCAGTGGGGAGACGGCGCGAACTCACGCTGCATTGCGCACTATGTATTGCCGAAGCAAAGCTTAc                              <  2:2637549/148‑1 (MQ=255)
                                                                                                                     tACCGGAATAACCAGTATGGTGGGTGATGACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACTGGTTCAGTGGGGAGACGGCGCGAACTCACGCTGCATTGCGCACTATGTATTGCCGAAGCAAAGCTTACAGCaa                         >  1:214748/1‑148 (MQ=255)
                                                                                                                      aCCGGAATAACCAGTATGGTGGGTGATTACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACTGGTTCAGTGGGGAGACGGCGCGAACTCACGCTGCATTGCGCACTATGTATTGCCGAAGCAAAGCTTACAGCAAg                        <  1:1521817/148‑1 (MQ=255)
                                                                                                                      aCCGGAATAACCAGTATGGTGGGTGATTACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACTGGTTCAGTGGGGAGACGGCGCGAACTCACGCTGCATTGCGCACTATGTATTGCCGAAGCAAAGCTTACAGCAAg                        <  1:2696754/148‑1 (MQ=255)
                                                                                                                       ccGGAATAACCAGTATGGTGGGTGATTACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACTGGTTCAGTGGGGAGACGGCGCGAACTCACGCTGCATTGCGCACTATGTATTGCCGAAGCAAAGCTTACAGCAAGc                       <  2:1445298/148‑1 (MQ=255)
                                                                                                                        cGGAATAACCAGTATGGTGGGTGATTACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACTGGTTCAGTGGGGAGACGGCGCGAACTCACGCTGCATTGCGCACTATGTATTGCCGAAGCAAAGCTTACAGCAAGcc                      >  2:2932569/1‑148 (MQ=255)
                                                                                                                          gAATAACCAGTATGGTGGGTGATTACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACTGGTTCAGTGGGGAGACGGCGCGAACTCACGCTGCATTGCGCACTATGTATTGCCGAAGCAAAGCTTACAGCAAGCCGt                    <  2:1217855/148‑1 (MQ=255)
                                                                                                                          gAATAACCAGTATGGTGGGTGATTACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACTGGTTCAGTGGGGAGACGGCGCGAACTCACGCTGCATTGCGCACTATGTATTGCCGAAGCAAAGCTTACAGCAAGCCGt                    >  1:1606332/1‑148 (MQ=255)
                                                                                                                          gAATAACCAGTATGGTGGGTGATTACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACTGGTTCAGTGGGGAGACGGCGCGAACTCACGCTGCATTGCGCACTATGTATTGCCGAAGCAAAGCTTACAGCAAGCCGt                    >  1:1095612/1‑148 (MQ=255)
                                                                                                                          gAATAACCAGTATGGTGGGTGATTACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACTGGTTCAGTGGGGAGACGGCGCGAACTCACGCTGCATTGCGCACTATGTATTGCCGAAGCAAAGCTTACAGCAAGCCGt                    <  1:2983027/148‑1 (MQ=255)
                                                                                                                             tAACCAGTATGGTGGGTGATTACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACTGGTTCAGTGGGGAGACGGCGCGAACTCACGCTGCATTGCGCACTATGTATTGCCGAAGCAAAGATTACAGCAAGCCGTCAc                 >  2:2106006/1‑148 (MQ=255)
                                                                                                                               aCCAGTATGGTGGGTGATTACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACTGGTTCAGTGGGGAGACGGCGCGAACTCACGCTGCATTGCGCACTATGTATTGCCGAAGCAAAGCTTACAGCAAGCCGTCACTg               >  2:2880838/1‑148 (MQ=255)
                                                                                                                                 cAGTATGGTGGGTGATTACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACTGGTTCAGTGGGGAGACGGCGCGAACTCACGCTGCATTGCGCACTATGTATTGCCGAAGCAAAGCTTACAGCAAGCCGTCACTGtt             <  1:324746/148‑1 (MQ=255)
                                                                                                                                 cAGTATGGTGGGTGATTACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACTGGTTCAGTGGGGAGACGGCGCGAACTCACGCTGCATTGCGCACTATGTATTGCCGAAGCAAAGCTTACAGCAAGCCGTCACTGtt             <  1:582977/148‑1 (MQ=255)
                                                                                                                                 cAGTATGGTGGGTGATTACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACTGGTTCAGTGGGGAGACGGCGCGAACTCACGCTGCATTGCGCACTATGTATTGCCGAAGCAAAGCTTACAGCAAGCCGTCACTGtt             <  1:772156/148‑1 (MQ=255)
                                                                                                                                            gtgATTACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACTGGTTCAGTGGGGAGACGGCGCGAACTCACGCTGCATTGCGCACTATGTATTGCCGAAGCAAAGCTTACAGCAAGCCGTCACTGTTATTTCGGCAGt  >  2:58637/1‑148 (MQ=255)
                                                                                                                                            gtgATGACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACTGGTTCAGTGGGGAGACGGCGCGAACTCACGCTGCATTGCGCACTATGTATTGCCGAAGCAAAGCTTACAGCAAGCCGTCACTGTTATTTCGGCAGt  >  1:3178345/1‑148 (MQ=255)
                                                                                                                                                 |                                                                                                                                              
GACGCAAGGCGCGTTGGTTCGTGCCAATTTTGATACCCGCATAGGCGTGCGGGCGCTCATTACCGTTACCCAGGGCGGAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATGACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACTGGTTCAGTGGGGAGACGGCGCGAACTCACGCTGCATTGCGCACTATGTATTGCCGAAGCAAAGCTTACAGCAAGCCGTCACTGTTATTTCGGCAGT  >  CP009273/557447‑557734

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: