Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,647,850	T→G	7.3%	D430A (GAC→GCC) ‡	PP_4112 ←	sulfate transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,647,850	0	T	G	7.3%	206.7 / 6.6	69	D430A (GAC→GCC) ‡	PP_4112	sulfate transporter
Reads supporting (aligned to +/- strand):  ref base T (30/34);  new base G (3/2);  total (33/36)
Fisher's exact test for biased strand distribution p-value = 6.65e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.45e-01

AACAGCTTCTTGCACCGCTCCGACAAATGCACCACGCGCAAATGCTTGCCCGCCTTGGCATAGC‑GCTCGCGCAGGGTTTTCAGCGCCGCAACGGCCGAATAATCGACAAAGCTCAAATGCTGGCAA‑TCCAGCGTGACCTGGGCCGGGTCCCCGGCCGTGTCGAACTGGTTGAGGAACGAGGTGGTCGAGGCAAAGAACAGCGTGCCATGCAGTTGATAACACTTGCCGCCCTCCGCGTCCTCATGGGTGTCGGCATACAGCTCCCGCGCATGCTGCCAGGCAAAGTTGAGCGCCG  >  NC_002947/4647703‑4647997                                                                                                                                                      |                                                                                                                                                    aaCAGCTTCTTGCACCGCTCCGACAAATGCACCACGCGCAAATGCTTGCCCGCCTTGGCATAGC‑GCTCGCGCAGGGTTTTCAGCGCCGCAACGGCCGAATAATCGACAAAGCTCAAATGCTGGCAA‑TCCAGCGTGACCTGGGCCGGGTc                                                                                                                                                    >  8:59846/1‑149 (MQ=255)  aCAGCTTCTTGCACCGCTCCGACAAATGCACCACGCGCAAATGCTTGCCCGCCTTGGCATAGC‑GCTCGCGCAGGGTTTTCAGCGCCGCAACGGCCGAATAATCGACAAAGCTCAAATGCTGGCAA‑TCCAGCGTGACCTGGGCCGGGTcc                                                                                                                                                   >  1:220837/1‑149 (MQ=255)    aGCTTCTTGCACCGCTCCGACAAATGCACCACGCGCAAATGCTTGCCCGCCTTGGCATAGC‑GCTCGCGCAGGGTTTTCAGCGCCGCAACGGCCGAATAATCGACAAAGCTCAAATGCTGGCAA‑TCCAGCGTGACCTGGGCCGGGTcccc                                                                                                                                                 <  8:140456/149‑1 (MQ=255)     gCTTCTTGCACCGCTCCGACAAATGCACCACGCGCAAATGCTTGCCCGCCTTGGCATAGC‑GCTCGCGCAGGGTTTTCAGCGCCGCAACGGCCGAATAATCGACAAAGCTCAAATGCTGGCAA‑TCCAGCGTGACCTGGGCCGGGTCCCCg                                                                                                                                                <  5:211871/149‑1 (MQ=255)      cttcttGCACCGCTCCGACAAATGCACCACGCGCAAATGCTTGCCCGCCTTGTCATAGC‑GCTCGCGCAGGGTTTTCAGCGCCGCAACGGCCGAATAATCGACAAAGCTCAAATGCTGGCAA‑TCCAGCGTGACCTGGGCCGGGTCCCCgg                                                                                                                                               <  5:105750/149‑1 (MQ=255)      cttcttGCACCGCTCCGACAAATGCACCACGCGCAAATGCTTGCCCGCCTTGGCATAGC‑GCTCGCGCAGGGTTTTCAGCGCCGCAACGGCCGAATAATCGACAAAGCTCAAATGCTGGCAA‑TCCAGCGTGACCTGGGCCGGGTCCCCgg                                                                                                                                               <  5:127275/149‑1 (MQ=255)         cttGCACCGCTCCGACAAATGCACCACGCGCAAATGCTTGCCCGCCTTGGCATAGC‑GCTCGCGCAGGGTTTTCAGCGCCGCAACGGCCGAATAATCGACAAAGCTCAAATGCTGGCAA‑TCCAGCGTGACCTGGGCCGGGTCCCCGGcc                                                                                                                                             <  2:164109/148‑1 (MQ=255)          ttGCACCGCTCCGACAAATGCACCACGCGCAAATGCTTGCCCGCCTTGGCATAGC‑GCTCGCGCAGGGTTTTCAGCGCCGCAACGGCCGAATAATCGACAAAGCTCAAATGCTGGCAA‑TCCAGCGTGACCTGGGCCGGGGACCAGGccca                                                                                                                                           >  3:62874/1‑147 (MQ=255)           tGCACCGCTCCGACAAATGCACCACGCGCAAATGCTTGCCCGCCTTGGCATAGC‑GCTCGCGCAGGGTTTTCAGCGCCGCAACGGCCGAATAATCGACAAAGCTCAAATGCTGGCAA‑TCCAGCGTGACCTGGGCCGGGGACCCGGccca                                                                                                                                           <  4:62874/148‑3 (MQ=255)            gCACCGCTCCGACAAATGCACCACGCGCAAATGCTTGCCCGCCTTGGCATAGC‑GCTCGCGCAGGGTTTTCAGCGCCGCAACGGCCGAATAATCGACAAAGCTCAAATGCTGGCAACTCCAGCGTGACCTGGGCCGGGTCCCCGGcc                                                                                                                                             <  1:183497/146‑1 (MQ=255)             cACCGCTCCGACAAATGCACCACGCGCAAATGCTTGCCCGCCTTGGCATAGC‑GCTCGCGCAGGGATTTCAGCGCCGCAACGGCCGAATAATCGACAAAGCTCAAATGCTGGCAA‑TCCAGCGTGACCTGGGCCGGGTCCCCGGCCGTGTc                                                                                                                                        >  1:90808/1‑149 (MQ=255)             cACCGCTACGACAAATGCACCACTCGCAAATGATTGCCCGCCTTTGCATAGC‑GCTCGCGCAGGCTTTTCAGCGCCGCAACGGCCGAATAATCGACAAAGCTCAAATGCTGGCAA‑TCCAGCGTGACCTGGGCCGGGTCCCCGGCCGTGTc                                                                                                                                        >  5:306179/1‑149 (MQ=255)              aCCGCTCCGACAAATGCACCACGCGCAAATGCTTGCCCGCCTTGGCATAGCAGCTCGCGCAGGGTTTTCAGCGCCGCAACGGCCGAATAATCGACAAAGCTCAAATGCTGGCAA‑TCCAGCGTGACCTGGGCCGGGTCCCCGGCCGTGTc                                                                                                                                        >  1:75392/1‑149 (MQ=255)              aCCGCTCCGACAAATGCACCACGCGCAAATGCTTGCCCGCCTTGGCATAGC‑GCTCGCGCAGGGTTTTCAGCGCCGCAACGGCCGAATAATCGACAAAGCTCAAATGCTGGCAA‑TCCAGCGTGACCTGGGCCGGGTCCCCGGCCGTGTc                                                                                                                                        >  8:161300/1‑148 (MQ=255)                     cGACAAATGCACCACGCGCAAATGCTTGCCCGCCTTGGCATAGC‑GCTCGCGCAGGGTTTTCAGCGCCGCAACGGCCGAATAATCGACAAAGCTCAAATGCTGGCAA‑TCCAGCGTGACCTGGGCCGGGTCCCCGGCCGTGTCGAACTGGt                                                                                                                                <  2:75392/149‑1 (MQ=255)                         aaaTGCACCACGCGCAAATGCTTGCCCGCCTTGGCATAGC‑GCTCGCGCAGGGTTTTCAGCGCCGCAACGGCCGAATAATCGACAAAGCTCAAATGCTGGCAA‑TCCAGCGTGACCTGGGCCGGGTCCCCGGCCGTGTCGAACTGGTTGAg                                                                                                                            <  8:5857/149‑1 (MQ=255)                              caccacGCGCAAATGCTTGCCCGCCTTGGCATAGC‑GCTCGCGCAGGGTTTTCAGCGCCGCAACGGCCGAATAATCGACAAAGCTCAAATGCTGGCAA‑TCCAGCGTGACCTGGGCCGGGTCCCCGGCCGTGTCGAACTGGTTGAGGAACg                                                                                                                       <  7:59846/149‑1 (MQ=255)                              caccacGCGCAAATGCTTGCCCGCCTTGGCATAGC‑GCTCGCGCAGGGTTTTCAGCGCCGCAACGGCCGAATAATCGACAAAGCTCAAATGCTGGCAA‑TCCAGCGTGACCTGGGCCGGGTCCCCGGCCGTGTCGAACTGGTTGAGGAACg                                                                                                                       <  2:90808/149‑1 (MQ=255)                                ccacGCGCAAATGCTTGCCCGCCTTGGCATAGC‑GCTCGCGCAGGGTTTTCAGCGCCGCAACGGCCGAATAATCGACAAAGCTCAAATGCTGGCAA‑TCCAGCGTGACCTGGGCCGGGTCCCCGGCCGTGTCGAACTGGTTGAGGAACGAg                                                                                                                     >  2:265979/1‑149 (MQ=255)                                 cacGCGCAAATGCTTGCCCGCCTTGGCATAGC‑GCTCGCGCAGGGTTTTCAGCGCCGCAACGGCCGAATAATCGACAAAGCTCAAATGCTGGCAA‑TCCAGCGTGACCTGGGCCGGGTCCCCGGCCGTGTCGAACTGGTTGAGGAACGAgg                                                                                                                    >  6:257963/1‑149 (MQ=255)                                  acGCGCAAATGCTTGCCCGCCTTGGCATGGC‑GCTCGCGCAGGGTTTTCAGCGCCGCAACGGCCGAATAATCGACAAAGCTCAAATGCTGGCAA‑TCCAGCGTGACCTGGGCCGGGGACCCGGcccaggt                                                                                                                                        >  6:93039/1‑123 (MQ=255)                                  acGCGCAAATGCTTGCCCGCCTTGGCATGGC‑GCTCGCGCAGGGTTTTCAGCGCCGCAACGGCCGAATAATCGACAAAGCTCAAATGCTGGCAA‑TCCAGCGTGACCTGGGCCGGGGACCCGGcccaggt                                                                                                                                        <  5:93039/128‑6 (MQ=255)                                  acGCGCAAATGCTTGCCCGCCTTGGCATAGC‑GCTCGCGCAGGGTTTTCAGCGCCGCAACGGCCGAATAATCGACAAAGCTCAAATGCTGGCAA‑TCCAGCGTGACCTGGGCCGGGTCCCCGGCCGTGTCGAACTGGTTGAGGAACGAggt                                                                                                                   >  7:128854/1‑149 (MQ=255)                                         aaTGCTTGCCCGCCTTGGCATAGC‑GCTCGCGCAGGGTTTTCAGCGCCGCAACGGCCGAATAATCGACAAAGCTCAAATGCTGGCAA‑TCCAGCGTGACCTGGGCCGGGTCCCCGGCCGTGTCGAACTGGTTGAGGAACGAGGTGGTCGAg                                                                                                            <  2:277696/149‑1 (MQ=255)                                          aTGCTTGCCCGCCTTGGCATAGC‑GCTCGCGCAGGGTTTTCAGCGCCGCAACGGCCGAATAATCGACAAAGCTCAAATGCTGGCAA‑TCCAGCGTGACCTGGGCCGGGTCCCCGGCCGTGTCGAACTGGTTGAGGAACGAGGTGGTCGAgg                                                                                                           <  2:220837/149‑1 (MQ=255)                                          aTGCTTGCCCGCCTTGGCATAGC‑GCTCGCGCAGGGTTTTCAGCGCCGCAACGGCCGAATAATCGACAAAGCTCAAATGCTGGCAA‑TCCAGCGTGACCTGGGCCGGGTCCCCGGCCGTGTCGAACTGGTTGAGGAACGAGGTGGTCGAgg                                                                                                           >  5:22837/1‑149 (MQ=255)                                          aTGCTTGCCCGCCTTGGCATAGC‑GCTCGCGCAGGGTTTTCAGCGCCGCAACGGCCGAATAATCGACAAAGCTCAAATGCTGGCAA‑TCCAGCGTGACCTGGGCCGGGTCCCCGGCCGTGTCGAACTGGTTGAGGAACGAGGTGGTCGAgg                                                                                                           <  1:265979/149‑1 (MQ=255)                                           tGCTTGCCCGCCTTGGCATAGC‑GCTCGCGCAGGGTTTTCAGCGCCGCAACGGCCGAATAATCGACAAAGCTCAAATGCTGGCAA‑TCCAGCGTGACCTGGGCCGGGTCCCCGGCCGTGTCGAACTGGTTGAGGAACGAGGTGGTCGAGGc                                                                                                          <  6:52682/149‑1 (MQ=255)                                                 cccGCCTTGGCATAGC‑GCTCGCGCAGGGTTTTCAGCGCCGCAACGGCCGAATAATCGACAAAGCTCAAATGCTGGCAA‑TCCAGCGTGACCTGGGCCGGGTCCCCGGCCGTGTCGAACTGGTTGAGGAACGAGGTGGTCGAGGCAAAGaa                                                                                                    <  2:104880/149‑1 (MQ=255)                                                     ccTTGGCATAGC‑GCTCGCGCAGGGTTTTCAGCGCCGCAACGGCCGAATAATCGACAAAGCTCAAATGCTGGCAA‑TCCAGCGTGACCTGGGCCGGGTCCCCGGCCGTGTCGAACTGGTTGAGGAACGAGGTGGTCGAGGCAAAGAACAGc                                                                                                <  5:257963/149‑1 (MQ=255)                                                      cTTGGCATAGC‑GCTCGCGCAGGGTTTTCAGCGCCGCAACGGCCGAATAATCGACAAAGCTCAAATGCTGGCAA‑TCCAGCGTGACCTGGGCCGGGTCCCCGGCCGTGTCGAACTGGTTGAGGAACGAGGTGGTCGAGGCAAAGAACAGCg                                                                                               >  1:204026/1‑149 (MQ=255)                                                             tAGC‑GCTCGCGCACGGCTTTCAGCGCCGCAACGGCCGAATAATCGACAAAGCTCAAATGCTGGCAA‑TCCAGCGTGACCTGGGCCGGGTCCCCGGCCGTGTCGAACTGGTTGAGGAACGAGGTGGTCGAGGCAAAGAACAGCGTGCCATg                                                                                        <  1:106666/149‑1 (MQ=255)                                                               gc‑gcTCGCGCAGGGTTTTCAGCGCCGCAACGGCCGAATAATCGACAAAGCTCAAATGCTGGCAA‑TCCAGCGTGACCTGGGCCGGGTCCCCGGCCGTGTCGAACTGGTTGAGGAACGAGGTGGTCGAGGCAAAGAACAGCGTGCCATGc                                                                                       <  2:204026/148‑1 (MQ=255)                                                                   cTCGCGCAGGGTTTTCAGCGCCGCAACGGCCGAATAATCGACAAAGCTCAAATGCTGGCAA‑TCCAGCGTGACCTGGGCCGGGTCCCCGGCCGTGTCGAACTGGTTGAGGAACGAGGTGGTCGAGGCAAAGAACAGc                                                                                                >  2:180061/1‑136 (MQ=255)                                                                   cTCGCGCAGGGTTTTCAGCGCCGCAACGGCCGAATAATCGACAAAGCTCAAATGCTGGCAA‑TCCAGCGTGACCTGGGCCGGGTCCCCGGCCGTGTCGAACTGGTTGAGGAACGAGGTGGTCGAGGCAAAGAACAGc                                                                                                <  1:180061/136‑1 (MQ=255)                                                                      gcgcAGGGTTTTCAGCGCCGCAACGGCCGAATAATCGACAAAGCTCAAATGCTGGCAA‑TCCAGCGTGACCTGGGCCGGGTCCCCGGCCGTGTCGAACTGGTTGAGGAACGAGGTGGTCGAGGCAAAGAACAGCGTGCCATGCAGTTGAt                                                                                <  8:128854/149‑1 (MQ=255)                                                                       cgcAGGGTTTTCAGCGCCGCAACGGCCGAATAATCGACAAAGCTCAAATGCTGGCAA‑TCCAGCGTGACCTGGGCCGGGTCCCCGGCCGTGTCGAACTGGTTGAGGAACGAGGTGGTCGAGGCAAAGAACAGCGTGCCATGCAGTTGAt                                                                                >  8:84429/1‑148 (MQ=255)                                                                        gcAGGGTTTTCAGCGCCGCAACGGCCGAATAATCGACAAAGCTCAAATGCTGGCAA‑TCCAGCGTGACCTGGGCCGGGTCCCCGGCCGTGTCGAACTGGTTGAGGAACGAGGTGGTCGAGGCAAAGAACAGCGTGcc                                                                                           <  4:219478/136‑1 (MQ=255)                                                                        gcAGGGTTTTCAGCGCCGCAACGGCCGAATAATCGACAAAGCTCAAATGCTGGCAA‑TCCAGCGTGACCTGGGCCGGGTCCCCGGCCGTGTCGAACTGGTTGAGGAACGAGGTGGTCGAGGCAAAGAACAGCGTGCCATGCAGTTGATaa                                                                              <  4:148813/149‑1 (MQ=255)                                                                          aGGGTTTTCAGCGCCGCAACGGCCGAATAATCGACAAAGCTCAAATGCTGGCAA‑TCCAGCGTGACCTGGGCCGGGTCCCCGGCCGTGTCGAACTGGTTGAGGAACGAGGTGGTCGAGGCAAAGAACAGCGTGCCATGCAGTTGATAac                                                                             >  4:263023/1‑148 (MQ=255)                                                                                ttCAGCGCCGCAACGGCCGAATAATCGACAAAGCTCAAATGCTGGCAA‑TCCAGCGTGACCTGGGCCGGGTCCCCGGCCGTGTCGAACTGGTTTAGGAACGAGGTGGTCGAGGCAAAGAACAGCGTGCCATGCAGTTGATAACACTTgc                                                                       >  5:176841/1‑148 (MQ=255)                                                                                ttCAGCGCCGCAACGGCCGAATAATCGACAAAGCTCAAATGCTGGCAA‑TCCAGCGTGACCTGGGCCGGGTCCCCGGCCGTGTCGAACTGGTTGAGGAACGAGGTGGTCGAGGCAAAGAACAGCGTGCCATGCAGTTGATAACACTTgc                                                                       <  6:176841/148‑1 (MQ=255)                                                                                    gcgcCGCAACGGCCGAATAATCGACAAAGCTCAAATGCTGGCAA‑TCCAGCGTGACCTGGGCCGGGTCCCCGGCCGTGTCGAACTGGTTGAGGAACGAGGTGGTCGAGGCAAAGAACAGCGTGCCATGCAGTTGATAACACTTGCCGccc                                                                  >  3:152042/1‑149 (MQ=255)                                                                                    gcgcCGCAACGGCCGAATAATCGACAAAGCTCAAATGCTGGCAA‑TCCAGCGTGACCTGGGCCGGGTCCCCGGCCGTGTCGAACTGGTTGAGGAACGAGGTGGTCGAGGCAAAGAACAGCGTGCCATGCAGTTGATAACACTTGCCGccc                                                                  >  7:141412/1‑149 (MQ=255)                                                                                     cgccgcAACGGCCGAATAATCGACAAAGCTCAAATGCTGGCAA‑TCCAGCGTGACCTGGGCCGGGTCCCCGGCCGTGTCGAACTGGTTGAGGAACGAGGTGGTCGAGGCAAAGAACAGCGTGCCATGCAGTTGATAACACTTGCCGCCCt                                                                 >  7:8606/1‑149 (MQ=255)                                                                                        cgcAACGGCCGAATAATCGACAAAGCTCAAATGCTGGCAA‑TCCAGCGTGACCTGGGCCGGGTCCCCGGCCGTGTCGAACTGGTTGAGGAACGAGGTGGTCGAGGCAAAGAACAGCGTGCCATGCAGTTGATAACACTTGCCGCCCTCcg                                                              <  7:161300/149‑1 (MQ=255)                                                                                        cgcAACGGCCGAATAATCGACAAAGCTCAAATGCTGGCAA‑TCCAGCGTGACCTGGGCCGGGTCCCCGGCCGTGTCGAACTGGTTGAGGAACGAGGTGGTCGAGGCAAAGAACAGCGTGCCATGCAGTTGATAACACTTGCCGCCCTCcg                                                              >  7:188705/1‑149 (MQ=255)                                                                                          cAACGGCCGAATAATCGACAAAGCTCAAATGCTGGCAA‑TCCAGCGTGACCTGGGCCGGGTCCCCGGCCGTGTCGAACTGGTTGAGGAACGAGGTGGTCGAGGCAAAGAACAGCGTGCCATGCAGTTGATAACACTTGCCGCCCTCcgcg                                                            <  8:80667/149‑1 (MQ=255)                                                                                                    ataatCGACAAAGCTCAAATGCTGGAAA‑TCCAGCGTTACCTGGGCTGGGTCCCAGGCCGTGTAGAAATGTTTGATGAACGAGGTGGTCGAGGCAAAGAACTGAGTGCCATGCAGTTGATAAGACTTGCCGCCCTCCGCGTACTTATgg                                                   <  6:306179/148‑1 (MQ=255)                                                                                                     taatCGACAAAGCTCAAATGCTGGCAA‑TCCAGCGTGACCTGGGCCGGGTCCCCGGCCGTGTCGAACTGGTTGAGGAACGAGGTGGTCGAGGCAAAGAACAGCGTGCCATGCAGTTGATAACACTTGCCGCCCTCCGCGTCCTCATGGgt                                                 >  5:69616/1‑149 (MQ=255)                                                                                                       atCGACAAAGCTCAAATGCTGGCAA‑TCCAGCGTGACCTGGGCCGGGTCCCCGGCCGTGTCGAACTGGTTGAGGAACGAGGTGGTCGAGGCAAAGAACAGCGTGCCATGCAGTTGATAACACTTGCCGCCCTCCGCGTCCTCATGGgtgt                                               >  4:99441/1‑149 (MQ=255)                                                                                                       atCGACAAAGCTCAAATGCTGGCAA‑TCCAGCGTGACCTGGGCCGGGTCCCCGGCCGTGTCGAACTGGTTGAGGAACGAGGTGGTCGAGGCAAAGAACAGCGTGCCATGCAGTTGATAACACTTGCCGCCCTCCGCGTCCTCATGGgtg                                                <  6:22837/148‑1 (MQ=255)                                                                                                          gACAAAGCTCAAATGCTGGCAA‑TCCAGCGTGACCTGGGCCGGGTCCCCGGCCGTGTCGAACTGGTTGAGGAACGAGGTGGTCGAGGCAAAGAACAGCGt                                                                                              >  8:53569/1‑99 (MQ=255)                                                                                                          gACAAAGCTCAAATGCTGGCAA‑TCCAGCGTGACCTGGGCCGGGTCCCCGGCCGTGTCGAACTGGTTGAGGAACGAGGTGGTCGAGGCAAAGAACAGCGt                                                                                              <  7:53569/99‑1 (MQ=255)                                                                                                            cAAAGCTCAAATGCTGGCAA‑TCCAGCGTGACCTGGGCCGGGTCCCCGGCCGTGTCGAACTGGTTGAGGAACGAGGTGGTCGAGGCAAAGAACAGCGTGCCATGCAGTTGATAACACTTGCCGCCCTCCGCGTCCTCATGGGTGTCGGCa                                          <  8:8606/149‑1 (MQ=255)                                                                                                             aaaGCTCAAATGCTGGCAA‑TGCAGCGTGACCTGGGCCGGGTCCCCGGCCGTGTCGAACTGGTTGAGGAACGAGGTGGTCGAGGCAAAGAACAGCGTGCCATGCAGTTGATAACACTTGCCGCCCTCCGCGTCCTCATGGGTGTCGGCa                                          <  5:289731/148‑1 (MQ=255)                                                                                                                  tCAAATGCTGGCAA‑TCCAGCGTGACCTGGGCCGGGTCCCCGGCCGTGTCGAACTGGTTGAGGAACGAGGTGGTCGAGGCAAAGAACAGCGTGCCATGCAGTTGATAACACTTGCCGCCCTCCGCGTCCTCATGGGTGTCGGCATACAGc                                    >  8:188124/1‑149 (MQ=255)                                                                                                                   cAAATGCTGGCAA‑TCCAGCGTGACCTGGGCCGGGTCCCCGGCCGTGTCGAACTGGTTGAGGAACGAGGTGGTCGAGGCAAAGAACAGCGTGCCATGCAGTTGATAACACTTGCCGCCCTCCGCGTCCTCATGGGTGTCGGCATACAGCt                                   >  5:108293/1‑149 (MQ=255)                                                                                                                        gCTGGCAA‑TCCAGCGTGACCTGGGCCGGGTCCCCGGCCGTGTCGAACTGGTTGAGGAACGAGGTGGTCGAGGCAAAGAACAGCGTGCCATGCAGTTGATAACACTTGCCGCCCTCCGCGTCCTCATGGGTGTCGGCATACAGCTCCcgc                              >  1:286810/1‑149 (MQ=255)                                                                                                                            gCAA‑TCCAGCGTGACCTGGGCCGGGTCCCCGGCCGTGTCGAACTGGTTGAGGAACGAGGTGGTCGAGGCAAAGAACAGCGTGCCATGCAGTTGATAACACTTGCCGCCCTCCGCGTCCTCATGGGTGTCGGCATACAGCTCCCGCGCAt                          <  3:263023/149‑1 (MQ=255)                                                                                                                            gCAA‑TCCAGCGTGACCTGGGCCGGGTCCCCGGCCGTGTCGAACTGGTTGAGGAACGAGGTGGTCGAGGCAAAGAACAGCGTGCCATGCAGTTGATAACACTTGCCGCCCTCCGCGTCCTCATGGGTGTCGGCATACAGCTCCCGCGCAt                          >  3:30879/1‑149 (MQ=255)                                                                                                                              aa‑TCCAGCGTGACCTGGGCCGGGTCCCCGGCCGTGTCGAAGTGGTTGAGGAACGAGGTGGTAGAGGCAATGTACAGCGTGGCATGCAGTTGATTCGACTTCACGCCCGCCGCGACATCAAGGGTGTCGGCAGACGGCTCCCGCGCAggc                        >  8:214663/1‑149 (MQ=255)                                                                                                                               a‑TCCAGCGTGACCTGGGCCGGGTCCCCGGCCGTGTCGAACTGGTTGAGGAACGAGGTGGTCGAGGCAAAGAACAGCGTGCCATGCAGTTGATAACACTTGCCGCCCTCCGCGTCCTCATGGgtgt                                               >  2:259534/1‑125 (MQ=255)                                                                                                                               a‑TCCAGCGTGACCTGGGCCGGGTCCCCGGCCGTGTCGAACTGGTTGAGGAACGAGGTGGTCGAGGCAAAGAACAGCGTGCCATGCAGTTGATAACACTTGCCGCCCTCCGCGTCCTCATGGgtgt                                               <  1:259534/125‑1 (MQ=255)                                                                                                                                    aGCGTGACCTGGGCCGGGTCCCCGGCCGTGTCGAACTGGTTGAGGAACGAGGTGGTCGAGGCAAAGAACAGCGTGCCATGCAGTTGATAACACTTGCCGCCCTCCGCGTCCTCATGGGTGTCGGCATACAGCTCCCGCGCAtgctgc                     >  2:248216/1‑147 (MQ=255)                                                                                                                                     gCGTGACCTGGGCCGGGTCCCCGGCCGTGTCGAACTGGTTGAGGAACGAGGTGGTCGAGGCAAAGAACAGCGTGCCATGCAGTTGATAACACTTGCCGCCCTCCGCGTCCTCATGGGTGTCGGCATACAGCTCCCGCGCATGCTGCCAg                  <  6:69616/149‑1 (MQ=255)                                                                                                                                       gTGACCTGGGCCGGGTCCCCGGCCGTGTCGAACTGGTTGAGGAACGAGGTGGTCGAGGCAAAGAACAGCGTGCCATGCAGTTGATAACACTTGCCGCCCTCCGCGTCCTCATGGGTGTCGGCATACAGCTCCCGCGCATGCTGCCAGGc                <  8:141412/149‑1 (MQ=255)                                                                                                                                             acggCCGGGGACCCGGCCGTGTCGAACTGGTTGAGGAACGAGGTGGTCGAGGCAAAGAACAGCGTGCCATGCAGTTGATAACACTTGCCGCCCTCCGCGTCCTCATGGGTGTCGGCATACAGCTCCCGCGCATGCTGCCAGGCAAAGtt          >  2:125509/3‑149 (MQ=255)                                                                                                                                                     gTCCCCGGCCGTGTCGAACTGGTTGAGGAACGAGGTGGTCGAGGCAAAGCACAGCGTGCCATGCAGTTGATAACACTTGCCGCCCTCCGCGTCCTCATGGGTGTCGGCATACAGCTCCCGCGCATGCTGCCAGCCAAAGTTCAGcgccg  >  6:242867/1‑149 (MQ=255)                                                                                                                                                      |                                                                                                                                                    AACAGCTTCTTGCACCGCTCCGACAAATGCACCACGCGCAAATGCTTGCCCGCCTTGGCATAGC‑GCTCGCGCAGGGTTTTCAGCGCCGCAACGGCCGAATAATCGACAAAGCTCAAATGCTGGCAA‑TCCAGCGTGACCTGGGCCGGGTCCCCGGCCGTGTCGAACTGGTTGAGGAACGAGGTGGTCGAGGCAAAGAACAGCGTGCCATGCAGTTGATAACACTTGCCGCCCTCCGCGTCCTCATGGGTGTCGGCATACAGCTCCCGCGCATGCTGCCAGGCAAAGTTGAGCGCCG  >  NC_002947/4647703‑4647997

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