Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,129,816 G→T 7.0% A187S (GCC→TCC)  flgF → flagellar component of cell‑proximal portion of basal‑body rod

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,129,8160GT 7.0% 1106.9 / 40.9 301A187S (GCC→TCC) flgFflagellar component of cell‑proximal portion of basal‑body rod
Reads supporting (aligned to +/- strand):  ref base G (135/144);  new base T (6/15);  total (142/159)
Fisher's exact test for biased strand distribution p-value = 1.11e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 308 total aligned reads displayed.

TGCTGCCGATGGCACAATCTCGGCGCTCAATCCGGGCGATCCGGCAAATACGGTTGCGCCAGTAGGGCGTCTTAAACTGGTGAAAGCCACGGGCAGCGAAGTGCAGCGCGGTGACGACGGCATTTTTCGTTTAAGCGCAGAAACCCAGGCCACGCGTGGGCCGGTACTGCAGGCAGATCCAACCTTGCGTGTGATGTCGGGGGTTCTGGAAGGCAGTAACGTCAATGCCGTTGCGGCAATGAGCGACATGATTGCCAGCGCGCGGCGTTTTGAAATGCAGATGAAGGTGATCAGCAGC  >  CP009273/1129668‑1129965
                                                                                                                                                    |                                                                                                                                                     
 gctgcCGATGGCACAATCTCGGCGCTCAATCCGGGCGATCCGGCAAATACGGTTGCGCCAGTAGGGCGTCTTAAACTGGTGAAAGCCACGGGCAGCGAAGTGCAGCGCGGTGACGACGGCATTTTTCGTTTAAGCGCAGAAACCCAGGcc                                                                                                                                                     >  1:3865932/1‑150 (MQ=255)
  ctgcCGATGGCACAATCTCGGCGCTCAATCCGGGCGATCCGGCAAATACGGTTGCGCCAGTAGGGCGTCTTAAACTGGTGAAAGCCACGGGCAGCGAAGTGCAGCGCGGTGACGACGGCATTTTTCGTTTAAGCGCAGAAACCCAGTCCa                                                                                                                                                    >  2:391317/1‑150 (MQ=255)
  ctgcCGATGGCACAATCTCGGCGCTCAATCCGGGCGATCCGGCAAATACGGTTGCGCCAGTAGGGCGTCTTAAACTGGTGAAAGCCACGGGCAGCGAAGTGCAGCGCGGTGACGACGGCATTTTTCGTTTAAGCGCAGAAACCCAGTCCa                                                                                                                                                    >  2:3821364/1‑150 (MQ=255)
  ctgcCGATGGCACAATCTCGGCGCTCAATCCGGGCGATCCGGCAAATACGGTTGCGCCAGTAGGGCGTCTTAAACTGGTGAAAGCCACGGGCAGCGAAGTGCAGCGCGGTGACGACGGCATTTTTCGTTTAAGCGCAGAAACCCAGGCCa                                                                                                                                                    >  2:1312738/1‑150 (MQ=255)
   tgcCGATGGCACAATCTCGGCGCTCAATCCGGGCGATCCGGCAAATACGGTTGCGCCAGTAGGGCGTCTTAAACTGGTGAAAGCCACGGGCAGCGAAGTGCAGCGCGGTGACGACGGCATTTTTCGTTTAAGCGCAGAAACCCAGGCCAc                                                                                                                                                   <  1:3169872/150‑1 (MQ=255)
       gATGGCACAATCTCGGCGCTCAATCCGGGCGATCCGGCAAATACGGTTGCGCCAGTAGGGCGTCTTAAACTGGTGAAAGCCACGGGCAGCGAAGTGCAGCGCGGTGACGACGGCATTTTTCGTTTAAGCGCAGAAACCCAGGCCACGCGt                                                                                                                                               >  1:2303742/1‑150 (MQ=255)
        aTGGCACAATCTCGGCGCTCAATCCGGGCGATCCGGCAAATACGGTTGCGCCAGTAGGGCGTCTTAAACTGGTGAAAGCCACGGGCAGCGAAGTGCAGCGCGGTGACGACGGCATTTTTCGTTTAAGCGCAGAAACCCAGGCCACGCGTg                                                                                                                                              >  2:5005126/1‑150 (MQ=255)
         tGGCACAATCTCGGCGCTCAATCCGGGCGATCCGGCAAATACGGTTGCGCCAGTAGGGCGTCTTAAACTGGTGAAAGCCACGGGCAGCGAAGTGCAGCGCGGTGACGACGGCATTTTTCGTTTAAGCGCAGAAACCCAGGCCACGCGTgg                                                                                                                                             <  1:1383574/150‑1 (MQ=255)
          ggCACAATCTCGGCGCTCAATCCGGGCGATCCGGCAAATACGGTTGCGCCAGTAGGGCGTCTTAAACTGGTGAAAGCCACGGGCAGCGAAGTGCAGCGCGGTGACGACGGCATTTTTCGTTTAAGCGCAGAAACCCAGGCCACGCGTggg                                                                                                                                            <  1:3588078/150‑1 (MQ=255)
            cacaATCTCGGCGCTCAATCCGGGCGATCCGGCAAATACGGTTGCGCCAGTAGGGCGTCTTAAACTGGTGAAAGCCACGGGCAGCGAAGTGCAGCGCGGTGACGACGGCATTTTTCGTTTAAGCGCAGAAACCCAGGCCACGCGTGGGcc                                                                                                                                          >  2:2298089/1‑150 (MQ=255)
              caATCTCGGCGCTCAATCCGGGCGATCCGGCAAATACGGTTGCGCCAGTAGGGCGTCTTAAACTGGTGAAAGCCACGGGCAGCGAAGTGCAGCGCGGTGACGACGGCATTTTTCGTTTAAGCGCAGAAACCCAGGCCACGCGTGGGCCgg                                                                                                                                        <  2:3590436/150‑1 (MQ=255)
               aaTCTCGGCGCTCAATCCGGGCGATCCGGCAAATACGGTTGCGCCAGTAGGGCGTCTTAAACTGGTGAAAGCCACGGGCAGCGAAGTGCAGCGCGGTGACGACGGCATTTTTCGTTTAAGCGCAGAAACCCAGGCCACGCGTGGGCCGGt                                                                                                                                       <  1:428685/150‑1 (MQ=255)
                   tcGGCGCTCAATCCGGGCGATCCGGCAAATACGGTTGCGCCAGTAGGGCGTCTTAAACTGGTGAAAGCCACGGGCAGCGAAGTGCAGCGCGGTGACGACGGCATTTTTCGTTTAAGCGCAGAAACCCAGACCACGCGTGGGCCGGTACTg                                                                                                                                   >  2:1673638/1‑150 (MQ=255)
                    cGGCGCTCAATCCGGGCGATCCGGCAAATACGGTTGCGCCAGTAGGGCGTCTTAAACTGGTGAAAGCCACGGGCAGCGAAGTGCAGCGCGGTGACGACGGCATTTTTCGTTTAAGCGCAGAAACCCAGGCCACGCGTGGGCCGGTACTgc                                                                                                                                  <  1:3595704/150‑1 (MQ=255)
                    cGGCGCTCAATCCGGGCGATACGGCAAATACGGTTGCGCCAGTAGGGCGTCTTATACTGGTGAAAGCCACGGGCAGCGAAGTGCAGCGCGGTGACGACGGCATTTTTCGTTTAAGCGCAGAAACCCAGGCCACGCGTGGGCCGGTACTgc                                                                                                                                  <  1:2148235/150‑1 (MQ=255)
                     ggCGCTCAATCCGGGCGATCCGGCAAATACGGTTGCGCCAGTAGGGCGTCTTAAACTGGTGAAAGCCACGGGCAGCGAAGTGCAGCGCGGTGACGACGGCATTTTTCGTTTAAGCGCAGAAACCCAGTCCACGCGTGGGCCGGTACTgca                                                                                                                                 <  2:1214731/150‑1 (MQ=255)
                     ggCGCTCAATCCGGGCGATCCGGCAAATACGGTTGCGCCAGTAGGGCGTCTTAAACTGGTGAAAGCCACGGGCAGCGAAGTGCAGCGCGGTGACGACGGCATTTTTCGTTTAAGCGCAGAAACCCAGTCCACGCGTGGGCCGGTACTgca                                                                                                                                 <  1:4207144/150‑1 (MQ=255)
                      gcgcTCAATCCGGGCGATCCGGCAAATACGGTTGCGCCAGTAGGGCGTCTTAAACTGGTGAAAGCCACGGGCAGCGAAGTGCAGCGCGGTGACGACGGCATTTTTCGTTTAAGCGCAGAAACCCAGGCCACGCGTGGGCCGGTACTgcag                                                                                                                                <  1:1254674/150‑1 (MQ=255)
                       cgcTCAATCCGGGCGATCCGGCAAATACGGGTGCGCCAGTAGGGGGTCTTAAACTGGTGAAAGCCACGGGCAGCGAAGTGCAGCGCGGTGACGACGGCATTTTTCGTTTAAGCGCAGAAACCCAGGCCACGCGTGGGCCGGTACTgcagg                                                                                                                               <  2:3143752/150‑1 (MQ=255)
                         cTCAATCCGGGCGATCCGGCAAATACGGTTGCGCCAGTAGGGCGTCTTAAACTGGTGAAAGCCACGGGCAGCGAAGTGCAGCGCGGTGACGACGGCATTTTTCGTTTAAGCGCAGAAACCCAGGCCACGCGTGGGCCGGTACTgcaggca                                                                                                                             <  1:2638859/150‑1 (MQ=255)
                           cAATCCGGGCGATCCGGCAAATACGGTTGCGCCAGTAGGGCGTCTTAAACTGGTGAAAGCCACGGGCAGCGAAGTGCAGCGCGGTGACGACGGCATTTTTCGTTTAAGCGCAGAAACCCAGGCCACGCGTGGGCCGGTACTGCAGGCAGa                                                                                                                           >  1:4097638/1‑150 (MQ=255)
                           cAATCCGGGCGATCCGGCAAATACGGTTGCGCCAGTAGGGCGTCTTAAACTGGTGAAAGCCACGGGCAGCGAAGTGCAGCGCGGTGACGACGGCATTTTTCGTTTAAGCGCAGAAACCCAGGCCACGCGTGGGCCGGTACTGCAGGCAGa                                                                                                                           <  2:3569355/150‑1 (MQ=255)
                           cAATCCGGGCGATCCGGCAAATACGGTTGCGCCAGTAGGGCGTCTTAAACTGGTGAAAGCCACGGGCAGCGAAGTGCAGCGCGGTGACGACGGCATTTTTCGTTTAAGCGCAGAAACCCAGGCCACGCGTGGGCCGGTACTGCAGGCAGa                                                                                                                           >  2:3068558/1‑150 (MQ=255)
                            aaTCCGGGCGATCCGGCAAATACGGTTGCGCCAGTAGGGCGTCTTAAACTGGTGAAAGCCACGGGCAGCGAAGTGCAGCGCGGTGACGACGGCATTTTTCGTTTAAGCGCAGAAACCCAGGCCACGCGTGGGCCGGTACTGCAGGCAGAt                                                                                                                          >  2:572868/1‑150 (MQ=255)
                                         cGGCAAATACGGTTGCGCCAGTAGGGCGTCTTAAACTGGTGAAAGCCACGGGCAGCGAAGTGCAGCGCGGTGACGACGGCATTTTTCGTTTAAGCGCAGAAACCCAGGCCACGCGTGGGCCGGTACTGCAGGCAGATCCAACCTTGCgtg                                                                                                             >  2:1725859/1‑150 (MQ=255)
                                         cGGCAAATACGGTTGCGCCAGTAGGGCGTCTTAAACTGGTGAAAGCCACGGGCAGCGAAGTGCAGCGCGGTGACGACGGCATTTTTCGTTTAAGCGCAGAAACCCAGGCCACGCGTGGGCCGGTACTGCAGGCAGATCCAACCTTGCgtg                                                                                                             >  2:3888749/1‑150 (MQ=255)
                                           gCAAATACGGTTGCGCCAGTAGGGCGTCTTAAACTGGTGAAAGCCACGGGCAGCGAAGTGCAGCGCGGTGACGACGGCATTTTTCGTTTAAGCGCAGAAACCCAGGCCACGCGTGGGCCGGTACTGCAGGCAGATCCAACCTTGCGtgtg                                                                                                           <  1:2407570/150‑1 (MQ=255)
                                             aaaTACGGTTGCGCCAGTAGGGCGTCTTAAACTGGTGAAAGCCACGGGCAGCGAAGTGCAGCGCGGTGACGACGGCATTTTTCGTTTAAGCGCAGAAACCCAGGCCACGCGTGGGCCGGTACTGCAGGCAGATCCAACCTTGCGTGTGAt                                                                                                         <  1:3243405/150‑1 (MQ=255)
                                              aaTACGGTTGCGCCAGTAGGGCGTCTTAAACTGGTGAAAGCCACGGGCAGCGAAGTGCAGCGCGGTGACGACGGCATTTTTCGTTTAAGCGCAGAAACCCAGGCCACGCGTGGGCCGGTACTGCAGGCAGATCCAACCTTGCGTGTGATg                                                                                                        >  1:1906888/1‑150 (MQ=255)
                                              aaTACGGTTGCGCCAGTAGGGCGTCTTAAACTGGTGAAAGCCACGGGCAGCGAAGTGCAGCGCGGTGACGACGGCATTTTTCGTTTAAGCGCAGAAACCCAGGCCACGCGTGGGCCGGTACTGCAGGCAGATCCAACCTTGCGTGTGATg                                                                                                        <  2:2162476/150‑1 (MQ=255)
                                              aaTACGGTTGCGCCAGTAGGGCGTCTTAAACTGGTGAAAGCCACGGGCAGCGAAGTGCAGCGCGGTGACGACGGCATTTTTCGTTTAAGCGCAGAAACCCAGGCCACGCGTGGGCCGGTACTGCAGGCAGATCCAACCTTGCGTGTGATg                                                                                                        <  1:3451662/150‑1 (MQ=255)
                                               aTACGGTTGCGCCAGTAGGGCGTCTTAAACTGGTGAAAGCCACGGGCAGCGAAGTGCAGCGCGGTGACGACGGCATTTTTCGTTTAAGCGCAGAAACCCAGGCCACGCGTGGGCCGGTACTGCAGGCAGATCCAACCTTGCGTGTGATGt                                                                                                       >  2:1087714/1‑150 (MQ=255)
                                                 aCGGTTGCGCCAGTAGGGCGTCTTAAACTGGTGAAAGCCACGGGCAGCGAAGTGCAGCGCGGTGACGACGGCATTTTTCGTTTAAGCGCAGAAACCCAGGCCACGCGTGGGCCGGTACTGCAGGCAGATCCAACCTTGCGTGTGATGTCg                                                                                                     <  2:245853/150‑1 (MQ=255)
                                                    gTTGCGCCAGTAGGGCGTCTTAAACTGGTGAAAGCCACGGGCAGCGAAGTGCAGCGCGGTGACGACGGCATTTTTCGTTTAAGCGCAGAAACCCAGTCCACGCGTGGGCCGGTACTGCAGGCAGATCCAACCTTGCGTGTGATGTCgggg                                                                                                  <  2:4528295/150‑1 (MQ=255)
                                                    gTTGCGCCAGTAGGGCGTCTTAAACTGGTGAAAGCCACGGGCAGCGAAGTGCAGCGCGGTGACGACGGCATTTTTCGTTTAAGCGCAGAAACCCAGGCCACGCGTGGGCCGGTACTGCAGGCAGATCCAACCTTGCGTGTGATGTCgggg                                                                                                  >  2:5053986/1‑150 (MQ=255)
                                                           cAGTAGGGCGTCTTAAACTGGTGAAAGCCACGGGCAGCGAAGTGCAGCGCGGTGACGACGGCATTTTTCGTTTAAGCGCAGAAACCCAGGCCACGCGTGGGCCGGTACTGCAGGCAGATCCAACCTTGCGTGTGATGTCGGGGGTTCTgg                                                                                           >  2:1553209/1‑150 (MQ=255)
                                                            aGTAGGGCGTCTTACCCTGGAGAAAGCCACGGGCAGCGAAGTGCAGCGCGGTGACGACGGCATTTTTCGTTTAAGCGCAGAAACCCAGGCCACGCGTGGGCCGGTACTGCAGGCAGATCCAACCTTGCGTGTGATGTCGGGGGTTCTGGa                                                                                          <  2:2675943/150‑1 (MQ=255)
                                                            aGTAGGGCGTCTTAAACTGGTGAAAGCCACGGGCAGCGAAGTGCAGCGCGGTGACGACGGCATTTTTCGTTTAAGCGCAGAAACCCAGGCCACGCGTGGGCCGGTACTGCAGGCAGATCCAACCTTGCGTGTGATGTCGGGGGTTCTGGa                                                                                          <  1:3784946/150‑1 (MQ=255)
                                                            aGTAGGGCGTCTTAAACTGGTGAAAGCCACGGGCAGCGAAGTGCAGCGCGGTGACGACGGCATTTTTCGTTTAAGCGCAGAAACCCAGGCCACGCGTGGGCCGGTACTGCAGGCAGATCCAACCTTGCGTGTGATGTCGGGGGTTCTGGa                                                                                          >  1:1568153/1‑150 (MQ=255)
                                                               aGGGCGTCTTAAACTGGTGAAAGCCACGGGCAGCGAAGTGCAGCGCGGTGACGACGGCATTTTTCGTTTAAGCGCAGAAACCCAGGCCACGCGTGGGCCGGTACTGCAGGCAGATCCAACCTTGCGTGTGATGTCGGGGGTTCTGGAAgg                                                                                       >  2:886111/1‑150 (MQ=255)
                                                                gggCGTCTTAAACTGGTGAAAGCCACGGGCAGCGAAGTGCAGCGCGGTGACGACGGCATTTTTCGTTTAAGCGCAGAAACCCAGGCCACGCGTGGGCCGGTACTGCAGGCAGATCCAACCTTGCGTGTGATGTCGGGGGTTCTGGAAGGc                                                                                      <  1:5015383/150‑1 (MQ=255)
                                                                gggCGTCTTAAACTGGTGAAAGCCACGGGCAGCGAAGTGCAGCGCGGTGACGACGGCATTTTTCGTTTAAGCGCAGAAACCAAGTCCACGCGTGGGCCGGTACTGCAGGCAGATCCAACCTTGCGTGTGATGTCGGGGGTTCTGGAAGGc                                                                                      <  1:5109318/150‑1 (MQ=255)
                                                                   cGTCTTAAACTGGTGAAAGCCACGGGCAGCGAAGTGCAGCGCGGTGACGACGGCATTTTTCGTTTAAGCGCAGAAACCCAGGCCACGCGTGGGCCGGTACTGCAGGCAGATCCAACCTTGCGTGTGATGTCGGGGGTTCTGGAAGGCAGt                                                                                   >  1:2111422/1‑150 (MQ=255)
                                                                    gTCTTAAACTGGTGAAAGCCACGGGCAGCGAAGTGCAGCGCGGTGACGACGGCATTTTTCGTTTAAGCGCAGAAACCCAGGCCACGCGTGGGCCGGTACTGCAGGCAGATCCAACCTTGCGTGTGATGTCGGGGGTTCTGGAAGGCAGTa                                                                                  >  2:725007/1‑150 (MQ=255)
                                                                      cTTAAACTGGTGAAAGCCACGGGCAGCGAAGTGCAGCGCGGTGACGACGGCATTTTTCGTTTAAGCGCAGAAACCCAGGCCACGCGTGGGCCGGTACTGCAGGCAGATCCAACCTTGCGTGTGATGTCGGGGGTTCTGGAAGGCAGTCAc                                                                                >  2:4057678/1‑150 (MQ=255)
                                                                      cTTAAACTGGTGAAAGCCACGGGCAGCGAAGTGCAGCGCGGTGACGACGGCATTTTTCGTTTAAGCGCAGAAACCCAGGCCACGCGTGGGCCGGTACTGCAGGCAGATCCAACCTTGCGTGTGATGTCGGGGGTTCTGGAAGGCAGTAAc                                                                                <  1:2778927/150‑1 (MQ=255)
                                                                        tAAACTGGTGAAAGCCACGGGCAGCGAAGTGCAGCGCGGTGACGACGGCATTTTTCGTTTAAGCGCAGAAACCCAGGCCACGCGTGGGCCGGTACTGCAGGCAGATCCAACCTTGCGTGTGATGTCGGGGGTTCTGGAAGGCAGTAACGt                                                                              <  1:223579/150‑1 (MQ=255)
                                                                         aaaCTGGTGAAAGCCACGGGCAGCGAAGTGCAGCGCGGTGACGACGGCATTTTTCGTTTAAGCGCAGAAACCCAGTCCACGCGTGGGCCGGTACTGCAGGCAGATCCAACCTTGCGTGTGATGTCGGGGGTTCTGGAAGGCAGTAACGTc                                                                             <  1:3027694/150‑1 (MQ=255)
                                                                         aaaCTGGTGAAAGCCACGGGCAGCGAAGTGCAGCGCGGTGACGACGGCATTTTTCGTTTAAGCGCAGAAACCCAGTCCACGCGTGGGCCGGTACTGCAGGCAGATCCAACCTTGCGTGTGATGTCGGGGGTTCTGGAAGGCAGTAACGTc                                                                             <  1:391317/150‑1 (MQ=255)
                                                                            cTGGTGAAAGCCACGGGCAGCGAAGTGCAGCGCGGTGACGACGGCATTTTTCGTTTAAGCGCAGAAACCCAGGCCACGCGTGGGCCGGTACTGCAGGCAGATCCAACCTTGCGTGTGATGTCGGGGGTTCTGGAAGGCAGTAACGTCAAt                                                                          <  1:54771/150‑1 (MQ=255)
                                                                            cTGGTGAAAGCCACGGGCAGCGAAGTGCAGCGCGGTGACGACGGCATTTTTCGTTTAAGCGCAGAAACCCAGGCCACGCGTGGGCCGGTACTGCAGGCAGATCCAACCTTGCGTGTGATGTCGGGGGTTCTGGAAGGCAGTAACGTCAAt                                                                          <  2:2809113/150‑1 (MQ=255)
                                                                              ggTGAAAGCCACGGGCAGCGAAGTGCAGCGCGGTGACGACGGCATTTTTCGTTTAAGCGCAGAAACCCAGGCCACGCGTGGGCCGGTACTGCAGGCAGATCCAACCTTGCGTGTGATGTCGGGGGTTCTGGAAGGCAGTAACGTCAATGc                                                                        <  2:2753551/150‑1 (MQ=255)
                                                                              ggTGAAAGCCACGGGCAGCGAAGTGCAGCGCGGTGACGACGGCATTTTTCGTTTAAGCGCAGAAACCCAGGCCACGCGTGGGCCGGTACTGCAGGCAGATCCAACCTTGCGTGTGATGTCGGGGGTTCTGGAAGGCAGTAACGTCAATGc                                                                        <  2:4815936/150‑1 (MQ=255)
                                                                              ggTGAAAGCCACGGGCAGCGAAGTGCAGCGCGGTGACGACGGCATTTTTCGTTTAAGCGCAGAAACCCAGGCCACGCGTGGGCCGGTACTGCAGGCAGATCCAACCGTGCGTGTGATGTCGGGGGTTCTGGAAGGCAGTAACGTCAATGc                                                                        >  1:3662213/1‑150 (MQ=255)
                                                                               gTGAAAGCCACGGGCAGCGAAGTGCAGCGCGGTGACGACGGCATTTTTCGTTTAAGCGCAGAAACCCAGGCCACGCGTGGGCCGGTACTGCAGGCAGATCCAACCTTGCGTGTGATGTCGGGGGTTCTGGAAGGCAGTAACGTCAATGcc                                                                       >  2:209260/1‑150 (MQ=255)
                                                                                 gAAAGCCACGGGCAGCGAAGTGCAGCGCGGTGACGACGGCATTTTTCGTTTAAGCGCAGAAACCCAGGCCACGCGTGGGCCGGTACTGCAGGCAGATCCAACCTTGCGTGTGATGTCGGGGGTTCTGGAAGGCAGTAACGTCAATGCCGt                                                                     >  2:5166057/1‑150 (MQ=255)
                                                                                   aaGCCACGGGCAGCGAAGTGCAGCGCGGTGACGACGGCATTTTTCGATTAAGCGCAGAAACCCAGGCCACGCGTGGGCCGGTACTGCAGGCAGATCCAACCTTGCGTGTGATGTCGGGGGTTCTGGAAGGCAGTAACGTCAATGCCGTTg                                                                   <  2:1554040/150‑1 (MQ=255)
                                                                                    aGCCACGGGCAGCGAAGTGCAGCGCGGTGACGACGGCATTTTTCGTTTAAGCGCAGAAACCCAGGCCACGCGTGGGCCGGTACTGCAGGCAGATCCAACCTTGCGTGTGATGTCGGGGGTTCTGGAAGGCAGTAACGTCAATGCCGTTGc                                                                  <  2:539110/150‑1 (MQ=255)
                                                                                     gCCACGGGCAGCGAAGTGCAGCGCGGTGACGACGGCATTTTTCGTTTAAGCGCAGAAACCCAGGCCACGCGTGGGCCGGTACTGCAGGCAGATCCAACCTTGCGTGTGATGTCGGGGGTTCTGGAAGGCAGTAACGTCAATGCCGTTGCg                                                                 >  2:4690977/1‑150 (MQ=255)
                                                                                      ccACGGGCAGCGAAGTGCAGCGCGGTGACGACGGCATTTTTCGTTTAAGCGCAGAAACCCAGGCCACGCGTGGGCCGGTACTGCAGGCAGATCCAACCTTGCGTGTGATGTCGGGGGTTCTGGAAGGCAGTAACGTCAATGCCGTTGCgg                                                                <  1:4420223/150‑1 (MQ=255)
                                                                                       cACGGGCAGCGAAGTGCAGCGCGGTGACGACGGCATTTTTCGTTTAAGCGCAGAAACCCAGGCCACGCGTGGGCCGGTACTGCAGGCAGATCCAACCTTGCGTGTGATGTCGGGGGTTCTGGAAGGCAGTAACGTCAATGCCGTTGCGGc                                                               >  2:3989522/1‑150 (MQ=255)
                                                                                       cACGGGCAGCGAAGTGCAGCGCGGTGACGACGGCATTTTTCGTTTAAGCGCAGAAACCCAGGCCACGCGTGGGCCGGTACTGCAGGCAGATCCAACCTTGCGTGTGATGTCCGGGGTTCTGGAAGGCAGTAACGTCAATGCCGTTGCGGc                                                               >  2:41261/1‑150 (MQ=255)
                                                                                           ggCAGCGAAGTGCAGCGCGGTGACGACGGCATTTTTCGTTTAAGCGCAGAAACCCAGGCCACGCGTGGGCCGGTACTGCAGGCAGATCCAACCTTGCGTGTGATGTCGGGGGTTCTGGAAGGCAGTAACGTCAATGCCGTTGCGGCAATg                                                           <  1:5238003/150‑1 (MQ=255)
                                                                                            gCAGCGAAGTGCAGCGCGGTGACGACGGCATTTTTCGTTTAAGCGCAGAAACCCAGGCCACGCGTGGGCCGGTACTGCAGGCAGATCCAACCTTGCGTGTGATGTCGGGGGTTCTGGAAGGCAGTAACGTCAATGCCGTTGCGGCAATGa                                                          <  1:4781056/150‑1 (MQ=255)
                                                                                             cAGCGAAGTGCAGCGCGGTGACGACGGCATTTTTCGTTTAAGCGCAGAAACCCAGGCCACGCGTGGGCCGGTACTGCAGGCAGATCCAACCTTGCGTGTGATGTCGGGGGTTCTGGAAGGCAGTAACGTCAATGCCGTTGCGGCAATGAg                                                         <  1:4679613/150‑1 (MQ=255)
                                                                                              aGCGAAGTGCAGCGCGGTGACGACGGCATTTTTCGTTTAAGCGCAGAAACCCAGGCCACGCGTGGGCCGGTACTGCAGGCAGATCCAACCTTGCGTGTGATGTCGGGGGTTCTGGAAGGCAGTAACGTCAATGCCGTTGCGGCAATGAGc                                                        <  1:1312738/150‑1 (MQ=255)
                                                                                               gCGAAGTGCAGCGCGGTGACGACGGCATTTTTCGTTTAAGCGCAGAAACCCAGGCCACGCGTGGGCCGGTACTGCAGGCAGATCCAACCTTGCGTGTGATGTCGGGGGTTCTGGAAGGCAGTAACGTCAATGCCGTTGCGGCAATGAGCg                                                       <  1:1412160/150‑1 (MQ=255)
                                                                                                cGAAGTGCAGCGCGGTGACGACGGCATTTTTCGTTTAAGCGCAGAAACCCAGGCCACGCGTGGGCCGGTACTGCAGGCAGATCCAACCTTGCGGGTGATGTCGGGGGTTCTGGAAGGCAGTAACGTCAATGCCGTTGCGGCAATGAGCGa                                                      >  1:3104626/1‑150 (MQ=255)
                                                                                                 gAAGTGCAGCGCGGTGACGACGGCATTTTTCGTTTAAGCGCAGAAACCCAGGCCACGCGTGGGCCGGTACTGCAGGCAGATCCAACCTTGCGTGTGATGTCGGGGGTTCTGGAAGGCAGTAACGTCAATGCCGTTGCGGCAATGAGCGAc                                                     >  1:4014223/1‑150 (MQ=255)
                                                                                                 gAAGTGCAGCGCGGTGACGACGGCATTTTTCGTTTAAGCGCAGAAACCCAGGCCACGCGTGGGCCGGTACTGCAGGCAGATCCAACCTTGCGTGTGATGTCGGGGGTTCTGGAAGGCAGTAACGTCAATGCCGTTGCGGCAATGAGCGAc                                                     >  1:843760/1‑150 (MQ=255)
                                                                                                 gAAGTGCAGCGCGGTGACGACGGCATTTTTCGTTTAAGCGCAGAAACCCAGGCCACGCGTGGGCCGGTACTGCAGGCAGATCCAACCTTGCGTGTGATGTCGGGGGTTCTGGAAGGCAGTAACGTCAATGCCGTTGCGGCAATGAGCGAc                                                     >  2:4449000/1‑150 (MQ=255)
                                                                                                 gAAGTGCAGCGCGGTGACGACGGCATTTTTCGTTTAAGCGCAGAAACCCAGGCCACGCGTGGGCCGGTACTGCAGGCAGATCCAACCTTGCGTGTGATGTCGGGGGTTCTGGAAGGCAGTAACGTCAATGCCGTTGCGGCAATGAGCGAc                                                     <  2:3556718/150‑1 (MQ=255)
                                                                                                  aaGTGCAGCGCGGTGACGACGGCATTTTTCGTTTAAGCGCAGAAACCCAGGCCACGCGTGGGCCGGTACTGCAGGCAGATCCAACCTTGCGTGTGATGTCGGGGGTTCTGGAAGGCAGTAACGTCAATGCCGTTGCGGCAATGAGCGACa                                                    >  1:3525404/1‑150 (MQ=255)
                                                                                                   aGTGCAGCGCGGTGACGACGGCATTTTTCGTTTAAGCGCAGAAACCCAGGCCACGCGTGGGCCGGTACTGCAGGCAGATCCAACCTTGCGTGTGATGTCGGGGGTTCTGGAAGGCAGTAACGTCAATGCCGTTGCGGCAATGAGCGACAt                                                   >  1:4530209/1‑150 (MQ=255)
                                                                                                      gCAGCGCGGTGACGTCGGCATTTTTCGTTTAAGCGCAGAAACCCAGGCCACGCGTGGGCCGGTACTGCAGGCAGATCCAACCTTGCGTGTGATGTCGGGGGTTCTGGAAGGCAGTAACGTCAATGCCGTTGCGGCAATGAGCGACATGAt                                                <  2:5184562/150‑1 (MQ=255)
                                                                                                      gCAGCGCGGTGACGACGGCATTTTTCGTTTAAGCGCAGAAACCCAGGCCACGCGTGGGCCGGTACTGCAGGCAGATCCAACCTTGCGTGTGATGTCGGGGGTTCTGGAAGGCAGTAACGTCAATGCCGTTGCGGCAATGAGCGACATGAt                                                <  1:117791/150‑1 (MQ=255)
                                                                                                      gCAGCGCGGTGACGACGGCATTTTTAGTTTAAGCGCAGAAACCCAGGCCACGCGTGGGCCGGTACTGCAGGCAGATCCAACCTTGCGTGTGATGTCGGGGGTTCTGGAAGGCAGTAACGTCAATGCCGTTGCGGCAATGAGCGACATGAt                                                <  2:2595541/150‑1 (MQ=255)
                                                                                                       cAGCGCGGTGACGACGGCATTTTTCGTTTAAGCGCAGAAACCCAGGCCACGCGTGGGCCGGTACTGCAGGCAGATCCAACCTTGCGTGTGATGTCGGGGGTTCTGGAAGGCAGTAACGTCAATGCCGTTGCGGCAATGAGCGACATGAtt                                               >  2:358087/1‑150 (MQ=255)
                                                                                                        aGCGCGGTGACGACGGCATTTTTCGTTTAAGCGCAGAAACCCAGGCCACGCGTGGGCCGGTACTGCAGGAAGATCCAACCTTGCGTGTGATGTCGGGGGTTCTGGAAGGCAGTAACGTCAATGCCGTTGCGGCAATGAGCGACATGATTg                                              >  2:204925/1‑150 (MQ=255)
                                                                                                          cgcgGTGACGACGGCATTTTTCGTTTAAGCGCAGAAACCCAGGCCACGCGTGGGCCGGTACTGCAGGCAGATCCAACCTTGCGTGTGATGTCGGGGGTTCTGGAAGGCAGTAACGTCAATGCCGTTGCGGCAATGAGCGACATGATTGcc                                            <  2:1595402/150‑1 (MQ=255)
                                                                                                           gcgGTGACGACGGCATTTTTCGTTTAAGCGCAGAAACCCAGGCCACGCGTGGGCCGGTACTGCAGGCAGATCCAACCTTGCGTGTGATGTCGGGGGTTCTGGAAGGCAGTAACGTCAATGCCGTTGCGGCAATGAGCGACATGATTGCCa                                           <  1:2054557/150‑1 (MQ=255)
                                                                                                                gacgacGGCATTTTTCGTTTAAGCGCAGAAACCCAGGCCACGCGTGGGCCGGTACTGCAGGCAGATCCAACCTTGCGTGTGATGTCGGGGGTTCTGGAAGGCAGTAACGTCAATGCCGTTGCGGCAATGAGCGACATGATTGCCAgcgcg                                      >  1:699210/1‑150 (MQ=255)
                                                                                                                gacgacGGCATTTTTCGTTTAAGCGCAGAAACCCAGGCCACGCGTGGGCCGGTACTGCAGGCAGATCCAACCTTGCGTGTGATGTCGGGGGTTCTGGAAGGCAGTAACGTCAATGCCGTTGCGGCAATGAGCGACATGATTGCCAgcgcg                                      <  1:2250870/150‑1 (MQ=255)
                                                                                                                  cgacgGCATTTTTCGTTTAAGCGCAGAAACCCAGGCCACGCGTGGGCCGGTACTGCAGGCAGATCCAACCTTGCGTGTGATGTCGGGGGTTCTGGAAGGCAGTAACGTCAATGCCGTTGCGGCACTGAGCGACATGATTGCCAGcgcgcg                                    >  2:2725374/1‑150 (MQ=255)
                                                                                                                  cgacgGCATTTTTCGTTTAAGCGCAGAAACCCAGGCCACGCGTGGGCCGGTACTGCAGGCAGATCCAACCTTGCGTGTGATGTCGGGGGTTCTGGAAGGCAGTAACGTCAATGCCGTTGCGGCAATGAGCGACATGATTGCCAGcgcgcg                                    >  1:1699090/1‑150 (MQ=255)
                                                                                                                     cgGCATTTTTCGTTTAAGCGCAGAAACCCAGGCCACGCGTGGGCCGGTACTGCAGGCAGATCCAACCTTGCGTGTGATGTCGGGGGTTCTGGAAGGCAGTAACGTCAATGCCGTTGCGGCAATGAGCGACATGATTGCCAGCGCgcggcg                                 >  2:4303672/1‑150 (MQ=255)
                                                                                                                     cgGCATTTTTCGTGTAAGCGCAGAAACCCAGGCCACGCGTGGGCCGGTACTGCAGGCAGATCCAACCTTGCGTGTGATGTCGGGGGTTCTGGAAGGCAGTAACGTCAATGCCGTTGCGGCAATGAGCGACATGATTGCCAGCGCgcggcg                                 <  2:1024262/150‑1 (MQ=255)
                                                                                                                         aTTTTTCGTTTAAGCGCAGAAACCCAGGCCACGCGTGGGCCGGTACTGCAGGCAGATCCAACCTTGCGTGTGATGTCGGGGGTTCTGGAAGGCAGTAACGTCAATGCCGTTGCGGCAATGAGCGACATGATTGCCAGCGCGCGGCGtttt                             >  1:3267047/1‑150 (MQ=255)
                                                                                                                          tttttCGTTTAAGCGCAGAAACCCAGGCCACGCGTGGGCCGGTACTGCAGGCAGATCCAACCTTGCGTGTGATGTCGGGGGTTCTGGAAGGCAGTAACGTCAATGCCGTTGCGGCAATGAGCGACATGATTGCCAGCGCGCGGCGTTTTg                            >  1:1811990/1‑150 (MQ=255)
                                                                                                                            tttCGTTTAAGCGCAGAAACCCAGGCCACGCGTGGGCCGGTACTGCAGGCAGATCCAACCTTGCGTGTGATGTCGGGGGTTCTGGAAGGCAGTAACGTCAATGCCGTTGCGGCAATGAGCGACATGATTGCCAGCGCGCGGCGTTTTGaa                          <  1:4079515/150‑1 (MQ=255)
                                                                                                                              tCGTTTAAGCGCAGAAACCCAGGCCACGCGTGGGCCGGTACTGCAGGCAGATCCAACCTTGCGTGTGATGTCGGGGGTTCTGGAAGGCAGTAACGTCAATGCCGTTGCGGCAATGAGCGACATGATTGCCAGCGCGCGGCGTTTTGAAAt                        <  1:1808831/150‑1 (MQ=255)
                                                                                                                               cGTTTAAGCGCAGAAACCCAGGCCACGCGTGGGCCGGTACTGCAGGCAGATCCAACCTTGCGTGTGATGTCGGGGGTTCTGGAAGGCAGTAACGTCAATGCCGTTGCGGCAATGAGCGACATGATTGCCAGCGCGCGGCGTTTTGAAATg                       >  2:957279/1‑150 (MQ=255)
                                                                                                                                gTTTAAGCGCAGAAACCCAGGCCACGCGTGGGCCGGTACTGCAGGCAGATCCAACCTTGCGTGTGATGTCGGGGGTTCTGGAAGGCAGTAACGTCAATGCCGTTGCGGCAATGAGCGACATGATTGCCAGCGCGCGGCGTTTTGAAATGc                      <  1:4389555/150‑1 (MQ=255)
                                                                                                                                     aGCGCAGAAACCCAGGCCACGCGTGGGCCGGTACTGCAGGCAGATCCAACCTTGCGTGTGATGTCGGGGGTTCTGGAAGGCAGTCACGTCAATGCCGTTGCGGCAATGAGCGACATGATTGCCAGCGCGCGGCGTTTTGAAATGCAGATg                 >  2:4251726/1‑150 (MQ=255)
                                                                                                                                               cccAGGCCACGCGTGGGCCGGTACTGCAGGCAGATCCAACCTTGCGTGTGATGTCGGGGGTTCTGGAAGGCAGTAACGTCAATGCCGTTGCGGCAATGAGCGACATGATTGCCAGCGCGCGGCGTTTTGAAATGCAGATGAAGGTGATca       <  1:1087714/150‑1 (MQ=255)
                                                                                                                                                ccAGGCCACGCGTGGGCCGGTACTGCAGGCAGATCCAACCTTGCGTGTGATGTCGGGGGTTCTGGAAGGCAGTAACGTCAATGCCGTTGCGGCAATGAGCGACATGATTGCCAGCGCGCGGCGTTTTGAAATGCAGATGAAGGTGATcag      >  2:1835779/1‑150 (MQ=255)
                                                                                                                                                  aGGCCACGCGTGGGCCGGTACTGCAGGCAGATCCAACCTTGCGTGTGATGTCGGGGGTTCTGGAAGGCAGTAACGTCAATGCCGTTGCGGCAATGAGCGACATGATTGCCAGCGCGCGGCGTTTTGAAATGCAGATGAAGGTGATcagca    <  1:1857667/150‑1 (MQ=255)
                                                                                                                                                  aGGCCACGCGTGGGCCGGTACTGCAGGCAGATCCAACCTTGCGTGTGATGTCGGGGGTTCTGGAAGGCAGTAACGTCAATGCCGTTGCGGCAATGAGCGACATGATTGCCAGCGCGCGGCGTTTTGAAATGCAGATGAAGGTGATcagca    >  2:1129867/1‑150 (MQ=255)
                                                                                                                                                  aGGCCACGCGTGGGCCGGTACTGCAGGCAGATCCAACCTTGCGTGTGATGTCGGGGGTTCTGGAAGGCAGTAACGTCAATGCCGTTGCGGCAATGAGCGACATGATTGCCAGCGCGCGGCGTTTTGAAATGCAGATGAAGGTGATcagca    >  2:3171008/1‑150 (MQ=255)
                                                                                                                                                   ggCCACGCGTGGGCCGGTACTGCAGGCAGATCCAACCTTGCGTGTGATGTCGGGGGTTCTGGAAGGCAGTAACGTCAATGCCGTTGCGGCAATGAGCGACATGATTGCCAGCGCGCGGCGTTTTGAAATGCAGATGAAGGTGATcagcag   >  1:292404/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
TGCTGCCGATGGCACAATCTCGGCGCTCAATCCGGGCGATCCGGCAAATACGGTTGCGCCAGTAGGGCGTCTTAAACTGGTGAAAGCCACGGGCAGCGAAGTGCAGCGCGGTGACGACGGCATTTTTCGTTTAAGCGCAGAAACCCAGGCCACGCGTGGGCCGGTACTGCAGGCAGATCCAACCTTGCGTGTGATGTCGGGGGTTCTGGAAGGCAGTAACGTCAATGCCGTTGCGGCAATGAGCGACATGATTGCCAGCGCGCGGCGTTTTGAAATGCAGATGAAGGTGATCAGCAGC  >  CP009273/1129668‑1129965

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: