Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,362,134 C→A 8.9% A12A (GCG→GCT pspF ← psp operon transcriptional activator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,362,1340CA 8.9% 874.6 / 41.5 243A12A (GCG→GCTpspFpsp operon transcriptional activator
Reads supporting (aligned to +/- strand):  ref base C (80/141);  new base A (6/16);  total (86/157)
Fisher's exact test for biased strand distribution p-value = 4.88e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 245 total aligned reads displayed.

GAAATAAACGGCCCTTGCCAACGGGAGGAGAGATAATGCAGGCGGCTGGCAATCAGCTCTTTACCGGTGCCGCGTTCGCCGATGATGAGCACCGGTTTGTCCAGCGGTGCGAGATGCGAAACCTGTTCCAGCACTTCGAGAAAGCTGTTCGCCTCACCAAGTAAATTATCTTTGTATTCTGCCATGATGAAATTCGCCACTTGTTAGTGTAATTCGCTAACTCATCCTGGCATGTTGCTGTTGATTCTTCAATCAGATCTTTATAAATCAAAAAGATAAAAAATTGGCACGCAAATT  >  CP009273/1361985‑1362281
                                                                                                                                                     |                                                                                                                                                   
gAAATAAACGGCCCTTGCCAACGGGAGGAGAGATAATGCAGGCGGCTGGCAATCAGCTCTTTACCGGTGCCGCGTTCGCCGATGATGAGCACCGGTTTGTCCAGCGGTGCGAGATGCGAAACCTGTTCCAGCACTTCGAGAAAGCTGTTc                                                                                                                                                     <  1:3885278/150‑1 (MQ=255)
gAAATAAACGGCCCTTGCCAACGGGAGGAGAGATAATGCAGGCGGCTGGCAATCAGCTCTTTACCGGTGCCGCGTTCGCCGATGATGAGCACCGGTTTGTCCAGCGGTGCGAGATGCGAAACCTGTTCCAGCACTTCGAGAAAGCTGTTc                                                                                                                                                     <  2:1763524/150‑1 (MQ=255)
  aaTAAACGGCCCTTGCCAACGGGAGGAGAGATAATGCAGGCGGCTGGCAATCAGCTCTTTACCGGTGCCGCGTTCGCCGATGATGAGCACCGGTTTGTCCAGCGGTGCGAGATGCGAAACCTGTTCCAGCACTTCGAGAAAGCTGTTCGc                                                                                                                                                   <  2:2867103/150‑1 (MQ=255)
   aTAAACGGCCCTTGCCAACGGGAGGAGAGATAATGCAGGCGGCTGGCAATCAGCTCTTTACCGGTGCCGCGTTCGCCGATGATGAGCACCGGTTTGTCCAGCGGTGCGAGATGCGAAACCTGTTCCAGCACTTCGAGAAAGCTGTTCGcc                                                                                                                                                  >  1:1799059/1‑150 (MQ=255)
     aaaCGGCCCTTGCCAACGGGAGGAGAGATAATGCAGGCGGCTGGCAATCAGCTCTTTACCGGTGCCGCGTTCGCCGATGATGAGCACCGGTTTGTCCAGCGGTGCGAGATGCGAAACCTGTTCCAGCACTTCGAGAAAGCTGTTCGCCTc                                                                                                                                                <  1:103609/150‑1 (MQ=255)
     aaaCGGCCCTTGCCAACGGGAGGAGAGATAATGCAGGCGGCTGGCAATCAGCTCTTTACCGGTGCCGCGTTCGCCGATGATGAGCACCGGTTTGTCCAGCGGTGCGAGATGCGAAACCTGTTCCAGCACTTCGAGAAAGCTGTTCGCCTc                                                                                                                                                >  2:1655176/1‑150 (MQ=255)
     aaaCGGCCCTTGCCAACGGGAGGAGAGATAATGCAGGCGGCTGGCAATCAGCTCTTTACCGGTGCCGCGTTCGCCGATGATGAGCACCGGTTTGTCCAGCGGTGCGAGATGCGAAACCTGTTCCAGCACTTCGAGAAAGCTGTTCGCCTc                                                                                                                                                >  2:3273328/1‑150 (MQ=255)
       aCGGCCCTTGCCAACGGGAGGAGAGATAATGCAGGCGGCTGGCAATCAGCTCTTTACCGGTGCCGCGTTCGCCGATGATGAGCACCGGTTTGTCCAGCGGTGCGAGATGCGAAACCTGTTCCAGCACTTCGAGAAAGCTGTTCGCCTCAc                                                                                                                                              <  2:1607232/150‑1 (MQ=255)
           cccTTGCCAACGGGAGGAGAGATAATGCAGGCGGCTGGCAATCAGCTCTTTACCGGTGCCGCGTTCGCCGATGATGAGCACCGGTTTGTCCAGCGGTGCGAGATGCGAAACCTGTTCCAGCACTTCGAGAAAGCTGTTCGCCTCACCAAg                                                                                                                                          <  1:888044/150‑1 (MQ=255)
           cccTTGCCAACGGGAGGAGAGATAATGCAGGCGGCTGGCAATCAGCTCTTTACCGGTGCCGCGTTCGCCGATGATGAGCACCGGTTTGTCCAGCGGTGCGAGATGCGAAACCTGTTCCAGCACTTCGAGAAAGCTGTTCGCCTCACCAAg                                                                                                                                          >  1:2075312/1‑150 (MQ=255)
             cTTGCCAACGGGAGGAGAGATAATGCAGGCGGCTGGCTATCAGCTCTTTACCGGTGCCGCGTTCGCCGATGATGAGCACCGGTTTGTCCAGCGGTGCGAGATGCGAAACCTGTTCCAGCACTTCGAGAAAGATGTTCGCATCACCAAGTa                                                                                                                                        <  2:832122/150‑1 (MQ=255)
             cTTGCCAACGGGAGGAGAGATAATGCAGGCGGCTGGCAATCAGCTCTTTACCGGTGCCGCGTTCGCCGATGATGAGCACCGGTTTGTCCAGCGGTGCGAGATGCGAAACCTGTTCCAGCACTTCGAGAAAGCTGTTCGCCTCACCAAGTa                                                                                                                                        >  1:2189929/1‑150 (MQ=255)
             cTTGCCAACGGGAGGAGAGATAATGCAGGCGGCTGGCAATCAGCTCTTTACCGGTGCCGCGTTCGCCGATGATGAGCACCGGTTTGTCCAGCGGTGCGAGATGCGAAACCTGTTCCAGCACTTCGAGAAAGCTGTTCGCCTCACCAAGTa                                                                                                                                        <  2:2313973/150‑1 (MQ=255)
              ttGCCAACGGGAGGAGAGATAATGCAGGCGGCTGGCAATCAGCTCTTTACCGGTGCCGCGTTCGCCGATGATGAGCACCGGTTTGTCCAGCGGTGCGAGATGCGAAACCTGTTCCAGCACTTCGAGAAAGCTGTTCGCCTCACCAAGTaa                                                                                                                                       <  1:487882/150‑1 (MQ=255)
              ttGCCAACGGGAGGAGAGATAATGCAGGCGGCTGGCAATCAGCTCTTTACCGGTGCCGCGTTCGCCGATGATGAGCACCGGTTTGTCCAGCGGTGCGAGATGCGAAACCTGTTCCAGCACTTCGAGAAAGCTGTTCGCCTCACCAAGTaa                                                                                                                                       <  1:2824156/150‑1 (MQ=255)
              ttGCCAACGGGAGGAGAGATAATGCAGGCGGCTGGCAATCAGCTCTTTACCGGTGCCGCGTTCGCCGATGATGAGCACCGGTTTGTCCAGCGGTGCGAGATGCGAAACCTGTTCCAGCACTTCGAGAAAGCTGTTCGCCTCACCAAGTaa                                                                                                                                       <  1:2772463/150‑1 (MQ=255)
                    aCGGGAGGAGAGATAATGCAGGCGGCTGGCAATCAGCTCTTTACCGGTGCCGCGTTCGCCGATGATGAGCACCGGTTTGTCCAGCGGTGCGAGATGCGAAACCTGTTCCAGCACTTCGAGAAAGCTGTTCGCCTCACCAAGTAAATTATc                                                                                                                                 >  1:2043224/1‑150 (MQ=255)
                    aCGGGAGGAGAGATAATGCAGGCGGCTGGCAATCAGCTCTTTACCGGTGCCGCGTTCGCCGATGATGAGCACCGGTTTGTCCAGCGGTGCGAGATGCGAAACCTGTTCCAGCACTTCGAGAAAGCTGTTAGCCTCACCAAGTAAATTATc                                                                                                                                 <  2:2363838/150‑1 (MQ=255)
                    aCGGGAGGAGAGATAATGCAGGCGGCTGGCAATCAGCTCTTTACCGGTGCCGCGTTCGCCGATGATGAGCACCGGTTTGTCCAGCGGTGCGAGATGCGAAACCTGTTCCAGCACTTCGAGAAAGCTGTTAGCCTCACCAAGTAAATTATc                                                                                                                                 <  2:2363732/150‑1 (MQ=255)
                     cGGGAGGAGAGATAATGCAGGCGGCTGGCAATCAGCTCTTTACCGGTGCCGCGTTCGCCGATGATGAGCACCGGTTTGTCCAGCGGTGCGAGATGCGAAACCTGTTCCAGCACTTCGAGAAAGCTGTTAGCCTCACCAAGTAAATTATCt                                                                                                                                >  1:2574729/1‑150 (MQ=255)
                       ggaggaGAGATAATGCAGGCGGCTGGCAATCAGCTCTTTACCGGTGCCGCGTTCGCCGATGATGAGCACCGGTTTGTCCAGCGGTGCGAGATGCGAAACCTGTTCCAGCACTTCGAGAAAGCTGTTCGCCTCACCAAGTAAATTATCttt                                                                                                                              <  2:1011286/150‑1 (MQ=255)
                       ggaggaGAGATAATGCAGGCGGCTGGCAATCAGCTCTTTACCGGTGCCGCGTTCGCCGATGATGAGCACCGGTTTGTCCAGCGGTGCGAGATGCGAAACCTGTTCCAGCACTTCGAGAAAGCTGTTCGCCTCACCAAGTAAATTATCttt                                                                                                                              <  2:1436545/150‑1 (MQ=255)
                        gaggagAGATAATGCAGGCGGCTGTCAATCAGCTCTTTACCGGTGCCGCGTTCGCCGATGATGAGCACCGGTTTGTCCAGCGGTGCGAGATGCGAAACCTGTTCCAGCACTTCGAGAAAGCTGTTCGCCTCACCAAGTAAATTATCTTTg                                                                                                                             <  2:3339500/150‑1 (MQ=255)
                             gagaTAATGCAGGCGGCTGGCAATCAGCTCTTTACCGGTGCCGCGTTCGCCGATGATGAGCACCGGTTTGTCCAGCGGTGCGAGATGCGAAACCTGTTCCAGCACTTCGAGAAAGCTGTTCGCCTCACCAAGTAAATTATCTTTGTATTc                                                                                                                        >  1:1013728/1‑150 (MQ=255)
                                  aaTGCAGGCGGCTGGCAATCAGCTCTTTACCGGTGCCGCGTTCGCCGATGATGAGCACCGGTTTGTCCAGCGGTGCGAGATGCGAAACCTGTTCCAGCACTTCGAGAAAGCTGTTCGCCTCACCAAGTAAATTATCTTTGTATTCTGCCa                                                                                                                   <  2:345604/150‑1 (MQ=255)
                                  aaTGCAGGCGGCTGGCAATCAGCTCTTTACCGGTGCCGCGTTCGCCGATGATGAGCACCGGTTTGTCCAGCGGTGCGAGATGCGAAACCTGTTCCAGCACTTCGAGAAAGCTGTTCGCCTCACCAAGTAAATTATCTTTGTATTCTGCCa                                                                                                                   <  1:344522/150‑1 (MQ=255)
                                   aTGCAGGCGGCTGGCAATCAGCTCTTTACCGGTGCCGCGTTCGCCGATGATGAGCACCGGTTTGTCCAGCGGTGCGAGATGCGAAACCTGTTCCAGCACTTCGAGAAAGCTGTTCGCCTCACCAAGTAAATTATCTTTGTATTCTGCCat                                                                                                                  <  1:2215321/150‑1 (MQ=255)
                                      cAGGCGGCTGGCAATCAGCTCTTTACCGGTGCCGCGTTCGCCGATGATGAGCACCGGTTTGTCCAGCGGTGCGAGATGCGAAACCTGTTCCAGCACTTCGAGAAAGCTGTTCGCCTCACCAAGTAAATTATCTTTGTATTCTGCCatgat                                                                                                               >  2:2133065/1‑150 (MQ=255)
                                      cAGGCGGCTGGCAATCAGCTCTTTACCGGTGCCGCGTTCGCCGATGATGAGCACCGGTTTGTCCAGCGGTGCGAGATGCGAAACCTGTTCCAGCACTTCGAGAAAGCTGTTCGCCTCACCAAGTAAATTATCTTTGTATTCTGCCatgat                                                                                                               <  1:809517/150‑1 (MQ=255)
                                       aGGCGGCTGGCAATCAGCTCTTTACCGGTGCCGCGTTAGCCGATGATGAGCACCGGTTTGTCCAGCGGTGCGAGATGCGAAACCTTTTCCAGCACTTCGAGAAAGCTGTTCGCCTCACCAAGTAAATTATCTTTGTATTCTGCCatgatg                                                                                                              >  1:4437760/1‑150 (MQ=255)
                                         gcggcTGGCAATCAGCTCTTTACCGGTGCCGCGTTCGCCGATGATGAGCACCGGTTTGTCCAGCGGTGCGAGATGCGAAACCTGTTCCAGCACTTCGAGAAAGCTGTTCGCCTCACCAAGTAAATTATCTTTGTATTCTGCCATGATGaa                                                                                                            <  1:163172/150‑1 (MQ=255)
                                         gcggcTGGCAATCAGCTCTTTACCGGTGCCGCGTTCGCCGATGATGAGCACCGGTTTGTCCAGCGGTGCGAGATGCGAAACCTGTTCCAGCACTTCGAGAAAGCTGTTCGCCTCACCAAGTAAATTATCTTTGTATTCTGCCATGATGaa                                                                                                            <  2:2629002/150‑1 (MQ=255)
                                              tGGCAATCAGCTCTTTACCGGTGCCGCGTTCGCCGATGATGAGCACCGGTTTGTCCAGCGGTGCGAGATGCGAAACCTGTTCCAGCACTTCGAGAAAGCTGTTCGCCTCACCAAGTAAATTATCTTTGTATTCTGCCATGATGAAATTCg                                                                                                       >  2:4532171/1‑150 (MQ=255)
                                                    tccgCTCTTGCCCGGTGCCGCGTTCGCCGCTGATGAGCACCGGTTTGTCCAGCGGTGCGAGATGCGAAACCTGTACAAGCACTTCGAGACAGCTGTTCGCCTCACCAAGTAAATTATCTTTGTATTCTGCCATGATGAAATTCGCCACtt                                                                                                 <  2:880603/147‑1 (MQ=255)
                                                    tCAGCTCTTTACCGGTGCCGCGTTCGCCGATGATGAGCACCGGTTTGTCCAGCGGTGCGAGATGCGAAACCTGTTCCAGCACTTCGAGAAAGCTGTTCGCCTCACCAAGTAAATTATCTTTGTATTCTGCCATGATGAAATTCGCCACtt                                                                                                 <  1:3605778/150‑1 (MQ=255)
                                                     cAGCTCTTTACCGGTGCCGCGTTCGCCGATGATGAGCACCGGTTTGTCCAGCGGTGCGAGATGCGAAACCTGTTCCAGCACTTCGAGAAAGCTGTTCGCCTCACCAAGTAAATTATCTTTGTATTCTGCCATGATGAAATTCGCCACTTg                                                                                                <  1:1954004/150‑1 (MQ=255)
                                                     cAGCTCTTTACCGGTGCCGCGTTCGCCGATGATGAGCACCGGTTTGTCCAGCGGTGCGAGATGCGAAACCTGTTCCAGCACTTCGAGAAAGCTGTTCGCCTCACCAAGTAAATTATCTTTGTATTCTGCCATGATGAAATTCGCCACTTg                                                                                                <  1:4308811/150‑1 (MQ=255)
                                                      aGCTCTTTACCGGTGCCGCGTTCGCCGATGATGAGCACCGGTTTGTCCAGCGGTGCGAGATGCGAAACCTGTTCCAGCACTTCGAGAAAGCTGTTCGCCTCACCAAGTAAATTATCTTTGTATTCTGCCATGATGAAATTCGCCACTTGt                                                                                               <  2:3685840/150‑1 (MQ=255)
                                                       gCTCTTTACCGGTGCCGCGTTCGCCGATGATGAGCACCGGTTTGTCCAGCGGTGCGAGATGCGAAACCTGTTCCAGCACTTCGAGAAAGCTGTTCGCCTCACCAAGTAAATTATCTTTGTATTCTGCCATGATGAAATTCGCCACTTGtt                                                                                              <  1:4317665/150‑1 (MQ=255)
                                                        ctctTTACCGGTGCCGCGTTCGCCGATGATGAGCACCGGTTTGTCCAGCGGTGCGAGATGCGAAACCTGTTCCAGCACTTCGAGAAAGCTGTTCGCCTCACCAAGTAAATTATCTTTGTATTCTGCCATGATGAAATTCGCCACTTGTTa                                                                                             <  1:3622325/150‑1 (MQ=255)
                                                         tctTTACCGGTGCCGCGTTCGCCGATGATGAGCACCGGTTTGTCCAGCGGTGCGAGATGCGAAACCTGTTCCAGCACTTCGAGAAAGCTGTTCGCCTCACCAAGTAAATTATCTTTGTATTCTGCCATGATGAAATTCGCCACTTGTTAg                                                                                            >  1:3290265/1‑150 (MQ=255)
                                                          ctTTACCGGTGCCGCGTTCGCCGATGATGAGCACCGGTTTGTCCAGCGGTGCGAGATGCGAAACCTGTTCCAGCACTTCGAGAAAGCTGTTCGCCTCACCAAGTAAATTATCTTTGTATTCTGCCATGATGAAATTCGCCACTTGTTAgt                                                                                           <  1:3272118/150‑1 (MQ=255)
                                                          ctTTACCGGTGCCGCGTTCGCCGATGATGAGCACCGGTTTGTCCAGCGGTGCGAGATGCGAAACCTGTTCCAGCACTTCGAGAAAGCTGTTCGCCTCACCAAGTAAATTATCTTTGTATTCTGCCATGATGAAATTCGCCACTTGTTAgt                                                                                           <  1:3274725/150‑1 (MQ=255)
                                                           tttACCGGTGCCGCGTTCGCCGATGATGAGCACCGGTTTGTCCAGCGGTGCGAGATGCGAAACCTGTTCCAGCACTTCGAGAAAGCTGTTCGCCTCACCAAGTAAATTATCTTTGTATTCTGCCATGATGAAATTCGCCACTTGTTAgtg                                                                                          <  1:2496821/150‑1 (MQ=255)
                                                            ttACCGGTGCCGCGTTCGCCGATGATGAGCACCTGTTTGTCCAGCGGTGCGAGATGCGAAACCTGTTCCAGCACTTCGAGAAAGCTGTTCGCCTCACCAAGTAAATTATCTTTGTATTCTGCCATGATGAAATTCGCCACTTGTTAgtgt                                                                                         <  2:4222775/150‑1 (MQ=255)
                                                             tACCGGTGCCGCGTTCGCCGATGATGAGCACCGGTTTGTCCAGCGGTGCGAGATGCGAAACCTGTTCCAGCACTTCGAGAAAGCTGTTCGCCTCACCAAGTAAATTATCTTTGTATTCTGCCATGATGAAATTCGCCACTTGTTAGTGTa                                                                                        >  1:722517/1‑150 (MQ=255)
                                                             tACCGGTGACGCGTTCGCCGATGATGAGCACCGGTTTGTCCAGCGGTGCGAGATGCGAAACCTGTTCCAGCACTTCGAGAAAGCTGTTCGCCTCACCAAGTAAATTATCTTTGTATTCTGCCATGATGAAATTCGCCACTTGTTAGTGTa                                                                                        <  2:1801336/150‑1 (MQ=255)
                                                               ccGGTGCCGCGTTCGCCGATGATGAGCACCGGTTTGTCCAGCGGTGCGAGATGCGAAACCTGTTCCAGCACTTCGAGAAAGCTGTTCGCCTCACCAAGTAAATTATCTTTGTATTCTGCCATGATGAAATTCGCCACTTGTTAGTGTAAt                                                                                      <  1:2068494/150‑1 (MQ=255)
                                                                        cgTTCGCCGATGATGAGCACCGGTTTGTCCAGCGGTGCGAGATGCGAAACCTGTTCCAGCACTTCGAGAAAGCTGTTCGCCTCACCAAGTAAATTATCTTTGTATTCTGCCATGATGAAATTCGCCACTTGTTAGTGTAATTCGCTAACt                                                                             <  2:1963807/150‑1 (MQ=255)
                                                                        cgTTCGCCGATGATGAGCACCGGTTTGTCCAGCGGTGCGAGATGCGAAACCTGTTCCAGCACTTCGAGAAAGCTGTTCGCCTCACCAAGTAAATTATCTTTGTATTCTGCCATGATGAAATTCGCCACTTGTTAGTGTAATTCGCTAACt                                                                             >  1:2798190/1‑150 (MQ=255)
                                                                             gCCGATGATGAGCACCGGTTTGTCCAGCGGTGCGAGATGCGAAACCTGTTCCAGCACTTCGAGAAAGCTGTTCGCCTCACCAAGTAAATTATCTTTGTATTCTGCCATGATGAAATTCGCCACTTGTTAGTGTAATTCGCTAACTCATcc                                                                        >  2:3821764/1‑150 (MQ=255)
                                                                              ccGATGATGAGCACCGGTTTGTCCAGCGGTGCGAGATGCGAAACCTGTTCCAGCACTTCGAGAAAGCTGTTCGCCTCACCAAGTAAATTATCTTTGTATTCTGCCATGATGAAATTCGCCACTTGTTAGTGTAATTCGCTAACTCATCCt                                                                       >  2:3106939/1‑150 (MQ=255)
                                                                               cGATGATGAGCACCGGTTTGTCCAGCGGTGCGAGATGCGAAACCTGTTCCAGCACTTCGAGAAAGCTGGTCGCCTCACCAAGTAAATTATCTTTGTATTCTGCCATGATGACATTCGCCACTTGTTAGGGTAATTCGCTCACTCATCCTg                                                                      <  2:3967390/150‑1 (MQ=255)
                                                                                gatgatGAGCACCGGTTTGTCCAGCGGTGCGAGATGCGAAACCTGTTCCAGCACTTCGAGAAAGCTGTTCGCCTCACCAAGTAAATTATCTTTGTATTCTGCCATGATGAAATTCGCCACTTGTTAGTGTAATTCGCTAACTCATCCTgg                                                                     >  1:4250975/1‑150 (MQ=255)
                                                                                  tgatgaGCACCGGTTTGTCCAGCGGTGCGAGATGCGAAACCTGTTCCAGCACTTCGAGAAAACTGTTCGCCTCACCAAGTAAATTATCTTTGTATTCTGCCATGATGAAATTCGCCACTTGTTAGTGTAATTCGCTAACTCATCCTGGCa                                                                   >  2:2746930/1‑150 (MQ=255)
                                                                                    atgaGCACCGGTTTGTCCAGCGGTGCGAGATGCGAAACCTGTTCCAGCACTTCGAGAAAGCTGTTCGCCTCACCAAGTAAATTATCTTTGTATTCTGCCATGATGAAATTCGCCACTTGTTAGTGTAATTCGCTAACTCATCCTGGCATg                                                                 <  2:657642/150‑1 (MQ=255)
                                                                                     tgaGCACCGGTTTGTCCAGCGGTGCGAGATGCGAAACCTGTTCCAGCACTTCGAGAAAGCTGTTCGCCTCACCAAGTAAATTATCTTTGTATTCTGCCATGATGAAATTCGCCACTTGTTAGTGTAATTCGCTAACTCATCCTGGCATGt                                                                >  1:3713667/1‑150 (MQ=255)
                                                                                     tgaGCACCGGTTTGTCCAGCGGTGCGAGATGCGAAACCTGTTCCAGCACTTCGAGAAAGCTGTTCGCCTCACCAAGTAAATTATCTTTGTATTCTGCCATGATGAAATTCGCCACTTGTTAGTGTAATTCGCTAACTCATCCTGGCATGt                                                                >  1:3946232/1‑150 (MQ=255)
                                                                                      gaGCACCGGTTTGTCCAGCGGTGCGAGATGCGAAACCTGTTCCAGCACTTCGAGAAAGCTGTTCGCCTCACCAAGTAAATTATCTTTGTATTCTGCCATGATGAAATTCGCCACTTGTTAGTGTAATTCGCTAACTCATCCTGGCATGtt                                                               <  1:610431/150‑1 (MQ=255)
                                                                                      gaGCACCGGTTTGTCCAGCGGTGCGAGATGCGAAACCTGTTCCAGCACTTCGAGAAAGCTGTTCGCCTCACCAAGTAAATTATCTTTGTATTCTGCCATGATGAAATTCGCCACTTGTTAGTGTAATTCGCTAACTCATCCTGGCATGtt                                                               <  1:586371/150‑1 (MQ=255)
                                                                                         cACCGGTTTGTCCAGCGGTGCGAGATGCGAAACCTGTTCCAGCACTTCGAGAAAGCTGTTCGCCTCACCAAGTAAATTATCTTTGTATTCTGCCATGATGAAATTCGCCACTTGTTAGTGTAATTCGCTAACTCATCCTGGCATGTTGCt                                                            >  1:544673/1‑150 (MQ=255)
                                                                                          aCCGGTTTGTCCAGCGGTGCGAGATGCGAAACCTGTTCCAGCACTTCGAGAAAGCTGTTCGCCTCACCAAGTAAATTATCTTTGTATTCTGCCATGATGAAATTCGCCACTTGTTAGTGTAATTCGCTAACTCATCCTGGCATGTTGCTg                                                           <  1:2729288/150‑1 (MQ=255)
                                                                                          aCCGGTTTGTCCAGCGGTGCGAGATGCGAAACCTGTTCCAGCACTTCGAGAAAGCTGTTAGCCTCACCAAGTAAATTATCTTTGTATTCTGCCATGATGAAATTCGCCACTTGTTAGTGTAATTCGCTAACTCATCCTGGCATGTTGCTg                                                           >  1:310315/1‑150 (MQ=255)
                                                                                              gTTTGTCCAGCGGTGCGAGATGCGAAACCTGTTCCAGCACTTCGAGAAAGCTGTTCGCCTCACCAAGTAAATTATCTTTGTATTCTGCCATGATGAAATTCGCCACTTGTTAGTGTAATTCGCTAACTCATCCTGGCATGTTGCTGTTGa                                                       >  1:2254461/1‑150 (MQ=255)
                                                                                                 tGTCCAGCGGTGCGAGATGCGAAACCTGTTCCAGCACTTCGAGAAAGCTGTTCGCCTCACCAAGTAAATTATCTTTGTATTCTGCCATGATGAAATTCGCCACTTGTTAGTGTAATTCGCTAACTCATCCTGGCATGTTGCTGTTGAttc                                                    <  2:3083873/150‑1 (MQ=255)
                                                                                                  gTCCAGCGGTGCGAGATGCGAAACCTGTTCCAGCACTTCGAGAAAGCTGTTCGCCTCACCAAGTAAATTATCTTTGTATTCTGCCATGATGAAATTCGCCACTTGTTAGTGTAATTCGCTAACTCATCCTGGCATGTTGCTGTTGAttct                                                   <  1:2384125/150‑1 (MQ=255)
                                                                                                     cAGCGGTGCGAGATGCGAAACCTGTTCCAGCACTTCGAGAAAGCTGTTCGCCTCACCAAGTAAATTATCTTTGTATTCTGCCATGATGAAATTCGCCACTTGTTAGTGTAATTCGCTAACTCATCCTGGCATGTTGCTGTTGATTCTTCa                                                >  2:2562295/1‑150 (MQ=255)
                                                                                                        cGGTGCGAGATGCGAAACCTGTTCCAGCACTTCGAGAAAGCTGTTCGCCTCACCAAGTAAATTATCTTTGTATTCTGCCATGATGAAATTCGCCACTTGTTAGTGTAATTCGCTAACTCATCCTGGCATGTTGCTGTTGATTCTTCAATc                                             <  2:3409133/150‑1 (MQ=255)
                                                                                                         ggTGCGAGATGCGAAACCTGTTCCAGCACTTCGAGAAAGCTGTTCGCCTCACCAAGTAAATTATCTTTGTATTCTGCCATGATGAAATTCGCCACTTGTTAGTGTAATTCGCTAACTCATCCTGGCATGTTGCTGTTGATTCTTCAATCa                                            <  1:173820/150‑1 (MQ=255)
                                                                                                              gagaTGCGAAACCTGTTCCAGCACTTCGAGAAAGCTGTTCGCCTCACCAAGTAAATTATCTTTGTATTCTGCCATGATGAAATTCGCCACTTGTTAGTGTAATTCGCTAACTCATCCTGGCATGTTGCTGTTGATTCTTCAATCAGATCt                                       >  1:426305/1‑150 (MQ=255)
                                                                                                               agaTGCGAAACCTGTTCCAGCACTTCGAGAAAGCTGTTCGCCTCACCAAGTAAATTATCTTTGTATTCTGCCATGATGAAATTCGCCACTTGTTAGTGTAATTCGCTAACTCATCCTGGCATGTTGCTGTTGATTCTTCAATCAGATCtt                                      >  2:1878884/1‑150 (MQ=255)
                                                                                                                gaTGCGAAACCTGTTCCAGCACTTCGAGAAAGCTGTTCGCCTCACCAAGTAAATTATCTTTGTATTCTGCCATGATGAAATTCGCCACTTGTTAGTGTAATTCGCTAACTCATCCTGGCATGTTGCTGTTGATTCTTCAATCAGATCttt                                     <  1:3782698/150‑1 (MQ=255)
                                                                                                                 aTGCGAAACCTGTTCCAGCACTTCGAGAAAGCTGTTCGCCTCACCAAGTAAATTATCTTTGTATTCTGCCATGATGAAATTCGCCACTTGTTAGTGTAATTCGCTAACTCATCCTGGCATGTTGCTGTTGATTCTTCAATCAGATCTTta                                    <  1:2990958/150‑1 (MQ=255)
                                                                                                                  tGCGAAACCTGTTCCAGCACTTCGAGAAAGCTGTTAGCCTCACCAAGTAAATTATCTTTGTATTCTGCCATGATGAAATTCGCCACTTGTTAGTGTAATTCGCTAACTCATCCTGGCATGTTGCTGTTGATTCTTCAATCAGATCTTtat                                   >  2:3789470/1‑150 (MQ=255)
                                                                                                                   gCGAAACCTGTTCCAGCACTTCGAGAAAGCTGTTCGCCTCACCAAGTAAATTATCTTTGTATTCTGCCATGATGAAATTCGCCACTTGTTAGTGTAATTCGCTAACTCATCCTGGCATGTTGCTGTTGATTCTTCAATCAGATCTTtata                                  >  2:2934919/1‑150 (MQ=255)
                                                                                                                     gAAACCTGTTCCAGCACTTCGAGAAAGCTGTTCGCCTCACCAAGTAAATTATCTTTGTATTCTGCCATGATGAAATTCGCCACTTGTTAGTGTAATTCGCTAACTCATCCTGGCATGTTGCTGTTGATTCTTCAATCAGATCTTTATaaa                                >  2:2272519/1‑150 (MQ=255)
                                                                                                                     gAAACCTGTTCCAGCACTTCGAGAAAGCTGTTCGCCTCACCAAGTAAATTATCTTTGTATTCTGCCATGATGAAATTCGCCACTTGTTAGTGTAATTCGCTAACTCATCCTGGCATGTTGCTGTTGATTCTTCAATCAGATCTTTATaaa                                >  1:381893/1‑150 (MQ=255)
                                                                                                                     gAAACCTGTTCCAGCACTTCGAGAAAGCTGTTCGCCTCACCAAGTAAATTATCTTTGTATTCTGCCATGATGAAATTCGCCACTTGTTAGTGTAATTCGCTAACTCATCCTGGCATGTTGCTGTTGATTCTTCAATCAGATCTTTATaaa                                >  1:2677275/1‑150 (MQ=255)
                                                                                                                      aaaCCTGTTCCAGCACTTCGAGAAAGCTGTTCGCCTCACCAAGTAAATTATCTTTGTATTCTGCCATGATGAAATTCGCCACTTGTTAGTGTAATTCGCTAACTCATCCTGGCATGTTGCTGTTGATTCTTCAATCAGATCTTTATAAAt                               <  1:1587106/150‑1 (MQ=255)
                                                                                                                       aaCCTGTTCCAGCACTTCGAGAAAGCTGTTCGCCTCACCAAGTAAATTATCTTTGTATTCTGCCATGATGAAATTCGCCACTTGTTAGTGTAATTCGCTAACTCATCCTGGCATGTTGCTGTTGATTCTTCAATCAGATCTTTATAAATc                              <  2:1066104/150‑1 (MQ=255)
                                                                                                                            gTTCCAGCACTTCGAGAAAGCTGTTCGCCTCACCAAGTAAATTATCTTTGTATTCTGCCATGATGAAATTCGCCACTTGTTAGTGTAATTCGCTAACTCATCCTGGCATGTTGCTGTTGATTCTTCAATCAGATCTTTATAAATCaaaaa                         <  1:1769673/150‑1 (MQ=255)
                                                                                                                            gTTCCAGCACTTCGAGAAAGCTGTTAGCCTCACCAAGTAAATTATCTTTGTATTCTGCCATGATGAAATTCGCCACTTGTTAGTGTAATTCGCTAACTCATCCTGGCATGTTGCTGTTGATTCTTCAATCAGATCTTTATAACTCaaaaa                         <  2:2782666/150‑1 (MQ=255)
                                                                                                                              tCCAGCACTTCGAGAAAGCTGTTAGCCTCACCTATTAAATTATCTTTGTTTTCTGCCATGTTGAAATTCGCCACTTGTTAGTGTAATTCGCTAACTCATCCTGGCATGTTGCTGTTGATTCTTCAATCAGATCTTTATAAATCAAAAAGa                       <  1:495819/150‑1 (MQ=255)
                                                                                                                               ccAGCACTTCGAGAAAGCTGTTCGCCTCACCAAGTAAATTATCTTTGTATTCTGCCATGATGAAATTCGCCACTTGTTAGTGTAATTCGCTAACTCATCCTGGCATGTTGCTGTTGATTCTTCAATCAGATCTTTATAAATCAAAAAGAt                      <  2:426329/150‑1 (MQ=255)
                                                                                                                               ccAGCACTTCGAGAAAGCTGTTCGCCTCACCAAGTAAATTATCTTTGTATTCTGCCATGATGAAATTCGCCACTTGTTAGTGTAATTCGCTAACTCATCCTGGCATGTTGCTGTTGATTCTTCAATCAGATCTTTATAAATCAAAAAGAt                      <  2:426292/150‑1 (MQ=255)
                                                                                                                                cAGCACTTCGAGAAAGCTGTTCGCCTCACCAAGTAAATTATCTTTGTATTCTGCCATGATGAAATTCGCCACTTGTTAGTGTAATTCGCTAACTCATCCTGGCATGTTGCTGTTGATTCTTCAATCAGATCTTTATAAATCAAAAAGATa                     >  1:2174745/1‑150 (MQ=255)
                                                                                                                                cAGCACTTCGAGAAAGCTGTTCGCCTCACCAAGTAAATTATCTTTGTATTCTGCCATGATGAAATTCGCCACTTGTTAGTGTAATTCGCTAACTCATCCTGGCATGTTGCTGTTGATTCTTCAATCAGATCTTTATAAATCAAAAAGATa                     >  1:3809443/1‑150 (MQ=255)
                                                                                                                                cAGCACTTCGAGAAAGCTGTTCGCCTCACCAAGTAAATTATCTTTGTATTCTGCCATGATGAAATTCGCCACTTGTTAGTGTAATTCGCTAACTCATCCTGGCATGTTGCTGTTGATTCTTCAATCAGATCTTTATAAATCAAAAAGATa                     <  1:279652/150‑1 (MQ=255)
                                                                                                                                 aGCACTTCGAGAAAGCTGTTCGCCTCACCAAGTAAATTATCTTTGTATTCTGCCATGATGAAATTCGCCACTTGTTAGTGTAATTCGCTAACTCATCCTGGCATGTTGCTGTTGATTCTTCAATCAGATCTTTATAAATCAAAAAGATaa                    >  1:867683/1‑150 (MQ=255)
                                                                                                                                    aCTTCGAGAAAGCTGTTCGCCTCACCAAGTAAATTATCTTTGTATTCTGCCATGATGAAATTCGCCACTTGTTAGTGTAATTCGCTAACTCATCCTGGCATGTTGCTGTTGATTCTTCAATCAGATCTTTATAAATCAAAAAGATaaaaa                 >  1:22917/1‑150 (MQ=255)
                                                                                                                                    aCTTCGAGAAAGCTGTTCGCCTCACCAAGTAAATTATCTTTGTATTCTGCCATGATGAAATTCGCCACTTGTTAGTGTAATTCGCTAACTCATCCTGGCATGTTGCTGTTGATTCTTCAATCAGATCTTTATAAATCAAAAAGATaaaaa                 >  1:2326307/1‑150 (MQ=255)
                                                                                                                                     cTTCGAGAAAGCTGTTCGCCTCACCAAGTAAATTATCTTTGTATTCTGCCATGATGAAATTCGCCACTTGTTAGTGTAATTCGCTAACTCATCCTGGCATGTTGCTGTTGATTCTTCAATCAGATCTTTATAAATCAAAAAGATaaaaaa                <  2:1476489/150‑1 (MQ=255)
                                                                                                                                         gagaAAGCTGTTCGCCTCACCAAGTAAATTATCTTTGTATTCTGCCATGATGAAATTCGCCACTTGTTAGTGTAATTCGCTAACTCATCCTGGCATGTTGCTGTTGATTCTTCAATCAGATCTTTATAAATCAAAAAGATAAAAAATTgg            <  1:1848461/150‑1 (MQ=255)
                                                                                                                                         gagaAAGCTGTTCGCCTCACCAAGTAAATTATCTTTGTATTCTGCCATGATGAAATTCGCCACTTGTTAGTGTAATTCGCTAACTCATCCTGGCATGTTGCTGTTGATTCTTCAATCAGATCTTTATAAATCAAAAAGATAAAAAATTgg            >  2:1704299/1‑150 (MQ=255)
                                                                                                                                         gagaAAGCTGTTCGCCTCACCAAGTAAATTATCTTTGTATTCTGCCATGATGAAATTCGCCACTTGTTAGTGTAATTCGCTAACTCATCCTGGCATGTTGCTGTTGATTCTTCAATCAGATCTTTATAAATCAAAAAGATAAAAAATTgg            >  2:1272415/1‑150 (MQ=255)
                                                                                                                                         gagaAAGCTGTTCGCCTCACCAAGTAAATTATCTTTGTATTCTGCCATGATGAAATTCGCCACTTGTTAGTGTAATTCGCTAACTCATCCTGGCATGTTGCTGTTGATTCTTCAATCAGATCTTTATAAATCAAAAAGATAAAAAATTgg            <  2:2075312/150‑1 (MQ=255)
                                                                                                                                          agaAAGCTGTTCGCCTCACCAAGTAAATTATCTTTGTATTCTGCCATGATGAAATTCGCCACTTGTTAGTGTAATTCGCTAACTCATCCTGGCATGTTGCTGTTGATTCTTCAATCAGATCTTTATAAATCAAAAAGATAAAAAATTGGc           <  2:2208594/150‑1 (MQ=255)
                                                                                                                                               gCTGTTCGCCTCACCAAGTAAATTATCTTTGTATTCTGCCATGATGAAATTCGCCACTTGTTAGTGTAATTCGCTAACTCATCCTGGCATGTTGCTGTTGATTCTTCAATCAGATCTTTATAAATCAAAAAGATAAAAAATTGGCACGCa      <  1:3832624/150‑1 (MQ=255)
                                                                                                                                                 tGTTCGCCTCACCAAGTAAATTATCTTTGTATTCTGCCATGATGAAATTCGCCACTTGTTAGTGTAATTCGCTAACTCATCCTGGCATGTTGCTGTTGATTCTTCAATCAGATCTTTATAAATCAAAAAGATAAAAAATTGGCACGCaaa    <  2:4310240/150‑1 (MQ=255)
                                                                                                                                                  gTTAGCCTCACCAAGTAAATTATCTTTGTATTCTGCCATGATGAAATTCGCCACTTGTTAGTGTAATTCGCTAACTCATCCTGGCATGTTGCTGTTGATTCTTCAATCAGATCTTTATAAATCAAAAAGATAAAAAATTGGCACGCAAAt   <  2:1668727/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
GAAATAAACGGCCCTTGCCAACGGGAGGAGAGATAATGCAGGCGGCTGGCAATCAGCTCTTTACCGGTGCCGCGTTCGCCGATGATGAGCACCGGTTTGTCCAGCGGTGCGAGATGCGAAACCTGTTCCAGCACTTCGAGAAAGCTGTTCGCCTCACCAAGTAAATTATCTTTGTATTCTGCCATGATGAAATTCGCCACTTGTTAGTGTAATTCGCTAACTCATCCTGGCATGTTGCTGTTGATTCTTCAATCAGATCTTTATAAATCAAAAAGATAAAAAATTGGCACGCAAATT  >  CP009273/1361985‑1362281

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: