Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 244,782:1 +G 8.9% intergenic (+161/+63) rayT → / ← lfhA RAYT REP element‑mobilizing transposase; TnpA(REP)/pseudogene, flagellar system protein, promoterless fragment; flagellar biosynthesis

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273244,7821.G 8.9% 165.9 / 14.0 45intergenic (+161/+63)rayT/lfhARAYT REP element‑mobilizing transposase; TnpA(REP)/pseudogene, flagellar system protein, promoterless fragment; flagellar biosynthesis
Reads supporting (aligned to +/- strand):  ref base . (21/20);  new base G (2/2);  total (23/22)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.04e-01

GTGAGTGATTGGCCATTCTCAACGTTCCATCGCGATGTCGCGCGAGGGTTATATCCCATCGATTGGGCGGGGGACGTAACGGATTTTAGTGCCGGGGAGCGTATCATTTCATAATTGTGCGCAGATGCCTGATGCGACGCTAGCGCGTCTTATCATGCCTACAAACTTGTGCCGGATCGGTAGGCCGGATAAGGCGTTTATGCCGCATCCGGCAGTTATGCGCAGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACAAACTTGTGCCTGAACGGTAGGACGGATAAGGCGTTCACGCCGCATCCGGCAGTTGTGCACCAATGCTGGATACTTACATCAACGCCCCCATAATCTTCAACTCCAGCTCATCCGGCACTTCGTTATACGACAGCACATGCAGCCCCGGCGCAAACAACCTTGCATAACGCGCCAGCAAAGGGCGCAGCTGCGGTGGCACCAGCAGCACCGGGTCTTTCCCCGCCGCTTTCATCTGCTCTTTCACCTGTGGCATCGTACTCTGGAACTGGTTGAGCAT  >  CP009273/244508‑245048
                                                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                                          
gtgagtgaTTGGCCATTCTCAACGTTCCATCGCGATGTCGCGCGAGGGTTATATCCCATCGATTGGGCGGGGGACGTAACGGATTTTAGTGCCGGGGAGCGTATCATTTCATAATTGTGCGCAGATGCCTGATGCGACGCTAGCGCGTCTTATCATGCCTACAAACTTGTGCCGGATCGGTAGGCCGGATAAGGCGTTTATGCCGCATCCGGCAGTTATGCGCAGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACAAACTTGTGCCt                                                                                                                                                                                                                                                                             >  1:121549/1‑275 (MQ=255)
                     aCGTTCCATCGCGATGTCGCGCGAGGGTTATATCCCATCGATTGGGCGGGGGACGTAACGGATTTTAGTGCCGGGGAGCGTATCATTTCATAATTGTGCGCAGATGCCTGATGCGACGCTAGCGCGTCTTATCATGCCTACAAACTTGTGCCGGATCGGTAGGCCGGCTAAGGCGTTTATGCCGCATCCGGCAGTTATGCGCAGAAGCCTGATGGGACGCTGGCGCGGCTTAGCATGCCTACCAACTTGGGCCTGAACGGTAGGAGGGCTAAgg                                                                                                                                                                                                                                                        >  2:304803/1‑274 (MQ=255)
                          ccATCGCGATGTCGCGCGAGGGTTATATCCCATCGCTTGGGCGGGGGACGTAACGGATTTTAGTGCCGGGGAGCGTATCATTTCATAATTGTGCGCAGATGCCTGATGCGACGCTAGCGCGTCTTATCATGCCTACAAACTTGTGCCGGATCGGTAGGCCGGATAAGGCGTTTATGCCCCATCCGGCAGTTATGCTCAGATGCCTGATGCGACGCTCGCGCGTCTTAACATTCCCACAACCTTGTGCCTGACCGGTAGGCCGGAGNAGGCggcct                                                                                                                                                                                                                                                  >  2:277354/1‑271 (MQ=38)
                          ccATCGCGATATCGCGCGAGGGTTATATCCCATCGATTGGGCGGGGGACGTAACGGATTTTAGTGCCGGGGAGCGTATCATTTCATAATTGTGCGCAGATGCCTGATGCGACGCTAGCGCGTCTTATCATGCCTACAAACTTGTGCCGGATCGGTAGGCCGGATAAGGCGTTTATGCCGCATCCGGCCGTTATGCGCATATGCCTGATGCGATGCTGGCCCGTCTTCTTATGCCCACAAACTTGGGCCTGAACGGTAGGAccttaaaggtgagg                                                                                                                                                                                                                                                   >  2:105129/1‑261 (MQ=38)
                                        ctcggacttgatatCCCATCGATTGGGCGGGGGACGTAACGGATTTTAGTGCCGGGGATCGTATCATTTAATAATTGTGCGCAGATGCCTGATGCGAAGCTAGCGCGTCTTATCATGCCTACAAAATTGTGCCGGATCGGTAGGCCGGATAAGGCGTTTATGCCGCATCCGGCAGTTATGCGCAGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACAAACTTGTGCCTGAACGGTAGGACGGATAAGGCGTTCACGCCGCATCCGGc                                                                                                                                                                                                                                     <  2:121548/264‑1 (MQ=255)
                                                       ccgtagagtggaagggggAAGTANGGGGTTTTGGTGCCGGGGGGCAGATCATTTCATAATTGTGCGCAGATGCCTGATGCGACGCTAGCGGGTCTTTTCATGCCTACAAACTTGTGCCGGAGCGGTAGGCCGGATAAGGCGTTTATGCCGCATCCGGCAGTTATGCGCAGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACAAACTTGTGCCTGAACGGTAGGACGGATAAGGCGTTCACGCCGCATCCGGCAGTTGTGCACCAATg                                                                                                                                                                                                                      <  2:432286/261‑1 (MQ=25)
                                                                 ggCGGGGGACGTAACGGATTTTAGTGCCGGGGAGCGTATCATTTCATAATTGTGCGCAGATGCCTGATGCGACGCTAGCGCGTCTTATCATGCCTACAAACTTGTGCCGGATCGGTAGGCCGGATAAGGCGTTTATGCCGCATCCGGCAGTTATGCGCAGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACAAACTTGTTCCTGAACGGTAGGACGGATAAGGCGTTCACGCCGCAACCGGCCGTTGTTCACCaagtcttgattcct                                                                                                                                                                                                            >  2:642420/1‑262 (MQ=255)
                                                                 ggCGGGGGACGTAACGGATTTTAGTGCCGGGGAGCGTATCATTTCATAATTGTGCGCAGATGCCTGATGCGACGCTAGCGCGTCTTATCATGCCTACAAACTTGTGCCGGATCGGTAGGCCGGATAAGGCGTTTATGCCGCATCCGGCAGTTATGCGCAGATGCATGATGCGACGCTGGCGCGTCTTATAATGCCTACACACGTGTGCCTGAACGGTAGGCCGGATAAGGCGGCAACGACGCACCCGGCCGTTGCGCCCCAATCCTGGCTaatg                                                                                                                                                                                                            >  2:299911/1‑271 (MQ=37)
                                                                        gACGTAACTGATTTTAGTGCCGGGGAGCGTATCATTTCATAATTGTGCGCAGATGCCTGATGCGACGCTAGCGCGTCTTATCATGCCTACAAACTTGTGCCGGATCGGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAGTTATGCGCAGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACAAACTTGTGCCTGAACGGTAGGACGGATAAGGCGTTCACGCCGCATCCGGCAGTTGTGCACCAATGCTGGATACTTACATCAAc                                                                                                                                                                                                    <  1:642434/275‑1 (MQ=255)
                                                                                ggcttttcatGCCGGGGAGCGTATACTTTCATGATGGTGCGCAGATGCCTGGTGCGGCGCTAGCGCGTCTTATCATGCCTACAACCTTGTGCCGGATCGGTAGGCCGGATAAGGCGTTTATGCCGCATCCGGCAGTTATGCGCAGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACAAACTTGTGCCTGAACGGTAGGACGGATAAGGCGTTCACGCCGCATCCGGCAGTTGTGCACCAATGCTGGATACTTACATCAACGCCCCCa                                                                                                                                                                                             <  2:520894/265‑1 (MQ=37)
                                                                                  aTTTTAGTGCCGGGGAGCGTATCATTTCATAATTGTGCGCAGATGCCTGATGCGACGCTAGCGCGTCTTATCATGCCTACAAACTTGTGCCGGATCGGTAGGCCGGATAAGGCGTTTATGCCGCATCCGGCAGTTATGCGCAGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACAAACTTGTGCCTGAACGGTAGGACGGATAAGGCGTTCACGCCGCATCCGGCAGTTGTGCACCAATGCTGGATACTTACATCAACGCCCCCATaa                                                                                                                                                                                          <  1:304808/275‑1 (MQ=255)
                                                                                           ccGGGGAGCGTATCAGTTCATAATGGTGCGCAGATGCCTGATGCGACGCTAGCGCGTCTTATCATGCCTACAAACTTGTGCCGGATCGGTAGGCCGGATAAGGCGTTTATGCCGCATCCGGCAGTTATGCGCAGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACAAACTTGTGCCNGAACGTTAGGACGGATAAGGCGTTCCCGCCGCATCCGGCAGTTGTGCACCAATGCTGGATACTTACATCAACGcc                                                                                                                                                                                                 >  1:520903/1‑259 (MQ=255)
                                                                                                              attaTTGTGCGCAGATGCCTGATGCGACGCTAGCGCGTCTTATCATGCCTACAAACTTGTGCCGGATCGGTAGGCCGGATAAGGAGTTTATGCCGCATCCGGCAGTTATGCGCAGATGCCTGATGCGACGCTGGCGCGTCTTAGCATGCCTACAAACTTGTGCCTGAACGGTAGGACGGATAAGGCGTTCACGCCGCATCCGGCAGTTGTGCACCAATGCTGGATACTTACATCAACGCCCCCATAATCTTCAACTCCAGCTCATCCGGCACtt                                                                                                                                                               <  1:277358/271‑1 (MQ=255)
                                                                                                              atcaTTGTGTGCAGATGCCTGATGCGTCGCTAGCGCGTCTTATCATGCCTACAAACTTGTGCCGGATCGGTAGGCCGGATAAGGCGTTTATGCCGCATCCGGCAGTTATGCGCAGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACAAACTTGTGCCTGAACGGTAGGACGGATAAGGCGTTCACGCCGCATCCGGCAGTTGTGCACCAATGCTGGATACTTACATCAACGCCCCCATAATCTTCAACTCCAGCTCATCCGGCACtt                                                                                                                                                               <  1:310190/271‑1 (MQ=255)
                                                                                                                                           cTAGCGCGTCTTATCATGCCTACAAACTTGGGCCGGATCGGTAGGCCGGATAAGGCGTTTATGCCGCATCCGGCAGTTATGCGCAGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACAAACTTGTGCCTGAACGGTAGGACGGATAAGGCGTTCACGCCGCATCCGGCAGTTGTGCACCAATGCTGGATACTTACATCAACGCCCCCATAATCTTCAACTCCAGCTCATCCGGCACTTCGTTATACGACAGCACATGCAGCCCCg                                                                                                                                    <  1:105130/272‑1 (MQ=255)
                                                                                                                                                             ccTACAAACTTGTGCCGGATCGGTAGGCCGGATAAGGCGTTTATGCCGCATCCGGCAGTTATGCGCAGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACAAACTTGTGCCTGAACGGTAGGACGGATAAGGCGTTCAcgc                                                                                                                                                                                                                                               >  1:596805/1‑147 (MQ=255)
                                                                                                                                                             ccTACAAACTTGTGCCGGATCGGTAGGCCGGATAAGGCGTTTATGCCGCATCCGGCAGTTATGCGCAGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACAAACTTGTGCCTGAACGGTAGGACGGATAAGGCGTTCACGCCGCATCCGGCAGTtgt                                                                                                                                                                                                                               >  1:475098/1‑163 (MQ=255)
                                                                                                                                                               tACAAACTTGTGCCGGATCGGTAGGCCGGATAAGGCGTTTATGCCGCATCCGGCAGTTATGCGCAGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACAAACTTGTGCCTGAACGGTAGGACGGATAAGGCg                                                                                                                                                                                                                                                      >  1:499715/1‑138 (MQ=255)
                                                                                                                                                                 cAAACTTGTGCCGGATCGGTAGGCCGGATAAGGCGTTTATGCCGCATCCGGCAGTTATGCGCAGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACAAACTTGTGCCTGAACGGTAGGACGGATAAGGCGTTCACGCCGCa                                                                                                                                                                                                                                           <  2:596792/147‑1 (MQ=255)
                                                                                                                                                                 cAAACTTGTGCCGGATCGGTAGGCCGGATAAGGCGTTTATGCCGCATCCGGCAGTTATGCGCAGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACAAACTTGTGCCTGAACGGTAGGACGGATAAGGCGTTCACGCCGCATCCGGCAGTTGTGCAc                                                                                                                                                                                                                           <  2:475090/163‑1 (MQ=255)
                                                                                                                                                                   aaCTTGTGCCGGATCGGTAGGCCGGATAAGGCGTTTATGCCGCATCCGGCAGTTATGCGCAGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACAAACTTGTGCCTGAACGGTAGGACGGATAAGGCGTTCa                                                                                                                                                                                                                                                  <  2:499706/138‑1 (MQ=255)
                                                                                                                                                                              gATCGGTAGGCCGGATAAGGCGTTTATGCCGCATCCGGCAGTTATGCGCAGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACAAACTTGTGCCTGAACGGTAGGACGGATAAGGCGTTCACGCCGCATCCGGCAGTTGTGCACCAATGCTGGATACTTACATCAACGCCCGCATAATCTTCAACTCCAGCTCATCCGGCACTTCGTTATACGACAGCACATGCAGCCCCGGCGCAAACAACCTGGCATAACGCGCCAGCAAAGGgcg                                                                                               >  1:405672/1‑274 (MQ=255)
                                                                                                                                                                                                  cGTTTATGCCGCATCCGGCAGTTATGCGCAGATGCCTGATGCGACGCGGGCGCGTCTTATCATGCCTACAAACTTGTGCCTGAACGGTAGGACGGATAAGGCGTTCACGCCGCATCCGGCAGTTGTGCACCAATGCTGGATACTTACATCAATGCCCCCATAATCTTCAACTCCAGCTCATCCGGCACTTCGTTAGACGACAGCACATGCAGCCCCGGCGCAAACAACCTTACATAACGCGCCAGCAAAGAGCACAGCTGCGAAGGCACcagca                                                                           >  1:703376/1‑274 (MQ=255)
                                                                                                                                                                                                                   gcattgaTGCGCAGTTGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACAAACTTGTGCCTGAACGGTAGGACGGATAAGGCGTTCACGCCGCATCCGGCAGTTGTGCACCAATGCTGGATACTTACTTCAACGCCCCCATAATCTTCAACTCCAGCTCATCCGGCACTTCGTTATACGACAGCACATGCAGCCCCGGCGCAAACAACCTTGCATAACGCGCCAGCAAAGGGCGCAGCTGCGGTGGCACCAGCAGCACCGGGTCTTTCCccg                                                         <  1:299916/269‑1 (MQ=255)
                                                                                                                                                                                                                       ttATGCGCAGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACAAACTTGTGCCTGAACGGTAGGACGGATAAGGCGTTCACGCCGCATCCGGCAGTTGTGCACCAATGCTGGATACTTACATCAACGCCCCCATAATCTTCAACTCCAGCTCATCCGGCACTTCGTTATACGACAGCACATGCAGCCCCGGCGCa                                                                                                                               >  1:142017/1‑201 (MQ=255)
                                                                                                                                                                                                                           gcgcAGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACAAACTTGTGCCTGAACGGTAGGACGGATAAGGCGTTCACGCCGCATCCGGCAGTTGTGCACCAATGCTGGATACTTACATCAACGCCCCCATAATCTTCAACTCCAGCTCATCCGGCACTTCGTTATACGACAGCACATGCAGCCCCGGCGCAaaca                                                                                                                           <  2:142016/201‑1 (MQ=255)
                                                                                                                                                                                                                                 aTGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACAAACTTGTGCCTGAACGGt                                                                                                                                                                                                                                                                     >  2:339457/1‑57 (MQ=255)
                                                                                                                                                                                                                                     cTGATGCGACGCTGGCGCGTCTTATCATGCCTACAAACTTGTGCCTGAACGGTAGGa                                                                                                                                                                                                                                                                 <  1:339463/57‑1 (MQ=255)
                                                                                                                                                                                                                                           cGACGCTGGCGCGTCTTATCATGCCTACAAACTTGTGCCTGAACGGTAGGACGGATAAGGCGTTCACGCCGCATCCGGCAGTTGTGCACCAATGCTg                                                                                                                                                                                                                   >  1:471977/1‑97 (MQ=255)
                                                                                                                                                                                                                                               gCTGGCGCGTCTTATCATGCCTACAAACTTGTGCCTGAACGGTAGGACGGATAAGGCGTTCACGCCGCATCCGGCAGTTGTGCACCAATGCTGGATa                                                                                                                                                                                                               <  2:471969/97‑1 (MQ=255)
                                                                                                                                                                                                                                                  ggCGCGTCTTATCATGCCTACAAACTTGTGCCTGAACGGTAGGACGGATAAGGCGTTCACGCCGCATCCGGCAGTTGTGCAcc                                                                                                                                                                                                                          >  1:90006/1‑83 (MQ=255)
                                                                                                                                                                                                                                                  ggCGCGTCTTATCATGCCTACAAACTTGTGCCTGAACGGTAGGACGGATAAGGCGTTCACGCCGCATCCGGCAGTTGTGCACCAAGGCTGGATACTTACAGCAACGCCCCCATAATCGTCAACTCCAGCTCATCCGGCACTTCGTTATACGACAGCACATGCAGCCCCGGCGCAAACAACCTTGCATAACGCGCCAGCAAAGGGCGCAGCTGCGGTGGCACCAGCAGCACCGGGTCTGTCCCCGCCGCTTTCATCTGCTCTGTCACCTGTGGCa                           <  1:134588/274‑1 (MQ=255)
                                                                                                                                                                                                                                                    cgcgTCTTATTATGCCTACAAACTTGTGCCTGAACGGTATGACGGATAAGGCGTTCACGCCGCATCCGGCAGTTGTGCACCAATGCTGGATACTTACATCAACGCCCCCATAATCTTCAACTCCAGCTCATCCGGCACTTCGTTATACGACAGCACATGCAGCCCCGGCGCAAACAACCTTGCATAACGCGCCAGCAAAGGGCGCAGCTGCGGTGGCACCAGCAGCACCGGGTCTTTCCCCGCCGCTTTCATCTGCTCTTTCACCTGTGGCATCg                        <  1:441025/275‑1 (MQ=255)
                                                                                                                                                                                                                                                      cgTCTTATCATGCCTACAAACTTGTGCCTGAACGGTAGGACGGATAAGGCGTTCACGCCGCATCCGGCAGTTGTGCACCAATg                                                                                                                                                                                                                      <  2:90005/83‑1 (MQ=255)
                                                                                                                                                                                                                                                        tCTTATCATGCCTACAAACTTGTGCCTGGAACGGTAGGACGGATAAGGCGTTCACGCCGCATCCGGCAGTTGTGTACCAATGCTGGATACTTACATCAACGCCCCCATAATCTTCAACTCCAGCTCATCCGGCACTTCGTtat                                                                                                                                                         >  2:466596/1‑143 (MQ=255)
                                                                                                                                                                                                                                                        tCTTATCATGCCTACAAACTTGTGCCTGGAACGGTAGGACGGATAAGGCGTTCAAGCCGCATCCGGCAGTTGTGCACCAATGCTGGATACTTACCTCACCGCCCCCATAATCTTCAACTCCAGCTCATCCGTCACTTCGTtat                                                                                                                                                         >  2:710976/1‑143 (MQ=37)
                                                                                                                                                                                                                                                            aTCATGCCTACAAACTTGTGCCTGGAACGGTAGGACGGATAAGGCGTTCACGCCGCATCCGGCAGTTGTGTACCAATGCTGGATACTTACATCAACGCCCCCATAATCTTCAACTCCAGCTCATCCGGCACTTCGTTATACGa                                                                                                                                                     <  1:466603/143‑1 (MQ=255)
                                                                                                                                                                                                                                                            aTCATGCCTACAAACTTGTGCCTGGAACGGTAGGACGGATAAGGCGTTCACGCCGCATCCGGCAGTTGTGCACCAATGCTGGATACTTACATCAACGCCCCCATAATCTTCAACTCCAGCTCATCCGGCACTTCGTTATACGa                                                                                                                                                     <  1:710990/143‑1 (MQ=255)
                                                                                                                                                                                                                                                                  ccTACAAACTTGTGCCTGAACGGTAGGACGGAAAAGCCGTCCACGCCGCATCCGGCAGTTGTGCACCAATGCTGGATACTTACATCAAAGCCCCCATAATCTTCAACTCCAGCTCATCCGGCACTTCGTTATACGACAGCAAATGCAGCCCCGGCGCAAACAACCTTGCCTAACGCGCCAGCAAAGGGCGCAGCGGCGGCGGCAACAGCAGCACAGGGTCTTTCCCCACCGCTTTAATCAGCCCTTGAACCCGTGGCATCGAACTCTGGTACggg          >  2:670810/1‑275 (MQ=255)
                                                                                                                                                                                                                                                                   cTACAAACTTGTGCCTGAACGGTAGGTCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGTGCACCAATGCTGGATACTTACATCAACGCCCCCATAATCTTCAACCCCAGCTCATCCGGCACTTCGTTATACGACAGCACATGCAGCCCCGGCGCAAACAACCTTGCATAACGCGCCAGCAAAGGGCGCAGCTGCGGTGGCACCAGCAGCACCGGGTCTTTCCCCGCCGCTTTCATCTGCTCTTTCACCTGTGGCATCGTACTCTGGAACTGGt         <  1:199159/275‑1 (MQ=255)
                                                                                                                                                                                                                                                                      cAAACTTGTGCCTGAACGGTAGGACGGATAAGGCGTTCACGCCGCATCCGGCAGTTGTGCACCAATGCTGGATACTTACATCAACGCCCCCATCATCTTCAACTCCAGCTCATCCGGCACTTCGTTATACGACAGCACATGCAGCCCCGGCGCAAACAACCTTCCATAACGCGCCAGCCAAGGGCGCAGCTGCGGTGGCACCAGCAGCACCCGGTCTTTCCCCCCCGCTCTCACTCGCTCTCTCCCCTTGGGCATCGGActctcgcaccgtcttg      >  2:622393/1‑263 (MQ=255)
                                                                                                                                                                                                                                                                         aCTTGTGCCTGAACGGTAGGACGGATAAGGCGTTCACGCCGCATCCGGCAGTTGTGCACCAATg                                                                                                                                                                                                                      >  2:339952/1‑64 (MQ=255)
                                                                                                                                                                                                                                                                            tgtgCCTGAACGGTAGGACGGATAAGGCGTTCACGCCGCATCCGGCAGTTGTGCACCAATGCTGGATACTTACATCAGCGCCCCCATAATATTCAACTACAGCACATCCGGCACTTCGTTATACGACAGCACAAGCAGCCCCGGCGCAAACAACCTTGCCTAACGCGCCAGCAAAGGGAGCAGCTGCGCGGGCAACAACACCACAGAGTCTTGCCCCGCCACCTCCATCCGCTCTTTCACCTGTGCCACCGTACTCTGGaannnngcgcgnct  >  2:52945/1‑261 (MQ=255)
                                                                                                                                                                                                                                                                             gtgCCTGAACGGTAGGACGGATAAGGCGTTCACGCCGCATCCGGCAGTTGTGCACCAATGCTgg                                                                                                                                                                                                                  <  1:339958/64‑1 (MQ=255)
                                                                                                                                                                                                                                                                              tgCCTGAACGGTAGGACGGATAAGGCGTTCACGCCGCATCCGGCAGt                                                                                                                                                                                                                                  >  1:320644/1‑47 (MQ=255)
                                                                                                                                                                                                                                                                              tgCCTGAACGGTAGGACGGATAAGGCGTTCACGCCGCATCCGGCAGTTGTGCACCAATGCTGGATACTTACATCAACGCCCCCATAATCTTCAACTCCAGCTCATCCGGCACTTCGTTATACGACAGCCCATGCAGCCCCGGCGCAAACAACCTTGCATAACGCGCCAGCAAAGGGCGCAGCTGCGGTGGCACCAGCAGCACCGGGTCTTTCCCCGCCGCTTTc                                                 >  1:153913/1‑224 (MQ=255)
                                                                                                                                                                                                                                                                                  tGAACGGTAGGACGGATAAGGCGTTCACGCCGCGTCCGGCAGTtgtg                                                                                                                                                                                                                              <  2:320638/47‑1 (MQ=255)
                                                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                                          
GTGAGTGATTGGCCATTCTCAACGTTCCATCGCGATGTCGCGCGAGGGTTATATCCCATCGATTGGGCGGGGGACGTAACGGATTTTAGTGCCGGGGAGCGTATCATTTCATAATTGTGCGCAGATGCCTGATGCGACGCTAGCGCGTCTTATCATGCCTACAAACTTGTGCCGGATCGGTAGGCCGGATAAGGCGTTTATGCCGCATCCGGCAGTTATGCGCAGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACAAACTTGTGCCTGAACGGTAGGACGGATAAGGCGTTCACGCCGCATCCGGCAGTTGTGCACCAATGCTGGATACTTACATCAACGCCCCCATAATCTTCAACTCCAGCTCATCCGGCACTTCGTTATACGACAGCACATGCAGCCCCGGCGCAAACAACCTTGCATAACGCGCCAGCAAAGGGCGCAGCTGCGGTGGCACCAGCAGCACCGGGTCTTTCCCCGCCGCTTTCATCTGCTCTTTCACCTGTGGCATCGTACTCTGGAACTGGTTGAGCAT  >  CP009273/244508‑245048

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