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.A90.0% 647.0 / 66.9 190coding (1159/1395 nt)wcaMcolanic acid biosynthesis protein
Reads supporting (aligned to +/- strand):  ref base . (11/8);  new base A (77/94);  total (88/102)
Fisher's exact test for biased strand distribution p-value = 3.37e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.71e-01
Rejected as polymorphism: Polymorphic indel expands or contracts a homopolymer stretch.

Only 100 of 196 total aligned reads displayed.

TTTTCATTCACCGCATGAACATTTGCAAGAGACAGCAGTGTTTCTTTTTTCGCCATAAAAACACCACGAACGTCTTTGCGCATGTCAAAGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGC‑‑ATAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATGTTAGTTAGCGCCACAAAGGAAACAGCATTCCCGGCGGAGATTTGAATGCCGCGCAATTTATAAGCAAGATGAGTATTATCCAGTTGAATA  >  CP009273/2108074‑2108364
                                                                                                                                                    |                                                                                                                                                 
  ttcattcaCCGCATGAACATTTGCAAGAGACAGCAGTGTTTCTTTTTTCGCCATAAAAACACCACGAACGTCTTTGCGCATGTCAAAGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGC‑‑ATaaaaaa                                                                                                                                             <  2:3268732/150‑1 (MQ=255)
   tcattcaCCGCATGAACATTTGCAAGAGACAGCAGTGTTTCTTTTTTCGCCATAAAAACACCACGAACGTCTTTGCGCATGTCAAAGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGC‑‑ATaaaaaaa                                                                                                                                            >  2:408543/1‑150 (MQ=255)
       tcaCCGCATGAACATTTGCAAGAGACAGCAGTGTTTCTTTTTTCGCCATAAAAACACCACGAACGTCTTTGCGCATGTCAAAGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGC‑‑ATaaaaaaaatgt                                                                                                                                        >  1:3550697/1‑146 (MQ=255)
         aCCGCATGAACATTTGCAAGAGACAGCAGTGTTTCTTTTTTCGCCATAAAAACACCACGAACGTCTTTGCGCATGTCAAAGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGC‑‑ATAAAAAAAATGTTg                                                                                                                                       <  2:2106508/150‑1 (MQ=255)
         aCCGCATGAACATTTGCAAGAGACAGCAGTGTTTCTTTTTTCGCCATAAAAACACCACGAACGTCTTTGCGCATGTCAAAGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGC‑‑ATAAAAAAAATGTTg                                                                                                                                       <  2:2032300/150‑1 (MQ=255)
            gCATGAACATTTGCAAGAGACAGCAGTGTTTCTTTTTTCGCCATAAAAACACCACGAACGTCTTTGCGCATGTCAAAGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGC‑‑ATAAAAAAAATGTTGCgg                                                                                                                                    >  1:3809402/1‑150 (MQ=255)
             cATGAACATTTGCAAGAGACAGCAGTGTTTCTTTTTTCGCCATAAAAACACCACGAACGTCTTTGCGCATGTCAAAGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGC‑‑ATAAAAAAAATGTTGCGGt                                                                                                                                   >  2:1329160/1‑150 (MQ=255)
              aTGAACATTTGCAAGAGACAGCAGTGTTTCTTTTTTCGCCATAAAAACACCACGAACGTCTTTGCGCATGTCAAAGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGC‑‑ATAAAAAAAATGTTGCGGtt                                                                                                                                  >  1:333929/1‑150 (MQ=255)
              aTGAACATTTGCAAGAGACAGCAGTGTTTCTTTTTTCGCCATAAAAACACCACGAACGTCTTTGCGCATGTCAAAGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGC‑‑ATAAAAAAAATGTTGCGGtt                                                                                                                                  <  1:3912195/150‑1 (MQ=255)
              aTGAACATTTGCAAGAGACAGCAGTGTTTCTTTTTTCGCCATAAAAACACCACGAACGTCTTTGCGCATGTCAAAGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGC‑‑ATAAAAAAAATGTTGCGGtt                                                                                                                                  >  1:585078/1‑150 (MQ=255)
               tGAACATTTGCAAGAGACAGCAGTGTTTCTTTTTTCGCCATAAAAACACCACGAACGTCTTTGCGAATGTCAAAGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGC‑‑ATAAAAAAAATGTTGCGGttt                                                                                                                                 <  1:1630936/150‑1 (MQ=255)
                 aaCATTTGCAAGAGACAGCAGTGTTTCTTTTTTCGCCATAAAAACACCACGAACGTCTTTGCGCATGTCAAAGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGC‑‑ATAAAAAAAATGTTGCGgtgtgt                                                                                                                               >  2:282695/1‑150 (MQ=255)
                  aCATTTGCAAGAGACAGCAGTGTTTCTTTTTTCGCCATAAAAACACCACGAACGTCTTTGCGCATGTCAAAGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGC‑‑ATAAAAAAAATGTTGCGGTttgtt                                                                                                                              <  1:1215928/150‑1 (MQ=255)
                   cATTTGCAAGAGACAGCAGTGTTTCTTTTTTCGCCATAAAAACACCACGAACGTCTTTGCGCATGTCAAAGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGC‑‑ATAAAAAAAATGTTGCGGTttgttg                                                                                                                             <  1:91085/150‑1 (MQ=255)
                    aTTTGCAAGAGACAGCAGTGTTTCTTTTTTCTCCATAAAAACACCACGAACGTCTTTGCGCATGTCAAAGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCAAATTGATATTTCGC‑‑ATAAAAAAAATGTTGCGGTTtgttgt                                                                                                                            >  2:4264072/1‑150 (MQ=255)
                    aTTTGCAAGAGACAGCAGTGTTTCTTTTTTCGCCATAAAAACACCACGAACGTCTTTGCGCATGTCAAAGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGC‑‑ATAAAAAAAATGTTGCGGTTtgttgt                                                                                                                            >  2:4049414/1‑150 (MQ=255)
                    aTTTGCAAGAGACAGCAGTGTTTCTTTTTTCGCCATAAAAACACCACGAACGTCTTTGCGCATGTCAAAGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGC‑‑ATAAAAAAATGTTGCGGTTTGTtgtg                                                                                                                           <  2:970243/150‑1 (MQ=255)
                     tttGCAAGAGACAGCAGTGTTTCTTTTTTCGCCATAAAAACACCACGAACGTCTTTGCGCATGTCAAAGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGC‑‑ATAAAAAAAATGTTGCGGTTTGTtgtg                                                                                                                           <  2:727634/150‑1 (MQ=255)
                     tttGCAAGAGACAGCAGTGTTTCTTTTTTCGCCATAAAAACACCACGAACGTCTTTGCGCATGTCAAAGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGC‑‑ATAAAAAAAATGTTGCGGTTTGTtgtg                                                                                                                           <  2:4147425/150‑1 (MQ=255)
                         cAAGAGACAGCAGTGTTTCTTTTTTCGCCATAAAAACACCACGAACGTCTTTGCGCATGTCAAAGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGC‑‑ATAAAAAAAATGTTGCGGTTTGTTGTGTAAc                                                                                                                       <  2:61967/150‑1 (MQ=255)
                         cAAGAGACAGCAGTGTTTCTTTTTTCGCCATAAAAACACCACGAACGTCTTTGCGCATGTCAAAGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGC‑‑ATAAAAAAATGTTGCGGTTTGTTGTGTAACt                                                                                                                      <  2:2385695/150‑1 (MQ=255)
                             agaCAGCAGTGTTTCTTTTTTCGCCATAAAAACACCACGAACGTCTTTGCGCATGTCAAAGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGC‑‑ATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCa                                                                                                                   <  1:156695/150‑1 (MQ=255)
                              gaCAGCAGTGTTTCTTTTTTCGCCATAAAAACACCACGAACGTCTTTGCGCATGTCAAAGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGC‑‑ATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAg                                                                                                                  >  2:3380188/1‑150 (MQ=255)
                                 agcagTGTTTCTTTTTTCGCCATAAAAACACCACGAACGTCTTTGCGCATGTCAAAGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGC‑‑ATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGa                                                                                                               >  1:2081116/1‑150 (MQ=255)
                                  gcagTGTTTCTTTTTTCGCCATAAAAACACCACGAACGTCTTTGCGCATGTCAAAGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGC‑‑ATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGac                                                                                                              >  1:3119475/1‑150 (MQ=255)
                                    agTGTTTCTTTTTTCGCCATAAAAACACCACGAACGTCTTTGCGCATGTCAAAGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGC‑‑ATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGacgc                                                                                                            >  1:3658413/1‑150 (MQ=255)
                                    agTGTTTCTTTTTTCGCCATAAAAACACCACGAACGTCTTTGCGCATGTCAAAGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGC‑‑ATAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGacgca                                                                                                           >  1:1902581/1‑150 (MQ=255)
                                     gtgtTTCTTTTTTCGCCATAAAAACACCACGAACGTCTTTGCGCATGTCAAAGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGC‑‑ATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGacgca                                                                                                           <  2:4181483/150‑1 (MQ=255)
                                           cTTTTTTCGCCATAAAAACACCACGAACGTCTTTGCGCATGTCAAAGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGC‑‑ATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTc                                                                                                     >  2:3199436/1‑150 (MQ=255)
                                           cTTTTTTCGCCATAAAAACACCACGAACGTCTTTGCGCATGTCAAAGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGC‑‑ATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTc                                                                                                     <  2:3560592/150‑1 (MQ=255)
                                           cTTTTTTCGCCATAAAAACACCACGAACGTCTTTGCGCATGTCAAAGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGC‑‑ATAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCa                                                                                                    >  2:2617660/1‑150 (MQ=255)
                                             tttttCGCCATAAAAACACCACGAACGTCTTTGCGCATGTCAAAGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGC‑‑ATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCAt                                                                                                   >  2:1274432/1‑150 (MQ=255)
                                             tttttCGCCATAAAAACACCACGAACGTCTTTGCGCATGTCAAAGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGC‑‑ATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCAt                                                                                                   >  1:733043/1‑150 (MQ=255)
                                                  cGCCATAAAAACACCACGAAAGTCTTTGCGCATGTCAAAGTTCATGCTCAATTCGGGTCCAACTGAGGATTCCTGAATCACATTGATATTTCGA‑‑ATAAAAAAAATGTTGCTGTTGGTTTTGTAACTCCAGCGACGCACGCTTCTTATCaa                                                                                              >  2:1004646/1‑150 (MQ=255)
                                                   gCCATAAAAACACCACGAACGTCTTTGCGCATGTCAAAGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGC‑‑ATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAAt                                                                                             <  2:3277115/150‑1 (MQ=255)
                                                   gCCATAAAAACACCACGAACGTCTTTGCGCATGTCAAAGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGC‑‑ATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAAt                                                                                             <  1:2940784/150‑1 (MQ=255)
                                                      aTAAAAACACCACGAACGTCTTTGCGCATGTCAAAGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGC‑‑ATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATgtt                                                                                          >  2:1290270/1‑150 (MQ=255)
                                                       tAAAAACACCACGAACGTCTTTGCGCATGTCAAAGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGC‑‑ATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATgtta                                                                                         <  2:271142/150‑1 (MQ=255)
                                                        aaaaaCACCACGAACGTCTTTGCGCATGTCAAAGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGC‑‑ATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATgttag                                                                                        <  1:205912/150‑1 (MQ=255)
                                                        aaaaaCACCACGAACGTCTTTGCGCATGTCAAAGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGC‑‑ATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATgttag                                                                                        >  2:1930510/1‑150 (MQ=255)
                                                         aaaaCACCACGGACCTCGTTGGGCCTTTAAAAAGTCCTTGTCAATGCGGGTTCAACTGTGGTTTCCGGGCTCACCATGATTTTTCGCAAAAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATgttag                                                                                        <  2:2376632/150‑1 (MQ=255)
                                                         aaaaCACCACGAACGTCTTTGCGCGTGTCAAAGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGC‑‑ATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATgttagt                                                                                       <  2:3758595/150‑1 (MQ=255)
                                                         aaaaCACCACGAACGTCTTTGCGCATGTCAAAGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGC‑‑ATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATgttagt                                                                                       >  2:4113906/1‑150 (MQ=255)
                                                            acacCACGAACGTCTTTGCGCATGTCAAAGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGC‑‑ATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATGttagttag                                                                                    >  1:2384341/1‑150 (MQ=255)
                                                                   gAACGTCTTTGCGCATGTCAAAGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGC‑‑ATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATGTTAGTTAGCGCcaca                                                                             >  1:502379/1‑150 (MQ=255)
                                                                    aaCGTCTTTGCGCATGTCAAAGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGC‑‑ATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATGTTAGTTAGCGCCACaa                                                                            <  1:26393/150‑1 (MQ=255)
                                                                    aaCGTCTTTGCGCATGTCAAAGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGC‑‑ATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATGTTAGTTAGCGCCACaa                                                                            <  1:27593/150‑1 (MQ=255)
                                                                     aCGTCTTTGCGCATGTCAAAGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGC‑‑ATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATGTTAGTTAGCGCCACaaa                                                                           >  1:217875/1‑150 (MQ=255)
                                                                         cTTTGCGCATGTCAAAGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGC‑‑ATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATGTTAGTTAGCGCCACAAAGGaa                                                                       <  1:2012148/150‑1 (MQ=255)
                                                                         cTTTGCGCATGTCAAAGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGC‑‑ATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATGTTAGTTAGCGCCACAAAGGaa                                                                       <  2:2271411/150‑1 (MQ=255)
                                                                         cTTTGCGCATGTCAAAGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGC‑‑ATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATGTTAGTTAGCGCCACAAAGGaa                                                                       >  1:3440839/1‑150 (MQ=255)
                                                                                cATGTCAAAGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGC‑‑ATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATGTTAGTTAGCGCCACAAAGGAAACAGCAt                                                                <  2:2935690/150‑1 (MQ=255)
                                                                                cATGTCAAAGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGC‑‑ATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATGTTAGTTAGCGCCACAAAGGAAACAGCAt                                                                <  2:2935073/150‑1 (MQ=255)
                                                                                cATGTCAAAGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGC‑‑ATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATGTTAGTTAGCGCCACAAAGGAAACAGCAt                                                                >  2:2330571/1‑150 (MQ=255)
                                                                                cATGTCAAAGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGC‑‑ATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATGTTAGTTAGCGCCACAAAGGAAACAGCAt                                                                >  2:1340533/1‑150 (MQ=255)
                                                                                cATGTCAAAGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGC‑‑ATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATGTTAGTTAGCGCCACAAAGGAAACAGCAt                                                                >  1:4066762/1‑150 (MQ=255)
                                                                                cATGTCAAAGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGC‑‑ATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATGTTAGTTAGCGCCACAAAGGAAAAAGCAt                                                                >  2:2577534/1‑150 (MQ=255)
                                                                                 aTGTCAAAGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGC‑‑ATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATGTTAGTTAGCGCCACAAAGGAAACAGCAtt                                                               >  2:742374/1‑150 (MQ=255)
                                                                                 aTGTCAAAGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGC‑‑ATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATGTTAGTTAGCGCCACAAAGGAAACAGCAtt                                                               >  2:742350/1‑150 (MQ=255)
                                                                                  tGTCAAAGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGC‑‑ATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATGTTAGTTAGCGCCACAAAGGAAACAGCATTc                                                              >  1:395358/1‑150 (MQ=255)
                                                                                     cAAAGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGC‑‑ATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATGTTAGTTAGCGCCACAAAGGAAACAGCATTCCcg                                                           <  1:1813552/150‑1 (MQ=255)
                                                                                      aaaGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGC‑‑ATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATGTTAGTTAGCGCCACAAAGGAAACAGCATTCCcgg                                                          >  2:1445182/1‑150 (MQ=255)
                                                                                           tCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGC‑‑ATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATGTTAGTTAGCGCCACAAAGGAAACAGCATTCCCGGCGgag                                                     <  2:518590/150‑1 (MQ=255)
                                                                                             aTGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGC‑‑ATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATGTTAGTTAGCGCCACAAAGGAAACAGCATTCCCGGCGGAGAt                                                   <  2:2170525/150‑1 (MQ=255)
                                                                                             aTGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGC‑‑ATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATGTTAGTTAGCGCCACAAAGGAAACAGCATTCCCGGCGGAGAt                                                   >  1:399007/1‑150 (MQ=255)
                                                                                             aTGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGC‑‑ATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATGTTAGTTAGCGCCACAAAGGAAACAGCATTCCCGGCGGAGAt                                                   <  2:3484999/150‑1 (MQ=255)
                                                                                             aTGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGC‑‑ATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATGTTAGTTAGCGCCACAAAGGAAACAGCATTCCCGGCGGAGAt                                                   <  2:3070085/150‑1 (MQ=255)
                                                                                               gCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGC‑‑ATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATGTTAGTTAGCGCCACAAAGGAAACAGCATTCCCGGCGGAGAttt                                                 >  1:3541115/1‑150 (MQ=255)
                                                                                               gCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGC‑‑ATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATGTTAGTTAGCGCCACAAAGGAAACAGCATTCCCGGCGGAGAttt                                                 <  1:2913039/150‑1 (MQ=255)
                                                                                               gCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGC‑‑ATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATGTTAGTTAGCGCCACAAAGGAAACAGCATTCCCGGCGGAGAttt                                                 <  1:1922409/150‑1 (MQ=255)
                                                                                                cTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGC‑‑ATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATGTTAGTTAGCGCCACAAAGGAAACAGCATTCCCGGCGGAGATTTg                                                <  2:579728/150‑1 (MQ=255)
                                                                                                cTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGC‑‑ATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATGTTAGTTAGCGCCACAAAGGAAACAGCATTCCCGGCGGAGATTTg                                                >  2:2854843/1‑150 (MQ=255)
                                                                                                 tCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGC‑‑ATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATGTTAGTTAGCGCCACAAAGGAAACAGCATTCCCGGCGGAGATTTGa                                               >  1:10664/1‑150 (MQ=255)
                                                                                                  cAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGC‑‑ATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATGTTAGTTAGCGCCACAAAGGAAACAGCATTCCCGGCGGAGATTTGaa                                              >  2:3251/1‑150 (MQ=255)
                                                                                                   aaTGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGC‑‑ATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATGTTAGTTAGCGCCACAAAGGAAACAGCATTCCCGGCGGAGATTTGAAt                                             >  1:320508/1‑150 (MQ=255)
                                                                                                   aaTGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGC‑‑ATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATGTTAGTTAGCGCCACAAAGGAAACAGCATTCCCGGCGGAGATTTGAAt                                             <  1:1290344/150‑1 (MQ=255)
                                                                                                       cGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGC‑‑ATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATGTTAGTTAGCGCCACAAAGGAAACAGCATTCCCGGCGGAGATTTGAATGCcg                                         <  1:690120/150‑1 (MQ=255)
                                                                                                           tCCAACTGAGGATTCCTGCATCACATTGATATTTCGC‑‑ATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATGTTAGTTAGCGCCACAAAGGAAACAGCATTCCCGGCGGAGATTTGAATGCCGCGCa                                     <  2:1062630/150‑1 (MQ=255)
                                                                                                            ccAACTGAGGATTCCTGCATCACATTGATATTTCGC‑‑ATAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATGTTAGTTAGCGCCACAAAGGAAACAGCATTCCCGGCGGAGATTTGAATGCCGCGCAAt                                   <  1:3685703/150‑1 (MQ=255)
                                                                                                             cAACTGAGGATTCCTGCATCACATTGATATTTCGC‑‑ATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCCATGTTAGTTAGCGCCACAAAGGAAACAGCATTCCCGGCGGAGATTTGAATGCCGCGCAAt                                   <  1:3840007/150‑1 (MQ=255)
                                                                                                               aCTGAGGATTCCTGCATCACATTGATATTTCGC‑‑ATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATGTTAGTTAGCGCCACAAAGGAAACAGCATTCCCGGCGGAGATTTGAATGCCGCGCAAttt                                 <  2:2372790/150‑1 (MQ=255)
                                                                                                               aCTGAGGATTCCTGCATCACATTGATATTTCGC‑‑ATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATGTTAGTTAGCGCCACAAAGGAAACAGCATTCCCGGCGGAGATTTGAATGCCGCGCAAttt                                 <  1:1338401/150‑1 (MQ=255)
                                                                                                                  gAGGATTCCTGCATCACATTGATATTTCGC‑‑ATAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATGTTAGTTAGCGCCACAAAGGAAACAGCATTCCCGGCGGAGATTTGAATGCCGCGCAATTTATaa                             >  2:1012916/1‑150 (MQ=255)
                                                                                                                   aGGATTCCTGCATCACATTGATATTTCGC‑‑ATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATGTTAGTTAGCGCCACAAAGGAAACAGCATTCCCGGCGGAGATTTGAATGCCGCGCAATTTATaa                             <  2:3996582/150‑1 (MQ=255)
                                                                                                                   aGGATTCCTGCATCACATTGATATGTCGC‑‑ATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATGTTAGTTAGCGCCACAAAGGAAACAGCATTCCCGGCGGAGATTTGAATGCCGCGCAATTTATaa                             <  2:4297022/150‑1 (MQ=255)
                                                                                                                      aTTCCTGCATCACATTGATATTTCGC‑‑ATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATGTTAGTTAGCGCCACAAAGGAAACAGCATTCCCGGCGGAGATTTGAATGCCGCGCAATTTATAAGCa                          <  1:4240098/150‑1 (MQ=255)
                                                                                                                      aTTCCTGCATCACATTGATATTTCGC‑‑ATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATGTTAGTTAGCGCCACAAAGGAAACAGCATTCCCGGCGGAGATTTGAATGCCGCGCAATTTATAAGCa                          >  2:1569267/1‑150 (MQ=255)
                                                                                                                        tCCTGCATCACATTGATATTTCGC‑‑ATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATGTTAGTTAGCGCCACAAAGGAAACAGCATTCCCGGCGGAGATTTGAATGCCGCGCAATTTATAAGCAAg                        <  2:3289372/150‑1 (MQ=255)
                                                                                                                         ccTGCATCACATTGATATTTCGC‑‑ATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATGTTAGTTAGCGCCACAAAGGAAACAGCATTCCCGGCGGAGATTTGAATGCCGCGCAATTTATAAGCAAGa                       >  1:3253553/1‑150 (MQ=255)
                                                                                                                            gCATCACATTGATATTTCGC‑‑ATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATGTTAGTTAGCGCCACAAAGGAAACAGCATTCCCGGCGGAGATTTGAATGCCGCGCAATTTATAAGCAAGATGa                    >  1:417098/1‑150 (MQ=255)
                                                                                                                            gCATCACATTGATATTTCGC‑‑ATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATGTTAGTTAGCGCCACAAAGGAAACAGCATTCCCGGCGGAGATTTGAATGCCGCGCAATTTATAAGCAAGATGa                    <  2:383738/150‑1 (MQ=255)
                                                                                                                             cATCACATTGATATTTCGC‑‑ATAAAAAAAATGTTGAGGTTTGTTGTGTAACTCCAGAGACTCACGCTTCATCTCAATTTTAGTTAGCTCCACAAAGGAAACAGCATTCCCTTCGGAGATTTGAATGCCGCGCAATTTATAAGCAAGATGAg                   >  1:3773811/1‑150 (MQ=255)
                                                                                                                                 acaTTGATATTTCGC‑‑ATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATGTTAGTTAGCGCCACAAAGGAAACAGCATTCCCGGCGGAGATTTGAATGCCGCGCAATTTATAAGCAAGATGAGtatt               <  1:1378100/150‑1 (MQ=255)
                                                                                                                                  caTTGATATTTCGC‑‑ATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATGTTAGTTAGCGCCACAAAGGAAACAGCATTCCCGGCGGAGATTTGAATGCCGCGCAATTTATAAGCAAGATGAGtatta              <  2:2210690/150‑1 (MQ=255)
                                                                                                                                   aTTGATATTTCGC‑‑ATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCCATGTTAGTTAGCGCCACAAAGGAAACAGCATTCCCGGCGGAGATTTGAATGCCGCGCAATTTATAAGCAAGATGAGtattat             <  1:1252973/150‑1 (MQ=255)
                                                                                                                                      gATATTTCGC‑‑ATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATGTTAGTTAGCGCCACAAAGGAAACAGCATTCCCGGCGGAGATTTGAATGCCGCGCAATTTATAAGCAAGATGAGTATTATCCa          <  2:911014/150‑1 (MQ=255)
                                                                                                                                      gATAGTGCGC‑‑ATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCCTCTCAATGTTAGTTAGCGCCACAAAGGAAACAGCATTCCCGGCGGAGATTTGAATGCCGCGCAATTTATAAGCAAGATGAGTATTATCCa          <  2:1237037/150‑1 (MQ=255)
                                                                                                                                        tatTTCGC‑‑ATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATGTTAGTTAGCGCCACAAAGGAAACAGCATTCCCGGCGGAGATTTGAATGCCGCGCAATTTATAAGCAAGATGAGTATTATCCAGt        >  1:3737927/1‑150 (MQ=255)
                                                                                                                                        tatTTCGC‑‑ATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATGTTAGTTAGCGCCACAAAGGAAACAGCATTCCCGGCGGAGATTTGAATGCCGCGCAATTTATAAGCAAGATGAGTATTATCCAGt        <  1:1719366/150‑1 (MQ=255)
                                                                                                                                         atTTCGC‑‑ATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATGTTAGTTAGCGCCACAAAGGAAACAGCATTCCCGGCGGAGATTTGAATGCCGCGCAATTTATAAGCAAGATGAGTATTATCCAGtt       >  1:1390414/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                 
TTTTCATTCACCGCATGAACATTTGCAAGAGACAGCAGTGTTTCTTTTTTCGCCATAAAAACACCACGAACGTCTTTGCGCATGTCAAAGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGC‑‑ATAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATGTTAGTTAGCGCCACAAAGGAAACAGCATTCCCGGCGGAGATTTGAATGCCGCGCAATTTATAAGCAAGATGAGTATTATCCAGTTGAATA  >  CP009273/2108074‑2108364

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: