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% 976.8 / 24.4 280coding (1159/1395 nt)wcaMcolanic acid biosynthesis protein
Reads supporting (aligned to +/- strand):  ref base . (5/4);  new base A (122/149);  total (127/153)
Fisher's exact test for biased strand distribution p-value = 7.36e-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 309 total aligned reads displayed.

TCCTTTTTCATTCACCGCATGAACATTTGCAAGAGACAGCAGTGTTTCTTTTTTCGCCATAAAAACACCACGAACGTCTTTGCGCATGTCAAAGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGCAT‑‑AAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATGTTAGTTAGCGCCACAAAGGAAACAGCATTCCCGGCGGAGATTTGAATGCCGCGCAATTTATAAGCAAGATGAGTATTATCCAGTTGAAT  >  CP009273/2108070‑2108363
                                                                                                                                                      |                                                                                                                                                 
    ttttCATTCACCGCATGAACATTTGCAAGAGACAGCAGTGTTTCTTTTTTCGCCATAAAAACACCACGAACGTCTTTGCGCATGTCAAAGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGCAT‑‑aaaa                                                                                                                                              >  2:4058864/1‑150 (MQ=255)
     tttCATTCACCGCATGAACATTTGCAAGAGACAGCAGTGTTTCTTTTTTCGCCATAAAAACACCACGAACGTCTTTGCGCATGTCAAAGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGCAT‑‑aaaaa                                                                                                                                             >  2:1858066/1‑150 (MQ=255)
     tttCATTCACCGCATGAACATTTGCAAGAGACAGCAGTGTTTCTTTTTTCGCCATAAAAACACCACGAACGTCTTTGCGCATGTCAAAGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGCAT‑‑aaaaa                                                                                                                                             <  1:2073486/150‑1 (MQ=255)
        cattcaCCGCATGAACATTTGCAAGAGACAGCAGTGTTTCTTTTTTCGCCATAAAAACACCACGAACGTCTTTGCGCATGTCAAAGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTTATATTTCGCAT‑‑aaaaaaaa                                                                                                                                          >  1:2278095/1‑149 (MQ=255)
         attcaCCGCATGAACATTTGCAAGAGACAGCAGTGTTTCTTTTTTCGCCATAAAAACACCACGAACGTCTTTGCGCATGTCAAAGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGTATCACATTGATATTTCGCAT‑‑aaaaaaaat                                                                                                                                         <  2:3280418/150‑3 (MQ=255)
         attcaCCGCATGAACATTTGCAAGAGACAGCAGTGTTTCTTTTTTCGCCATAAAAACACCACGAACGTCTTTGCGCATGTCAAAGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGCAT‑‑aaaaaaaat                                                                                                                                         <  2:1969443/150‑3 (MQ=255)
         attcaCCGCATGAACATTTGCAAGAGACAGCAGTGTTTCTTTTTTCGCCATAAAAACACCACGAACGTCTTTGCGCATGTCAAAGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGCAT‑‑aaaaaaaat                                                                                                                                         >  1:11563/1‑148 (MQ=255)
           tcaCCGCATGAACATTTGCAAGAGACAGCAGTGTTTCTTTTTTCGCCATAAAAACACCACGAACGTCTTTGCGCATGTCAAAGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGCAT‑‑aaaaaaaatgt                                                                                                                                       <  1:461491/150‑5 (MQ=255)
            caCCGCATGAACATTTGCAAGAGACAGCAGTGTTTCTTTTTTCGCCATAAAAACACCACGAACGTCTTTGCGCATGTCAAAGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGCAT‑‑aaaaaaaatgtt                                                                                                                                      >  2:1838846/1‑145 (MQ=255)
             aCCGCATGAACATTTGCAAGAGACAGCAGTGTTTCTTTTTTCGCCATAAAAACACCACGAACGTCTTTGCGCATGTCAAAGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGCATAAAAAAAATGTTg                                                                                                                                      <  1:1234493/150‑1 (MQ=255)
              ccGCATGAACATTTGCAAGAGACAGCAGTGTTTCTTTTTTCGCCATAAAAACACCACGAACGTCTTTGCGCATGTCAAAGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGTATCACATTGATATTTCGCATAAAAAAAATGTTGc                                                                                                                                     <  1:2829411/150‑1 (MQ=255)
               cGCATGAACATTTGCAAGAGACAGCAGTGTTTCTTTTTTCGCCATAAAAACACCACGAACGTCTTTGCGCATGTCAAAGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGTATCACATTGATATTTCGCATAAAAAAAATGTTGCg                                                                                                                                    >  1:3194050/1‑150 (MQ=255)
               cGCATGAACATTTGCAAGAGACAGCAGTGTTTCTTTTTTCGCCATAAAAACACCACGAACGTCTTTGCGCATGTCAAAGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGCATAAAAAAAATGTTGCg                                                                                                                                    <  2:1805604/150‑1 (MQ=255)
                 cATGAACATTTGCAAGAGACAGCAGTGTTTCTTTTTTCGCCATAAAAACACCACGAACGTCTTTGCGCATGTCAAAGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGCATAAAAAAAATGTTGCGGt                                                                                                                                  <  2:648547/150‑1 (MQ=255)
                      aCATTTGCAAGAGACAGCAGTGTTTCTTTTTTCGCCATAAAAACACCACGAACGTCTTTGCGCATGTCAAAGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGCATAAAAAAAATGTTGCGGTttgtt                                                                                                                             <  1:3601259/150‑1 (MQ=255)
                       cATTTGCAAGAGACAGCAGTGTTTCTTTTTTCGCCATAAAAACACCACGAACGTCTTTGCGCATGTCAAAGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGCATAAAAAAAATGTTGCGGTttgttg                                                                                                                            >  1:3557605/1‑150 (MQ=255)
                       cATTTGCAAGAGACAGCAGTGTTTCTTTTTTCGCCATAAAAACACCACGAACGTCTTTGCGCATGTCAAAGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGCATAAAAAAAATGTTGCGGTttgttg                                                                                                                            <  1:845776/150‑1 (MQ=255)
                        aTTTGCAAGAGACAGCAGTGTTTCTTTTTTCGCCATAAAAACACCACGAACGTCTTTGCGCATGTCAAAGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGTATCACATTGATATTTCGCATAAAAAAAATGTTGCGGTTtgttgt                                                                                                                           >  2:3120938/1‑150 (MQ=255)
                        aTTTGCAAGAGACAGCAGTGTTTCTTTTTTCGCCATAAAAACACCACGAACGTCTTTGCGCATGTCAAAGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGCATAAAAAAAATGTTGCGGTGtgttgt                                                                                                                           >  2:4482961/1‑150 (MQ=255)
                        aTTTGCAAGAGACAGCAGTGTTTCTTTTTTCGCCATAAAAACACCACGAACGTCTTTGCGCATGTCAAAGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGCATAAAAAAAATGTTGAGGTTtgttgt                                                                                                                           <  2:1999707/150‑1 (MQ=255)
                          ttGCAAGAGACAGCAGTGTTTCTTTTTTCGCCATAAAAACACCACGAACGTCTTTGCGCATGTCAAAGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGCATAAAAAAAATGTTGCGGTTTGTTgtgt                                                                                                                         <  1:145508/150‑1 (MQ=255)
                             cAAGAGACAGCAGTGTTTCTTTTTTCGCCATAAAAACACCACGAACGTCTTTGCGCATGTCAAAGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGCATAAAAAAAATGTTGCGGTTTGTTGTGTAAc                                                                                                                      >  1:2601224/1‑150 (MQ=255)
                                   aCAGCAGTGTTTCTTTTTTCGCCATAAAAACACCACGAACGTCTTTGCGCATGTCAAAGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGCATAAAAAAAATGTTGCGGTCTGTTGTGTAACTCCAGc                                                                                                                <  1:946424/150‑1 (MQ=255)
                                    cagcagTGTTTCTTTTTTCGCCATAAAAACACCACGAACGTCTTTGCGCATGTCAAAGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGCATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCg                                                                                                               <  2:2570708/150‑1 (MQ=255)
                                      gcagTGTTTCTTTTTTCGCCATAAAAACACCACGAACGTCTTTGCGCATGTCAAAGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGTATCACATTGATATTTCGCATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGac                                                                                                             >  1:2853731/1‑150 (MQ=255)
                                      gcagTGTTTCTTTTTTCGCCATAAAAACACCACGAACGTCTTTGCGCATGTCAAAGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGCATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGac                                                                                                             <  1:1400815/150‑1 (MQ=255)
                                         gtgtTTCTTTTTTCGCCATAAAAACACCACGAACGTCTTTGCGCATGTCAAAGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGCATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGacgca                                                                                                          <  2:2090999/150‑1 (MQ=255)
                                         gtgtTTCTTTTTTCGCCATAAAAACACCACGAACGTCTTTGCGCATGTCAAAGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGCATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGacgca                                                                                                          >  2:4132257/1‑150 (MQ=255)
                                           gtTTCTTTTTTCGCCATAAAAACACCACGAACGTCTTTGCGCATGTCAAAGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGCATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGacgcacg                                                                                                        >  1:2296552/1‑150 (MQ=255)
                                            tttCTTTTTTCGCCATAAAAACACCACGAACGTCTTTGCGCATGTCAAAGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGTATCACATTGATATTTCGCATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGacgcacgc                                                                                                       >  1:4478826/1‑150 (MQ=255)
                                            tttCTTTTTTCGCCATAAAAACACCACGAACGTCTTTGCGCATGTCAAAGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGCATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGacgcacgc                                                                                                       >  2:3214055/1‑150 (MQ=255)
                                               cTTTTTTCGCCATAAAAACACCACGAACGTCTTTGCGCATGTCAAAGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGTATCACATTGATATTTCGCATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTc                                                                                                    >  1:671032/1‑150 (MQ=255)
                                               cTTTTTTCGCCATAAAAACACCACGAACGTCTTTGCGCATGTCAAAGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGCATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTc                                                                                                    >  2:2643834/1‑150 (MQ=255)
                                                 tttttCGCCATAAAAACACCACGAACGTCTTTGCGCATGTCAAAGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGCATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCAt                                                                                                  <  1:4433560/150‑1 (MQ=255)
                                                  ttttCGCCATAAAAACACCACGAACGTCTTTGCGCATGTCAAAGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGCATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCAtc                                                                                                 <  2:1808858/150‑1 (MQ=255)
                                                    ttCGCCATAAAAACACCACGAACGTCTTTGCGCATGTCAAAGTTCATGCTCAATGCGGGTCCAACGGAGGATTCCTGCATCACATTGATATTTCGCATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCAtctc                                                                                               >  1:2213627/1‑150 (MQ=255)
                                                      cGCCATAAAAACACCACGAACGTCTTTGCGCATGTCAAAGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGTATCACATTGATATTTCGCATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCaa                                                                                             >  1:3488716/1‑150 (MQ=255)
                                                       gCCATAAAAACACCACGAACGTCTTTGCGCATGTCAAAGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGCATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAAt                                                                                            <  1:1197616/150‑1 (MQ=255)
                                                        ccATAAAAACACCACGAACGTCTTTGCGCATGTCAAAGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGCATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATg                                                                                           >  2:690623/1‑150 (MQ=255)
                                                         cATAAAAACACCACGAACGTCTTTGCGCATGTCAAAGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGTATCACATTGATATTTCGCATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATgt                                                                                          >  2:636266/1‑150 (MQ=255)
                                                         cATAAAAACACCACGAACGTCTTTGCGCATGTCAAAGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGCATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATgt                                                                                          >  2:1057210/1‑150 (MQ=255)
                                                           tAAAAACACCACGAACGTCTTTGCGCATGTCAAAGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGCATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATgtta                                                                                        <  1:2982642/150‑1 (MQ=255)
                                                             aaaaCACCACGAACGTCTTTGCGCATGTCAAAGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGCATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATgttagt                                                                                      >  1:697548/1‑150 (MQ=255)
                                                              aaaCACCACGAACGTCTTTGCGCATGTCAAAGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGTATCACATTGATATTTCGCATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATgttagtt                                                                                     <  1:1064190/150‑1 (MQ=255)
                                                              aaaCACCACGAACGTCTTTGCGCATGTCAAAGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGCATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATgttagtt                                                                                     >  1:2418789/1‑150 (MQ=255)
                                                                   ccacGAACGTCTTTGCGCATGTCAAAGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGCATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATGTTAGTTAgcgc                                                                                >  2:344754/1‑150 (MQ=255)
                                                                    cacGAACGTCTTTGCGCATGTCAAAGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGCATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATGTTAGTTAGCGcc                                                                               <  1:2218939/150‑1 (MQ=255)
                                                                    cacGAACGTCTTTGCGCATGTCAAAGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGCATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATGTTAGTTAGCGcc                                                                               <  1:4546755/150‑1 (MQ=255)
                                                                     acGAACGTCTTTGCGCATGTCAAAGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGCATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATGTTAGTTAGCGCca                                                                              <  1:3943281/150‑1 (MQ=255)
                                                                        aaCGTCTTTGCGCATGTCAAAGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGCATAAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATGTTAGTTAGCGCcaca                                                                            >  1:1124104/1‑150 (MQ=255)
                                                                          cGTCTTTGCGCATGTCAAAGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGCATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATGTTAGTTAGCGCCACAAAg                                                                         >  2:3158172/1‑150 (MQ=255)
                                                                           gTCTTTGCGCATGTCAAAGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGTATCACATTGATATTTCGCATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATGTTAGTTAGCGCCACAAAgg                                                                        >  2:519496/1‑150 (MQ=255)
                                                                             cTTTGCGCATGTCAAAGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGCATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATGTTAGTTAGCGCCACAAAGGaa                                                                      >  1:431586/1‑150 (MQ=255)
                                                                             cTTTGCGCATGTCAAAGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGCATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATGTTAGTTAGCGCCACAAAGGaa                                                                      >  1:3065273/1‑150 (MQ=255)
                                                                              tttGCGCATGTCAAAGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGTATCACATTGATATTTCGCATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATGTTAGTTAGCGCCACAAAGGaaa                                                                     >  1:3700307/1‑150 (MQ=255)
                                                                              tttGCGCATGTCAAAGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGCATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATGTTAGTTAGCGCCACAAAGGaaa                                                                     >  2:3857310/1‑150 (MQ=255)
                                                                               ttGCGCATGTCAAAGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGCATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATGTTAGTTAGCGCCACAAAGGAAAc                                                                    <  1:769561/150‑1 (MQ=255)
                                                                                  cgcATGTCAAAGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGCATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATGTTAGTTAGCGCCACAAAGGAAACAGc                                                                 >  1:1255616/1‑150 (MQ=255)
                                                                                     aTGTCAAAGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGCATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATGTTAGTTAGCGCCACAAAGGAAACAGCAtt                                                              >  1:3880781/1‑150 (MQ=255)
                                                                                       gTCAAAGTTCATGCTCACTGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGCATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATGTTAGTTAGCGCCACAAAGGAAACAGCATTcc                                                            <  1:313563/150‑1 (MQ=255)
                                                                                        tCAAAGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGTATCACATTGATATTTCGCATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATGTTAGTTAGCGCCACAAAGGAAACAGCATTccc                                                           <  1:2770908/150‑1 (MQ=255)
                                                                                        tCAAAGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGCATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATGTTAGTTAGCGCCACAAAGGAAACAGCATTccc                                                           <  1:1887131/150‑1 (MQ=255)
                                                                                         cAAAGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGCATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATGTTAGTTAGCGCCACAAAGGAAACAGCATTCCcg                                                          <  1:4410906/150‑1 (MQ=255)
                                                                                         cAAAGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGCATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATGTTAGTTAGCGCCACAAAGGAAACAGCATTCCcg                                                          <  2:616259/150‑1 (MQ=255)
                                                                                            aGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGCATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATGTTAGTTAGCGCCACAAAGGAAACAGCATTCCcggcg                                                       <  2:4322642/150‑1 (MQ=255)
                                                                                              ttCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGCATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATGTTAGTTAGCGCCACAAAGGAAACAGCATTCCCGGCGga                                                     >  1:4267669/1‑150 (MQ=255)
                                                                                               tCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGCATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATGTTAGTTAGCGCCACAAAGGAAACAGCATTCCCGGCGgag                                                    >  2:1552635/1‑150 (MQ=255)
                                                                                               tCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGCATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATGTTAGTTAGCGCCACAAAGGAAACAGCATTCCCGGCGgag                                                    >  1:1280050/1‑150 (MQ=255)
                                                                                               tCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGCATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATGTTAGTTAGCGCCACAAAGGAAACAGCATTCCCGGCGgag                                                    <  1:1055164/150‑1 (MQ=255)
                                                                                                 aTGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGCATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATGTTAGTTAGCGCCACAAAGGAAACAGCATTCCCGGCGGAGAt                                                  >  1:1589321/1‑150 (MQ=255)
                                                                                                  tGCTCAATGCGGGTCCAACTGAGGATTCCTGTATCACATTGATATTTCGCATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATGTTAGTTAGCGCCACAAAGGAAACAGCATTCCCGGCGGAGAtt                                                 <  2:2305605/150‑1 (MQ=255)
                                                                                                  ggCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGCATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATGTTAGTTAGCGCCACAAAGGAAACAGCATTCCCGGCGGAGAtt                                                 <  1:1902581/149‑1 (MQ=255)
                                                                                                    cTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGCATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATGTTAGTTAGCGCCACAAAGGAAACAGCATTCCCGGCGGAGATTTg                                               <  1:506534/150‑1 (MQ=255)
                                                                                                       aaTGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGCATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATGTTAGTTAGCGCCACAAAGGAAACAGCATTCCCGGCGGAGATTTGAAt                                            <  2:1021387/150‑1 (MQ=255)
                                                                                                        aTGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGCATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATGTTAGTTAGCGCCACAAAGGAAACAGCATTCCCGGCGGAGATTTGAATg                                           <  1:3622462/150‑1 (MQ=255)
                                                                                                          gCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGCATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATGTTAGTTAGCGCCACAAAGGAAACAGCATTCCCGGCGGAGATTTGAATGcc                                         <  1:1655413/150‑1 (MQ=255)
                                                                                                           cGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGCATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATGTTAGTTAGCGCCACAAAGGAAACAGCATTCCCGGCGGAGATTTGAATGCcg                                        >  1:2696942/1‑150 (MQ=255)
                                                                                                            gggTCCAACTGAGGATTCCTGCATCACATTGATATTTCGCATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATGTTAGTTAGCGCCACAAAGGAAACAGCATTCCCGGCGGAGATTTGAATGCcgc                                       <  2:11767/150‑1 (MQ=255)
                                                                                                                 cAACTGAGGATTCCTGCATCACATTGATATTTCGCATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATGTTAGTTAGCGCCACAAAGGAAACAGCATTCCCGGCGGAGATTTGAATGCCGCGCAAt                                  >  1:1787043/1‑150 (MQ=255)
                                                                                                                   aCTGAGGATTCCTGCATCACATTGATATTTCGCATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATGTTAGTTAGCGCCACAAAGGAAACAGCATTCCCGGCGGAGATTTGAATGCCGCGCAAttt                                <  2:2208210/150‑1 (MQ=255)
                                                                                                                    cTGAGGATTCCTGCATCACATTGATATTTCGCATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATGTTAGTTAGCGCCACAAAGGAAACAGCATTCCCGGCGGAGATTTGAACGCCGCGCAATTta                               <  2:2667024/150‑1 (MQ=255)
                                                                                                                      gAGGATTCCTGCATCACATTGATATTTCGCATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATGTTAGTTAGCGCCACAAAGGAAACAGCATTCCCGGCGGAGATTTGAATGCCGCGCAATTtata                             <  2:2122074/150‑1 (MQ=255)
                                                                                                                        ggATTCCTGCATCACATTGATATTTCGCATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATGTTAGTTAGCGCCACAAAGGAAACAGCATTCCCGGCGGAGATTTGAATGCCGCGCAATTTATAAg                           <  1:1526988/150‑1 (MQ=255)
                                                                                                                        ggATTCCTGCATCACATTGATATTTCGCATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATGTTAGTTAGCGCCACAAAGGAAACAGCATTCCCGGCGGAGATTTGAATGCCGCGCAATTTATAAg                           <  1:3111982/150‑1 (MQ=255)
                                                                                                                               tGCATCACATTGATATTTCGCATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATGTTAGTTAGCGCCACAAAGGAAACAGCATTCCCGGCGGAGATTTGAATGCCGCGCAATTTATAAGCAAGATg                    <  1:2040250/150‑1 (MQ=255)
                                                                                                                                  aTCACATTGATATTTCGCATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATGTTAGTTAGCGCCACAAAGGAAACAGCATTCCCGGCGGAGATTTGAATGCCGCGCAATTTATAAGCAAGATGAGt                 <  2:238913/150‑1 (MQ=255)
                                                                                                                                  aTCACATTGATATTTCGCATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATGTTAGTTAGCGCCACAAAGGAAACAGCATTCCCGGCGGAGATTTGAATGCCGCGCAATTTATAAGCAAGATGAGt                 <  2:594848/150‑1 (MQ=255)
                                                                                                                                    cacaTTGATATTTCGCATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATGTTAGTTAGCGCCACAAAGGAAACAGCATTCCCGGCGGAGATTTGAATGCCGCGCAATTTATAAGCAAGATGAGtat               >  1:3719364/1‑150 (MQ=255)
                                                                                                                                    cacaTTGATATTTCGCATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATGTTAGTTAGCGCCACAAAGGAAACAGCATTCCCGGCGGAGATTTGAATGCCGCGCAATTTATAAGCAAGATGAGtat               >  1:920546/1‑150 (MQ=255)
                                                                                                                                     acaTTGATATTTCGCATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATGTTAGTTAGCGCCACAAAGGAAACAGCATTCCCGGCGGAGATTTGAATGCCGCGCAATTTATAAGCAAGATGAGtatt              <  1:2902607/150‑1 (MQ=255)
                                                                                                                                       aTTGATATTTCGCATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATGTTAGTTAGCGCCACAAAGGAAACAGCATTCCCGGCGGAGATTTGAATGCCGCGCAATTTATAAGCAAGATGAGtattat            <  1:2892160/150‑1 (MQ=255)
                                                                                                                                       aTTGATATTTCGCATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATGTTAGTTAGCGCCACAAAGGAAACAGCATTCCCGGCGGAGATTTGAATGCCGCGCAATTTATAAGCAAGATGAGtattat            <  1:4131814/150‑1 (MQ=255)
                                                                                                                                       aTTGATATTTCGCATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCAACTCAATGTAAGTTAGCGCCACAAAGGAAACAGCATTCCCGGCGGAGATTTGAATGCCGCGCAATTTATAAGCAAGATGAGtattat            >  1:3370488/1‑150 (MQ=255)
                                                                                                                                         tGATATTTCGCATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATGTTAGTTAGCGCCACAAAGGAAACAGCATTCCCGGCGGAGATTTGAATGCCGCGCAATTTATAAGCAAGATGAGTATTATcc          <  2:1922900/150‑1 (MQ=255)
                                                                                                                                         tGATATTTCGCATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATGTTAGTTAGCGCCACAAAGGAAACAGCATTCCCGGCGGAGATTTGAATGCCGCGCAATTTATAAGCAAGATGAGTATTATcc          >  2:129526/1‑150 (MQ=255)
                                                                                                                                            tatTTCGCATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATGTTAGTTAGCGCCACAAAGGAAACAGCATTCCCGGCGGAGATTTGAATGCCGCGCAATTTATAAGCAAGATGAGTATTATCCAGt       <  1:3811915/150‑1 (MQ=255)
                                                                                                                                            tatTTCGCATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATGTTAGTTAGCGCCACAAAGGAAACAGCATTCCCGGCGGAGATTTGAATGCCGCGCAATTTATAAGCAAGATGAGTATTATCCAGt       >  1:29558/1‑150 (MQ=255)
                                                                                                                                             atTTCGCATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATGTTAGTTAGCGCCACAAAGGAAACAGCATTCCCGGCGGAGATTTGAATGCCGCGCAATTTATAAGCAAGATGAGTATTATCCAGtt      >  1:3635105/1‑150 (MQ=255)
                                                                                                                                             atTTCGCATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATGTTAGTTAGCGCCACAAAGGAAACAGCATTCCCGGCGGAGAGTTGAATGCCGCGCAATTTATAAGCAAGATGAGTATTATCCAGtt      >  1:2645150/1‑150 (MQ=255)
                                                                                                                                                 cGCATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATGTTAGTTAGCGCCACAAAGGAAACAGCATTCCCGGCGGAGATTTGAATGCCGCGCAATTTATAAGCAAGATGAGTATTATCCAGTTGAat  <  2:1683635/150‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                                 
TCCTTTTTCATTCACCGCATGAACATTTGCAAGAGACAGCAGTGTTTCTTTTTTCGCCATAAAAACACCACGAACGTCTTTGCGCATGTCAAAGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGCAT‑‑AAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATGTTAGTTAGCGCCACAAAGGAAACAGCATTCCCGGCGGAGATTTGAATGCCGCGCAATTTATAAGCAAGATGAGTATTATCCAGTTGAAT  >  CP009273/2108070‑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: