Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,858,834	G→A	5.0%	A443T (GCT→ACT)	PP_5134 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,858,834	0	G	A	5.0%	236.9 / 2.2	79	A443T (GCT→ACT)	PP_5134	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (27/48);  new base A (2/2);  total (29/50)
Fisher's exact test for biased strand distribution p-value = 6.22e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.46e-01

CCGCGGGCATGCGGCGCAAGACAGCATTCGCCTGGGTGAGCTGGATGAAAAGGCCTGGCGTGAGGGCAAGTTGCCGGAGGCGTTGGTCAAAGCGCGGGCGCACCCGCAGTGGTCCACCCTGAATCCCGGGG‑CCAAGCTGAATCAGGCTGAGCGGCAGGAGCAGTTGGAGGCGCTGGTTTCACAGATTTGAGCCAGCAGGGTTGGCCTCTTCGCGGGTAAACCCGCTCCTACGGGCGGCCT‑TTGTACAACCACCGATGTGCGTCCACCACCAATCCGGCGTAGGAGCGGGTTTACC  >  NC_002947/5858689‑5858983                                                                                                                                                   |                                                                                                                                                       ccGCGGGCATGCGGCGCAAGACAGCATTCGCCTGGGTGAGCTGGATGAAAAGGCCTGGCGTGAGGGCAAGTTGCCGGAGGCGTTGGTCAAAGCGCGGGCGCACCCGCAGTGGTCCACCCTGAATCCCGGGG‑CCAAGCTGAATCAGGCTGa                                                                                                                                                    >  1:698249/1‑150 (MQ=255) ccGCGGGCATGCGGCGCAAGACAGCATTCGCCTGGGTGAGCTGGATGAAAAGGCCTGGCGTGAGGGCAAGTTGCCGGAGGCGTTGGTCAAAGCGCGGGCGCACCCGCAGTGGTCCACCCTGAATCCCGGGG‑CCAAGCTGAATCAGGCTGa                                                                                                                                                    >  1:1679197/1‑150 (MQ=255)   gcgGGCATGCGGCGCAAGACAGCATTCGCCTGGGTGAGCTGGATGAAAAGGCCTGGCGTGAGGGCAAGTTGCCGGAGGCGTTGGTCAAAGCGCGGGCGCACCCGCAGTGGTCCACCCTGAATCCCGGGG‑CCAAGCTGAATCAGGCTGAgg                                                                                                                                                  >  2:105251/1‑149 (MQ=255)       gCATGCGGCGCAAGACAGCATTCGCCTGGGTGAGCTGGATGAAAAGGCCTGGCGTGAGGGCAAGTTGCCGGAGGCGTTGGTCAAAGCGCGGGCGCACCCGCAGTGGTCCACCCTGAATCCCGGGG‑CCAAGCTGAATCAGGCTGAGCGGCa                                                                                                                                              >  1:1036211/1‑150 (MQ=255)          tGCGGCGCAAGACAGCATTCGCCTGGGTGAGCTGGATGAAAAGGCCTGGCGTGAGGGCAAGTTGCCGGAGGCGTTGGTCAAAGCGCGGGCGCACCCGCAGTGGTCCACCCTGAATCCCGGGG‑CCAAGCTGAATCAGGCTGAGCGGCAGGa                                                                                                                                           <  2:752070/150‑1 (MQ=255)             ggcgCAAGACAGCATTCGCCTGGGTGAGCTGGATGAAAAGGCCTGGCGTGAGGGCAAGTTGCCGGAGGCGTTGGTCAAAGCGCGGGCGCACCCGCAGTGGTCCACCCTGAATCCCGGGG‑CCAAGCTGAATCAGGCTGCGCGGCAGGAGCa                                                                                                                                        <  1:964930/150‑1 (MQ=255)               cgcAAGACAGCATTCGCCTGGGTGAGCTGGATGAAAAGGCCTGGCGTGAGGGCAAGTTGCCGGAGGCGTTGGTCAAAGCGCGGGCGCACCCGCAGTGGTCCACCCTGAATCCCGGGG‑CCAAGCTGAATCAGGCTGAGCGGCAGGAGCAGt                                                                                                                                      <  2:698249/150‑1 (MQ=255)                gcAAGACAGCATTCGCCTGGGTGAGCTGGATGAAAAGGCCTGGCGTGAGGGCAAGTTGCCGGAGGCGTTGGTCAAAGCGCGGGCGCACCCGCAGTGGTCCACCCTGAATCCCGGGG‑CCAAGCTGAATCAGGCTGAGCGGCAGGAGCAGtt                                                                                                                                     >  2:432028/1‑150 (MQ=255)                    gACAGCATTCGCCTGGGTGAGCTGGATGAAAAGGCCTGGCGTGAGGGCAAGTTGCCGGAGGCGTTGGTCAAAGCGCGGGCGCACCCGCAGTGGTCCACCCTGAATCCCGGGG‑CCAAGCTGAATCAGGCTGAGCGGCAGGAGCAGTTGGAg                                                                                                                                 >  1:1412947/1‑150 (MQ=255)                           ttCGCCTGGGTGAGCTGGATGAAAAGGCCTGGCGTGAGGGCAAGTTGCCGGAGGCGTTGGTCAAAGCGCGGGCGCACCCGCAGTGGTCCACCCTGAATCCCGGGG‑CCAAGCTGAATCAGGCTGAGCGGCAGGAGCAGTTGGAGGCGCTgg                                                                                                                          >  2:1714149/1‑150 (MQ=255)                           ttCGCCTGGGTGAGCTGGATGAAAAGGCCTGGCGTGAGGGCAAGTTGCCGGAGGCGTTGGTCAAAGCGCGGGCGCACCCGCAGTGGTCCACCCTGAATCCCGGGG‑CCAAGCTGAATCAGGCTGAGCGGCAGGAGCAGTTGGAGGCGCTgg                                                                                                                          <  1:1728376/150‑1 (MQ=255)                                     tGAGCTGGATGAAAAGGCCTGGCGTGAGGGCAAGTTGCCGGAGGCGTTGGTCAAAGCGCGGGCGCACCCGCAGTGGTCCACCCTGAATCCCGGGG‑CCAAGCTGAATCAGGCTGAGCGGCAGGAGCAGTTGGAGGCGCTGGTTTCACAGAt                                                                                                                <  1:977304/150‑1 (MQ=255)                                     tGAGCTGGATGAAAAGGCCTGGAGTGGCCGGTTCTCGCCGGAGGCGTTGGTCAAAGCGCGGGCGCACCCGCAGTGGTCCACCCTGAATCCCGGGG‑CCAAGCTGAATCAGGCTGAGCGGCAGGAGCAGTTGGGGGCGCTGGTTTCACAGAt                                                                                                                <  2:1679197/150‑1 (MQ=255)                                        gCTGGATGAAAAGGCCTGGCGTGAGGGCAAGTTGCCGGAGGCGTTGGTCAAAGCGCGGGCGCACCCGCAGTGGTCCACCCTGAATCCCGGGG‑CCAAGCTGAATCCGGCCACGCAGCATGAGCAGTTGGAGGCGCTGGTTTCACAGATTTg                                                                                                             >  2:1401477/1‑150 (MQ=255)                                         cTGGATGAAAAGGCCTGGCGTGAGGGCAAGTTGCCGGAGGCGTTGGTCAAAGCGCGGGCGCACCCGCAGTGGTCCACCCTGAATCCCGGGG‑CCAAGCTGAATCAGGCTGAGCGGCAGGAGCAGTTGGAGGCGCTGGTTTCACAGATTTGa                                                                                                            <  1:1761169/150‑1 (MQ=255)                                         cTGGATGAAAAGGCCGGGCGTGAGGGCAAGTTGCCGGAGGCGTTGGTCACAGAGCGGCTGCGCCAGCTCTGGTCTAGCTTGAATCCCGGGT‑CCACGCTGAATCAGGCTGAGCGGCAGGAGCAGTTGGAGGCGCTGGTTTCACAGATTTGa                                                                                                            <  1:432028/150‑1 (MQ=255)                                           ggATGAAAAGGCCTGGCGTGAGGGCAAGTTGCCGGAGGCGTTGGTCAAAGCGCGGGCGCACCCGCAGTGGTCCACCCTGAATCCCGGGG‑CCAAGCTGAATCAGGCTGAGCGGCAGGAGCAGTTGGAGGCGCTGGTTTCACAGATTTGAGc                                                                                                          >  2:741060/1‑150 (MQ=255)                                            gATGAAAAGGCCTGGCGTGAGGGCAAGTTGCCGGAGGCGTTGGTCAAAGCGCGGGCGCACCCGCAGTGGTCCACCCTGAATCCCGGGG‑CCAAGCTGAATCAGGCTGAGCGGCAGGAGCAGTTGGAGGCGCTGGTTTCACAGATTTGAGcc                                                                                                         >  2:295028/1‑150 (MQ=255)                                             aTGAAAAGGCCTGGCGTGAGGGCAAGTTGCCGGAGGCGTTGGTCAAAGCGCGGGCGCACCCGCAGTGGTCCACCCTGAATCCCGGGG‑CCAAGCTGAATCAGGCTGAGCGGCAGGAGCAGTTGGAGGCGCTGGTTTCACAGATTTGAGCca                                                                                                        >  2:1875772/1‑150 (MQ=255)                                                  aaGGCCTGGCGTGAGGGCAAGTTGCCGGAGGCGTTGGTCAAAGCGCGGGCGCATACGCAGTGGTCCACCCTGAATCCCGGGG‑CCAAGCTGAATCAGGCTGAGCGGCAGGAGCAGTTGGAGGCGCTGGTTTCACAGATTTGAGCCAGCAgg                                                                                                   >  2:1606648/1‑150 (MQ=255)                                                      ccTGGCGTGAGGGCAAGTTGCCGGAGGCGTTGGTCAAAGCGCGGGCGCACCCGCAGTGGTCCACCCTGAATCCCGGGG‑CCAAGCTGAATCAGGCTGAGCGGCAGGAGCAGTTGGAGGCGCTGGTTTCACAGATTTGAGCCAGCAGGGTTg                                                                                               <  1:295028/150‑1 (MQ=255)                                                      ccTGGCGTGAGGGCAAGTTGCCGGAGGCGTTGGTCAAAGCGCGGGCGCACCCGCAGTGGTCCACCCTGAATCCCGGGG‑CCAAGCTGAATCAGGCTGAGCGGCAGGAGCAGTTGGAGGCGCTGGTTTCACAGATTTGAGCCAGCAGGGTTg                                                                                               <  1:915880/150‑1 (MQ=255)                                                       cTGGCGTGAGGGCAAGTTGCCGGAGGCGTTGGTCAAAGCGCGGGCGCACCCGCAGTGGTCCACCCTGAATCCCGGGG‑CCAAGCTGAATCAGGCTGAGCGGCAGGAGCAGTTGGAGGCGCTGGTTTCACAGATTTGAGCCAGCAGGGTTgg                                                                                              >  1:223233/1‑150 (MQ=255)                                                        tGGCGTGAGGGCAAGTTGCCGGAGGCGTTGGTCAAAGCGCGGGCGCACCCGCAGTGGTCCACCCTGAATCCCGGGG‑CCAAGCTGAATCAGGCTGAGCGGCAGGAGCAGTTGGAGGCGCTGGTTTCACAGATTTGAGCCAGCAGGGTTGGc                                                                                             <  1:632418/150‑1 (MQ=255)                                                        tGGCGTGAGGGCAAGTTGCCGGAGGCGTTGGTCAAAGCGCGGGCGCACCCGCAGTGGTCCACCCTGAATCCCGGGG‑CCAAGCTGAATCAGGCTGAGCGGCAGGAGCAGTTGGAGGCGCTGGTTTCACAGATTTGAGCCAGCAGGGTTGGc                                                                                             <  2:715429/150‑1 (MQ=255)                                                        tGGCGTGAGGGCAAGTTGCCGGAGGCGTTGGTCAAAGCGCGGGCGCACCCGCAGTGGTCCACCCTGAATCCCGGGG‑CCAAGCTGAATCAGGCTGAGCGGCAGGAGCAGTTGGAGGCGCTGGTTTCACAGATTTGAGCCAGCAGGGTTGGc                                                                                             <  2:435893/150‑1 (MQ=255)                                                        tGGCGTGAGGGCAAGTTGCCGGAGGCGTTGGTCAAAGCGCGGGCGCACCCGCAGTGGTCCACCCTGAATCCCGGGG‑CCAAGCTGAATCAGGCTGAGCGGCAGGAGCAGTTGGAGGCGCTGGTTTCACAGATTTGAGCCAGCAGGGTTGGc                                                                                             >  1:534535/1‑150 (MQ=255)                                                        tGGCGTGAGGGCAAGTTGCCGGAGGCGTTGGTCAAAGCGCGGGCGCACCCGCAGTGGTCCACCCTGAATCCCGGGG‑CCAAGCTGAATCAGGCTGAGCGGCAGGAGCAGTTGGAGGCGCTGGTTTCACAGATTTGAGCCAGCAGGGTTGGc                                                                                             <  2:184206/150‑1 (MQ=255)                                                        tGGCGTGAGGGCAAGTTGCCGGAGGCGTTGGTCAAAGCGCGGGCGCACCCGAAGTGGTCCCCCCTGTCTCCCGGGG‑CCAAGCTGAATCTGGATGAGCGGCAGGAGGAGGTGGAGGCGCTGGTTTCACAGATTTGAGCCAGCAGGGTTGGc                                                                                             <  2:1655662/150‑1 (MQ=255)                                                                gggCAAGTTGCCGGAGGCGTTGGTCAAAGCGCGGGCGCACCCGCAGTGGTCCACCCTGAATCCCGGGG‑CCAAGCTGAATCAGGCTGAGCGGCAGGAGCAGTTGGAGGCGCTGGTTCCACAGATTTGAGCCAGCAGGGTTGGCCTCTTcgc                                                                                     <  1:1401477/150‑1 (MQ=255)                                                                 ggCAAGTTGCCGGAGGCGTTGGTCAAAGCGCGGGCGCACCCGCAGTGGTCCACCCTGAATCCCGGGG‑CCAAGCTGAATCAGGCTGAGCGGCAGGAGCAGTTGGAGGCGCTGGTTTCACAGATTTGAGCCAGCAGGGTTGGCCTCTTcgcg                                                                                    <  1:1182048/150‑1 (MQ=255)                                                                  gCAAGTTGCCGGAGGCGTTGGTCAAAGCGCGGGCGCACCCGCAGTGGTCCACCCTGAATCCCGGGG‑CCAAGCTGAATCAGGCTGAGCGGCAGGAGCAGTTGGAGGCGCTGGTTTCACAGATTTGAGCCAGCAGGGTTGGCCTCTTCGCgg                                                                                   <  2:145789/150‑1 (MQ=255)                                                                  gCAAGTTGCCGGAGGCGTTGGTCAAAGCGCGGGCGCACCCGCAGTGGTCCACCCTGAATCCCGGGG‑CCAAGCTGAATCAGGCTGAGCGGCAGGAGCAGTTGGAGGCGCTGGTTTCACAGATTTGAGCCAGCAGGGTTGGCCTCTTCGCgg                                                                                   <  2:1082761/150‑1 (MQ=255)                                                                  gCAAGTTGCCGGAGGCGTTGGTCAAAGCGCGGGCGCACCCGCAGTGGTCCACCCTGAATCCCGGGG‑CCAAGCTGAATCAGGCTGAGCGGCAGGAGCAGTTGGAGGCGCTGGTTTCACAGATTTGAGCCAGCAGGGTTGGCCTCTTCGCgg                                                                                   <  1:1668611/150‑1 (MQ=255)                                                                  gCAAGTTGCCGGAGGCGTTGGTCAAAGCGCGGGCGCACCCGCAGTGGTCCACCCTGAATCCCGGGG‑CCAAGCTGAATCAGGCTGAGCGGCAGGAGCAGTTGGAGGCGCTGGTTTCACAGATTTGAGCCAGCAGGGTTGGCCTCTTCGCgg                                                                                   <  1:1393511/150‑1 (MQ=255)                                                                  gCAAGTTGCCGGAGGCGTTGGTCAAAGCGCGGGCGCACCCGCAGTGGTCCACCCTGAAGCCCGGGG‑CCAAGCTGAATCAGGCTGAGCGGCAGGAGCAGTTGGAGGCGCTGGTTTCACAGATTTGAGCCAGCAGGGTTGGCCTCTTCCCgg                                                                                   <  1:274559/150‑1 (MQ=255)                                                                      gTTGCCGGAGGCGTTGGTCAAAGCGCGGGCGCACCCGCAGTGGTCCACCCTGAATCCCGGGG‑CCAAGCTGAATCAGGCTGAGCGGCAGGAGCAGTTGGAGGCGCTGGTTTCACAGATTTGAGCCAGCAGGGtt                                                                                                >  1:413113/1‑133 (MQ=255)                                                                      gTTGCCGGAGGCGTTGGTCAAAGCGCGGGCGCACCCGCAGTGGTCCACCCTGAATCCCGGGG‑CCAAGCTGAATCAGGCTGAGCGGCAGGAGCAGTTGGAGGCGCTGGTTTCACAGATTTGAGCCAGCAGGGtt                                                                                                <  2:413113/133‑1 (MQ=255)                                                                            ggAGGCGTTGGTCAAAGCGCGGGCGCACCCGCAGTGGTCCACCCTGAATCCCGGGG‑CCAAGCTGAATCAGGCTGAGCGGCAGGAGCAGTTGGAGGCGCTGGTTTCACAGATTTGAGCCAGCAGGGTTGGCCTCTTCGCGGGTAAACCCGc                                                                         >  2:1226376/1‑150 (MQ=255)                                                                            ggAGGCGTTGGTCAAAGCGCGCGCGCACCCGCAGTGGTCCACCCTGAATCCCGGGG‑CCAAGCTGAATCAGGCTGCGCGGCAGGAGCAGTTGGAGGCGCTGGTTTCACAGATTTGAGCCAGCAGGGTTGGCCTCTTCGCGGGTAAACCCGc                                                                         <  1:836269/150‑1 (MQ=255)                                                                             gAGGCGTTGGTCAAAGCGCGGGCGCACCCGCAGTGGTCCACCCTGAATCCCGGGG‑CCAAGCTGAATCAGGCTGAGCGGCAGGAGCa                                                                                                                                        <  2:597577/86‑1 (MQ=255)                                                                             gAGGCGTTGGTCAAAGCGCGGGCGCACCCGCAGTGGTCCACCCTGAATCCCGGGG‑CCAAGCTGAATCAGGCTGAGCGGCAGGAGCa                                                                                                                                        >  1:597577/1‑86 (MQ=255)                                                                              aGGCGTTGGTCAAAGCGCGGGCGCACCCGCAGTGGTCCACCCTGAATCCCGGGG‑CCAAGCTGAATCAGGCTGAGCGGCAGGAGCAGTTGGAGGCGCTGGTTTCACAGATTTGAGCCAGCAGGGTTGGCCTCTTCGCGGGTAAACCCGCTc                                                                       <  1:741060/150‑1 (MQ=255)                                                                              aGGCGTTGGTCAAAACGCGGCCGCACCCGCTGTAGTCTACCCTAAATCCCGAAGCCCCTGATTTGTCA‑GCGCAAGGGCAGGAGCAGTTGGAGGCGCTGGTTTCACAGATTTGAGCCAGCAGGGTTGGCCTCTTCGCGGGTAAACCCGCTc                                                                       <  1:974667/150‑1 (MQ=255)                                                                               ggCGTTGGTCAAAGCGCGGGCGCACCCGCAGTGGTCCACCCTGAATCCCGGGG‑CCAAGCTGAATCAGGCTGAGCGGCAGGAGCAGTTGGAGGCGCTGGTTTCACAGATTTGAGCCAGCAGGGTTGGCCTCTTCGCGGGTAAACCCGCTcc                                                                      <  1:1714149/150‑1 (MQ=255)                                                                                    tGGTCAAAGCGCGGGCGCACCCGCAGTGGTCCACCCTGAATCCCGGGG‑CCAAGCTGAATCAGGCTGAGCGGCAGGAGCAGTTGGAGGCGCTGGTTTCACAGATTTGAGCCAGCAGGGTTGGCCTCTTCGCGGGTAAACCCGCTCCTACgg                                                                 >  1:429458/1‑150 (MQ=255)                                                                                     ggTCAAAGCGCGGGCGCACACGCAGTGGTCCACCCTGAATCCCGGGG‑CCAAGCTGAATCAGGCTGAGCGGCAGGAGCAGTTGGAGGCGCTGGTTTCACAGATTTGAGCCAGCAGGGTTGGCCTCTTCGCGGGTAAACCCGCTCCTaaaaa                                                                <  1:1006032/150‑5 (MQ=255)                                                                                             cgcgGGCGCACCCGCAGTGGTCCACCCTGAATCCCGGGG‑CCAAGCTGAATCAGGCTGAGCGGCAGGAGCAGTTGGAGGCGCTGGTTTCACAGATTTGAGCCAGCAGGGTTGGCCTCTTCGCGGGTAAACCCGCTCCTACGGGCGGCCt‑tt                                                       <  1:1226376/150‑1 (MQ=255)                                                                                                   cgcACCCGCAGTGGTCCACCCTGAATCCCGGGG‑CCAGGCTGAATCAGGCTGAGCGGCAGGAGCAGTTGGAGGCGCTGGTTTCACAGATTTGAGCCAGCAGGGTTGGCCTCTTCGCGGGTAAACCCGCTCCTACGGGCGGCCT‑TTGTACaa                                                 <  1:956874/150‑1 (MQ=255)                                                                                                   cgcACCCGCAGTGGTCCACCCTGAATCCCGGGG‑CCAAGCTGAATCAGGCTGAGCGGCAGGAGCAGTTGGAGGCGCTGGTTTCACAGATTTGAGCCAGCAGGGTTGGCCTCTTCGCGGGTAAACCCGCTCCTACGGGCGGCCT‑TTGTACaa                                                 <  1:957155/150‑1 (MQ=255)                                                                                                    gcACCCGCAGTGGTCCACCCTGAATCCCGGGG‑CCAAGCTGAATCAGGCTGAGCGGCAGGAGCAGTTGGAGGCGCTGGTTTCACAGATTTGAGCCAGCAGGGTTGGCCTCTTCGCGGGTAAACCCGCTCCTACGGGCGGCCC‑TTGTACAac                                                <  2:1713644/150‑1 (MQ=255)                                                                                                     cACCCGCAGTGGTCCACCCTGAATCCCGGGG‑CCAAGCTGAATCAGGCTGAGCGGCAGGAGCAGTTGGAGGCGCTGGTTTCACAGATTTGAGCCAGCAGGGTTGGCCTCTTCGCGGGTAAACCCGCTCCTACGGGCGGCCT‑TTGTACAacc                                               <  2:534535/150‑1 (MQ=255)                                                                                                       cccGCAGTGGTCCACCCTGAATCCCGGGG‑CCAAGCTGAATCAGGCTGAGCGGCAGGAGCAGTTGGAGGCGCTGGTTTCACAGATTTGAGCCAGCAGGGTTGGCCTCTTCGCGGGTAAACCCGCTCCTACGGGCGGCCT‑TTGTACAaccac                                             >  1:556503/1‑150 (MQ=255)                                                                                                       cccGCAGTGGTCCACCCTGAATCCCGGGG‑CCAAGCTGAATCAGGCTGAGCGGCAGGAGCAGTTGGAGGCGCTGGTTTCACAGATTTGAGCCAGCAGGGTTGGCCTCTTCGCGGGTAAACCCGCTCCTACGGGCGGCCT‑TTGTACAaccac                                             <  2:556503/150‑1 (MQ=255)                                                                                                        ccGCAGTGGTCCACCCTGAATCCCGGGG‑CCAAGCTGAATCAGGCTGAGCGGCAGGAGCAGTTGGAGGCGCTGGTTTCACAGATTTGAGCCAGCAGGGTTGGCCTCTTCGCGGGTAAACCCGCTCCTACGGGCGGCCT‑TTGTACAaccacc                                            >  2:890224/1‑150 (MQ=255)                                                                                                          gCAGTGGTCCACCCTGAATCCCGGGG‑CCAAGCTGAATCAGGCTGAGCGGCAGGAGCAGTTGGAGGCGCTGGTTTCACAGATTTGAGCCAGCAGGGTTGGCCTCTTCGCGGGTAAACCCGCTCCTACGGGCGGCCT‑TTGTACAACCACCGa                                          <  1:1606648/150‑1 (MQ=255)                                                                                                                  ccACCCTGAATCCCGGGG‑CCAAGCTGAATCAGGCTGAGCGGCAGGAGCAGTTGGAGGCGCTGGTTTCACAGATTTGAGCCAGCAGGGTTGGCCTCTTCGCGGGTAAACCCGCTCCTACGGGCGGCCT‑TTGTACAACCACCGATGTGCGTc                                  >  2:550655/1‑150 (MQ=255)                                                                                                                     cccTGAATCCCGGGG‑CCAAGCTGAATCAGGCTGAGCGGCAGGAGCAGTTGGAGGCGCTGGTTTCACAGATTTGAGCCAGCAGGGTTGGCCTCTTCGCGGGTAAACCCGCTCCTACGGGCGGCCT‑TTGTACaa                                                 >  2:601919/1‑132 (MQ=255)                                                                                                                     cccTGAATCCCGGGG‑CCAAGCTGAATCAGGCTGAGCGGCAGGAGCAGTTGGAGGCGCTGGTTTCACAGATTTGAGCCAGCAGGGTTGGCCTCTTCGCGGGTAAACCCGCTCCTACGGGCGGCCT‑TTGTACaa                                                 <  1:601919/132‑1 (MQ=255)                                                                                                                     cccTGAATCCCGGGG‑CCAAGCTGAATCAGGCTGAGCGGCAGGAGCAGTTGGAGGCGCTGGTTTCACAGATTTGAGCCAGCAGGGTTGGCCTCTTCGCGGGTAAACCCGCTCCTACGGGCGGCCT‑TTGTACAACCACCGATGTGCGTccac                               <  2:1743253/150‑1 (MQ=255)                                                                                                                        tGAATCCCGGGG‑CCAAGCTGAATCAGGCTGAGCGGCAGGAGCAGTTGGAGGCGCTGGTTTCACAGATTTGAGCCAGCAGGGTTGGCCTCTTCGCGGGTAAACCCGCTCCTACGGGCGGCCT‑TTGTACAACCACCGATGTGCGTccaccac                            <  2:1412947/150‑1 (MQ=255)                                                                                                                           aTCCCGGGG‑CCAAGCTGAATCAGGCTGAGCGGCAGGAGCAGTTGGAGGCGCTGGTTCCGCAGTTCAGACGTATCATGGTTGGCCTCTTCTCGGGTAAACCAGCTCTTACGGGCGGCCT‑TTGTACAACCACCGATGTGCGTCCACCACCaa                         <  1:1875772/150‑1 (MQ=255)                                                                                                                                 ggg‑CCAAGCTGAATCAGGCTGAGCGGCAGGAGCAGTTGGAGGCGCTGGTTTCACAGATTTGAGCCAGCAGGGTTGGCCTCTTCGCGGGTAAACCCGCTCCTACGGGCGGCCT‑TTGTACAACCACCGATGTGCGTCCACCACCAATCCGGc                   <  2:322363/150‑1 (MQ=255)                                                                                                                                  gg‑CCAAGCTGAATCAGGCTGAGCGGCAGGAGCAGtt                                                                                                                                     <  2:741917/36‑1 (MQ=255)                                                                                                                                  gg‑CCAAGCTGAATCAGGCTGAGCGGCAGGAGCAGtt                                                                                                                                     >  1:741917/1‑36 (MQ=255)                                                                                                                                      cAAGCTGAATCAGGCTGAGCGGCAGGAGCAGTTGGAGGCGCTGGTTTCACAGATTTGAGCCAGCAGGGTTGGCCTCTTCGCGGGTAAACCCGCCCCTACGGGCGGCCT‑TTGTACAACCACCGATGTGCGTcc                                 >  1:428567/1‑132 (MQ=255)                                                                                                                                      cAAGCTGAATCAGGCTGAGCGGCAGGAGCAGTTGGAGGCGCTGGTTTCACAGATTTGAGCCAGCAGGGTTGGCCTCTTCGCGGGTAAACCCGCCCCTACGGGCGGCCT‑TTGTACAACCACCGATGTGCGTcc                                 <  2:428567/132‑1 (MQ=255)                                                                                                                                         gCTGAATCAGGCTGAGCGGCAGGAGCAGTTGGAGGCGCTGGTTTCACAGATTTGAGCCAGCAGGGTTGGCCTCTTCGCgg                                                                                   >  2:954893/1‑80 (MQ=255)                                                                                                                                         gCTGAATCAGGCTGAGCGGCAGGAGCAGTTGGAGGCGCTGGTTTCACAGATTTGAGCCAGCAGGGTTGGCCTCTTCGCgg                                                                                   <  1:954893/80‑1 (MQ=255)                                                                                                                                          cTGAATCAGGCTGAGCGGCAGGAGCAGTTGGAGGCGCTGGTTTCACAGATTTGAGCCAGCAGGGTTGGCCTCTTCGCGGGTAAACCCGCTCCTACGGGCGGCCT‑TTGTCCAACCACCGAAGTGCGTCCGCGACGAATCCTGCGTAGTAGc           <  2:223233/150‑1 (MQ=38)                                                                                                                                          cTGAATCAGACTGAGCGGCAGGAGCAGTTGGAGGCGCTGGTTTCACAGATTTGAGACAGCAGGGTTGGCCTCTTCGCGGGTAAACCCGCTCCTACGGGCGGCCT‑TTGTACAacc                                               >  2:1033176/1‑114 (MQ=255)                                                                                                                                          cTGAATCAGACTGAGCGGCAGGAGCAGTTGGAGGCGCTGGTTTCACAGATTTGAGACAGCAGGGTTGGCCTCTTCGCGGGTAAACCCGCTCCTACGGGCGGCCT‑TTGTACAacc                                               <  1:92068/114‑1 (MQ=255)                                                                                                                                          cTGAATCAGACTGAGCGGCAGGAGCAGTTGGAGGCGCTGGTTTCACAGATTTGAGACAGCAGGGTTGGCCTCTTCGCGGGTAAACCCGCTCCTACGGGCGGCCT‑TTGTACAacc                                               >  2:92068/1‑114 (MQ=255)                                                                                                                                          cTGAATCAGACTGAGCGGCAGGAGCAGTTGGAGGCGCTGGTTTCACAGATTTGAGACAGCAGGGTTGGCCTCTTCGCGGGTAAACCCGCTCCTACGGGCGGCCT‑TTGTACAacc                                               <  1:1033176/114‑1 (MQ=255)                                                                                                                                           tGAATCAGGCTGAGCGGCAGGAGCAGTTGGAGGCGCTGGTTTCACAGATTTGAGCCAGCAGGGTTGGCCTCTTCGCGGGTAAACCCGCTCCTACGGGCGGCCT‑TTGTACAACCACCGATGTGCGTCCACCACCAATCCGGCGTAGGAGCg          <  1:105251/150‑1 (MQ=255)                                                                                                                                            gAATCAGGCTGAGCGGCAGGAGCAGTTGGAGGCGCTGGTTTCACAGATTTGAGCCAGCAGGGTTGGCCTCTTCGCGGGTAAACCCGCTCCTACGGGCGGCCT‑TTGTACAACCACCGATGTGCGTCCACCACCAATCCGGCGTAGGAGCgg         >  2:1490269/1‑150 (MQ=255)                                                                                                                                            gAATCAGGCTGAGCGGCAGGAGCAGTTGGAGGCGCTGGTTTCACAGATTTGAGCCAGCAGGGTTGGCCTCTTCGCGGGTAAACCCGCTCCTACGGGCGGCCT‑TTGTACAACCACCGATGTGCGTCCACCACCAATCCGGCGTAGGAGCgg         <  1:550655/150‑1 (MQ=255)                                                                                                                                            gAATCAGGCTGAGCGGCAGGAGCAGTTGGAGGCGCTGGTTTCACAGATTTGAGCCAGCAGGGTTGGCCTCTTCGCGGGTAAACCCGCTCATACGAAGAGCGTCGTGT‑CAACCACCGATGTGCGTCCACCACCAATCCGGCGTAGGAGCgg         >  2:529067/1‑150 (MQ=255)                                                                                                                                             aaTCAGGCTGAGCGGCAGGAGCAGTTGGAGGCGCTGGTTTCACAGATTTGAGCAAGCAGGGTTGGCCTCTTCGCGGGTAAACCCGCTCCTACGGGCGGCCT‑TTGTACAACCACCGATGTGCGTCCACCACCAATCCGGCGTAGGAGCggg        <  2:1036211/150‑1 (MQ=255)                                                                                                                                                  ggCTGAGCGGCAGGAGCAGTTGGAGGCGCTGGTTTCACAGATTTGAGCCAGCAGGGTTGGCCTCTTCGCGGGTAAACCCGCTCCTACGGGCGGCCT‑TTGTACAACCACCGATGTGCGTCCACCACCAATCCGGCGTAGGAGCGGGTTTAc   <  1:1490269/150‑1 (MQ=255)                                                                                                                                                  ggCTGAGCGGCAGGAGCAGTTGGAGGCGCTGGTTTCACAGATTTGAGCCAGCAGGGTTGGCCTCTTCGCGGGTAAACCCGCTCCTACGGGCGGCCT‑TTGTACAACCACCGATGTGCGTACACCTCCAATCCGGCGTAGGAGCGGGTTTAc   <  1:529067/150‑1 (MQ=255)                                                                                                                                                   gCTGAGCGGCAGGAGCAGTTGGAGGCGCTGGTTTCACAGATTCGAGCCAGCAGTGTGGGCCTCTTCGCGGGTCCCCCCGCTCCTACGGGCGGCCT‑TTGTACAACCACCGATGTGCGTCCACCACCAATCCGGCGTAGGAGCGGGTTTAcc  <  1:890224/150‑1 (MQ=255)                                                                                                                                                   |                                                                                                                                                       CCGCGGGCATGCGGCGCAAGACAGCATTCGCCTGGGTGAGCTGGATGAAAAGGCCTGGCGTGAGGGCAAGTTGCCGGAGGCGTTGGTCAAAGCGCGGGCGCACCCGCAGTGGTCCACCCTGAATCCCGGGG‑CCAAGCTGAATCAGGCTGAGCGGCAGGAGCAGTTGGAGGCGCTGGTTTCACAGATTTGAGCCAGCAGGGTTGGCCTCTTCGCGGGTAAACCCGCTCCTACGGGCGGCCT‑TTGTACAACCACCGATGTGCGTCCACCACCAATCCGGCGTAGGAGCGGGTTTACC  >  NC_002947/5858689‑5858983

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