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.A99.0% 638.6 / 0.9 193coding (1159/1395 nt)wcaMcolanic acid biosynthesis protein
Reads supporting (aligned to +/- strand):  ref base . (2/0);  new base A (87/104);  total (89/104)
Fisher's exact test for biased strand distribution p-value = 2.11e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.95e-01
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.
Rejected as polymorphism: Polymorphic indel expands or contracts a homopolymer stretch.

Only 100 of 208 total aligned reads displayed.

CTTTTTCATTCACCGCATGAACATTTGCAAGAGACAGCAGTGTTTCTTTTTTCGCCATAAAAACACCACGAACGTCTTTGCGCATGTCAAAGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGCATAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATGTTAGTTAGCGCCACAAAGGAAACAGCATTCCCGGCGGAGATTTGAATGCCGCGCAATTTATAAGCAAGATGAGTATTATCCAGTTGA  >  CP009273/2108072‑2108361
                                                                                                                                                    |                                                                                                                                               
cTTTTTCATTCACCGCATGAACATTTGCAAGAGACAGCAGTGTTTCTTTTTTCGCCATAAAAACACCACGAACGTCTTTGCGCATGTCAAAGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGCAt                                                                                                                                                  >  1:788918/1‑148 (MQ=255)
cTTTTTCATTCACCGCATGAACATTTGCAAGAGACAGCAGTGTTTCTTTTTTCGCCATAAAAACACCACGAACGTCTTTGAGCATGTCAAAGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGCAt                                                                                                                                                  >  2:1211264/1‑148 (MQ=255)
     tcattcaCCGCATGAACATTTGCAAGAGACAGCAGTGTTTCTTTTTTCGCCATAAAAACACCACGAACGTCTTTGCGCATGTCAAAGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGCATaaaaa                                                                                                                                           >  1:2652378/1‑148 (MQ=255)
        ttcaCCGCATGAACATTTGCAAGAGACAGCAGTGTTTCTTTTTTCGCCATAAAAACACCACGAACGTCTTTGCGCATGTCAAAGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGCATaaaaaaaa                                                                                                                                        >  2:724136/1‑147 (MQ=255)
         tcaCCGCATGAACATTTGCAAGAGACAGCAGTGTTTCTTTTTTCGCCATAAAAACACCACGAACGTCTTTGCGCATGTCAAAGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGCATaaaaaaaat                                                                                                                                       <  2:1386119/148‑3 (MQ=255)
         tcaCCGCATGAACATTTGCAAGAGACAGCAGTGTTTCTTTTTTCGCCATAAAAACACCACGAACGTCTTTGCGCATGTCAAAGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGCATaaaaaaaat                                                                                                                                       <  2:3338985/148‑3 (MQ=255)
         tcaCCGCATGAACATTTGCAAGAGACAGCAGTGTTTCTTTTTTCGCCATAAAAACACCACGAACGTCTTTGCGCATGTCAAAGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGCATaaaaaaaat                                                                                                                                       <  2:60381/148‑3 (MQ=255)
          caCCGCATGAACATTTGCAAGAGACAGCAGTGTTTCTTTTTTCGCCATAAAAACACCACGAACGTCTTTGCGCATGTCAAAGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGCATaaaaaaaatg                                                                                                                                      <  2:554965/148‑4 (MQ=255)
            ccGCATGAACATTTGCAAGAGACAGCAGTGTTTCTTTTTTCGCCATAAAAACACCACGAACGTCTTTGCGCATGTCAAAGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGCATaaaaaaaatgtt                                                                                                                                    >  1:698089/1‑143 (MQ=255)
                aTGAACATTTGCAAGAGACAGCAGTGTTTCTTTTTTCGCCATAAAAACACCACGAACGTCTTTGCGCATGTCAAAGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGCATAAAAAAAATGTTGCgg                                                                                                                                 >  2:1176516/1‑148 (MQ=255)
                 tGAACATTTGCAAGAGACAGCAGTGTTTCTTTTTTCGCCATAAAAACACCACGAACGTCTTTGCGCATGTCAAAGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGCATAAAAAAAATGTTGCGGt                                                                                                                                <  1:1870312/148‑1 (MQ=255)
                 tGAACATTTGCAAGAGACAGCAGTGTTTCTTTTTTCGCCATAAAAACACCACGAACGTCTTTGCGCATGTCAAAGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGCATAAAAAAAATGTTGCGGt                                                                                                                                >  1:2016089/1‑148 (MQ=255)
                 tGAACATTTGCAAGAGACAGCAGTGTTTCTTTTTTCGCCATAAAAACACCACGAACGTCTTTGCGCATGTCAAAGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGCATAAAAAAAATGTTGCGGt                                                                                                                                <  1:500562/148‑1 (MQ=255)
                    aCATTTGCAAGAGACAGCAGTGTTTCTTTTTTCGCCATAAAAACACCACGAACGTCTTTGCGCATGTCAAAGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGCATAAAAAAAATGTTGCGGTttg                                                                                                                             <  2:832640/148‑1 (MQ=255)
                        ttGCAAGAGACAGCAGTGTTTCTTTTTTCGCCATAAAAACACCACGAACGTCTTTGCGCATGTCAAAGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGCATAAAAAAAATGTTGCGGTTtgttgt                                                                                                                         >  1:1818762/1‑148 (MQ=255)
                        ttGCAAGAGACAGCAGTGTTTCTTTTTTCGCCATAAAAACACCACGAACGTCTTTGCGCATGTCAAAGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGCATAAAAAAAATGTTGCGGTTtgttgt                                                                                                                         <  2:8077/148‑1 (MQ=255)
                         tGCAAGAGACAGCAGTGTTTCTTTTTTCGCCATAAAAACACCACGAACGTCTTTGCGCATGTCAAAGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGCATAAAAAAAATGTTGCGGTTTGTtgtg                                                                                                                        <  2:2373917/148‑1 (MQ=255)
                         tGCAAGAGACAGCAGTGTTTCTTTTTTCGCCATAAAAACACCACGAACGTCTTTGCGCATGTCAAAGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGCATAAAAAAAATGTTGCGGTTTGTtgtg                                                                                                                        <  1:936448/148‑1 (MQ=255)
                            aaGAGACAGCAGTGTTTCTTTTTTCGCCATAAAAACACCACGAACGTCTTTGCGCATGTCAAAGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGCATAAAAAAAATGTTGCGGTTTGTTGTGTaa                                                                                                                     <  2:672957/148‑1 (MQ=255)
                            aaGAGACAGCAGTGTTTCTTTTTTCGCCATAAAAACACCACGAACGTCTTTGCGCATGTCAAAGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGCATAAAAAAAATGTTGCGGTTTGTTGTGTaa                                                                                                                     <  2:1892497/148‑1 (MQ=255)
                            aaGAGACAGCAGTGTTTCTTTTTTCGCCATAAAAACACCACGAACGTCTTTGCGCATGTCAAAGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGCATAAAAAAAATGTTGCGGTTTGTTGTGTaa                                                                                                                     <  2:772917/148‑1 (MQ=255)
                            aaGAGACAGCAGTGTTTCATTTTTCGCCATAAAAACACCACGAACGTCTTTCCGCATGTCAAACTTCATGCTCAATGCGGGTCCAACAGAGGATTCCTGCATCAAATTGATATTTCGCATAAAAAAAATGTTGCGGTTTGTTGTGTaa                                                                                                                     >  2:900068/1‑148 (MQ=255)
                              gagaCAGCAGTGTTTCTTTTTTCGCCATAAAAACACCACGAACGTCTTTGCGCATGTCAAAGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGCATAAAAAAAATGTTGCGGTTTGTTGTGTACCt                                                                                                                   <  1:284436/148‑1 (MQ=255)
                               agaCAGCAGTGTTTCTTTTTTCGCCATAAAAACACCACGAACGTCTTTGCGCATGTCAAAGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGCATAAAAAAAATGTTGCGGTTTGTTGTGTAACTc                                                                                                                  <  2:2683421/148‑1 (MQ=255)
                                gaCAGCAGTGTTTCTTTTTTCGCCATAAAAACACCACGAACGTCTTTGCGCATGTCAAAGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGCATAAAAAAAATGTTGCGGTTTGTTGTGTAACTcc                                                                                                                 <  2:1027212/148‑1 (MQ=255)
                                  cagcagTGTTTCTTTTTTCGCCATAAAAACACCACGAACGTCTTTGCGCATGTCAAAGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGCATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAg                                                                                                               <  1:1452488/148‑1 (MQ=255)
                                        tgtTTCTTTTTTCGCCATAAAAACACCACGAACGTCTTTGCGCATGTCAAAGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGCATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGacgc                                                                                                         >  2:1124576/1‑148 (MQ=255)
                                        tgtTTCTTTTTTCGCCATAAAAACACCACGAACGTCTTTGCGCATGTCAAAGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGCATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGacgc                                                                                                         >  2:2416078/1‑148 (MQ=255)
                                            tCTTTTTTCGCCATAAAAACACCACGAACGTCTTTGCGCATGTCAAAGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGCATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGacgcacgc                                                                                                     <  2:3442997/148‑1 (MQ=255)
                                             cTTTTTTCGCCATAAAAACACCACGAACGTCTTTGCGCATGTCAAAGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGCATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCt                                                                                                    >  1:671356/1‑148 (MQ=255)
                                                 tttCGCCATAAAAACACCACGAACGTCTTTGCGCATGTCAAAGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGCATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCAt                                                                                                >  1:2250313/1‑148 (MQ=255)
                                                 tttCGCCATAAAAACACCACGAACGTCTTTGCGCATGTCAAAGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGCATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCAt                                                                                                <  1:1505806/148‑1 (MQ=255)
                                                  ttCGCCATAAAAACACCACGAACGTCTTTGCGCATGTCAAAGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGCATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCAtc                                                                                               <  1:3119272/148‑1 (MQ=255)
                                                       cATAAAAACACCACGAACGTCTTTGCGCATGTCAAAGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGCATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAAt                                                                                          <  1:746164/148‑1 (MQ=255)
                                                       cATAAAAACACCACGAACGTCTTTGCGCATGTCAAAGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGCATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAAt                                                                                          >  1:3478046/1‑148 (MQ=255)
                                                        aTAAAAACACCACGAACGTCTTTGCGCATGTCAAAGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGCATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATg                                                                                         <  2:317166/148‑1 (MQ=255)
                                                        aTAAAAACACCACGAACGTCTTTGCGCATGTCAAAGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGCATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATg                                                                                         <  2:1628925/148‑1 (MQ=255)
                                                        aTAAAAACACCACGAACGTCTTTGCGCATGTCAAAGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGCATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATg                                                                                         >  1:1257677/1‑148 (MQ=255)
                                                        aTAAAAACACCACGAACGTCTTTGCGCATGTCAAAGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGCATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATg                                                                                         <  1:2915930/148‑1 (MQ=255)
                                                        aTAAAAACACCACGAACGTCTTTGCGCATGTCAAAGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGCATAAGAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATg                                                                                         >  1:508996/1‑148 (MQ=255)
                                                         tAAAAACACCACGAACGTCTTTGCGCATGTCAAAGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGCATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATgt                                                                                        <  2:2494296/148‑1 (MQ=255)
                                                         tAAAAACACCACGAACGTCTTTGCGCATGTCAAAGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGCATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATgt                                                                                        <  2:3253497/148‑1 (MQ=255)
                                                          aaaaaCACCACGAACGTCTTTGCGCATGTCAAAGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGCATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATgtt                                                                                       <  1:3035945/148‑1 (MQ=255)
                                                              acacCACGAACGTCTTTGCGCATGTCAAAGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGCATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATgttagtt                                                                                   <  1:3414335/148‑1 (MQ=255)
                                                               caccacGAACGTCTTTGCGCATGTCAAAGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGCATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATgttagtta                                                                                  >  1:3173975/1‑148 (MQ=255)
                                                               caccacGAACGTCTTTGCGCATGTCAAAGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGCATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATgttagtta                                                                                  >  1:1629185/1‑148 (MQ=255)
                                                                  cacGAACGTCTTTGCGCATGTCAAAGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGCATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATGTTAGTTAgcg                                                                               <  2:1839023/148‑1 (MQ=255)
                                                                      aaCGTCTTTGCGCATGTCAAAGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGCATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATGTTAGTTAGCGCcac                                                                           >  2:2089067/1‑148 (MQ=255)
                                                                       aCGTCTTTGCGCATGTCAAAGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGCATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATGTTAGTTAGCGCcaca                                                                          <  2:3308234/148‑1 (MQ=255)
                                                                           cTTTGCGCATGTCAAAGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGCATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATGTTAGTTAGCGCCACAAAgg                                                                      >  2:496802/1‑148 (MQ=255)
                                                                           cTTTGCGCATGTCAAAGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGCATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATGTTAGTTAGCGCCACAAAgg                                                                      >  1:1053037/1‑148 (MQ=255)
                                                                                  cATGTCAAAGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGCATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATGTTAGTTAGCGCCACAAAGGAAACAGc                                                               >  2:3467844/1‑148 (MQ=255)
                                                                                  cATGTCAAAGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGCATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATGTTAGTTAGCGCCACAAAGGAAACAGc                                                               <  2:3059256/148‑1 (MQ=255)
                                                                                  cATGTCAAAGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGCATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATGTTAGTTAGCGCCACAAAGGAAACAGc                                                               <  2:728252/148‑1 (MQ=255)
                                                                                  cATGTCAAAGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGCATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATGTTAGTTAGCGCCACAAAGGAAACAGc                                                               <  1:3313094/148‑1 (MQ=255)
                                                                                   aTGTCAAAGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGCATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATGTTAGTTAGCGCCACAAAGGAAACAGCa                                                              >  1:2519555/1‑148 (MQ=255)
                                                                                     gTCAAAGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGCATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATGTTAGTTAGCGCCACAAAGGAAACAGCAtt                                                            >  1:461785/1‑148 (MQ=255)
                                                                                      tCAAAGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGCATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATGTTAGTTAGCGCCACAAAGGAAACAGCATTc                                                           <  1:841865/148‑1 (MQ=255)
                                                                                       cAAAGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGCATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATGTTAGTTAGCGCCACAAAGGAAACAGCATTcc                                                          >  2:2773997/1‑148 (MQ=255)
                                                                                            ttCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGCATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATGTTAGTTAGCGCCACAAAGGAAACAGCATTCCcggcg                                                     >  2:1948355/1‑148 (MQ=255)
                                                                                             tCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGCATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATGTTAGTTAGCGCCACAAAGGAAACAGCATTCCcggcgg                                                    >  2:1991015/1‑148 (MQ=255)
                                                                                              cATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGCATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCTACGCACGCTTCATCTCAATGTTAGTTAGCGCCACAAAGGAAACAGCATTCCCGGCGga                                                   >  2:3444459/1‑148 (MQ=255)
                                                                                               aTGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCTCATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATGTTAGTTAGCGCCACAAAGGAAACAGCATTCCCGGCGgag                                                  <  1:179720/148‑1 (MQ=255)
                                                                                                tGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGCATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATGTTAGTTAGCGCCACAAAGGAAACAGCATTCCCGGCGgaga                                                 >  1:1396667/1‑148 (MQ=255)
                                                                                                tGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGCATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATGTTAGTTAGCGCCACAAAGGAAACAGCATTCCCGGCGgaga                                                 <  2:1159263/148‑1 (MQ=255)
                                                                                                 gCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGCATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATGTTAGTTAGCGCCACAAAGGAAACAGCATTCCCGGCGGAGAt                                                >  2:352379/1‑148 (MQ=255)
                                                                                                  cTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGCATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATGTTAGTTAGCGCCACAAAGGAAACAGCATTCCCGGCGGAGAtt                                               <  2:1818762/148‑1 (MQ=255)
                                                                                                    cAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGCATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATGTTAGTTAGCGCCACAAAGGAAACAGCATTCCCGGCGGAGATTTg                                             >  2:3229909/1‑148 (MQ=255)
                                                                                                      aTGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGCATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATGTTAGTTAGCGCCACAAAGGAAACAGCATTCCCGGCGGAGATTTGaa                                           >  1:1291187/1‑148 (MQ=255)
                                                                                                        gCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGCATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATGTTAGTTAGCGCCACAAAGGAAACAGCATTCCCGGCGGAGATTTGAATg                                         <  1:1532424/148‑1 (MQ=255)
                                                                                                         cGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGCATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATGTTAGTTAGCGCCACAAAGGAAACAGCATTCCCGGCGGAGATTTGAATGc                                        <  1:2151233/148‑1 (MQ=255)
                                                                                                              ccAACTGAGGATTCCTGCATCACATTGATATTTCGCATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATGTTAGTTAGCGCCACAAAGGAAACAGCATTCCCGGCGGAGATTTGAATGCCgcgc                                   >  1:470053/1‑148 (MQ=255)
                                                                                                                aaCTGAGGATTCCTGCATCACATTGATATTTCGCATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATGTTAGTTAGCGCCACAAAGGAAACAGCATTCCCGGCGGAGATTTGAATGCCGCGCaa                                 >  2:2133315/1‑148 (MQ=255)
                                                                                                                  cTGAGGATTCCTGCATCACATTGATATTTCGCATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATGTTAGTTAGCGCCACAAAGGAAACAGCATTCCCGGCGGAGATTTGAATGCCGCGCAAtt                               >  1:2427239/1‑148 (MQ=255)
                                                                                                                    gAGGATTCCTGCATCACATTGATATTTCGCATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATGTTAGTTAGCGCCACAAAGGAAACAGCATTCCCGGCGGAGATTTGAATGCCGCGCAATTta                             <  1:317593/148‑1 (MQ=255)
                                                                                                                     aGGATTCCTGCATCACATTGATATTTCGCATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATGTTAGTTAGCGCCACAAAGGAAACAGCATTCCCGGCGGAGATTTGAATGCCGCGCAATTtat                            >  2:3027934/1‑148 (MQ=255)
                                                                                                                     aGGAGTCCTGCATCACCTTGATATTTCGCATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAAGGTTAGTTAGCGCCACAAAGGAAACAGCATTCCCGGCGGAGATTTGAATGCCGCGCAATTtat                            <  1:1914698/148‑1 (MQ=255)
                                                                                                                      ggATTCCTGCATCACATTGATATTTCGCATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATGTTAGTTAGCGCCACAAAGGAAACAGCATTCCCGGCGGAGATTTGAATGCCGCGCAATTtata                           <  1:1176516/148‑1 (MQ=255)
                                                                                                                       gATTCCTGCATCACATTGATATTTCGCATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATGTTAGTTAGCGCCACAAAGGAAACAGCATTCCCGGCGGAGATTTGAATGCCGCGCAATTTATaa                          <  2:3293818/148‑1 (MQ=255)
                                                                                                                        aTTCCTGCATCACATTGATATTTCGCATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATGTTAGTTAGCGCCACAAAGGAAACAGCATTCCCGGCGGAGATTTGAATGCCGCGCAATTTATAAg                         <  2:1799479/148‑1 (MQ=255)
                                                                                                                        aTTCCTGCATCACATTGATATTTCGCATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATGTTAGTTAGCGCCACAAAGGAAACAGCATTCCCGGCGGAGATTTGAATGCCGCGCAATTTATAAg                         <  1:3211595/148‑1 (MQ=255)
                                                                                                                        aTTCCTGCATCACATTGATATTTCGCATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATGTTAGTTAGCGCCACAAAGGAAACAGCATTCCCGGCGGAGATTTGAATGCCGCGCAATTTATAAg                         <  1:3211487/148‑1 (MQ=255)
                                                                                                                         ttCCTGCATCACATTGATATTTCGCATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATGTTAGTTAGCGCCACAAAGGAAACAGCATTCCCGGCGGAGATTTGAATGCCGCGCAATTTATAAGc                        >  1:3297823/1‑148 (MQ=255)
                                                                                                                           ccTGCATCACATTGATATTTCGCATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATGTTAGTTAGCGCCACAAAGGAAACAGCATTCCCGGCGGAGATTTGAATGCCGCGCAATTTATAAGCaa                      <  2:1256024/148‑1 (MQ=255)
                                                                                                                           ccTGCATCACATTGATATTTCGCATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATGTTAGTTAGCGCCACAAAGGAAACAGCATTCCCGGCGGAGATTTGAATGCCGCGCAATTTATAAGCaa                      <  1:1255814/148‑1 (MQ=255)
                                                                                                                             tGCATCACATTGATATTTCGCATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATGTTAGTTAGCGCCACAAAGGAAACAGCATTCCCGGCGGAGATTTGAATGCCGCGCAATTTATAAGCAAGa                    <  1:3206139/148‑1 (MQ=255)
                                                                                                                               cATCACATTGATATTTCGCATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATGTTAGTTAGCGCCACAAAGGAAACAGCATTCCCGGCGGAGATTTGAATGCCGCGCAATTTATAAGCAAGATg                  <  2:1365796/148‑1 (MQ=255)
                                                                                                                               cATCACAGTGATATTTCGCATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATGTTAGTTAGCGCCACAAAGGAAACAGCATTCCCGGCGGAGATTTGAATGCCGCGCAATTTATAAGCAAGATg                  <  2:2908369/148‑1 (MQ=255)
                                                                                                                                aTCACATTGATATTTCGCATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATGTTAGTTAGCGCCACAAAGGAAACAGCATTCCCGGCGGAGATTTGAATGCCGCGCAATTTATAAGCAAGATGa                 >  2:2977390/1‑148 (MQ=255)
                                                                                                                                 tCACATTGATATTTCGCATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATGTTAGTTAGCGCCACAAAGGAAACAGCATTCCCGGCGGAGATTTGAATGCCGCGCAATTTATAAGCAAGATGAg                <  1:1002788/148‑1 (MQ=255)
                                                                                                                                 tCACATTGATATTTCGCATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATGTTAGTTAGCGCCACAAAGGAAACAGCATTCCCGGCGGAGATTTGAATGCCGCGCAATTTATAAGCAAGATGAg                <  1:2757316/148‑1 (MQ=255)
                                                                                                                                    caTTGATATTTCGCATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATGTTAGTTAGCGCCACAAAGGAAACAGCATTCCCGGCGGAGATTTGAATGCCGCGCAATTTATAAGCAAGATGAGtat             <  1:1948360/148‑1 (MQ=255)
                                                                                                                                      ttGATATTTCGCATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATGTTAGTTAGCGCCACAAAGGAAACAGCATTCCCGGCGGAGATTTGAATGCCGCGCAATTTATAAGCAAGATGAGtatta           <  1:550197/148‑1 (MQ=255)
                                                                                                                                      ttGATATTTCGCATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATGTTAGTTAGCGCCACAAAGGAAACAGCATTCCCGGCGGAGATTTGAATGCCGCGCAATTTATAAGCAAGATGAGtatta           <  1:2315313/148‑1 (MQ=255)
                                                                                                                                      ttGATATTTCGCATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATGTTAGTTAGCGCCACAAAGGAAACAGCATTCCCGGCGGAGATTTGAATGCCGCGCAATTTATAAGCAAGATGAGtatta           >  2:2583205/1‑148 (MQ=255)
                                                                                                                                         atatTTCGCATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATGTTAGTTAGCGCCACAAAGGAAACAGCATTCCCGGCGGAGATTTGAATGCCGCGCAATTTATAAGCAAGATGAGTATTATcc        >  2:1188358/1‑148 (MQ=255)
                                                                                                                                         atatTTCGCATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATGTTAGTTAGCGCCACAAAGGAAACAGCATTCCCGGCGGAGATTTGAATGCCGCGCAATTTATAAGCAAGATGAGTATTATcc        >  2:2938555/1‑148 (MQ=255)
                                                                                                                                           atTTCGCATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATGTTAGTTAGCGCCACAAAGGAAACAGCATTCCCGGCGGAGATTTGAATGCCGCGCAATTTATAAGCAAGATGAGTATTATCCAg      >  1:1760779/1‑148 (MQ=255)
                                                                                                                                               cGCATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATGTTAGTTAGCGCCACAAAGGAAACAGCATTCCCGGCGGAGATTTGAATGCCGCGCAATTTATAAGCAAGATGAGTATTATCCAGTTGa  <  1:2416078/148‑1 (MQ=255)
                                                                                                                                               cGCATAAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATGTTAGTTAGCGCCACAAAGGAAACAGCATTCCCGGCGGAGATTTGAATGCCGCGCAATTTATAAGCAAGATGAGTATTATCCAGTTGa  >  1:1152505/1‑148 (MQ=255)
                                                                                                                                                    |                                                                                                                                               
CTTTTTCATTCACCGCATGAACATTTGCAAGAGACAGCAGTGTTTCTTTTTTCGCCATAAAAACACCACGAACGTCTTTGCGCATGTCAAAGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGCATAAAAAAATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATGTTAGTTAGCGCCACAAAGGAAACAGCATTCCCGGCGGAGATTTGAATGCCGCGCAATTTATAAGCAAGATGAGTATTATCCAGTTGA  >  CP009273/2108072‑2108361

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: