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.A98.6% 1158.9 / 8.6 350coding (1159/1395 nt)wcaMcolanic acid biosynthesis protein
Reads supporting (aligned to +/- strand):  ref base . (3/2);  new base A (175/170);  total (178/172)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Polymorphic indel expands or contracts a homopolymer stretch.

Only 100 of 381 total aligned reads displayed.

CCTTTTTCATTCACCGCATGAACATTTGCAAGAGACAGCAGTGTTTCTTTTTTCGCCATAAAAACACCACGAACGTCTTTGCGCATGTCAAAGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGCAT‑‑AAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATGTTAGTTAGCGCCACAAAGGAAACAGCATTCCCGGCGGAGATTTGAATGCCGCGCAATTTATAAGCAAGATGAGTATTATCCAGTTGAAT  >  CP009273/2108071‑2108363
                                                                                                                                                     |                                                                                                                                                 
 cTTTTTCATTCACCGCATGAACATTTGCAAGAGACAGCAGTGTTTCTTTTTTCGCCATAAAAACACCACGAACGTCTTTGCGCATGTCAAAGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGCAT‑‑aa                                                                                                                                                <  2:3875882/150‑1 (MQ=255)
   ttttCATTCACCGCATGAACATTTGCAAGAGACAGCAGTGTTTCTTTTTTCGCCATAAAAACACCACGAACGTCTTTGCGCATGTCAAAGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGCAT‑‑aaaa                                                                                                                                              <  1:3934215/150‑1 (MQ=255)
   ttttCATTCACCGCATGAACATTTGCAAGAGACAGCAGTGTTTCTTTTTTCGCCATAAAAACACCACGAACGTCTTTGCGCATGTCAAAGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGCAT‑‑aaaa                                                                                                                                              <  1:229368/150‑1 (MQ=255)
        attcaCCGCATGAACATTTGCAAGAGACAGCAGTGTTTCTTTTTTCGCCATAAAAACACCACGAACGTCTTTGCGCATGTCAAAGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGCAT‑‑aaaaaaaat                                                                                                                                         >  1:3368352/1‑148 (MQ=255)
        attcaCCGCATGAACATTTGCAAGAGACAGCAGTGTTTCTTTTTTCGCCATAAAAACACCACGAACGTCTTTGCGCATGTCAAAGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGCAT‑‑aaaaaaaat                                                                                                                                         <  2:1828634/150‑3 (MQ=255)
        attcaCCGCATGAACATTTGCAAGAGACAGCAGTGTTTCTTTTTTCGCCATAAAAACACCACGAACGTCTTTGCGCATGTCAAAGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGCAT‑‑aaaaaaaat                                                                                                                                         <  2:461123/150‑3 (MQ=255)
        attcaCCGCATGAACATTTGCAAGAGACAGCAGTGTTTCTTTTTTCGCCATAAAAACACCACGAACGTCTTTGCGCATGTCAAAGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGCAT‑‑aaaaaaaat                                                                                                                                         <  2:4280203/150‑3 (MQ=255)
          tcaCCGCATGAACATTTGCAAGAGACAGCAGTGTTTCTTTTTTCGCCATAAAAACACCACGAACGTCTTTGCGCATGTCAAAGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGCAT‑‑aaaaaaaatgt                                                                                                                                       <  2:1057064/150‑5 (MQ=255)
          tcaCCGCATGAACATTTGCAAGAGACAGCAGTGGTTCTTTTTTCGCCATAAAAACACCACGATCGTCTTTGCGCATGTCAAAGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGCAT‑‑aaaaaaaatgt                                                                                                                                       <  2:2885505/150‑5 (MQ=255)
            aCCGCATGAACATTTGCAAGAGACAGCAGTGTTTCTTTTTTCGCCATAAAAACACCACGAACGTCTTTGCGCATGTCAAAGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGCATAAAAAAAATGTTg                                                                                                                                      <  1:4794401/150‑1 (MQ=255)
                 aTGAACATTTGCAAGAGACAGCAGTGTTTCTTTTTTCGCCATAAAAACACCACGAACGTCTTTGCGCATGTCAAAGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGCATAAAAAAAATGTTGCGGtt                                                                                                                                 <  1:2708166/150‑1 (MQ=255)
                  tGAACATTTGCAAGAGACAGCAGTGTTTCTTTTTTCGCCATAAAAACACCACGAACGTCTTTGCGCATGTCAAAGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGCATAAAAAAAATGTTGCGGttt                                                                                                                                >  2:569653/1‑150 (MQ=255)
                   gAACATTTGCAAGAGACAGCAGTGTTTCTTTTTTCGCCATAAAAACACCACGAACGTCTTTGCGCATGTCAAAGTTCATGCTCAATGCGTGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGCATAAAAAAAATGTTGCGGTttg                                                                                                                               >  1:1258450/1‑150 (MQ=255)
                     aCATTTGCAAGAGACAGCAGTGTTTCTTTTTTCGCCATAAAAACACCACGAACGTCTTTGCGCATGTCAAAGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGCATAAAAAAAATGTTGCGGTttgtt                                                                                                                             <  1:5189769/150‑1 (MQ=255)
                      cATTTGCAAGAGACAGCAGTGTTTCTTTTTTCGCCATAAAAACACCACGAACGTCTTTGCGCATGTCAAAGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGCATAAAAAAAATGTTGCGGTttgttg                                                                                                                            <  1:1546255/150‑1 (MQ=255)
                       aTTTGCAAGAGACAGCAGTGTTTCTTTTTTCGCCATAAAAACACCACGAACGTCTTTGCGCATGTCAAAGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGCATAAAAAAAATGTTGCGGTTtgttgt                                                                                                                           <  1:3531938/150‑1 (MQ=255)
                         ttGCAAGAGACAGCAGTGTTTCTTTTTTCGCCATAAAAACACCACGAACGTCTTTGCGCATGTCAAAGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGCATAAAAAAAAATGTTGCGGTTTGTtgtg                                                                                                                          >  1:331682/1‑150 (MQ=255)
                          tGCAAGAGACAGCAGTGTTTCTTTTTTCGCCATAAAAACACCACGAACGTCTTTGCGCATGTCAAAGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGCATAAAAAAAATGTTGCGGTTTGTTGTGTa                                                                                                                        >  1:2022609/1‑150 (MQ=255)
                             aaGAGACAGCAGTGTTTCTTTTTTCGCCATAAAAACACCACGAACGTCTTTGCGCATGTCAAAGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGCATAAAAAAAATGTTGCGGTTTGTTGTGTAACt                                                                                                                     >  1:621603/1‑150 (MQ=255)
                               gagaCAGCAGTGTTTCTTTTTTCGCCATAAAAACACCACGAACGTCTTTGCGCATGTCAAAGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGCATAAAAAAAATGTTGCGGTTTGTTGTGTAACTcc                                                                                                                   >  1:376599/1‑150 (MQ=255)
                               gagaCAGCAGTGTTTCTTTTTTCGCCATAAAAACACCACGAACGTCTTTGCGCATGTCAAAGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGCATAAAAAAAATGTTGCGGTTTGTTGTGTAACTcc                                                                                                                   <  1:558782/150‑1 (MQ=255)
                                 gaCAGCAGTGTTTCTTTTTTCGCCATAAAAACACCACGAACGTCTTTGCGCATGTCAAAGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGCATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAg                                                                                                                 <  1:1368891/150‑1 (MQ=255)
                                  aCAGCAGTGTTTCTTTTTTCGCCATAAAAACACCACGAACGTCTTTGCGCATGTCAAAGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATACTTCGCATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGc                                                                                                                >  1:1780445/1‑150 (MQ=255)
                                     gcagTGTTTCTTTTTTCGCCATAAAAACACCACGAACGTCTTTGCGCATGTCAAAGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGCATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGac                                                                                                             >  1:3116460/1‑150 (MQ=255)
                                      cagTGTTTCTTTTTTCGCCATAAAAACACCACGAACGTCTTTGCGCATGTCAAAGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGCATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGacg                                                                                                            <  1:638732/150‑1 (MQ=255)
                                       agTGTTTCTTTTTTCGCCATAAAAACACCACGAACGTCTTTGCGCATGTCAAAGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGCATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGacgc                                                                                                           >  2:3796167/1‑150 (MQ=255)
                                             tCTTTTTTCGCCATAAAAACACCACGAACGTCTTTGCGCATGTCAAAGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGCATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCtt                                                                                                     >  1:965341/1‑150 (MQ=255)
                                             tCTTTTTTCGCCATAAAAACACCACGAACGTCTTTGCGCATGTCAAAGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGCATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCtt                                                                                                     <  1:658247/150‑1 (MQ=255)
                                              cTTTTTTCGCCATAAAAACACCACGAACGTCTTTGCGCATGTCAAAGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGCATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTc                                                                                                    >  2:2135052/1‑150 (MQ=255)
                                                  tttCGCCATAAAAACACCACGAACGTCTTTGCGCATGTCAAAGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGCATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCAtct                                                                                                >  1:4682763/1‑150 (MQ=255)
                                                  tttCGCCATAAAAACACCACGAACGTCTTTGCGCATGTCAAAGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGCATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCAtct                                                                                                >  2:1688266/1‑150 (MQ=255)
                                                  tttCGCCATAAAAACACCACGAACGTCTTTGCGCATGTCAAAGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGCATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCAtct                                                                                                <  2:2651076/150‑1 (MQ=255)
                                                    tCGCCATAAAAACACCACGAACGTCTTTGCGCATGTCAAAGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGCATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCa                                                                                              <  1:254948/150‑1 (MQ=255)
                                                      gCCATAAAAACACCACGAACGTCTTTGCGCATGTCAAAGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGCATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAAt                                                                                            <  2:2196818/150‑1 (MQ=255)
                                                       ccATAAAAACACCACGAACGTCTTTGCGCATGTCAAAGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGCATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATg                                                                                           >  1:3837017/1‑150 (MQ=255)
                                                            aaaaCACCACGAACGTCTTTGCGCATGTCAAAGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGCATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATgttagt                                                                                      >  2:774591/1‑150 (MQ=255)
                                                            aaaaCACCACGAACGTCTTTGCGCATGTCAAAGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGCATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATgttagt                                                                                      <  2:818410/150‑1 (MQ=255)
                                                             aaaCACCACGAACGTCTTTGCGCATGTCAAAGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGCATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATgttagtt                                                                                     >  1:1044206/1‑150 (MQ=255)
                                                               acacCACGAACGTCTTTGCGCATGTCAAAGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGCATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATGttagttag                                                                                   <  1:2671175/150‑1 (MQ=255)
                                                                 accacGAACGTCTTTGCGCATGTCAAAGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGCATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATGTTAGTTAgcg                                                                                 <  2:655285/150‑1 (MQ=255)
                                                                   cacGAACGTCTTTGCGCATGTCAAAGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGCATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATGTTAGTTAGCGcc                                                                               <  2:5110875/150‑1 (MQ=255)
                                                                    acGAACTTATTTGCGCATGTCAAAGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGCATAAAAAAAATGTTGCGGTTGGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATGTTAGTTAGCGCca                                                                              <  2:4769888/150‑1 (MQ=255)
                                                                       aaCGTCTTTGCGCATGTCAAAGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGCATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATGTTAGTTAGCGCCACaa                                                                           <  2:2586863/150‑1 (MQ=255)
                                                                         cGTCTTTGCGCATGTCAAAGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGCATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATGTTAGTTAGCGCCACAAAg                                                                         <  2:4467754/150‑1 (MQ=255)
                                                                         cGTCTTTGCGCATGTCAAAGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGCATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATGTTAGTTAGCGCCACAAAg                                                                         >  1:2481430/1‑150 (MQ=255)
                                                                             tttGCGCATGTCAAAGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGCATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATGTTAGTTAGCGCCACAAAGGaaa                                                                     >  1:5238497/1‑150 (MQ=255)
                                                                               tGCGCATGTCAAAGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGCATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATGTTAGTTAGCGCCACAAAGGAAACa                                                                   >  2:2725286/1‑150 (MQ=255)
                                                                               tGCGCATGTCAAAGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACAGTGATATTTCGCATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATGTTAGTTAGCGCCACAAAGGAAACa                                                                   >  2:35062/1‑150 (MQ=255)
                                                                                 cgcATGTCAAAGTTCATGCTCTATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGCATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATGTTAGTTTGCGCCACAAAGGAAACAGc                                                                 <  1:10099/150‑1 (MQ=255)
                                                                                   cATGTCAAAGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGCATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATGTTAGTTAGCGCCACAAAGGAAACAGCAt                                                               >  2:3619749/1‑150 (MQ=255)
                                                                                    aTGTCAAAGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGCATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATGTTAGTTAGCGCCACAAAGGAAACAGCAtt                                                              >  1:5302110/1‑150 (MQ=255)
                                                                                     tGTCAAAGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGCATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATGTTAGTTAGCGCCACAAAGGAAACAGCATTc                                                             >  2:4056903/1‑150 (MQ=255)
                                                                                       tCAAAGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGCATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATGTTAGTTAGCGCCACAAAGGAAACAGCATTccc                                                           <  2:2863991/150‑1 (MQ=255)
                                                                                       tCAAAGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGCATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATGTTAGTTAGCGCCACAAAGGAAACAGCATTccc                                                           <  1:4562844/150‑1 (MQ=255)
                                                                                        cAAAGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGCATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATGTTAGTTAGCGCCACAAAGGAAACAGCATTCCcg                                                          >  1:2951459/1‑150 (MQ=255)
                                                                                        cAAAGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGCATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATGTTAGTTAGCGCCACAAAGGAAACAGCATTCCcg                                                          >  1:2837413/1‑150 (MQ=255)
                                                                                        cAAAGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGCATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATGTTAGTTAGCGCCACAAAGGAAACAGCATTCCcg                                                          >  2:456265/1‑150 (MQ=255)
                                                                                         aaaGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGCATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATGTTAGTTAGCGCCACAAAGGAAACAGCATTCCcgg                                                         <  1:2420593/150‑1 (MQ=255)
                                                                                           aGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGCATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATGTTAGTTAGCGCCACAAAGGAAACAGCATTCCcggcg                                                       >  1:5028825/1‑150 (MQ=255)
                                                                                              tCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGCATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATGTTAGTTAGCGCCACAAAGGAAACAGCATTCCCGGCGgag                                                    >  1:3979961/1‑150 (MQ=255)
                                                                                               cATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGCATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATGTTAGTTAGCGCCACAAAGGAAACAGCATTCCCGGAGgaga                                                   >  1:1085182/1‑150 (MQ=255)
                                                                                                aTGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGCATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATGTTAGTAAGCGCCACAAAGGAAACAGCATTCCCGGCGGAGAt                                                  >  1:4352221/1‑150 (MQ=255)
                                                                                                  gCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGCATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATGTTAGTTAGCGCCACAAAGGAAACAGCATTCCCGGCGGAGAttt                                                >  2:3156740/1‑150 (MQ=255)
                                                                                                    tCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGCATAAAAAAAAATGTTGCGGTGTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATGTTAGTTAGCGCCACAAAGGACACAGCATTCCCGGCGGAGATTTg                                               >  2:19370/1‑150 (MQ=255)
                                                                                                    tCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGCATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATGTTAGTTAGCGCCACAAAGGAAACAGCATTCCCGGCGGAGATTTGa                                              >  1:4905827/1‑150 (MQ=255)
                                                                                                      aaTGCGTGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGCATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATGTTAGTTAGCGCCACAAAGGAAACAGCATTCCCGGCGGAGATTTGAAt                                            <  2:1258450/150‑1 (MQ=255)
                                                                                                      aaTGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGCATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATGTTAGTTAGCGCCACAAAGGAAACAGCATTCCCGGCGGAGATTTGAAt                                            <  1:4238607/150‑1 (MQ=255)
                                                                                                      aaTGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGCATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATGTTAGTTAGCGCCACAAAGGAAACAGCATTCCCGGCGGAGATTTGAAt                                            <  1:4163481/150‑1 (MQ=255)
                                                                                                          cGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGCATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATGTTAGTTAGCGCCACAAAGGAAACAGCATTCCCGGCGGAGATTTGAATGCcg                                        <  2:4201306/150‑1 (MQ=255)
                                                                                                                 aaCTGAGTATTCCTGCATCACATTGATATTTCGCATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACTCTTCATCTCAATGGTAGTTAGCGCCACAAAGGACACAGCATTCCCGGCGGAGATTTGTATGCCGCGCAAtt                                 >  1:4070987/1‑150 (MQ=255)
                                                                                                                  aCTGAGGATTCCTGCATCACATTGATATTTCGCATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATGTTAGTTAGCGCCACAAAGGAAACAGCATTCCCGGCGGAGATTTGAATGCCGCGCAAttt                                >  2:2531212/1‑150 (MQ=255)
                                                                                                                    tGAGGATTCCTGCATCACATTGATATTTCGCATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATGTTAGTTAGCGCCACAAAGGAAACAGCATTCCCGGCGGAGATTTGAATGCCGCGCAATTtat                              <  2:1314880/150‑1 (MQ=255)
                                                                                                                     gAGGATTCCTGCATCACATTGATATTTCGCATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATGTTAGTTAGCGCCACAAAGGAAACAGCATTCCCGGCGGAGATTTGAATGCCGCGCAATTtata                             <  1:3200425/150‑1 (MQ=255)
                                                                                                                     gAGGATTCCTGCATCACATTGATATTTCGCATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATGTTAGTTAGCGCCACAAAGGAAACAGCATTCCCGGCGGAGATTTGAATGCCGCGCAATTtata                             <  1:5155790/150‑1 (MQ=255)
                                                                                                                       ggATTCCTGCATCACATTGATATTTCGCATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATGTTAGTTAGCGCCACAAAGGAAACAGCATTCCCGGCGGAGATTTGAATGCCGCGCAATTTATAAg                           <  1:757550/150‑1 (MQ=255)
                                                                                                                        gATTCCTGCATCACATTGATATTTCGCATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATGTTAGTTAGCGCCACAAAGGAAACAGCATTCCCGGCGGAGATTTGAATGCCGCGCAATTTATAAGc                          <  1:2398412/150‑1 (MQ=255)
                                                                                                                        gATTCCTGCATCACATTGATATTTCGCATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATGTTAGTTAGCGCCACAAAGGAAACAGCATTCCCGGCGGAGATTTGAATGCCGCGCAATTTATAAGc                          <  1:1709983/150‑1 (MQ=255)
                                                                                                                         aTTCCTGCATCACATTGATATTTCGCATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATGTTAGTTAGCGCCACAAAGGAAACAGCATTCCCGGCGGAGATTTGAATGCCGCGCAATTTATAAGCa                         >  1:1208897/1‑150 (MQ=255)
                                                                                                                             cTGCATCACATTGATATTTCGCATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATGTTAGTTAGCGCCACAAAGGAAACAGCATTCCCGGCGGAGATTTGAATGCCGCGCAATTTATAAGCAAGAt                     <  2:686602/150‑1 (MQ=255)
                                                                                                                             cTGCATCACATTGATATTTCGCATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATGTTAGTTAGCGCCACAAAGGAAACAGCATTCCCGGCGGAGATTTGAATGCCGCGCAATTTATAAGCAAGAt                     >  2:398036/1‑150 (MQ=255)
                                                                                                                             cTGCATCACATTGATATTTCGCATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATGTTAGTTAGCGCCACAAAGGAAACAGCATTCCCGGCGGAGATTTGAATGCCGCGCAATTTATAAGCAAGAt                     <  1:4498656/150‑1 (MQ=255)
                                                                                                                             cTGCATCACATTGATATTTCGCATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATGTTAGTTAGCGCCACAAAGGAAACAGCATTCCCGGCGGAGATTTGAATGCCGCGCAATTTATAAGCAAGAt                     >  1:1349045/1‑150 (MQ=255)
                                                                                                                                cATCACATTGATATTTCGCATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATGTTAGTTAGCGCCACAAAGGAAACAGCATTCCCGGCGGAGATTTGAATGCCGCGCAATTTATAAGCAAGATGAg                  <  2:2621575/150‑1 (MQ=255)
                                                                                                                                 aTCACATTGATATTTCGCATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATGTTAGTTAGCGCCACAAAGGAAACAGCATTCCCGGCGGAGATTTGAATGCCGCGCAATTTATAAGCAAGATGAGt                 >  1:3617731/1‑150 (MQ=255)
                                                                                                                                 aTCACATTGATATTTCGCATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATGTTAGTTAGCGCCACAAAGGAAACAGCATTCCCGGCGGAGATTTGAATGCCGCGCAATTTATAAGCAAGATGAGt                 >  1:2102937/1‑150 (MQ=255)
                                                                                                                                 aTCACATTGATATTTCGCATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATGTTAGTTAGCGCCACAAAGGAAACAGCATTCCCGGCGGAGATTTGAATGCCGCGCAATTTATAAGCAAGATGAGt                 >  1:4646894/1‑150 (MQ=255)
                                                                                                                                   cacaTTGATATTTCGCATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATGTTAGTTAGCGCCACAAAGGAAACAGCATTCCCGGCGGAGATTTGAATGCCGCGCAATTTATAAGCAAGATGAGtat               <  2:3378767/150‑1 (MQ=255)
                                                                                                                                    acaTTGATATTTCGCATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATGTTAGTTAGCGCCACAAAGGAAACAGCATTCCCGGCGGAGATTTGAATGCCGCGCAATTTATAAGCAAGATGAGtatt              <  2:3256527/150‑1 (MQ=255)
                                                                                                                                    acaTTGATATTTCGCATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATGTTAGTTAGCGCCACAAAGGAAACAGCATTCCCGGCGGAGATTTGAATGCCGCGCAATTTATAAGCAAGATGAGtatt              >  2:368476/1‑150 (MQ=255)
                                                                                                                                      aTTGATATTTCGCATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATGTTAGTTAGCGCCACAAAGGAAACAGCATTCCCGGCGGAGATTTGAATGCCGCGCAATTTATAAGCAAGATGAGtattat            <  1:3720803/150‑1 (MQ=255)
                                                                                                                                       ttGATATTTCGCATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATGTTAGTTAGCGCCACAAAGGAAACAGCATTCCCGGCGGAGATTTGAATGCCGCGCAATTTATAAGCAAGATGAGTATTATc           <  2:5221522/150‑1 (MQ=255)
                                                                                                                                        tGATATTTCGCATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATGTTAGTTAGCGCCACAAAGGAAACAGCATTCCCGGCGGAGATTTGAATGCCGCGCAATTTATAAGCAAGATGAGTATTATcc          >  1:3349197/1‑150 (MQ=255)
                                                                                                                                        tGATATTTCGCATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATGTTAGTTAGCGCCACAAAGGAAACAGCATTCCCGGCGGAGATTTGAATGCCGCGCAATTTATAAGCAAGATGAGTATTATcc          <  2:1877999/150‑1 (MQ=255)
                                                                                                                                         gATATTTCGCATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATGTTAGTTAGCGCCACAAAGGAAACAGCATTCCCGGCGGAGATTTGAATGCCGCGCAATTTATAAGCAAGATGAGTATTATCCa         <  2:5053654/150‑1 (MQ=255)
                                                                                                                                         gATATTTCGCATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATGTTAGTTAGCGCCACAAAGGAAACAGCATTCCCGGCGGAGATTTGAATGCCGCGCAATTTATAAGCAAGATGAGTATTATCCa         <  1:2135052/150‑1 (MQ=255)
                                                                                                                                         gATATTTCGCATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATGTTAGTTAGCGCCACAAAGGAAACAGCATTCCCGGCGGAGATTTGAATGCCGCGCAATTTATAAGCAAGATGAGTATTATCCa         <  1:3470292/150‑1 (MQ=255)
                                                                                                                                         gATATTTCGCATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATGTTAGTTAGCGCCACAAAGGAAACAGCATTCCCGGCGGAGATTTGAATGCCGCGCAATTTATAAGCAAGATGAGTATTATCCa         >  1:2201691/1‑150 (MQ=255)
                                                                                                                                          atatTTCGCATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATGTTAGTTAGCGCCACAAAGGAAACAGCATTCCCGGCGGAGATTTGAATGCCGCGCAATTTATAAGCAAGATGAGTATTATCCAg        >  1:737457/1‑150 (MQ=255)
                                                                                                                                              ttCGCATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATGTTAGTTAGCGCCACAAAGGAAACAGCATTCCCGGCGGAGATTTGAATGCCGCGCAATTTATAAGCAAGATGAGTATTATCCAGTTGa    >  2:2044567/1‑150 (MQ=255)
                                                                                                                                                cGCATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATGTTAGTTAGCGCCACAAAGGAAACAGCATTCCCGGCGGAGATTTGAATGCCGCGCAATTTATAAGCAAGATGAGTATTATCCAGTTGAat  <  2:4099643/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                 
CCTTTTTCATTCACCGCATGAACATTTGCAAGAGACAGCAGTGTTTCTTTTTTCGCCATAAAAACACCACGAACGTCTTTGCGCATGTCAAAGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGCAT‑‑AAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATGTTAGTTAGCGCCACAAAGGAAACAGCATTCCCGGCGGAGATTTGAATGCCGCGCAATTTATAAGCAAGATGAGTATTATCCAGTTGAAT  >  CP009273/2108071‑2108363

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: