Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,749,743 A→G 18.2% E171G (GAG→GGG)  sgbU → putative L‑xylulose 5‑phosphate 3‑epimerase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,749,7430AG18.2% 143.9 / 3.0 64E171G (GAG→GGG) sgbUputative L‑xylulose 5‑phosphate 3‑epimerase
Reads supporting (aligned to +/- strand):  ref base A (25/25);  new base G (5/8);  total (30/34)
Fisher's exact test for biased strand distribution p-value = 5.42e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGGGAAATCATGAGTAAAGCGATTCGCCTGGCGCGCGATCTCGGCATTCGCACCATCCAGCTGGCAGGTTATGACGTCTATTACGAAGACCACGACGAAGGCACCCGGCAACGTTTTGCTGAAGGGCTGGCGTGGGCAGTCGAACAGGCGGCGGCATCACAAGTAATGCTGGCGGTGGAGATTATGGATACCGCGTTTATGAACTCCATCAGCAAATGGAAAAAATGGGACGAGATGCTCGCCTCACCGTGGTTCACCGTCTACCCGGACGTCGGCAACCTCAGCGCCTGGGGCAATGATGTTCCCGCCGAACTGAAACTGGGCATTGACCGTATCGCCGCGATCCACCTGAAAGACACCCAGCCAGTCACCGGGCAAAGTCCCGGACAGTTCCGCGACGTGCCGTTTGGCGAAGGCTGCGTCGATTTCGTTGGCATCTTCAAAACGCTGCATAAACTGAACTATC  >  NC_000913/3749511‑3749976
                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                           
ttaagagccaggAGTTAAGCTATTCGCCTGGCTCCCGTTCTCGGCATTCGCCCCACCCAGCCTGCAGGTTAGGACGTCTATTACGAAGACCACGACGAAGGCACCCGGCAACGTTTTGCTGAAGGGCTGGCGTGGGCAGTCGAACAGGCGGCGGCATCACAAGTAATGCTGGCGGTGGAGATTATGGATACCGCGTTTATGAACTCCATCAGCAAATGGAAAAAATGGGACGAGATGCt                                                                                                                                                                                                                                      <  2:504293/228‑1 (MQ=255)
cgGGAAATCATGGGTAAAGCGTTTCGCCTGGCGTGCGATCTCGGCATTCGCACCCTCCAGGTGGCAGGTAATGACGTCTATTACGAAGACCACGACGAAGGCACCCGGCAACGTTTTGCTGAAGGGCTGGCGTGGGCAGTCGAACAGGCGGCGGCATCACAAGTAATGCTGGCGGTGGAGATTATGGATACCGCGTTTATGAACTCCATCAGCAAATGGAAAAAATGGGACGAGATGCt                                                                                                                                                                                                                                      <  2:956018/239‑1 (MQ=255)
       tCATGAGTAAAGCGATTCGCCTGGCGGGCGATCTCGGCATTCGCACCATCCAGCTGGCAGGTTATGACGTCTATTACGAAGACCACGACGAAGGCACCCGGCAACGTTTTGCTGAAGGGCTGGCGTGGGCAGTCGAACAGGCGGCGGCATCACAAGTAATGCTGGCGGTGGAGATTATGGATACCGCGTTTATGAACTCCATCAGCAAATGGAAAAAATGGGACGAGATGCTCGCCTCa                                                                                                                                                                                                                               <  1:92780/239‑1 (MQ=255)
       tCATGAGTAAAGCGATTCGCCTGGCGCGCGATCTCGGCATTCGCACCATCCAGCTGGCAGGTTATGACGTCTATTACGAAGACCACGACGAAGGCACCCGGCAACGTTTTGCTGAAGGGCTGGCGTGGGCAGTCGAACAGGCGGCGGCATCACAAGTAATGCTGGCGGTGGAGATTATGGATACCGCGTTTATGAACTCCATCAGCAAATGGAAAAAATGGGACGAGATGCTCGCCTCa                                                                                                                                                                                                                               <  1:10043/239‑1 (MQ=255)
       tCATGAGTAAAGCGATTCGCCTGGCGCGCGATCTCGGCATTCGCACCATCCAGCTGGCAGGTTATGACGTCTATTACGAAGACCACGACGAAGGCACCCGGCAACGTTTTGCTGAAGGGCTGGCGTGGGCAGTCGAACAGGCGGCGGCATCACAAGTAATGCTGGCGGTGGAGATTATGGATACCGCGTTTATGAACTCCATCAGCAAATGGAAAAAATGGGACGAGATGCTCGCCTCa                                                                                                                                                                                                                               <  1:1076400/239‑1 (MQ=255)
           gAGTAAAGCGATTCGCCTGGCGCGCGATCTCGGCATTCGCACCATCCAGCTGGCAGGTTATGACGTCTATTACGAAGACCACGACGAAGGCACCCGGCAACGTTTTGCTGAAGGGCTGGCGTGGGCAGTCGAACAGGCGGCGGCATCACAAGTAATGCTGGCGGTGGAGATTATGGATACCGCGTTTATGAACTCCATCAGCAAATGGAAAAAATGGGACGAGATGCTCGCCTCACCGt                                                                                                                                                                                                                           <  1:301646/239‑1 (MQ=255)
               aaaGCGATTCGGCTGGCGCGCGATCCGGGCATTCGCACCCACCGGGTGGCAGGTTATGGCGTCTATTACGAAGACCACGACGAAGGCACCCGGCAACGTTTTGCTGAAGGGCTGGCGTGGGCAGTCGAACAGGCGGCGGCATCACAAGTAATGCTGGCGGTGGAGATTATGGATACCGCGTTTATGAACTCCATCAGCAAATGGAAAAAATGGGACGAGATGCTCGCCTCACCGTGGtt                                                                                                                                                                                                                       <  2:758587/239‑1 (MQ=255)
                 caggtttCGCCTGGCGGCGGATCTCGGCTATCCCCCCATCCAGCGGGCAGGTTGTGACGTCTAATACGAAGACCACGACGAAGGCACCCGGCAACGTTGTGCTGAAGGGCTGGCGTGGGCAGTCGAACAGGCGGCGGCATCACAAGTAATGCTGGCGGTGGAGATTATGGATACCGCGTTTATGAACTCCATCAGCAAATGGAAAAAATGGGACGAGATGCTCGCCTCACCGTGGtt                                                                                                                                                                                                                       <  2:112504/232‑1 (MQ=255)
                 aGCGATTCGCCTGGCGCGCGATCTCGGCATTCGCACCATCCAGCTGGCAGGTTATGACGTCTATTACGAAGACCACGACGAAGGCACCCGGCAACGTTTTGCTGAAGGGCTGGCGTGGGCAGTCGAACAGGCGGCGGCATCACAAGTAATGCTGGCGGTGGAGATTATGGATACCGCGGTTATGAAATCCATCAGCAAATGGAAAAAATGGGACGAGATGACCCCCTCACCGGGGTTc                                                                                                                                                                                                                      >  1:454549/1‑238 (MQ=255)
                   ctaTTCGCCTGGCGCGCGATCTCGGCATTCGCACCATCCAGCTGGCAGGTTATGACGTCTATTACGAAGACCACGGCGAAGGCACACGGCAACGTTTTGCTGAAGGGCTGGCGTGGGCAGTCGAACAGGCGGCGGCATCACAAGTAATGCTGGCGGTGGAGATTATGGATACCGCGTTTATGAACTCCATCAGCAAATGGAAAAAATGGGACGAGATGCTCGCCTCACCGTGGTTCAcc                                                                                                                                                                                                                   <  1:554354/237‑1 (MQ=255)
                     aTTCGCCTGGCGCGCGATCTCGGCATTCGCACCATCCAGCTGGCAGGTTATGACGTCTATTACGAAGACCACGACGAAGGCACCCGGCAACGTTTTGCTGAAGGGCTGGCGTGGGCAGTCGAACAGGCGGCGGCATCACAAGTAATGCTGGCGGTGGGGATTATGGAAACCGCGGTTATGAACTCCATCCGCAAATGGAAAAAATGGGACAGGTAGCTCGCCCCACCGCGGTTTACCg                                                                                                                                                                                                                  >  1:416810/1‑238 (MQ=255)
                      ttCGCCTGGCGCGCGATCTCGGCATTCGCACCATCCAGCTGGCAGGTTATGACGTCTATTACGAAGACCACGACGAAGGCACCCGGCAACGTTTTGCTGAAGGGCTGGCGTGGGCAGTCGAACAGGCGGCGGCATCACAAGTAATGCTGGCGGTGGAGATTATGGGTACCGCGTTTATGAACTCCATCAGCAAATGGAAAAAATGGGACGGGATGCTCGCCTCACCCTGGTTCACCGt                                                                                                                                                                                                                 >  1:758587/1‑238 (MQ=255)
                      ttCGCCTGGCGCGCGATCTCGGCATTCGCACCATCCAGCTGGCAGGTTATGACGTCTATTACGAAGACCACGACGAAGGCACCCGGCAACGTTTTGCTGAAGGGCTGGCGTGGGCAGTCGAACAGGCGGCGGCATCACAAGTAATGCTGGCGGTGGAGATTATGGATACCGCGTTTATGAACTCCCTTCGCCAATGGGAAAAATGGGGCGGGATGCTCGCCTCACCg                                                                                                                                                                                                                            >  1:504293/1‑227 (MQ=255)
                             ggCGCGCGATCTCGGCATTCGCACCATCCAGCTGGCAGGTTATGACGTCTATTACGAAGACCACGACGAAGGCACCCGGCAACGTTTTGCTGAAGGGCTGGCGTGGGCAGTCGAACAGGCGGCGGCATCACAAGTAATGCTGGCGGTGGAGATTATGGATACCGCGGTTATGAACTCCATCAGCAAATGGAAAAAATGGGACTAGATGCTCGCCTCACCGGGGTTCACCGTCTACCCgg                                                                                                                                                                                                         >  1:716646/1‑239 (MQ=255)
                               gggcggATCTTGGCATTCGCGCCCACCAGCGGGCAGGTTATGAGGTTTATTACGAAGGCCACGGAGAAGGCCCCCGGCAACGTTTGGCTGAAGGGCTGGCGTGGGCAGTCGAACAGGCGGCGGCATCACAAGGAATGCTGGCGGTGGAGATTATGGATACCGCGTTTATGAACTCCATCAGCAAATGGAAAAAATGGGACGAGATGCTCGCCTCACCGTGGTTCACCGTCTACCCGGa                                                                                                                                                                                                        <  2:204087/233‑1 (MQ=255)
                               cgcgcgATCTCGGCATTCGCACCATCCAGCTGGCAGGTTATGACGTCTATTACGAAGACCACGACGAAGGCACCCGGCAACGTTTTGCTGAAGGGCTGGCGTGGGCAGTCGAACAGGCGGCGGCATCACAAGTAATGCTGGCGGTGGGGATTATGGATACCCCGTTTTTGAACTCCCTCAGCAAAAGGAAAAAATGGGGCGGGGTGCTCGCCTAACCGTGGTTCACCGTCTaa                                                                                                                                                                                                             >  1:112504/1‑232 (MQ=255)
                                             aTTCGCCCCAGCCCGCTGGCAGGTTATGGCGTTTATTACGAAGGCCACGACGAAGGCCCCCGGCAACGTTTTGCTGAAGGGCTGGCGTGGGCAGTCGAACAGGCGGCGGCATCACAAGTAATGCTGGCGGTGGAGATTATGGATACCGCGTTTATGAACTCCATCAGCAAATGGAAAAAATGGGACGAGATGCTCGCCTCACCGTGGTTCACCGTCTACCCGGACGTCGGCAACCTCa                                                                                                                                                                                          <  1:506448/238‑1 (MQ=255)
                                                cGCATCTTCCATCTGGCAGGTTATGACGTCTATTACGAAGCCCACGACGAAGGCACCCGGCAACGTTTTGCTGAAGGGCTGGCGTGGGCAGTCGAACAGGCGGCGGCATCACAAGTAATGCTGGCGGTGGAGATTTGGATACCGCGTTTATGAACTCCATCAGCAAATGGAAAAAATGGGACGAGATGCTCGCCTCACCGTGGTTCACCGTCTACCCGGACGTCGGCAACCTCAGCGcc                                                                                                                                                                                     <  1:556917/239‑1 (MQ=255)
                                                cGCACCATCCAGCTGGCAGGTTATGACGTCTATTACGAAGACCACGACGAAGGCACCCGGCAACGTTTTGCTGAAGGGCTGGCGTGGGCAGTCGAACAGGCGGCGGCATCACAAGTAATGCTGGCGGTGGAGATTATGGATACCGCGTTTATGAACTCCATCAGCAAATGGAAAAAATGGGACGAGATGCTCGCCTCACCGTGGTTCCCCGTCTAACCGGGCGTCGGCAACCCCAgcg                                                                                                                                                                                       >  1:969125/1‑238 (MQ=255)
                                                      ctCCAGCTGGCAGGTTATGGCGTTTATTAGGAAGACCCCGAAGAAGGCACCCGGCAACGTTTTGCTGAAGGGCTGGCGTGGGCAGTCGAACAGGCGGCGGTATAACAAGTAATGCTGGCGGTGGAGATTAGGGATACCGCGTTTATGAACTCCATCAGCAAATGGAAAAAATGGGACGAGATGCTCGCCTCACCGTGGTGTGCCGTCTCCCCGTACGTCGGCAACCTCAGCGCCTgggg                                                                                                                                                                                <  2:762105/238‑1 (MQ=255)
                                                               gCAGGTTATGACGTCTATTACGAAGACCACGACGAAGGCACCCGGCAACGTTTTGCTGAAGGGCTGGCGTGGGCAGTCGAACAGGCGGCGGCATCACAAGTAATGCTGGCGGGGGAGATTATGGAAACCGCGTTTATGAACTCCATCAGCAAATGGAAAAAATGGGGCGAGATGCTCGCCTCACCGTGGTTCACCGTCTACCCGGACGTCGGCCACCTCAGCGCCTGGGGGAATTATGt                                                                                                                                                                       >  1:877409/1‑239 (MQ=255)
                                                                    ttATGACGTCTATTACGAAGACCACGACGAAGGCACCCGGCAACGTTTTGCTGAAGGGCTGGCGTGGGCAGTCGAACAGGCGGCGGCATCACAAGTAATGCTGGCGGTGGAGATTATGGATACCGCGTTTATGAACTCCATCAGCAAATGGAAAAAATGGGACGAGATGCTCGCCTCACCGTGGtt                                                                                                                                                                                                                       <  1:579593/186‑1 (MQ=255)
                                                                             cTATTACGAAGACCACGACGAAGGCACCCGGCAACGTTTTGCTGAAGGGCTGGCGTGGGCAGTCGAACAGGCGGCGGCATCACAAGTAATGCTGGCGGTGGAGATTATGGATACCGCGTTTATGAACTCCATCAGCAAATGGGAAAAATGGGACGAGATGCTCGCCTCACCGTGGTTCACCGTCTACCCCGACGTCCGCAACCTCAGCGCCCGGGGCCATGAGGTTCCCGCCGAACTGa                                                                                                                                                         >  1:1104732/1‑239 (MQ=255)
                                                                             cTATTACGAAGACCACGACGAAGGCACCCGGCAACGTTTTGCTGAAGGGCTGGCGTGGGCAGTCGAACAGGCGGCGGCATCACAAGTAATGCTGGCGGTGGAGATTATGGATACCGCGTTTATGAACTCCATCAGCAAATGGAAAAAATGGGACGAGATGCTCGCCTCACCGTGGTTCACCGTCTACCCGGACGCCGGCAACCCCAGCGCCCGGGGCAATGAAGTTTCCCCCCAACTGa                                                                                                                                                         >  1:125145/1‑239 (MQ=255)
                                                                              tATTACGAAGACCACGACGAAGGCACCCGGCAACGTTTTGCTGAAGGGCTGGCGTGGGCAGTCGAACAGGCGGCGGCATCACAAGTAATGCTGGCGGTGGAGATTATGGATACCGCGTTTATGAACTCCATCAGCAAATGGGAAAAATGGGACGAGATGCTCGCCTCACCGTGGTTCACCGTCTACCCGGACGCCGGCAACCTAAACCCCTGGGGCAATGATGTTCCCGccgaaccgaa                                                                                                                                                        >  1:1017102/1‑239 (MQ=255)
                                                                                ttccGAAGCCCCCGCCGAAGGCCCCGGGCAACTTTTTGCGGAGGGGGTGGGGGGGGCAGTCGAAAAGGCGGCGGCATCAAAAGTAATGCTGGCGGTGGAGATTAGGGATACCGCGTTTATGAACCCCATCAGAAAAGGGAAAAAAGGGGACGAGATGCTCGCCTCACCGTGGTTCACCGTCTCCCCGGACGTCGGCAACCTCAGCGCCTGGGGCAATGATGTTCCCGCCGAACTGAAAc                                                                                                                                                      <  2:1085662/236‑1 (MQ=255)
                                                                                ttACGAAGGCCACGACGAAGGCACCCGGCAACGTTTTGCTGAAGGGCTGGCGTGGGCAGTCGAACAGGCGGCGGCATCACAAGTAATGCTGGCGGTGGAGATTATGGATACCGCGTTTATGAACTCCATCAGCAAATGGAAAAAATGGGACGAGATGCTCGCCTCACCGTGGTTCACCGt                                                                                                                                                                                                                 <  1:755737/180‑1 (MQ=255)
                                                                                        aCCACGACGAAGGCACCCGGCAACGTTTTGCTGAAGGGCTGGCGTGGGCAGTCGAACAGGCGGCGGCATCACAAGTAATGCTGGCGGTGGAGATTATGGATACCGCGTTTATGAACTCCATCAGCAAATGGAAAAAATTGGACGAGATGCTCGCCTCAcc                                                                                                                                                                                                                             >  2:458777/1‑160 (MQ=255)
                                                                                        aCCACGACGAAGGCACCCGGCAACGTTTTGCTGAAGGGCTGGCGTGGGCAGTCGAACAGGCGGCGGCATCACAAGTAATGCTGGCGGTGGAGATTATGGATACCGCGTTTATGAACTCCATCAGCAAATGGAAAAAATGGGACGAGATGCTCGCCTCAcc                                                                                                                                                                                                                             <  1:458777/160‑1 (MQ=255)
                                                                                        aCCACGACGAAGGCACCCGGCAACGTTTTGCTGAAGGGCTGGCGTGGGCAGTCGAACAGGCGGCGGCATCACAAGTAATGCTGGCGGTGGAGATTATGGATACCGCGGTTATGAACTCCCTCAGCAAATGGGAAAAATGGGACGGGATGCTCGCCTCACCGTGGTTCACCGTCTACCCGGACGTCGGCAACCTCAGCGCCCGGGGCAATGAGGTTCCCGCCCAACTGaaaacgggacat                                                                                                                                              >  1:164578/1‑230 (MQ=255)
                                                                                            cgacgaAGGCACCCGGCAACGTTTTGCTGAAGGGCTGGCGTGGGCAGTCGAACAGGCGGCGGCATCACAAGTAATGCTGGCGGTGGAGATTATGGATACCGCGTTTATGAACTCCATCAGCAAATGGAAAAAATGGGACGAGATGCTCGCCTCACCGTGGTTCACCGTCTACCCGGACGTCGGCAACCTCAGCCCCTGGGGCAATGATGTTCCCGCCGAACTGCAAATGGGCCTTGAcc                                                                                                                                          >  1:959940/1‑239 (MQ=255)
                                                                                                 aaGGCCCCCGGCAACGTTTTGCTGAAGGGCTGGCGTGGGCAGTCGAACAGGCGGCGGCATCACAAGTAATGCTGGCGGTGGAGATTATGGATACCGCGTTTATGAACTCCATCAGCAAATGGAAAAAATGGGACGAGATGCTCGCCTCACCGTGGTTCACCGTCTACCCGGACGTCGGCAACCTCAGCGcc                                                                                                                                                                                     <  1:924851/191‑1 (MQ=255)
                                                                                                      ccccGGCAAGTTTTTGCGGAGGGGCTGGCGTGGGCAGTCGAACAGGCGGCGGCATCAAAAGTAAGGCTGGCGGTGGAGATTAGGGATACCGCGTTTATGAACTCCATCAGCAAATGGAAAAAATGGGACGAGATGCTCGCCTCACCGTGGTTCACCGTCTACCCGGACGTCGGCAACCTCAGCGCCTGGGGCAATGATGTTCCCGCCGAACTGAAACTGGGCATTGACCGTATcgccg                                                                                                                                 <  1:875233/237‑1 (MQ=255)
                                                                                                           ggccctttttGCTGGGGGGTCGGCGGGGGGAGGCGAACGGGGGGCGGCATCACAATTAATGGGGGGGGTGGAGATTATGGATACCGCGTTTATCAACTCGATCATCAAATGGAAAAAGTGGGACGGGATGCTCCCCTCACCGTGGTTCACCGTCTCCCCGGAGGTCGGCACCCTCAGCGCCTGGGGCAATGATGTTCCCGCCGAACTGAAACTGGGCATTGACCGTATCGCCGCGATcc                                                                                                                          <  2:292301/233‑1 (MQ=255)
                                                                                                                        gagggggtggggtGGGGGTTCTAACGGGGGGCGGCATCACAAGTAGTGCGGGGGGTGGAGATTATGGACACCGCGTTTATGAACTCCATCAGAAAAGGGAAAAAGTGGGCCGAGATGCTCGCCTCACCGTGGTTCACCGTCTACCCGGACGTCGGCAACCTCAGCGCCTGGGGCAATGATGTTCCCGCCGAACTGAAACTGGGCATTGACCGTATCGCCGCGATCCACCTGAAAGaca                                                                                                              <  2:959940/235‑1 (MQ=255)
                                                                                                                          aGGGCGGGGGTGGGCAGTCGAACAGGCGGCGGCATCACAAGTAATGCTGGCGGTGGAGATGATGGATACCGCGTTTATGACCTCCATCAGCAAATGGAAAAAATGGGACGAGATGCTCGCCTCACCGTGGTTCACCGTCTACCCGGACGTCGGCAACCTCAGCGCCTGGGGAAATGATGTTCCCGCCGAACTGAAACTGGGCATTGACCGTATCGCCGCGATCCACCTGAAAGACAcc                                                                                                            <  1:196616/238‑1 (MQ=255)
                                                                                                                                      ggCAGTCGAACAGGCGGCGGCATCACAAGTAATGCTGGCGGGGGGGATTGTGGGTAACGCGGTTATGAAATCCCTCAGGCAATGGGAAAAATGGGGCGAGATGGTCGCCTCACCGTGGTTCACCGTCTCCCCGCACGTAGGCCACCCCCGACCCCGGGGCAATGATGTTCCCCCCGAACCGCAAATGGGGATTGACCCTAACGCCCCCAGCCCCCCCAAAAACACCCCACAAGTCACCg                                                                                               >  2:959829/1‑239 (MQ=255)
                                                                                                                                           tggAACGGGGGGGGGTATCAAAAGTAAGGGGGGGGGGGGAAATTGGGGACACCGTTTTTAAGACCCCCAGCGAAAATGGGAAAAAAGGGGACGAGAGCCTCGCCTCACCGGGTTTCACCGTCTCCCCGGACGCCGGCAACCTCAGCCCCTGGGGCAATGATGTTCCCGCCGAACTGAAACTGGGCATTGACCGTATCGCCGCGATCCACCTGAAAGACACCCAGCCAGTCACCGGGCa                                                                                           <  2:125145/236‑1 (MQ=255)
                                                                                                                                            caaaCGGGCGGGGGCAACACAAGTTTGGCTGGGTGGGGAAATTTTGGAACCCGCGTTTATCAACCCCCTCAGAAATGGGAAAAAAGGGGCGGGTGTGCCCCCCTCACCGTGGTTCCCCGTCTACCCGGACCTCGGCACCCTCAGCGCCTGGGGCATTGATGTCCCCGCCGAACTGAAACTGGGCATTGACCGTATCGCCGCGATCCACCTGAAAGACACCCAGCCAGTCACCGGGCa                                                                                           <  2:969125/235‑1 (MQ=255)
                                                                                                                                                             cacaAGTAATGCTGGCGGTGGAGATTATGGATACCGCGGTTATGAACTCCATCAGCAAATGGAAAAAATGGGACGAGATGCTCGCCTCACCGTGGTTCACCGTCTACCCGGACGTCGGCAACCTCCCCGCCTGGGGGAATGATGTTCCCGCCGAACTGAAACTGAGCATCGACAGGATCGCGCGCATCACCCCTAAAAACACCCCGCCAAGAACCCGGGAAAGTCCCGGAAAGATcgcg                                                                        >  2:113943/1‑236 (MQ=255)
                                                                                                                                                                          ggCGGTGGAGATTATGGATACCGCGTTTATGAACTCCATCAGCAAATGGAAAAAATGGGACGAGATGCTCGCCTCACCGTGGTTCACCGTCTACCCGGACGTCGGc                                                                                                                                                                                                 <  2:612156/106‑1 (MQ=255)
                                                                                                                                                                          ggCGGTGGAGATTATGGATACCGCGTTTATGAACTCCATCAGCAAATGGAAAAAATGGGACGAGATGCTCGCCTCACCGTGGTTCACCGTCTACCCGGACGTCGGc                                                                                                                                                                                                 >  1:612156/1‑106 (MQ=255)
                                                                                                                                                                                ggAGATTATGGATCCCGCGTTTTTGCACGCCATCCGAAAATGGAAAAAATGGGCTGAGATGCTCGCCTCCCCGTGGTTCACCGTCTACCCGGACGTCGGCCACCCCAGCGCCTGGGGCAAGGATGTTCCCGCCGAACTGAAACTGGGCATTGACCGTAGCGCCGCGATTCCCCTGAAAGACACCCAGCCAGTCACCGGGCAAAGTCCCGGACAGTTCCGCGACGTGCCGTTTGGCGAAg                                                     <  2:416810/239‑1 (MQ=255)
                                                                                                                                                                                ggAGATTATGGATACCGCGTTTATGAACTCCATCAGCAAATGGAAAAAATGGGACGAGATGCTCGCCTCACCGTGGTTCACCGTCTACCCGGACGTCGGCAACCCCAGCGCCTGGGGCAATGATGTTCCCGCCGAACTGAAACTGGGCATTGACCGCATCGCCGCGATCCCCCTGAAAGACACCCAGCCAGGCACCGGGCAAAGTCCCGGACCGTTCCGCGACGTGCCGTTTGGCGAAg                                                     >  1:48481/1‑239 (MQ=255)
                                                                                                                                                                                    aTTAGGGATCCCGCGTTTATGAACTCCATCAGCAAATGGAAAAAATGGGACGAGATGCTCGCCTCACCGTGGTTCACCGTCTACCCGGACGTCGGCAACCTCAGCGCCTGGGGCAATGATGTTCCCGCCGAACTGAAACTGGGCATTGACCGTATCGCCGCGATCCACCTGAAAGACACCCAGCCAGTCACCGGGCAAAGTCCCGGACAGTTCCGCGACGTGCCGTTTGGCg                                                        <  1:1070266/232‑1 (MQ=255)
                                                                                                                                                                                         ggATACCGCGTTTATGAACTCCATCAGCAAATGGAAAAAATGGGACGAGATGCTCGCCTCACCGTGGTTCACCGTCTACCCGGACGTCGGCAACCTCAGCGCCTGGGGCAATGATGTTCCCGCCGAACCGAAACTGGGCATTGCCCGTATCGCCGCGATCCCCCTGAAAGACACCCAGCCAGCCACCGGGCAAAGCCCCGGCCAGTTCCGCACCGAGCCGTTTGGAGAAgggcgcgcgg                                            >  2:253286/1‑231 (MQ=255)
                                                                                                                                                                                           aTACCGCGTTTATGAACTCCATCAGCAAATGGAAAAAATGGGACGAGATGCTCGCCTCACCGTGGTTCACCGTCTACCCGGACGTCGGCAACCTCAGCGCCTGGGGCAATGATGTTCCCGCCGAACTGAAACTGGGCATTGACCGTATCGCCGCGATCCACCTGAAAGACACCCAGCCAGTCACCGGGCAAAGTCCCGGACAGTTCCGCGACGTGCCGGTTGGCGAAGGCTGCGTCGAt                                          >  1:482343/1‑239 (MQ=255)
                                                                                                                                                                                             aCCGCGTTTATGAACTCCATCAGCAAATGGAAAAAATGGGACGAGATGCTCGCCTCACCGTGGTTCACCGTCTACCCGGACGTCGGCAACCTCAGCGCCTGGGGCAATGATGTTCCCGCCGAACTGAAACTGGGCATTGACCGTATCGCCGCGATCCACCTGAAAGACACCCAGCCAGTCACCGGGCAAAGTCCCGGa                                                                                 >  1:852122/1‑198 (MQ=255)
                                                                                                                                                                                              ccGCGTTTATGAACTCCATCAGCAAATGGAAAAAATGGGACGAGATGCTCGCCTCACCGTGGTTCACCGTCTACCCGGACGTCGGCAACCTCAGCGCCTGGGGCAATGATGTTCCCGCCGAACTGAAACTGGGCCTTGACCGTATCGCCCCGACCCACCTGAAAGACACCCAGCCCGTCACCGGGACAAGTCCCGGACAGTTCCGCGACGTGCCGGTTGGCGAAGGCTgggcccggg                                         >  2:1070266/1‑229 (MQ=255)
                                                                                                                                                                                               cgcgTTTATGAACTCCATCAGCAAATGGAAAAAATGGGACGAGATGCTCGCCTCACCGTGGTTCACCGTCTACCCGGACGTCGGCAACCTCAGCGCCTGGGGCAATGATGTTCCCGCCGAACTGAAACTGGGCATTGACCGTATCGCCGCGATCCACCTGAAAGACACCCAGCCAGTCACCGGGCAAAGTCCCGGACAGTTCCGCGCCGTGCCGTTTGGCGAAGGCTGCGTCGATTTCg                                      >  1:501321/1‑239 (MQ=255)
                                                                                                                                                                                               ccggTTTAAGAACCCCAACAGCAAAGTGAAAAAAGGGGACGGGATCCCGCCCTCCCCGTGGTTCACCGTCTACCCGGACGTCGGCAACCTCAGCGCCTGGGGCAATGATTTTCCCGCCGAACTGAAACTGGGCATTGACCGTATCGCCGCGATCCACCTGAAAGACACCCAGCCAGTCACCGGGCAAAGTCCCGGACAGTTCCGCGACGTGCCGTTTGGCGAAGGCTGCGTCGATTTCg                                      <  1:91548/236‑1 (MQ=255)
                                                                                                                                                                                                 cgTTTATGAACTCCATCAGCAAATGGAAAAAATGGGACGAGATGCTCGCCTCACCGTGGTTCACCGTCTACCCGGACGTCGGCAACCTCAGCGCCTGGGGCAATGATGTTCCCGCCGAACTGAAACTGGGCATTGACCGTATCGCCGCGATCCACCTGAAAGACACCCAGCCAGTCACCGGGCAAAGCCCCGGACAGTTCCGCGACGTGCCGTTTGGCGAAgg                                                    >  1:824637/1‑223 (MQ=255)
                                                                                                                                                                                                           ccccAACAGAAAATGGAAAAAATGGGAAGGTGTGCTCGCCTCACCGTGGTTCACGGTCTACCCGGGCGTCGGCAACCTCAGCGCCGGGGGCAATGATGTTCCCGCCGAACTGAAACTGGGCATTGACCGTATCGCCGCGATCCACCTGAAAGACACCCAGCCAGTCACCGGGCAAAGTCCCGGACAGTTCCGCGACGTGCCGTTTGGCGAAGGCTGCGTCGATTTCGTTGGCATCTTCa                          <  1:1069035/237‑1 (MQ=255)
                                                                                                                                                                                                           cTCCATCAGCAAATGGAAAAAATGGGACGAGATGCTCGCCTCACCGTGGTTCACCGTCTACCCGGACGTCGGCAACCTCAGCGCCTGGGGCAATGATGTTCCCGCCGAACTGAAACTGGGCATTGACCGTATCGCCGCGATCCACCTGAAAGACACCCAGCCAGTCACCGGGCAAAGTCCCGGACAGTTCCGCGACGTGCCGTTTGGCGAAGGCTGCGTCGATTTCGTTGGCATCTTCa                          >  1:1105253/1‑239 (MQ=255)
                                                                                                                                                                                                              cATCAGCAAATGGAAAAAATGGGACGAGATGCTCGCCTCACCGTGGTTCACCGTCTACCCGGACGTCGGCAACCTCAGCGCCTGGGGCAATGATGTTCCCGCCGAACTGAAACTGGGCATTGACCGTATCGCCGCGATCCACCTGAAAGACACCCAGCCAGTCACCGGGCAAAGTCCCGGACAGTTCCGCGACGTGCCGTTTGGCGAAGGCTGCGTCGATTTCGTTGGCATCTTCaaa                        >  1:564771/1‑238 (MQ=255)
                                                                                                                                                                                                               aTCAGCAAATGGAAAAAATGGGCCGGGAGGCTCGCCCCCCCGTGGTTCCCCGTCTACCCGGACGTCGGCAACCTCAGCGCCTGGGGCAATGATTTCCCCGCCGAACTGAAACTGGGCATTGACCGTATCGCCGCGATCCACCTGAAAGACACCCAGCCAGTCACCGGGCAAAGTCCCGGACAGTTCCGCGACGTGCCGTTTGGCGAAGGCTGCGTCGATTTCGTTGGCATCTTCAAAAc                      <  1:512456/239‑1 (MQ=255)
                                                                                                                                                                                                                 ctccccacatgAAAAAATGGGGACGGACTGCCCCCCCCCCCGGGGTCCGCCCTCCCCCCGGACGTGGGCACCCTCAGCGCTGGGGGCAATGTGGTCCCCGCCGAACTGAAACTGGGCATTGACCGTATCGCCGCGATCCACCTGAAAGACCCCCAGCCAGTCACCGGGCAAAGTCCCGGACAGTTCCGCGACGTGCCGTTTGGCGAAGGCTGCGTCGATTTCGTTGGCATCTTCAAAAc                      <  2:1079421/229‑1 (MQ=255)
                                                                                                                                                                                                                    cAAATGGAAAAAATGGGCCGGAATGCCCGCCCCACCGTGGTTCCCCGTCTACCGGGACGGCGGCAACCTCAGCGCCTGGGCAAATGATTTTCCCGCCGAACTGAAACTGGGCATTGACCGTATCGCCGCGATCCACCTGAAAGACACCCAGCCAGTCACCGGGCAAAGTCCCGGACAGTTCCGCGACGTGCCGTTTGGCGAAGGCTGCGTCGATTTCGTTGGCATCTTCAAAACGCTGc                 <  1:1070445/239‑1 (MQ=255)
                                                                                                                                                                                                                      aagagaaaaaaaTGGCGCGCGGTCCTCGCCTCCCGGTTGTTCCGCGTCTCCCCGGACGGCGGCACCCCCAGCGCCTGGGGCAATGATGCCCCCGCCGAACTGAAACTGGGCATTGACCGTATCGCCGCGATCCACCTGAAAGACACCCAGCCAGTCACCGGGCAAAGTCCCGGACAGTTCCGCGACGTGCCGTTTGGCGAAGGCTGCGTCGATTTCGTTGGCATCTTCAAAACGCTGCa                <  2:454549/233‑1 (MQ=255)
                                                                                                                                                                                                                      aaTGGAAAAAATGGGACGAGATGCTCGCCTCACCGTGGTTCACCGTCTACCCGGACGTCGGCAACCTCAGCGCCTGGGGCAATGATGTTCCCGCCGAACTGAAACTGGGCATTGACCGTATCGCCGCGATCCACCTGAAAGACACCCAGCCAGTCACCGGGCAAAGTCCCGGACAGTTCCGCGACGTGCCGTTTGGCGAAGGCTGCGTCGATTTCGTTGGCATCTTCAAAACGCTGCAt               >  1:1091935/1‑239 (MQ=255)
                                                                                                                                                                                                                        tGGAAAAAATGGGACGAGATGCTCGCCTCACCGTGGTTCACCGTCTACCCGGACGTCGGCAACCTCAGCGCCTGGGGCAATGATGTTCCCGCCGAACTGAAACTGGGCATTGACCGTATCGCCGCGATCCACCTGAAAGACACCCAGCCAGTCACCGGGCAAAGTCCCGGACAGTTCCGCGACGTGCCGTTTGGCGAAGGCTGCGTCGATTCCGTTGGCACCTTCAAAACGCTGCATaa             >  1:532623/1‑239 (MQ=255)
                                                                                                                                                                                                                           aaaaaaaGGGGACGGAATGCTCGCCTCCCCGTGTTTCCCGGTCTCCCCGGGCGGCGGCACCCTCAGCGCCTGGGGCAATGATGTTCCCGCCGAACTGAAACTGGGCATTGACCGTATCGCCGCGATCCACCTGAAAGACACCCAGCCAGTCACCGGGCAAAGTCCCGGACAGTTCCGCGACGTGCCGTTTGGCGAAGGCTGCGTCGATTTCGTTGGCATCTTCAAAACGCTGCATaaa            <  1:187870/237‑1 (MQ=255)
                                                                                                                                                                                                                                     ggACGAGATGCTCGCCTCACCGTGGTTCACCGTCTACCCGGACGTCGGCAACCTCAGCGCCTGGGGCAATGATGTTCCCGCCGAACTGAAACTGGGCATTGACCGTATCGCCGCGATCCACCTGAAAGACACCCAGCCAGTCACCGGGCAAAGTCCCGGACAGTTCCGCGACGTGCCGTTTGGCGAAGGCTGCGTCGATTTCGTTGGCATCTTCAAAACGCTGCCTAAACTGAACTAt   >  1:449911/1‑238 (MQ=255)
                                                                                                                                                                                                                                      gACGAGATGCTCGCCTCACCGTGGTTCACCGTCTACCCGGACGTCGGCAACCTCAGCGCCTGGGGCAATGATGTTCCCGCCGAACTGAAACTGGGCATTGACCGTATCGCCGCGATCCACCTGAAAGACACCCAGCCAGTCACCGGGCAAAGTCCCGGACAGTTCCGCGACGTGCCGTTTGGCGAAGGCTGCGTCGATTTCGTTGGCATCTTCAAAACGCTGCATAAACTGAACTATc  >  1:277048/1‑238 (MQ=255)
                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                           
CGGGAAATCATGAGTAAAGCGATTCGCCTGGCGCGCGATCTCGGCATTCGCACCATCCAGCTGGCAGGTTATGACGTCTATTACGAAGACCACGACGAAGGCACCCGGCAACGTTTTGCTGAAGGGCTGGCGTGGGCAGTCGAACAGGCGGCGGCATCACAAGTAATGCTGGCGGTGGAGATTATGGATACCGCGTTTATGAACTCCATCAGCAAATGGAAAAAATGGGACGAGATGCTCGCCTCACCGTGGTTCACCGTCTACCCGGACGTCGGCAACCTCAGCGCCTGGGGCAATGATGTTCCCGCCGAACTGAAACTGGGCATTGACCGTATCGCCGCGATCCACCTGAAAGACACCCAGCCAGTCACCGGGCAAAGTCCCGGACAGTTCCGCGACGTGCCGTTTGGCGAAGGCTGCGTCGATTTCGTTGGCATCTTCAAAACGCTGCATAAACTGAACTATC  >  NC_000913/3749511‑3749976

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 26 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: