Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,108,219:1 +A 100% coding (1159/1395 nt) wcaM ← colanic acid biosynthesis protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,108,2191.A96.8% 806.0 / 20.2 247coding (1159/1395 nt)wcaMcolanic acid biosynthesis protein
Reads supporting (aligned to +/- strand):  ref base . (3/5);  new base A (108/131);  total (111/136)
Fisher's exact test for biased strand distribution p-value = 7.34e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Polymorphic indel expands or contracts a homopolymer stretch.

Only 100 of 275 total aligned reads displayed.

TCCTTTTTCATTCACCGCATGAACATTTGCAAGAGACAGCAGTGTTTCTTTTTTCGCCATAAAAACACCACGAACGTCTTTGCGCATGTCAAAGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGCATAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATGTTAGTTAGCGCCACAAAGGAAACAGCATTCCCGGCGGAGATTTGAATGCCGCGCAATTTATAAGCAAGATGAGTATTATCCAGTTGAATATCA  >  CP009273/2108070‑2108367
                                                                                                                                                      |                                                                                                                                                     
tCCTTTTTCATTCACCGCATGAACATTTGCAAGAGACAGCAGTGTTTCTTTTTTCGCCATAAAAACACCACGAACGTCTTTGCGCATGTCAAAGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGCAt                                                                                                                                                        <  1:703920/150‑1 (MQ=255)
tCCTTTTTCATTCACCGCATGAACATTTGCAAGAGACAGCAGTGTTTCTTTTTTCGCCATAAAAACACCACGAACGTCTTTGCGCATGTCAAAGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGCAt                                                                                                                                                        <  1:3506861/150‑1 (MQ=255)
tCCTTTTTCATTCACCGCATGAACATTTGCAAGAGACAGCAGTGTTTCTTTTTTCGCCATAAAAACACCACGAACGTCTTTGCGCATGTCAAAGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGCAt                                                                                                                                                        <  1:2016025/150‑1 (MQ=255)
 ccTTTTTCATTCACCGCATGAACATTTGCAAGAGACAGCAGTGTTTCTTTTTTCGCCATAAAAACACCACGAACGTCTTTGCGCATGTCAAAGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGCATa                                                                                                                                                      <  2:3597793/150‑1 (MQ=255)
 ccTTTTTCATTCACCGCATGAACATTTGCAAGAGACAGCAGTGTTTCTTTTTTCGCCATAAAAACACCACGAACGTCTTTGCGCATGTCAAAGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGCATa                                                                                                                                                      <  2:2410368/150‑1 (MQ=255)
  cTTTTTCATTCACCGCATGAACATTTGCAAGAGACAGCAGTGTTTCTTTTTTCGCCATAAAAACACCACGAACGTCTTTGCGCATGTCAAAGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGCATaa                                                                                                                                                     <  2:4071041/150‑1 (MQ=255)
  cTTTTTCATTCACCGCATGAACATTTGCAAGAGACAGCAGTGTTTCTTTTTTCGCCATAAAAACACCACGAACGTCTTTGCGCATGTCAAAGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGCATaa                                                                                                                                                     <  1:265160/150‑1 (MQ=255)
    ttttCATTCACCGCATGAACATTTGCAAGAGACAGCAGTGTTTCTTTTTTCGCCATAAAAACACCACGAACGTCTTTGCGCATGTCAAAGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGCATaaaa                                                                                                                                                  <  2:4697522/150‑1 (MQ=255)
        cattcaCCGCATGAACATTTGCAAGAGACAGCAGTGTTTCTTTTTTCGCCATAAAAACACCACGAACGTCTTTGCGCATGTCAAAGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGCATaaaaaaaa                                                                                                                                              >  2:829365/1‑149 (MQ=255)
         attcaCCGCATGAACATTTGCAAGAGACAGCAGTGTTTCTTTTTTCGCCATAAAAACACCACGAACGTCTTTGCGCATGTCAATGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATAGTTCGCATaaaaaaaat                                                                                                                                             >  1:119269/1‑148 (MQ=255)
            caCCGCATGAACATTTGCAAGAGACAGCAGTGTTTCTTTTTTCGCCATAAAAACACCACGAACGTCTTTGCGCATGTCAAAGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGCATaaaaaaaatgtt                                                                                                                                          <  2:1235808/150‑6 (MQ=255)
            caCCGCATGAACATTTGCAAGAGACAGCAGTGTTTCTTTTTTCGCCATAAAAACAACACGAACGTCTTTGCGCATGTCAAAGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGCATaaaaaaaatgtt                                                                                                                                          >  1:3377044/1‑145 (MQ=255)
             aCCGCATGAACATTTGCAAGAGACAGCAGTGTTTCTTTTTTCGCCATAAAAACACCACGAACGTCTTTGCGCATGTCAAAGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGCATAAAAAAAATGTTg                                                                                                                                          <  2:1304043/150‑1 (MQ=255)
               cGCATGAACATTTGCAAGAGACAGCAGTGTTTCTTTTTTCGCCATAAAAACACCACGAACGTCTTTGCGCATGTCAAAGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGCATAAAAAAAATGTTGCg                                                                                                                                        <  2:2777605/150‑1 (MQ=255)
                gCATGAACATTTGCAAGAGACAGCAGTGTTTCTTTTTTCGCCATAAAAACACCACGAACGTCTTTGCGCATGTCAAAGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGCATAAAAAAAATGTTGCgg                                                                                                                                       <  2:4424830/150‑1 (MQ=255)
                  aTGAACATTTGCAAGAGACAGCAGTGTTTCTTTTTTCGCCATAAAAACACCACGAACGTCTTTGCGCATGTCAAAGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGCATAAAAAAAATGTTGCGGtt                                                                                                                                     >  1:1280901/1‑150 (MQ=255)
                   tGAACATTTGCAAGAGACAGCAGTGTTTCTTTTTTCGCCATAAAAACACCACGAACGTCTTTGCGCATGTCAAAGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGCATAAAAAAAATGTTGCGGttt                                                                                                                                    >  1:992098/1‑150 (MQ=255)
                    gAACATTTGCAAGAGACAGCAGTGTTTCTTTTTTCGCCATAAAAACACCACGAACGTCTTTGCGCATGTCAAAGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGCATAAAAAAAATGTTGCGGTttg                                                                                                                                   <  1:5051385/150‑1 (MQ=255)
                    gAACATTTGCAAGAGACAGCAGTGTTTCTTTTTTCGCCATAAAAACACCACGAACGTCTTTGCGCATGTCAAAGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGCATAAAAAAAATGTTGCGGTttg                                                                                                                                   <  2:4550314/150‑1 (MQ=255)
                     aaCATTTGCAAGAGACAGCAGTGTTTCTTTTTTCGCCATAAAAACACCACGAACGTCTTTGCGCATGTCAAAGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGCATAAAAAAAATGTTGCGGTttgt                                                                                                                                  <  2:145166/150‑1 (MQ=255)
                     aaCATTTGCAAGAGACAGCAGTGTTTCTTTTTTCGCCATAAAAACACCACGAACGTCTTTGCGCATGTCAAAGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGCATAAAAAAAATGTTGCGGTttgt                                                                                                                                  >  1:817786/1‑150 (MQ=255)
                        aTTTGCAAGAGACAGCAGTGTTTCTTTTTTCGCCATAAAAACACCACGAACGTCTTTGCGCATGTCAAAGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGCATAAAAAAAATGTTGCGGTTtgttgt                                                                                                                               >  2:3608897/1‑150 (MQ=255)
                         tttGCAAGAGACAGCAGTGTTTCTTTTTTCGCCATAAAAACACCACGAACGTCTTTGCGCATGTCAAAGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGCATAAAAAAAATGTTGCGGTTTGTtgtg                                                                                                                              <  1:4446030/150‑1 (MQ=255)
                            gCAAGAGACAGCAGTGTTTCTTTTTTCGCCATAAAAACACCACGAACGTCTTTGCGCATGTCAAAGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGCATAAAAAAAATGTTGCGGTTTGTTGTGTaa                                                                                                                           <  2:3908991/150‑1 (MQ=255)
                              aaGAGACAGCAGTGTTTCTTTTTTCGCCATAAAAACACCACGAACGTCTTTGCGCATGTCAAAGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGCATAAAAAAAATGTTGCGGTTTGTTGTGTAACt                                                                                                                         >  2:66293/1‑150 (MQ=255)
                              aaGAGACAGCAGTGTTTCTTTTTTCGCCATAAAAACACCACGAACGTCTTTGCGCATGTCAAAGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGCATAAAAAAAATGTTGCGGTTTGTTGTGTAACt                                                                                                                         >  2:1801936/1‑150 (MQ=255)
                                  gaCAGCAGTGTTTCTTTTTTCGCCATAAAAACACCACGAACGTCTTTGCGCATGTCAAAGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGCATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAg                                                                                                                     <  1:3637609/150‑1 (MQ=255)
                                   aCAGCAGTGTTTCTTTTTTCGCCATAAAAACACCACGAACGTCTTTGCGCATGTCAAAGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGCATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGc                                                                                                                    <  1:2869486/150‑1 (MQ=255)
                                    cagcagTGTTTCTTTTTTCGCCATAAAAACACCACGAACGTCTTTGCGCATGTCAAAGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGCATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCg                                                                                                                   <  2:3279448/150‑1 (MQ=255)
                                         gtgtTTCTTTTTTCGCCATAAAAACACCACGAACGTCTTTGCGCATGTCAAAGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGCATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGacgca                                                                                                              <  1:487821/150‑1 (MQ=255)
                                           gtTTCTTTTTTCGCCATAAAAACACCACGAACGTCTTTGCGCATGTCAAAGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGCATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGacgcacg                                                                                                            <  1:2219573/150‑1 (MQ=255)
                                            tttCTTTTTTCGCCATAAAAACACCACGAACGTCTTTGCGCATGTCAAAGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGCATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGacgcacgc                                                                                                           <  2:306168/150‑1 (MQ=255)
                                              tCTTTTTTCGCCATAAAAACACCACGAACGTCTTTGCGCATGTCAAAGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACCTTGATATTTCGCATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCtt                                                                                                         <  2:1123907/150‑1 (MQ=255)
                                              tCTTTTTTCGCCATAAAAACACCACGAACGTCTTTGCGCATGTCAAAGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGCATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCtt                                                                                                         >  1:1880695/1‑150 (MQ=255)
                                                 tttttCGCCATAAAAACACCACGAACGTCTTTGCGCATGTCAAAGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGCATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCAt                                                                                                      <  1:107510/150‑1 (MQ=255)
                                                    ttCGCCATAAAAACACCACGAACGTCTTTGCGCATGTCAAAGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGCATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCAtctc                                                                                                   >  1:3556004/1‑150 (MQ=255)
                                                         cATAAAAACACCACGAACGTCTTTGCGCATGTCAAAGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGCATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATgt                                                                                              >  1:4845725/1‑150 (MQ=255)
                                                           tAAAAACACCACGAACGTCTTTGCGCATGTCAAAGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGCATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATgtta                                                                                            <  2:2229507/150‑1 (MQ=255)
                                                             aaaaCACCACGAACGTCTTTGCGCATGTCAAAGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGCATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATgttagt                                                                                          >  1:1971861/1‑150 (MQ=255)
                                                             aaaaCACCACGAACGTCTTTGCGCATGTCAAAGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGCATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATgttagt                                                                                          <  1:8890/150‑1 (MQ=255)
                                                              aaaCACCACGAACGTCTTTGCGCATGTCAAAGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGCATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATgttagtt                                                                                         >  1:3258495/1‑150 (MQ=255)
                                                              aaaCACCACGAACGTCTTTGCGCATGTCAAAGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGCATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATgttagtt                                                                                         >  2:417587/1‑150 (MQ=255)
                                                              aaaCACCACGAACGTCTTTGCGCATGTCAAAGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGCATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATgttagtt                                                                                         <  2:1903834/150‑1 (MQ=255)
                                                              aaaCACCACGAACGTCTTTGCGCATGTCAAAGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGCATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATgttagtt                                                                                         <  1:4949101/150‑1 (MQ=255)
                                                                 caccacGAACGTCTTTGCGCATGTCAAAGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGCATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATGTTAGTTAgc                                                                                      >  1:2901752/1‑150 (MQ=255)
                                                                  accacGAACGTCTTTGCGCATGTCAAAGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGCATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATGTTAGTTAgcg                                                                                     <  1:4427526/150‑1 (MQ=255)
                                                                      cGAACGTCTTTGCGCATGTCAAAGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGCATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATGTTAGTTAGCGCcac                                                                                 >  1:4225962/1‑150 (MQ=255)
                                                                      cGAACGTCTTTGCGCATGTCAAAGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGCATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATGTTAGTTAGCGCcac                                                                                 >  2:1155919/1‑150 (MQ=255)
                                                                      cGAACGTCTTTGCGCATGTCAAAGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGCATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATGTTAGTTAGCGCcac                                                                                 >  2:4658091/1‑150 (MQ=255)
                                                                       gAACGTCTTTGCGCATGTCAAAGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGCATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATGTTAGTTAGCGCcaca                                                                                >  1:2551186/1‑150 (MQ=255)
                                                                        aaCGTCTTTGCGCATGTCAAAGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGCATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATGTTAGTTAGCGCCACaa                                                                               <  2:1716552/150‑1 (MQ=255)
                                                                          cGTCTTTGCGCATGTCAAAGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGCATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATGTTAGTTAGCGCCACAAAg                                                                             >  1:3527723/1‑150 (MQ=255)
                                                                           gTCTTTGCGCATGTCAAAGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGCATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATGTTAGTTAGCGCCACAAAgg                                                                            >  2:1077141/1‑150 (MQ=255)
                                                                            tCTTTGCGCATGTCAAAGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGCATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATGTTAGTTAGCGCCACAAAGGa                                                                           <  2:4228788/150‑1 (MQ=255)
                                                                               ttGCGCATGTCAAAGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGCATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATGTTAGTTAGCGCCACAAAGGAAAc                                                                        >  2:4303824/1‑150 (MQ=255)
                                                                                tGCGCATGTCAAAGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGCATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATGTTAGTTAGCGCCACAAAGGAAACa                                                                       <  2:3879968/150‑1 (MQ=255)
                                                                                  cgcATGTCAAAGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGCATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATGTTAGTTAGCGCCACAAAGGAAACAGc                                                                     >  2:3798004/1‑150 (MQ=255)
                                                                                  cgcATGTCAAAGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGCATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATGTTAGTTAGCGCCACAAAGGAAACAGc                                                                     >  1:587144/1‑150 (MQ=255)
                                                                                     aTGTCAAAGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGCATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATGTTAGTTAGCGCCACAAAGGAAACAGCAtt                                                                  <  2:3719841/150‑1 (MQ=255)
                                                                                     aTGTCAAAGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGCATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATGTTAGTTAGCGCCACAAAGGAAACAGCAtt                                                                  >  2:324682/1‑150 (MQ=255)
                                                                                        tCAAAGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGCATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATGTTAGTTAGCGCCACAAAGGAAACAGCATTccc                                                               <  1:3842853/150‑1 (MQ=255)
                                                                                         cAAAGTTCATGCTCAATGCGGGTCCAACTGAGTAGTCCTGCAACCCATTGATATTTCGCATAAACAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATGTTAGTTAGCGCCACAAAGGAAACAGCATTCCcg                                                              <  2:2955916/150‑1 (MQ=255)
                                                                                         cAAAGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGCATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATGTTAGTTAGCGCCACAAAGGAAACAGCATTCCcg                                                              >  1:2293959/1‑150 (MQ=255)
                                                                                          aaaGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGCATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATGTTAGTTAGCGCCACAAAGGAAACAGCATTCCcgg                                                             <  1:2176397/150‑1 (MQ=255)
                                                                                           aaGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGCATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATGTTAGTTAGCGCCACAAAGGAAACAGCATTCCcggc                                                            <  1:1426457/150‑1 (MQ=255)
                                                                                            aGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGCATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATGTTAGTTAGCGCCACAAAGGAAACAGCATTCCcggcg                                                           >  2:167647/1‑150 (MQ=255)
                                                                                            aGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGCATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATGTTAGTTAGCGCCACAAAGGAAACAGCATTCCcggcg                                                           <  1:2159005/150‑1 (MQ=255)
                                                                                              ttCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGCATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATGTTAGTTAGCGCCACAAAGGAAACAGCATTCCCGGCGga                                                         >  2:2092870/1‑150 (MQ=255)
                                                                                                    cTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGCATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATGTTAGTTAGCGCCACAAAGGAAACAGCATTCCCGGCGGAGATTTg                                                   <  2:234807/150‑1 (MQ=255)
                                                                                                     tCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGCATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATGTTAGTTAGCGCCACAAAGGAAACAGCATTCCCGGCGGAGATTTGa                                                  >  2:442042/1‑150 (MQ=255)
                                                                                                      cAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGCATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATGTTAGTTAGCGCCACAAAGGAAACAGCATTCCCGGCGGAGATTTGaa                                                 <  1:2041726/150‑1 (MQ=255)
                                                                                                          gCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGCATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATGTTAGTTAGCGCCACAAAGGAAACAGCATTCCCGGCGGAGATTTGAATGcc                                             <  1:3682303/150‑1 (MQ=255)
                                                                                                            gggTCCAACTGAGGATTCCTGCATCACATTGATATTTCGCATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATGTTAGTTAGCGCCACAAAGGAAACAGCATTCCCGGCGGAGATTTGAATGCcgc                                           >  2:3483916/1‑150 (MQ=255)
                                                                                                             ggTCCAACTGAGGATTCCTGCATCACATTGATATTTCGCATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATGTTAGTTAGCGCCACAAAGGAAACAGCATTCCCGGCGGAGATTTGAATGCcgcg                                          <  1:3919243/150‑1 (MQ=255)
                                                                                                             ggTCCAACTGAGGATTCCTGCATCACATTGATATTTCGCATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATGTTAGTTAGCGCCACAAAGGAAACAGCATTCCCGGCGGAGATTTGAATGCcgcg                                          <  2:1584873/150‑1 (MQ=255)
                                                                                                                 cAACTGAGGATTCCTGCATCACATTGATATTTCGCATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATGTTAGTTAGCGCCACAAAGGAAACAGCATTCCCGGCGGAGATTTGAATGCCGCGCAAt                                      >  2:3517580/1‑150 (MQ=255)
                                                                                                                  aaCTGAGGATTCCTGCATCACATTGATATTTCGCATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATGTTAGTTAGCGCCACAAAGGAAACAGCATTCCCGGCGGAGATTTGAATGCCGCGCAAtt                                     >  1:1399621/1‑150 (MQ=255)
                                                                                                                    cTGAGGATTCCTGCATCACATTGATATTTCGCATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGGCGCACGCTTCATCTCAATGTTAGTTAGCGCCACAAAGGAAACAGCATTCCCGGCGGAGATTTGAATGCCGCGCAATTta                                   <  2:3711504/150‑1 (MQ=255)
                                                                                                                    cTGAGGATTCCTGCATCACATTGATATTTCGCATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATGTTAGTTAGCGCCACAAAGGAAACAGCATTCCCGGCGGAGATTTGAATGCCGCGCAATTta                                   >  2:5002494/1‑150 (MQ=255)
                                                                                                                     tGAGGATTCCTGCATCACATTGATATTTCGCATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATGTTAGTTAGCGCCACAAAGGAAACAGCATTCCCGGCGGAGATTTGAATGCCGCGCAATTtat                                  >  1:3288355/1‑150 (MQ=255)
                                                                                                                      gAGGATTCCTGCATCACATTGATATTTCGCATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATGTTAGTTAGCGCCACAAAGGAAACAGCATTCCCGGCGGAGATTTGAATGCCGCGCAATTtata                                 >  2:45348/1‑150 (MQ=255)
                                                                                                                      gAGGATTCCTGCATCACATTGATATTTCGCATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATGTTAGTTAGCGCCACAAAGGAAACAGCATTCCCGGCGGAGATTTGAATGCCGCGCAATTtata                                 <  2:624526/150‑1 (MQ=255)
                                                                                                                       aGGATTCCTGCATCACATTGATATTTCGCATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATGTTAGTTAGCGCCACAAAGGAAACAGCATTCCCGGCGGAGATTTGAATGCCGCGCAATTTATaa                                <  2:3886785/150‑1 (MQ=255)
                                                                                                                        ggATTCCTGCATCACATTGATATTTCGCATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATGTTAGTTAGCGCCACAAAGGAAACAGCATTCCCGGCGGAGATTTGAATGCCGCGCAATTTATAAg                               <  1:5124185/150‑1 (MQ=255)
                                                                                                                         gATTCCTGCATCACATTGATATTTCGCATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATGTTAGTTAGCGCCACAAAGGAAACAGCATTCCCGGCGGAGATTTGAATGCCGCGCAATTTATAAGc                              <  1:973045/150‑1 (MQ=255)
                                                                                                                             ccTGCATCACATTGATATTTCGCATAAAAAAAATGTTGCCGGTTGGTGTGTAACGACAGCCACCCACGCATCAAAAAAATGGTAGGTAGCGCCACCAAAGAAAAACCAATACCCGCCGAATATTGAAATCCCCGCAAATTAAAATCAAGa                          >  2:212424/1‑150 (MQ=255)
                                                                                                                              cTGCATCACATTGATATTTCGCATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATGTTAGTTAGCGCCACAAAGGAAACAGCATTCCCGGCGGAGATTTGAATGCCGCGCAATTTATAAGCAAGAt                         <  1:148817/150‑1 (MQ=255)
                                                                                                                                gCATCACATTGATATTTCGCATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATGTTAGTTAGCGCCACAAAGGAAACAGCATTCCCGGCGGAGATTTGAATGCCGCGCAATTTATAAGCAAGATGa                       <  1:3353563/150‑1 (MQ=255)
                                                                                                                                  aTCACATTGATATTTCGCATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATGTTAGTTAGCGCCACAAAGGAAACAGCATTCCCGGCGGAGATTTGAATGCCGCGCAATTTATAAGCAAGATGAGt                     <  2:1411801/150‑1 (MQ=255)
                                                                                                                                  aTCACATTGATATTTCGCATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATGTTAGTTAGCGCCACAAAGGAAACAGCATTCCCGGCGGAGATTTGAATGCCGCGCAATTTATAAGCAAGATGAGt                     <  2:2912076/150‑1 (MQ=255)
                                                                                                                                   tCACATTGATATTTCGCATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATGTTAGTTAGCGCCACAAAGGAAACAGCATTCCCGGCGGAGATTTGAATGCCGCGCAATTTATAAGCAAGATGAGta                    <  1:1708859/150‑1 (MQ=255)
                                                                                                                                      caTTGATATTTCGCATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATGTTAGTTAGCGCCACAAAGGAAACAGCATTCCCGGCGGAGATTTGAATGCCGCGCAATTTATAAGCAAGATGAGtatta                 >  2:4865829/1‑150 (MQ=255)
                                                                                                                                         tGATATTTCGCATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATGTTAGTTAGCGCCACAAAGGAAACAGCATTCCCGGCGGAGATTTGAATGCCGCGCAATTTATAAGCAAGATGAGTATTATcc              <  2:1076928/150‑1 (MQ=255)
                                                                                                                                         tGATATTTCGCATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATGTTAGTTAGCGCCACAAAGGAAACAGCATTCCCGGCGGAGATTTGAATGCCGCGCAATTTATAAGCAAGATGAGTATTATcc              >  1:2404596/1‑150 (MQ=255)
                                                                                                                                          gATATTTCGCATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATGTTAGTTAGCGCCACAAAGGAAACAGCATTCCCGGCGGAGATTTGAATGCCGCGCAATTTATAAGCAAGATGAGTATTATCCa             <  1:66293/150‑1 (MQ=255)
                                                                                                                                            tatTTCGCATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATGTTAGTTAGCGCCACAAAGGAAACAGCATTCCCGGCGGAGATTTGAATGCCTCGCAATTTATAAGCAAGATGAGTATTATCCAGt           >  2:24739/1‑150 (MQ=255)
                                                                                                                                            tatTTCGCATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATGTTAGTTAGCGCCACAAAGGAAACAGCATTCCCGGCGGAGATTTGAATGCCGCGCAATTTATAAGCAAGATGAGTATTATCCAGt           >  2:92164/1‑150 (MQ=255)
                                                                                                                                             atTTCGCATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATGTTAGTTAGCGCCACAAAGGAAACAGCATTCCCGGCGGAGATTTGAATGCCGCGCAATTTATAAGCAAGATGAGTATTATCCAGtt          <  2:4036961/150‑1 (MQ=255)
                                                                                                                                             atTTCGCATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATGTTAGTTAGCGCCACAAAGGAAACAGCATTCCCGGCGGAGATTTGAATGCCGCGCAATTTATAAGCAAGATGAGTATTATCCAGtt          >  2:3816558/1‑150 (MQ=255)
                                                                                                                                                 cGCATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATGTTAGTTAGCGCCACAAAGGAAACAGCATTCCCGGCGGAGATTTGAATGCCGCGCAATTTATAAGCAAGATGAGTATTATCCAGTTGAat      <  2:2035748/150‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                                     
TCCTTTTTCATTCACCGCATGAACATTTGCAAGAGACAGCAGTGTTTCTTTTTTCGCCATAAAAACACCACGAACGTCTTTGCGCATGTCAAAGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGCATAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATGTTAGTTAGCGCCACAAAGGAAACAGCATTCCCGGCGGAGATTTGAATGCCGCGCAATTTATAAGCAAGATGAGTATTATCCAGTTGAATATCA  >  CP009273/2108070‑2108367

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: