Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,363,021 G→T 5.6% intergenic (+42/+161) yjeI → / ← yjeJ DUF4156 family lipoprotein/uncharacterized protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,363,0210GT 5.6% 256.6 / 2.8 77intergenic (+42/+161)yjeI/yjeJDUF4156 family lipoprotein/uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base G (34/38);  new base T (2/3);  total (36/41)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.77e-01

GTACTGCGACAGGTAAGCAAAGCAACTGGCTTTCCGGGCAACACGGAGAAGAAGGCGGTTCTATGCGCGGCGCAGCAAACGATCTGCGCAACCAGGCGGCTGCAATGGGCGGTAACGTGATTTATGGCATCAGTAGCCCGTCGCAGGGAATGTTGTCCAGTTTTGTCCCGACGGATAGCCAGATTATCGGTCAGGTTTATAAGTGCCCGAACTGATGTGATCCAGAGGCTGGGTGACAAACGTATAATTGCCTGATGCGCTTCGCTTATCAGGCCTACATAAGGCATTGCAACTGATTGAATGCTTGAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGACCAACGATCACGTTATCAGTCATTTCAGGCTGGATGGCTCCAGCCTCGCATTAGCAAAGATGAGATTATTTCGCCTGTGGTGCAGTTTTGCTGGTGGATTGAACACGTAAGTGATGCAGCGCGCGTAAGCATTGATCCTGGGTTAGCGGTTGGGCGTTATTCGCCGCTAC  >  CP009273/4362765‑4363286
                                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                                         
gTAGTGCGACAGGTAAGCAAAGGAACTGGCTTGCCGGGAAAAACGGAGAAGAAGGCGGTGCTATGGGCGGCGCAGCAAACGNTCTGCGCANCCAGGCGGCTGCCATGGGCGGTAACGTGATTTATGGCATCAGTAGCCCGTCGCAGGGAATGTTGTCCAGTTTTGTCCCGACGGATAGCCAGATTATCGGTCAGGTTTATAAGTGCCCGAACTGATGTGATCCAGAGGCTGGGTGACAACCGTATAATTGCCTGATGCGCTTCGCTTATCAGGc                                                                                                                                                                                                                                                          <  1:462557/274‑1 (MQ=255)
gTACTGCGACAGGTAAGCAAAGCAACTGGCTTTCCGGGCAACACGGAGAAGAAGGCGGTTCTATGCGCGGCGCAGCAATCGATCTGCGCAACCAGGCGGCTGCAATGGGCGGTAACGTGATGTATGGCATCAGTAGCCCGTCGCAGGGAATGTTGTCCAGTTTTGTCCCGACGGATAGCCAGATTATCGGTCAGGTTTATAAGTGCCCGAACTGATGTGATCCAGAGGCTGGGTGACAAACGTATAATTGCCTGATGCGCTTCGCTTATCAGGc                                                                                                                                                                                                                                                          <  1:690329/274‑1 (MQ=255)
    ggggagaGGCAAGAAAAGAAGGTGGCTTTCCGGGCAAAACGAGGAAGAAGGCGGTTCTATGCGCGGCGCAGCAAACGATCTGCGCAACCAGGCGGCTGCAATGGGCGGTAACGTGATTTATGGCATCAGTAGCCCGTCGCAGGGAATGTTGTCCAGTTTTGTCCCGACGGATAGCCAGATTATCGGTCAGGTTTATAAGTGCCCGAACTGATGTGATCCAGAGGCTGGGTGACAAACGTATAATTGCCTGATGCGCTTCGCTTATCAGGc                                                                                                                                                                                                                                                          <  2:379812/264‑1 (MQ=255)
              aagcaaagcaACTGGCTTTCCGGGCAACACGGAGAAGAAGGCGGTTCTATGCGCGGCGCAGCAAACGATCTGCGCAACCAGGCGGCTGCAATGGGCGGTAACGTGATTTATGGCATCAGTAGCCCGTCGCAGGGACTGTTGTCCAGTTTTGTCCCGACGGATAGCCAGCTTATCGGTCAGGTTTATAAGTGCCCGCACTGATGTGATCCAGCGGCTGGGTGACAAACGTATAATTGCCTGATGCGCTGCGCTTCTCAGGCCTTNTTAAGGCCtt                                                                                                                                                                                                                                            >  2:62364/1‑274 (MQ=255)
                              tagatccaatanactggccctctgcgcctggctgctcgtccccaatgcgtacactctgaccggcagcactcaattcgttgctggagctacaacccgccatcagtagcCCGTCGCTGGGAATGTTGTCCAGTTTTGTCCCGACGGATAGCCAGATTATCGGTCAGGTTTATAAGTGCCCGAACTGATGTGATCCAGAGGCTGGGTGACAAACGTATAATTGCCTGATGCGCTTCGCTTATCAGGCCTACATAAGGCATTGCAACTGATTGAATGCt                                                                                                                                                                                                                           <  2:65469‑M2/168‑1 (MQ=255)
                                                          ttCTATGCGCGGCGCAGCAAACGATCTGCGCAACCAGGCGGCTGCAATGGGCGGTAACGTGATTTATGGCATCAGTAGCCCGTCGCAGGGAATGTTGTCCAGTTTTGTCCCGACGGATAGCCAGATTATCGGTCAGGTTTATAAGTGCCCGAACTGATGTGATCCAGAGGCTGGGTGACAAACGTATAATTGCCTGATGcgcttcgc                                                                                                                                                                                                                                                                   >  1:197007/1‑207 (MQ=255)
                                                              aTGCGCGGCGCAGCAAACGATCTGCGCAACCAGGCGGCTGCAATGGGCGGTAACGTGATTTATGGCATCAGTAGCCCGTCGCAGGGAATGTTGTCCAGTTTTGTCCCGACGGATAGCCAGATTATCGGTCAGGTTTATAAGTGCCCGAACTGATGTGATCCAGAGGCTGGGTGACAAACGTATAATTGCCTGATGCGCTTCGCTTAt                                                                                                                                                                                                                                                               <  2:197004/207‑1 (MQ=255)
                                                                    ggcgCAGCAAACGATCTGCGCAACCAGGCGGCTGCAATGGGCGGTAACGTGATTTATGGCATCAGTAGCCCGTCGCAGGGAATGTTGTCCAGTTTTGTCCCGACGGATAGCCAGATTATCGGTCAGGTTTATAAGTGCCCGAACTGATGTGATCCAGAGGCTGGGTGACAAACGTATAATTGCCTGATGCGCTTCGCTTATc                                                                                                                                                                                                                                                              >  2:429169/1‑202 (MQ=255)
                                                                    ggcgCAGCAAACGATCTGCGCAACCAGGCGGCTGCAATGGGCGGTAACGTGATTTATGGCATCAGTAGCCCGTCGCAGGGAATGTTGTCCAGTTTTGTCCCGACGGATAGCCAGATTATCGGTCAGGTTTATAAGTGCCCGAACTGATGTGATCCAGAGGCTGGGTGACAAACGTATAATTGCCTGATGCGCTTCGCTTATCAGGCCTAc                                                                                                                                                                                                                                                      >  2:138224/1‑210 (MQ=255)
                                                                       gcagcaAACGATCTGCGCAACCAGGCGGCTGCAATGGGCGGTAACGTGATTTATGGCATCAGTAGCCCGTCGCAGGGAATGTTGTCCAGTTTTGTCCCGACGGATAGCCAGATTATCGGTCAGGTTTATAAGTGCCCGAACTGATGTGATCCAGAGGCTGGGTGACAAACGTATAATTGCCTGATGCGCTTCGCTTATCAGGCCTACATAAGGCATTGCAACTGATTGAATGCTTgaga                                                                                                                                                                                                                      <  1:404070/239‑1 (MQ=255)
                                                                        cagcaAACGATCTGCGCAACCAGGCGGCTGCAATGGGCGGTAACGTGATTTATGGCATCAGTAGCCCGTCGCAGGGAATGTTGTCCAGTTTTGTCCCGACGGATAGCCAGATTATCGGTCAGGTTTATAAGTGCCCGAACTGATGTGATCCAGAGGCTGGGTGACAAACGTATAATTGCCTGCTGCGCTTCGCTTATCAGGc                                                                                                                                                                                                                                                          <  1:587732/202‑1 (MQ=255)
                                                                        cagcaAACGATCTGCGCAACCAGGCGGCTGCAATGGGCGGTAACGTGATTTATGGCATCAGTAGCCCGTCGCAGGGAATGTTGTCCAGTTTTGTCCCGACGGATAGCCAGATTATCGGTCAGGTTTATAAGTGCCCGAACTGATGTGATCCAGAGGCTGGGTGACAAACGTATAATTGCCTGATGCGCTTCGCTTATCAGGc                                                                                                                                                                                                                                                          <  1:429176/202‑1 (MQ=255)
                                                                        cagcaAACGATCTGCGCAACCAGGCGGCTGCAATGGGCGGTAACGTGATTTATGGCATCAGTAGCCCGTCGCAGGGAATGTTGTCCAGTTTTGTCCCGACGGATAGCCAGATTATCGGTCAGGTTTATAAGTGCCCGAACTGATGTGATCCAGAGGCTGGGTGACAAACGTATAATTGCCTGATGCGCTTCGCTTATCAGGCCTACATaa                                                                                                                                                                                                                                                  <  1:138225/210‑1 (MQ=255)
                                                                         agcaAACGATCTGCGCAACCAGGCGGCTGCAATGGGCGGTAACGTGATTTATGGCATCAGTAGCCCGTCGCAGGGAATGTTGTCCAGTTTTGTCCCGACGGATAGCCAGATTATCGGTCAGGTTTATAAGTGCCCGAACTGATGTGATCCAGAGGCTGGGTGACAAACGTATAATTGCCTGATGCGCTTCGCTTATCAGGCCTACATAAGGCATTGCAACTGATTGAATGCTTGAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCgg                                                                                                                                                                                 >  1:224897/1‑274 (MQ=255)
                                                                                     gcgcAACCAGGCGGCTGCAATGGGCGGTAACGTGATTTATGGCATCAGTAGCCCGTCGCAGGGAATGTTGTCCAGTTTTGTCCCGACGGATAGCCAGATTATCGGTCAGGTTTATAAGTGCCCGAACTGATGTGATCCAGAGGCTGGGTGACAAACGTATAATTGCCTGATTCGCTACGCTTATCAGGCCTACATCGTGCATTGAAATTGAATgga                                                                                                                                                                                                                               >  1:102222/1‑214 (MQ=39)
                                                                                     gcgcAACCAGGCGGCTGCAATGGGCGGTAACGTGATTTATGGCATCAGTAGCCCGTCGCAGGGAATGTTGTCCAGTTTTGTCCCGACGGATAGCCAGATTATCGGTCAGGTTTATAAGTGCCCGAACTGATGTGATCCAGAGGCTGGGTGACAAACGTATAATTGCCTGATGCGCTTCGCTTATc                                                                                                                                                                                                                                                              >  1:198585/1‑185 (MQ=255)
                                                                                     gcgcAACCAGGCGGCTGCAATGGGCGGTAACGTGATTTATGGCATCAGTAGCCCGTCGCAGGGAATGTTGTCCAGTTTTGTCCCGACGGATAGCCAGATTATCGGTCAGGTTTATAAGTGCCCGAACTGATGTGATCCAGAGGCTGGGTGACAAACGTATAATTGCCTGATGCGCTTCGCTTATCAGGCCTACATAAGGCATTGCAACTGATTGCATGCTTGAGAATTTGTAGGACGGATAAg                                                                                                                                                                                                    >  1:176121/1‑243 (MQ=255)
                                                                                     gcgcAACCAGGCGGCTGCAATGGGCGGTAACGTGATTTATGGCATCAGTAGCCCGTCGCAGGGAATGTTGTCCAGTTTTGTCCCGACGGATAGCCAGATTATCGGTCAGGTTTATAAGTGCCCGAACTGATGTGATCCAGAGGCTGGGTGACAAACGTATAATTGCCTGATGCGCTTCGCTTATCAGGCCTACATAAGGCATTGCAACTGATTGAATGCTTGAGAATTTGTAAGCCGGCTAAGGCGTTTACGACGCATCCGGCATGACCAACGaa                                                                                                                                                                    >  1:132372/1‑274 (MQ=255)
                                                                                     gcgcAACCAGGCGGCTGCAATGGGCGGTAACGTGATTTATGGCATCAGTAGCCCGTCGAAGGGAATGTTGTCCAGTTTTGTCCCGACGGATAGCCAGATTATCGGTCAGGTTTATAAGTGCCCGAACTGATGTGATCCAGAGGCTGGGTGAAAAACGTAGAATTGCCTGATGCGCGTCGCTTATc                                                                                                                                                                                                                                                              >  2:11055/1‑185 (MQ=255)
                                                                                         aaCCAGGCGGCTGCAATGGGCGGTAACGTGATTTATGGCATCAGTAGCCCGTCGCAGGGAATGTTGTCCAGTTTTGTCCCGACGGATAGCCAGATTATCGGTCAGGTTTATAAGTGCCCGAACTGATGTGATCCAGAGGCTGGGTGACAAACGTATAATTGCCTGATTCGCTACGCTTATCAGGCCTACATCGTGCATTGAAATTGAATGGAtttg                                                                                                                                                                                                                           <  2:102221/216‑4 (MQ=38)
                                                                                         aaCCAGGCGGCTGCAATGGGCGGTAACGTGATTTATGGCATCAGTAGCCCGTCGCAGGGAATGTTGTCCAGTTTTGTCCCGACGGATAGCCAGATTATCGGTCAGGTTTATAAGTGCCCGAACTGATGTGATCCAGAGGCTGGGTGACAAACGTATAATTGCCTGATGCGCTTCGCTTATCAGGc                                                                                                                                                                                                                                                          <  2:198582/185‑1 (MQ=255)
                                                                                         aaCCAGGCGGCTGCAATGGGCGGTAACGTGATTTATGGCATCAGTAGCCCGTCGCAGGGAATGTTGTCCAGTTTTGTCCCGACGGATAGCCAGATTATCGGTCAGGTTTATAAGTGCCCGAACTGATGTGATCCAGAGGCTGGGTGACAAACGTATAATTGCCTGATGCGCTTCGCTTATCAGGc                                                                                                                                                                                                                                                          <  1:11055/185‑1 (MQ=255)
                                                                                          aCCAGGCGGCTGCAATGGGCGGTAACGTGATTTATGGCATCAGTAGCCCGTCGCAGGGAATGTTGTCCAGTTTTGTCCCGACGGATAGCCAGATTATCGGTCAGGTTTATAAGTGCCCGAACTGATGTGATCCAGAGGCTGGGTGACAAACGTATAATTGCCTGATGCGCTTCGCTTATCAGGCCTACATAAGGCATTGCAACTGATTGAATGCTTGAGAATTTGTAGGCCGGATAGGGCGTTTACGCCGCATGAGGCATGACAAACGAtaccgt                                                                                                                                                               >  2:127663/1‑270 (MQ=255)
                                                                                          aCCAGGCGGCTGCAATGGGCGGTAACGTGATTTATGGCATCAGTAGCCCGTCGCAGGGAATGTTGTCCAGTTNTGTCCCGACGGATAGCCAGATTATCGGTCAGGTTTATAAGTGCCCGAACTGATGTGATCCAGAGGCTGGGTGACAAACGTATAATTGCCTGATGCGCTTCGCTTATCAGGCCTACATAAGGCATTGCAACTGATTGAATGCTTGAGAATTTg                                                                                                                                                                                                                 <  1:620132/225‑1 (MQ=255)
                                                                                                cggcTGCAATGGGCGGTAACGTGATTTATGGCATCAGTAGCCCGTCGCAGGTAATGTTGTCCAGTTTTGTCCCGACGGATAGCCAGAGTATCGGTCAGGTTTATTAGTGCCCGAACTGATGTGATCCAGAGGCTGGGTGACAAACGTATAATTGCCTGATGCGCTTCGCTTATCAGGCCTACATAAGGCATTGCAACTGATTGAATGCTTGAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGACCAACGATCACGTTATCAg                                                                                                                                                         <  1:455640/275‑1 (MQ=255)
                                                                                                   cTGCAATGGGCGGTAACGTGATTTATGGCATCAGTAGCCCGTCGCAGGGAATGTTGTCCAGTTTTGTCCCGACGGATAGCCAGATTATCGGTCAGGTTTATAAGTGCCCGAACTGATGTGATCCAGAGGCTGGGTGACAAACGTATAATTGCCTGATGCGCTTCGCTTATc                                                                                                                                                                                                                                                              >  2:564879/1‑171 (MQ=255)
                                                                                                       aaTGGGCGGTAACGTGATTTATGGCATCAGTAGCCCGTCGCAGGGAATGTTGTCCAGTTTTGTCCCGACGGATAGCCAGATTATCGGTCAGGGTTATAAGTGCCCGAACTGATGTGATCCAGAGGCTGGGTGACAAACGTATAATTGCCTGATGCGCTTCGCTTATCAGGc                                                                                                                                                                                                                                                          <  1:564888/171‑1 (MQ=255)
                                                                                                              ggTTACTTGATTTATGGCATCAGTAGCCCGTCGTAGGGAATGTTGTCCAGTTTTGTCCCGACGGATAGCCAGATTATCGGTCAGGTTTATAAGTGCCCGAACTTCTGTGATCCAGAGGCNGTGTGACACACGTATAATTGCCTGATGCGCTTCGCTTATCAGGTCTACATAAGGCATTTCAACTGATTGAATTCTTGAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGACCAACGATCACGTTATCAGTCATTTCAGGc                                                                                                                                              <  1:596737/272‑1 (MQ=255)
                                                                                                                                 tcagtagcCCGTCGCAGGGAATGTTGTCCAGTTTTGTCCCGACGGATAGCCAGATTATCGGTCAGGTTTATAAGTGCCCGAACTGATGTGATCCAGAGGCTGGGTGACAAACGTATAATTGCCTGATGCGCTTCGCTTATCAGGCCTACATAAGGCATTGCAACTGATTGAATGCTTGAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGACCAACGATCACGTTATCAGTCATTTCAGGCTGGATGGCTCCAGCCTCGCAtt                                                                                                                        >  1:267019‑M2/9‑275 (MQ=255)
                                                                                                                                             ccgcgggAANGTTGTCCAGTTTTGTCCCGACGGATAGCCAGATTATCGGTCAGGTTTATAAGTGCCCGAACTGATGTGATCCAGAGGCTGGGTGACAAACGTATAATTGCCTGATGCGCTTCGCTTATCAGGCCTACATAAGGCATTGCAACTGATTGAATGCTTGAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGACCAACGATCACGTTATCAGTCATTTCAGGCTGGATGGCTCCAGCCTCGCATTAGCAAAGATgag                                                                                                            <  2:267014/271‑1 (MQ=255)
                                                                                                                                                      tgttgtCCAGTTTTGTCCCGACGGATAGCCAGATTATCGGTCAGGTTTATAAGTGCCCGAACTGATGTGATCCAGAGGCTGGGTGACAAACGTATAATTGCCTGATGCGCTTCGCTTATCAGGCCTACATAAGGCATTGCAACTGATTGAATGCTTGAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGACCaa                                                                                                                                                                        >  2:38257/1‑206 (MQ=255)
                                                                                                                                                          gtCCAGTTTTGTCCCGACGGATAGCCAGATTATCGGTCAGGTTTATAAGTGCCCGAACTGATGTGATCCAGAGGCTGGGTGACAAACGTATAATTGCCTGATGCGCTTCGCTTATCAGGCCTACATAAGGCATTGCAACTGATTGAATGCTTGAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGACCAACGAt                                                                                                                                                                    <  1:38257/206‑1 (MQ=255)
                                                                                                                                                            ccAGTTTTGTCCCGACGGATAGCCAGATTATCGGTCAGGTTTATAAGTGCCCGAACTGATGTGATCCAGAGGCTGGGTGACAAACGTATAATTGCCTGATGCGCTTCGCTTATCAGGCCTACATAAGGCATTGCAACTGATTGAATGCTTGAGAATTTGTAGGCCGGATAAgg                                                                                                                                                                                                   >  2:539796/1‑173 (MQ=255)
                                                                                                                                                              aGTTTTGTCCCGACGGATAGCCAGATTATCGGTCAGGTTTATAAGTGCCCGAACTGATGTGATCCAGAGGCTGGGTGACAAACGTATAATTGCCTGATGCGCTTCGCTTATCAGGCCTACATAAGGCATTGCAACTGATTGAATGCTTGAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGACCAACGAGCACGTTATCAGTCATTTCAGGCTGGATGGCTCCAGCCTCGCATTAGCAAAGAAAAAATAATTTCGACTgtggt                                                                                           >  1:274495/1‑275 (MQ=255)
                                                                                                                                                                ttttGTCCCGACGGATAGCCAGATTATCGGTCAGTTTTATAAGTGCCCGAACTGATGTGATCCAGAGGCTGGGTGACAAACGTATAATTGCCTGATGCGCTTCGCTTATCAGGCCTACATAAGGCAtt                                                                                                                                                                                                                                            >  1:466017/1‑128 (MQ=255)
                                                                                                                                                                ttttGTCCCGACGGATAGCCAGATTATCGGTCAGGTTTATAAGTGCCCGAACTGATGTGATCCAGAGGCTGGGTGACAAACGTATAATTGCCTGATGCGCTTCGCTTATCAGGCCTACATAAGGCATTGCAACTGATTGAATGCTTGAGAATTTGTAGGCCGGATAAGGCGtt                                                                                                                                                                                               <  1:539804/173‑1 (MQ=255)
                                                                                                                                                                ttttGTCCCGACGGATAGCCAGATTATCGGTCAGGTTTATAAGTGCCCGAACTGATGTGATCCAGAGGCTGGGTGACAAACGTATAATTGCCTGATGCGCTTCGCTTATCAGGCCTACATAAGGCATTGCAACTGATTGAATGCTTGAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGACCAACGATCACGTTATCAGTCATTTCAGGCTGGATGGCTCCa                                                                                                                                  >  2:706737/1‑234 (MQ=255)
                                                                                                                                                                    gTCCCGACGGATAGCCAGATTATCGGTCAGGTTTATAAGTGCCCGAACTGATGTGATCCAGAGGCTGGGTGACAAACGTATAATTGCCTGATGCGCTTCGCTTATCAGGCCTACATAAGGCATTGCaa                                                                                                                                                                                                                                        <  2:466010/128‑1 (MQ=255)
                                                                                                                                                                    gTCCCGACGGATAGCCAGATTATCGGTCAGGTTTATAAATGCCCGAACTGATGTGATCCAGAGGCTGGGTGACAAACGTATAATTGCCTGATGCGCTTCGCTTATCAGGCCTACATAAGGCATTGCAACTGATTGAATGCTTGAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGACCAACGATCACGTTATCAGTCATTTCAGGCTGGATAGCTCCAGCCt                                                                                                                              <  1:706747/234‑1 (MQ=255)
                                                                                                                                                                     tCCCGACGGATAGCCAGATTATCGGTCAGGTTTATAAGTGCCCGAACTGATGTGATCCAGAGGCTGGGTGACAAACGTATAATTGCCTGATGCGCTTCGCTTATCAGGCCTACATAAGGCATTGCAACTGATTGAATGCTTGAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGACCAACGATCACGTTATCAGTCAttt                                                                                                                                                   >  2:520/1‑212 (MQ=255)
                                                                                                                                                                         gACGGATAGCCAGATTATCGGTCAGGTTTATAAGTGCCCGAACTGATGTGATCCAGAGGCTGGGTGACAAACGTATAATTGCCTGATGCGCTTCGCTTATCAGGCCTACATAAGGCATTGCAACTGATTGAATGCTTGAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGACCAACGATCACGTTATCAGTCATTTCAgg                                                                                                                                               <  1:520/212‑1 (MQ=255)
                                                                                                                                                                               tggCCAGATNTTCGGTCAGGTTTATAAGGTCCCTAACTGATGTGTTCCAGAGGCTGGTTTATAAACGTATAATTGCCTGATGCGCTTCGCTTATCAGGCCTACATAAGGCATTGCAACTGATTGACTGCTTGAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGACCCACGATCACGTTATCAGTCATTTCAGGCTGGATGGCTCCAGCCTCGCATTATCAAAGATGAGATTATTTCGCCTGTGGTGCAGTTTTGCtggtgca                                                                          <  2:274490/273‑3 (MQ=255)
                                                                                                                                                                                  ccAGATTATCGGTCAGGTTTATAAGTGCCCGAACTGATGTGATCCAGAGGCTGGGTGACAAACGTATAATTGCTTGATGCGCTTCGCTTATCAGGCCTACATAAGGCATTGCAACTGATTGAATGATTGAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGACCAACGATCACGTTATCAGTCATTTCAGGCTGGATGGCTCCAGCCTCGCATTAGCAAAGATGAGATTATTTCGCCTGTGGTGCAGTTTTGCTGGTGGATTg                                                                       <  1:587960/275‑1 (MQ=255)
                                                                                                                                                                                       ttATCGGTCAGGTTTATAAGTGCCCGAACTGATGTGATCCAGAGGCTGGGTGACAAACGTATAATTGCCTGATGCGCTTCGCTTATCAGGCCTACATAAGGCATTGCAACTGATTGAATGCTTGAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGACCAACGATCACGTTATCAGTCATTTCAGGCTGGATGGCTCCAGCCTCGCATTAGCAAAGATGAGATTATTTCGCCTGTGGTGCAGTTTTGCTGGTGGATTGAacac                                                                  <  1:127664/275‑1 (MQ=255)
                                                                                                                                                                                          gcgagagggTGTATAAGGTCCTGAACTGATGTGATCCAGATGCTGGGTGACAATCGTATAATTGCCTGATGCGCTTCGCGTATTAGGACTACATAAGGCATTGCAACTGCTTGACTGCTTGGGAAGTTGGAGGCCGGATAAGGCTATTACGCCGCACACGGCATTAACACCGATCACGTTATCAGGCATTTCAGGCTGGATGGCTCCAGTCTCGCATTAGCAAAGATGAGACTATTTCGCCTGTGGTGCAGTGTTGCTGGTGGATTGAACACGt                                                                <  2:210381/267‑1 (MQ=255)
                                                                                                                                                                                           agnTCAGTNNNNTAAGTGCCCGAGCGGATGTTGTCCAGAGGCGGGGTGAAAAACGGATAATTGCCTGATGCGCTTCGCGTATCAGGCCTACATAAGGCATTGCAACTGATTGAATGCTTGAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGACCAACGATCACGTTATCAGTCATTTCAGGCTGGATGGCTCCAGCCTCGCATTAGCAAAGATGAGATTATTTCGCCTGTGGTGCAGTTTTGCTGGTGGATTGAACACGt                                                                <  2:51424/272‑1 (MQ=255)
                                                                                                                                                                                                aGGTTTATAAGTGCCCGAACTGATGTGATCCAGAGGCTGGGTGACAAACGTATAATTGCCTGATGCGCTTCGCTTATCAGGCCTACATAAGGCATTGCAACTGATTGAATGCTTGAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGACCAACGATCACGTTATCAGTCATTTCAGGCTGGATGGCTCCAGCCTCGCATTAGCAAAGATGAGATTATTTCGCCTGTGGTGCAGTTTTGCTGGTGGATTGAac                                                                    >  1:51424/1‑264 (MQ=255)
                                                                                                                                                                                                    ttATAAGTGCACGAACTGATGTGATCCAGAGGCTGGGTGACAAACGTATAATTGCCTGATGCGCTTCGCTTATCAGGCCTACATAAGGCATTGCAACTGATTGAATGCTTGAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGACCAACGATCACGTTATCAGTCATTTCAGGCTGGATGGCTCCAGCCTCGCATTAGCAAAGATGAGATTATTTCGCCTGTGGTGCAGTTTTGCTGGTGGATTGAACACGTAAGTGATGCa                                                      <  1:62365/274‑1 (MQ=255)
                                                                                                                                                                                                            gCCCGAACTGATGTGATCCAGAGGCTGGGTGACAAACGTATAATTGCCTGATGCGCTTCGCTTATCAGGCCTACATAAGGCATTGCAACTGTTTGAATGCTTGAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCCTCCGGCATGACCAACGATCACGTTATCAGTCATTTCAGGCTGGATGGCTCCCGCCTCGCATTAGCAAAGATGAGATTATTTCGCCTGTGGTGCAGTTTTTCTGGTGGATTGAac                                                                    >  1:210385/1‑252 (MQ=255)
                                                                                                                                                                                                             cccGAACTGATGTGATCCAGAGGCTGGGTGACAAACGTATAATTGCCTGATGCGCTTCGCTTATCAGGCCTACATAAGGCATTGCAACTGATTGAATGCTTGAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCa                                                                                                                                                                               >  2:707875/1‑144 (MQ=255)
                                                                                                                                                                                                                taaCTGATTTGAGTNTGATGCTGGGTGACAATCTTTTACTTTCCTGTTGTGGTTCGCTTATACGGCATACATAAAGCTTTGCAACTGAGTGAATGCTTGAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGACCAACGATCACGTTATCAGACATTTCATGCTGGATTGCTCCAGCCTCGCATTAGCAAAGATGAGATTATTTCGCCTGTGGTGCAGTTTTGCTGGTGGATTGAACACGTAAGTGATGCAGCGCGCGTAAGCa                                         <  2:518092/274‑1 (MQ=255)
                                                                                                                                                                                                                 tatggtgttgAGCCATATGCTGTGTGACAAACTTATAATTTCCTGTTTCGCTTCTCTTATCTGGCCTACATTAGGCATTGCAACTGATTGTATTCTTGAGAATTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGTCATGACCAACGATCACGTTATCAGTCATTTCAGGCTGGATGGCTCCAGCCTCGCATTAGCAAAGATGAGATTATTTCGCCTGTGGTGCAGTTTTGCTGGTGGATTGAACACGTAAGTGATGCAGCGCGCGTAAGCa                                         <  2:237486/266‑1 (MQ=255)
                                                                                                                                                                                                                 aaCTGATGTGATCCAGAGGCTGGGTGACAAACGTATAATTGCCTGATGCGCTTCGCTTATCAGGCCTACATAAGGCATTGCAACTGATTGAATGCTTGAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGAc                                                                                                                                                                           <  1:707885/144‑1 (MQ=255)
                                                                                                                                                                                                                            tCCAGAGGCTGGGTGACAAACGTATAATTGCCTGATGCGCTTCGCTTATc                                                                                                                                                                                                                                                              >  2:9349/1‑50 (MQ=255)
                                                                                                                                                                                                                                ggggCAGGGTTACAAAAGTATAATTGCCTGTGGCGCTTCTCTTATCTAGCCTACATTAGGCATTGCAACTGATTGAATGCTTGAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGACCAACGATCACGTTATCAGTCATTTCAGGCTGGATGGCTCCAGCCTCGCATTAGCAAAGATGAGATTATTTCGCCTGTGGTGCAGTTTTGCTGGTGGATTGAACACGTAAGTGATGCAGCGCGCGTAAGCATTGATCCTGGGTTAg                          <  2:224892/272‑1 (MQ=255)
                                                                                                                                                                                                                                gagGCTGGGTGACAAACGTATAATTGCCTGATTcgcta                                                                                                                                                                                                                                                                      <  2:320959/38‑2 (MQ=255)
                                                                                                                                                                                                                                gagGCTGGGTGACAAACGTATAATTGCCTGATGCGCTTCGCTTATc                                                                                                                                                                                                                                                              >  1:388109/1‑46 (MQ=255)
                                                                                                                                                                                                                                gagGCTGGGTGACAAACGTATAATTGCCTGATGCGCTTCGCTTATCAGGc                                                                                                                                                                                                                                                          <  1:9349/50‑1 (MQ=255)
                                                                                                                                                                                                                                gagGCTGGGTGACAAACGTATAATTGCCTGATGCGCTTCGCTTATCAGGCCTACATAAGGCATTGCAACTGATTGAATGCTTGAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCCTCCGGCATGACCAACGATCACGTTATCAGTCATTTCAGGCTGGATGGCTCCAGCCTCGCATTAGCNAAGATGAGNTTATTTCGCCTGTGGTGCAGTTTTGCTGGTGGATTGAACACGTAAGTGATGCAGCGCGCGTAAGCATTGATCCTGGGTTAg                          >  1:477367/1‑274 (MQ=255)
                                                                                                                                                                                                                                gagGCTGGGTGACAAACGTATAATTGCCTGATGCGCTTCGCTTATCAGGCCTACATAAGGCATTGCAACTGATTGAATGCTTGAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGACCaa                                                                                                                                                                        >  1:685407/1‑132 (MQ=255)
                                                                                                                                                                                                                                gagGCTGGGTGACAAACGTATAATTGCCTGATGCGCTTCGCTTATCAGGCCTACATAAGGCATTGCAACTGATTGAATGCTTGAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGACCAACGATCACGTTATCAGTCATTTCAGGCTGGATGGCTCCAGCCTCGCATTAGCAAAGATGAGATTAATTCGCCTGTGGTGCCAGTTTGCTGGGGGATTGAACAAGTAAGTGATGGAGCGCGCGTAAGCATTGATCCGGGGtgatc                         >  2:720712/1‑271 (MQ=255)
                                                                                                                                                                                                                                gagGCTGGGTGACAAACGTATAATTGCCTGATGCGCTTCGCCTATCAGGCCTACATAAGGCATTGCAACTGATTGAATGCTTGAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCa                                                                                                                                                                               >  2:633868/1‑125 (MQ=255)
                                                                                                                                                                                                                                    cTGGGTGACAAACGTATAATTGCCTGATGCGCTTCGCTTATCAGGc                                                                                                                                                                                                                                                          <  2:388103/46‑1 (MQ=255)
                                                                                                                                                                                                                                    cTGGGTGACAAACGTATAATTGCCTGATGCGCTTCGCTTATCAGGCCTACATAAGGCATTGCAACTGATTGAATGCTTGAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGACCAACGAt                                                                                                                                                                    <  2:685398/132‑1 (MQ=255)
                                                                                                                                                                                                                                    cTGGGTGACAAACGTATAATTGCCTGATGCGCTTCGCCTATCAGGCCTACATAAGGCATTGCAACTGATTGAATGCTTGAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGAc                                                                                                                                                                           <  1:633877/125‑1 (MQ=255)
                                                                                                                                                                                                                                     tGGGTGACAAACGTATAATCGCCTGATGCGCTTCGCTTATCAGGCCTACATAAGGCATTGCAACTGATTGAATGCTTGAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGACCAACGATCACGTTATCAGTCATTTCAGGCTGGATGGCTCCAGCCTCGCATTAGCAAAGATGAGATTATTTCGCCTGTGGTGCAGTTTTGCTGGTGGATTGAACACGTAAGTGATGAAGCGCGCGTAAGCATTGATCCTGGGTTAGCGGtt                     >  1:574586/1‑274 (MQ=255)
                                                                                                                                                                                                                                       ggTGACAAACGTATAATTGCCTGATGCGCTTCGCTTATCAGGCCTACATAAGGCATTGCAACTGATTGAATGCTTGAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGACCAACGATCACGTTATCAGTCATTTCAGGCTGGATGGCTCCAGCCTCGCATTAGCAAAGATGAGATTATTTCGCCTGTGGTGCAGTTTTGCTGGTGGAt                                                                         >  1:307102/1‑220 (MQ=255)
                                                                                                                                                                                                                                       ggTGACAAACGTATAATTGCCTGATGCGCTTCGCTTATCAGGCCTACATAAGGCATTGCAACTGATTGAATGCTTGAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGACCAACGATCACGTTATCAGTCATTTCAGGCTGAATGGCTCCAGCCTCGCATTAGCAAAGATGAGATTATTTCGCCTGTGGTTCAGTTTTGCTGGTTGATTGAACACGTAAGTTATGCAGCTCGCGTaa                                            >  1:232000/1‑249 (MQ=255)
                                                                                                                                                                                                                                        gTGACAAACGTATAATTGCCTGATGCGCTTCGCTTATCAGGCCTACATAAGGCATTGCAACTGATTGAATGCTTGAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGACCAACGATCACGTTATCAGTCATTTCAGGCTGGATGGCTCCAGCCTCGCATTAGCAAAGATGAGATTATTTCGCCTGTGGTGCAGTTTTGCTGGTGGATTGAACACGTAAGTGATGCAGCGCGCGTAAGCATTGATCCTGGGTTAGCGGTTggg                  >  1:299683/1‑274 (MQ=255)
                                                                                                                                                                                                                                           aaaaaaGTAAAATTGCCTGATGCGCTTCGCTTATCAGGCCTACATAAGGCATTGCAACTGATTGAATGCTTGAGAATTGGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGACCAACGATCACGTTATCAGTCATTTCAGGCTGGATGGCTCCAGCCTCGCATTAGCAAAGATGAGATTATTTCGCCTGTGGTGCAGTTTTGCTGGTGGATTGaa                                                                     <  2:307096/218‑1 (MQ=255)
                                                                                                                                                                                                                                                 gTATAATTGCCTGATGCTCTTCGCTTATCAGGCCTACATAAGGCATTGCAACTGATTGAATGCTTGAGAATTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGGCATGACCAACGATCACGTTATCAGTCATTTCAGGCTGGATGGCTCCAGCCTCGCATTAGCAAAGATGAGATTATTTCGCCTGTGGTGCAGTTTTGCTGGTGGATTGAACACGTAAGTGATGCAGCGCGCGTAAGCATTGATCCTGGGTTAGCGGTTGAGCGTTATTcg         <  1:720722/274‑1 (MQ=255)
                                                                                                                                                                                                                                                  tataATTGCCTGATGCGCTTCGCTTATCAGGCCTACATAAGGCATTGCAACTGATTGAATGCTTGAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGACCAACGATCACGTTATCAGTCATTTCAGGCTGGATGGCTACAGCCTCGCATTAGCAAAGATGAGATTATTTCGCCtgt                                                                                              >  1:581428/1‑188 (MQ=255)
                                                                                                                                                                                                                                                      aTTGCCTGATGCGCTTCGCTTATCAGGCCTACATAAGGCATTGCAACTGATTGAATGCTTGAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGACCAACGATCACGTTATCAGTCATTTCAGGCTGGATGGCTCCAGCCTCGCATTAGCAAAGATGAGATTATTTCGCCTgtggtg                                                                                          <  2:581419/188‑1 (MQ=255)
                                                                                                                                                                                                                                                          ccTGATTCGCTACGCTTATCAGGCCTACATAAGGCATTGCAACTGATTGAATGCTTGAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGACCAACGATCACGTTATCAGTCATTTCAGGCTGGATGGCTCCAGCCTCGCATTAGCCAAGATGAGATTATTTCGCCTGTGGTGCAGTTTTGCTGGTGGATTGAACACGTAAGTGATGCAGCGCGCGTACGCATTGATCCTGGGTTAGCGGTTGGGCGTTATTCGCCGCTAc  >  1:48475/1‑272 (MQ=255)
                                                                                                                                                                                                                                                           cTGATGCGCTTCGCTTATCAGGCCTACATAAGGCATTGCAACTGATTGAATGCTTGAGAATTTGTAGGCCGGATAAGGCGTTTAcgct                                                                                                                                                                                         >  1:724376/1‑87 (MQ=255)
                                                                                                                                                                                                                                                               tGCGCTTCGCTTATCAGGCCTACATAAGGCATTGCAACTGATTGAATGCTTTAGAATTTGTAGGCCGGATAAGGCGTTTACGCTGCAt                                                                                                                                                                                     <  2:724366/88‑1 (MQ=255)
                                                                                                                                                                                                                                                               tGCGCTACGCTTATCAGGCCTACATAAGGCATTGCAACTGATTGAATGCTTGAGAATTTGTAGGCCGGATAAGGCGtt                                                                                                                                                                                               >  2:708559/1‑78 (MQ=255)
                                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                                         
GTACTGCGACAGGTAAGCAAAGCAACTGGCTTTCCGGGCAACACGGAGAAGAAGGCGGTTCTATGCGCGGCGCAGCAAACGATCTGCGCAACCAGGCGGCTGCAATGGGCGGTAACGTGATTTATGGCATCAGTAGCCCGTCGCAGGGAATGTTGTCCAGTTTTGTCCCGACGGATAGCCAGATTATCGGTCAGGTTTATAAGTGCCCGAACTGATGTGATCCAGAGGCTGGGTGACAAACGTATAATTGCCTGATGCGCTTCGCTTATCAGGCCTACATAAGGCATTGCAACTGATTGAATGCTTGAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGACCAACGATCACGTTATCAGTCATTTCAGGCTGGATGGCTCCAGCCTCGCATTAGCAAAGATGAGATTATTTCGCCTGTGGTGCAGTTTTGCTGGTGGATTGAACACGTAAGTGATGCAGCGCGCGTAAGCATTGATCCTGGGTTAGCGGTTGGGCGTTATTCGCCGCTAC  >  CP009273/4362765‑4363286

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 10 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: