Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,518,539	A→T	5.9%	A14A (GCT→GCA)	PP_4850 ←	membrane protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,518,539	0	A	T	5.9%	258.0 / 3.3	84	A14A (GCT→GCA)	PP_4850	membrane protein
Reads supporting (aligned to +/- strand):  ref base A (44/35);  new base T (3/2);  total (47/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.44e-01

GTTGACCAGTTCGAACACCAGCAGACCGGCCAGCAGGCTGGGCAGCAACTTCAGGGGCAAAGCCAGCAGCAGCCCGGCCATGACGATGATCAGGCTGGCCAGGGTGATCTGACGGGGGGTAAAGGTCATACAGCCTCGACGGCAGACAGCGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAGGCGCA‑GGTCATTTCTTCTTCAGGCAATCGCTCATG‑AAGGCCTTGCGCTCATCCCCCTTCAGGGCTTTGGTGGTGGCGTCAGCATTGCAGGTCTTCATTTTTT  >  NC_002947/5518392‑5518684                                                                                                                                                    |                                                                                                                                                    gTTGACCAGTTCGAACACCAGCAGACCGGCCAGCAGGCTGGGCAGCAACTTCAGGGGCAAAGCCAGCAGCAGCCCGGCCATGACGATGATCAGGCTGGCCAGGGTGATCTGACGGGGGGTAAAGGTCATACAGCCTCGACGGCAGACAg                                                                                                                                                    >  8:409762/1‑149 (MQ=255)    gACCAGTTCGAACACCAGCAGACCGGCCAGCAGGCTGGGCAGCAACTTCAGGGGCAAAGCCAGCAGCAGCCCGGCCATGACGATGATCAGGCTGGCCAGGGTGATCTGACGGGGGGTAAAGGTCATACAGCCTCGACGGCAGACAGCg                                                                                                                                                  <  4:329104/148‑1 (MQ=255)     aCCAGTTCGAACACCAGCAGACCGGCCAGCAGGCTGGGCAGCAACTTCAGGGGCAAAGCCAGCAGCAGCCCGGCCATGACGATGATCAGGCTGGCCAGGGTGATCTGACGGGGGGTAAAGGTCATACAGCCTCGACGGCAGACAGCgg                                                                                                                                                 <  5:432298/148‑1 (MQ=255)      ccAGTTCGAACACCAGCAGACCGGCCAGCAGGCTGGGCCGCAACTTCAGGGGCAAAGCCAGCAGCAGCCCGGCCATGACGATGATCAGGCTGGCCAGGGTGATCTGACGGGGGGTAAAGGTCATACAGCCTCGACGGCAGACAGCggg                                                                                                                                                >  7:386020/1‑148 (MQ=255)         gTTCGAACACCAGCAGACCGGCCAGCAGGCTGGGCAGCAACTTCAGGGGCAAAGCCAGCAGCAGCCCGGCCATGACGATGATCAGGCTGGCCAGGGTGATCTGACGGGGGGTAAAGGTCATACAGCCTCGACGGCAGACAGCGGGAATg                                                                                                                                            >  4:298524/1‑149 (MQ=255)             gAACACCAGCAGACCGGCCAGCAGGCTGGGCAGCAACTTCAGGGGCAAAGCCAGCAGCAGCCCGGCCATGACGATGATCAGGCTGGCCAGGGTGATCTGACGGGGGGTAAAGGTCATACAGCCTCGACGGCAGACAGCGGGAATGCCg                                                                                                                                         >  1:35541/1‑148 (MQ=255)             ccgttctagcagACCGGCCAGCAGGCTGGTCAGCAACTTCTGGTGTAAAGCCAGCAGCAGCCCGGCCATGACGATGATCAGGCTGGCCAGGGTGATCTGACGGGGGGTAAAGGTCATACAGCCTCGACGGCAGACAGCGGGAATGCCg                                                                                                                                         <  2:293128/141‑1 (MQ=255)              aaCACCAGCAGACCGGCCAGCAGGCTGGGCAGCAACTTCAGGGGCAAAGCCAGCAGCAGCCCGGCCATGACGATGATCAGGCTGGCCAGGGTGATCTGACGGGGGGTAAAGGTCATACAGCCTCGACGGCAGACAGCGGGAATGCCGCa                                                                                                                                       <  1:381569/149‑1 (MQ=255)                   cagcagACCGGCCAGCAGGCTGGGCAGCAACTTCAGGGGCAAAGCCAGCAGCAGCCCGGCCATGACGATGATCAGGCTGGCCAGGGTGATCTGACGGGGGGTAAAGGTCATACAGCCTCGACGGCAGACAGCGGGAATGCCGCAGTCTg                                                                                                                                  >  6:53525/1‑149 (MQ=255)                   cagcagACCGGCCAGCAGGCTGGGCAGCAACTTCAGGGGCAAAGCCAGCAGCAGCCCGGCCATGACGATGATCAGGCTGGCCAGGGTGATCTGACGGGGGGTAAAGGTCATACAGCCTCGACGGCAGACAGCGGGAATGCCGCAGTCTg                                                                                                                                  >  2:217252/1‑149 (MQ=255)                    agcagACCGGCCAGCAGGCTGGGCAGCAACTTCAGGGGCAAAGCCAGCAGCAGCCCGGCCATGACGATGATCAGGCTGGCCAGGGTGATCTGACGGGGGGTAAAGGTCATACAGCCTCGACGGCAGACTGCGGCATTCCCGCTGTCTGc                                                                                                                                 <  4:334068/149‑1 (MQ=255)                     gcagACCGGCCAGCAGGCTGGGCAGAAACTTCAGGGGTAAAGCCAGCAGCAGCCCGGCCATGACGATGATCAGGCTGGCCAGGGTGATCTGACGGGGGGTAAAGGTCATACAGCCTCGACGGCAGACAGCGGGAATGCCg                                                                                                                                         >  1:293128/1‑140 (MQ=255)                      cagACCGGCCAGCAGGCTGGGCAGCAACTTCAGGGGCAAAGCCAGCAGCAGCCCGGCCATGACGATGATCAGGCTGGCCAGGGTGATCTGACGGGGGGTAAAGGTCATACAGCCTCGACGGCAGACAGCGGGAATGCCGCAGTCTGCCa                                                                                                                               <  2:35541/149‑1 (MQ=255)                         aCCGGCCAGCAGGCTGGGCAGCAACTTCAGGGGCAAAGCCAGCAGCAGCCCGGCCATGACGATGATCAGGCTGGCCAGGGTGATCTGACGGGGGGTAAAGGTCATACAGCCTCGACGGCAGACAGCGGGAATGCCGCAGTCTGcc                                                                                                                                >  4:231221/1‑145 (MQ=255)                         aCCGGCCAGCAGGCTGGGCAGCAACTTCAGGGGCAAAGCCAGCAGCAGCCCGGCCATGACGATGATCAGGCTGGCCAGGGTGATCTGACGGGGGGTAAAGGTCATACAGCCTCGACGGCAGACAGCGGGAATGCCGCAGTCTGcc                                                                                                                                <  3:231221/145‑1 (MQ=255)                         aCCGGCCAGCAGGCTGGGCAGCAACTTCAGGGGCAAAGCCAGCAGCAGCCCGGCCATGACGATGATCAGGCTGGCCAGGGTGATCTGACGGGGGGTAAAGGTCATACAGCCTCGACGGCAGACAGCGGGAATGCCGCAGTCTGCCAGcc                                                                                                                            >  8:232520/1‑149 (MQ=255)                           ccgCCAGCAGGCTGGGCAGCAACTTCAGGGGCAAAGCCAGCAGCAGCCCGGCCATGACGATGATCAGGCTGGCCAGGGTGATCTGACGGGGGGTAAAGGTCATACAGCCTCGACGGCAGACTGCGGCATTCCCGCTGTCTGccgtcg                                                                                                                            <  4:165977/145‑5 (MQ=255)                           cGGCCAGCAGGCTGGGCAGCAACTTCAGGGGCAAAGCCAGCAGCAGCCCGGCCATGACGATGATCAGGCTGGCCAGGGTGATCTGACGGGGGGTAAAGGTCATACAGCCTCGACCGCAGCCTGCGGCATTCCCGCTGTCTGccgtcgta                                                                                                                          >  3:165977/1‑143 (MQ=255)                            ggCCAGCAGGCTGGGCAGCAACTTCAGGGGCAAAGCCAGCAGCAGCCCGGCCATGACGATGATCAGGCTGGCCAGGGTGATCTGACGGGGGGTAAAGGTCATACAGCCTCGACGGCAGACAGCGGGAATGCCGCAGTCTGCCAGCCTTg                                                                                                                         <  3:358877/149‑1 (MQ=255)                                agcagGCTGGGCAGCAACTTCAGGGGCAAAGCCAGCAGCAGCCCGGCCATGACGATGATCAGGCTGGCCAGGGTGATCTGACGGGGGGTAAAGGTCATACAGCCTCGACGGCAGACAGCGGGAATGCCGCAGTCTGCCa                                                                                                                               <  3:425548/139‑1 (MQ=255)                                agcagGCTGGGCAGCAACTTCAGGGGCAAAGCCAGCAGCAGCCCGGCCATGACGATGATCAGGCTGGCCAGGGTGATCTGACGGGGGGTAAAGGTCATACAGCCTCGACGGCAGACAGCGGGAATGCCGCAGTCTGCCa                                                                                                                               >  4:425548/1‑139 (MQ=255)                                agcagGCTGGGCAGCAACTTCAGGGGCAAAGCCAGCAGCAGCCCGGCAATGACGATGATCAGGCTGGCCAGGGTGATCTGACGGGGGGTAAAGGTCATACAGCCTCGACGGCAGACAGCGGGAATGCCGCAGTCTGCCAGCCTTggcgg                                                                                                                     >  7:27287/1‑149 (MQ=255)                                     gCTGGGCAGCAACTTCAGGGGCAAAGCCAGCAGCAGCCCGGCCATGACGATGATCAGGCTGGCCAGGGTGATCTGACGGGGGGTAAAGGTCATACAGCCTCGACGGCAGACAGCGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAg                                                                                                                <  7:232520/149‑1 (MQ=255)                                     gCTGGGCAGCAACTTCAGGGGCAAAGCCAGCAGCAGCCCGGCCATGACGATGATCAGGCTGGCCAGGGTGATCTGACGGGGGGTAAAGGTCATACAGCCTCGACGGCAGACAGCGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAg                                                                                                                >  1:169965/1‑149 (MQ=255)                                              caACTTCAGGGGCAAAGCCAGCAGCAGCCCGGCCATGACGATGATCAGGCTGGCCAGGGTGATCTGACGGGGGGTAAAGGTCATACAGCCTCGACGGCAGACAGGGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAGgcg                                                                                                       >  2:209514/1‑149 (MQ=255)                                              caACTTCAGGGGCAAAGCCAGCAGCAGCCCGGCCATGACGATGATCAGGCTGGCCAGGGTGATCTGACGGGGGGTAAAGGTCATACAGCCTCGACGGCAGACAGCGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAGgcg                                                                                                       >  5:225343/1‑149 (MQ=255)                                                  ttCAGGGGCAAAGCCAGCAGCAGCCCGGCCATGACGATGATCAGGCTGGCCAGGGTGATCTGACGGGGGGTAAAGGTCATACAGCCTCGACGGCAGACAGCGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAGGCGCA‑gg                                                                                                  >  2:284688/1‑149 (MQ=255)                                                    cAGGGGCAAAGCCAGCAGCAGCCCGGCCATGACGATGATCAGGCTGGCCAGGGTGATCTGACGGGGGGTAAAGGTCATACAGCCTCGACGGCAGACAGCGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAGGCGCA‑GGt                                                                                                 <  3:298524/148‑1 (MQ=255)                                                       gggCAAAGCCAGCAGCAGCCCGGCCATGACGATGATCAGGCTGGCCAGGGTGATCTGACGGGGGGTAAAGGTCATACAGCCTCGACGGCAGACAGCGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAGGCGCA‑GGTCAtt                                                                                             <  2:169965/149‑1 (MQ=255)                                                         gCAAAGCCAGCAGCAGCCCGGCCATGACGATGATCAGGCTGGCCAGGGTGATCTGACGGGGGGTAAAGGTCATACAGCCTCGACGGCAGACAGCGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAGGCGCA‑GGTCATttc                                                                                           <  2:4107/149‑1 (MQ=255)                                                          cAAAGCCAGCAGCAGCCCGGCCATGACGATGATCAGGCTGGCCAGGGTGATCTGACGGGGGGTAAAGGTCATACAGCCTCGACGGCAGACAGCGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAGGCGCA‑GGTCATttct                                                                                          <  7:409762/149‑1 (MQ=255)                                                           aaaGCCAGCAGCAGCCCGGCCATGACGATGATCAGGCTGGCCAGGGTGATCTGACGGGGGGTAAAGGTCATACAGCCTCGACGGCAGACAGCGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAGGCGCA‑GGTCATttctt                                                                                         <  1:91141/149‑1 (MQ=255)                                                           aaaGCCAGCAGCAGCCCGGCCATGACGATGATCAGGCTGGCCAGGGTGATCTGACGGGGGGTAAAGGTCATACAGCCTCGACGGCAGACAGCGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAGGAGCA‑GGTCATttctt                                                                                         >  6:314968/1‑149 (MQ=255)                                                               ccAGCAGCAGCCCGGCCATGACGATGATCAGGCTGGCCAGGGTGATCTGACGGGGGGTAAAGGTCATACAGCCTCGACGGCAGACAGCGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAGGCGCA‑gg                                                                                                  >  7:151708/1‑136 (MQ=255)                                                               ccAGCAGCAGCCCGGCCATGACGATGATCAGGCTGGCCAGGGTGATCTGACGGGGGGTAAAGGTCATACAGCCTCGACGGCAGACAGCGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAGGCGCA‑gg                                                                                                  <  8:151708/136‑1 (MQ=255)                                                                cagcagcagCCCGGCCATGACGATGATCAGGCTGGCCAGGGTGATCTGACGGGGGGTAAAGGTCATACAGCCTCGACGGCAGACAGCGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAGGCGCA‑GGTCATttcttctt                                                                                      >  6:138048/1‑147 (MQ=255)                                                                cagcagcagCCCGGCCATGACGATGATCAGGCTGGCCAGGGTGATCTGACGGGGGGTAAAGGTCATACAGCCTCGACGGCAGACAGCGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAGGCGCA‑GGTCATTTCTTCTTCa                                                                                    <  2:93289/149‑1 (MQ=255)                                                                 agcagcagcCCGGCCATGACGATGATCAGGCTGGCCAGGGTGATCTGACGGGGGGTAAAGGTCATACAGCCTCGACGGCAGACAGCGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAGGCGCA‑GGTCATTTCTTCTTCAg                                                                                   >  3:390889/1‑149 (MQ=255)                                                                 agcagcagcCCGGCCATGACGATGATCAGGCTGGCCAGGGTGATCTGACGGGGGGTAAAGGTCATACAGCCTCGACGGCAGACAGCGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAGGCGCA‑GGTCATTTCTTCTTCAg                                                                                   <  5:314968/149‑1 (MQ=255)                                                                  ggagcagcCCGGCCATGACGATGATCAGGCTGGCCAGGGTGATCTGACGGGGGGTAAAGGTCATACAGCCTCGACGGCAGACAGCGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAGGCGCA‑GGTCATTTCTTCTTCAgg                                                                                  <  3:23355/147‑1 (MQ=255)                                                                  gcagcagcCCGGCCATGACGATGATCAGGCTGGCCAGGGTGATCTGACGGGGGGTAAAGGTCATACAGCCTCGACGGCAGACAGCGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAGGCGCA‑GGTCATTTCTTCTTCAgg                                                                                  >  8:20599/1‑149 (MQ=255)                                                                               cATGACGATGATCAGGCTGGCCAGGGTGATCTGACGGGGGGTAAAGGTCATACAGCCTCGACGGCAGACAGCGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAGGCGCA‑GGTCATTTCTTCTTCAGGCAATCGCTCATg                                                                      >  3:62705/1‑148 (MQ=255)                                                                                aTGACGATGATCAGGCTGGCCAGGGTGATCTGACGGGGGGTAAAGGTCATACAGCCTCGACGGCAGACAGCGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAGGCGCA‑GGTCATTTCTTCTTCAGGCAATCGCTCATG‑aa                                                                   <  5:53525/149‑1 (MQ=255)                                                                                   aCGATGATCAGGCTGGCCAGGGTGATCTGACGGGGGCTAAAGGTCATACAGCCACGACGGCAGACAGCTGGAATGCCGCAGTCTGCCAGCCTAGGCGGCGCACCCATAGGCGAA‑GGTCATTTCTTCTTCAGGCAATCCCTCATG‑AAgg                                                                 >  3:183299/1‑148 (MQ=255)                                                                                      atgatCAGGCTGGCCAGGGTGATCTGACGGGGGGTAAAGGTCATACAGCCTCGACGGCAGACAGCGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAGGCGCA‑GGTCATTTCTTCTTCAGGCAATCGCTCATG‑AAGGCCtt                                                             <  4:390889/149‑1 (MQ=255)                                                                                           cAGGCTGGCCAGGGTGATCTGACGGGGGGTAAAGGTCATACAGCCTCGACGGCAGACAGCGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAGGCGCA‑GGTCATTTCTTCTTCAGGCAATCGCTCATG‑AAGGCCTTGCGCt                                                        >  2:211611/1‑149 (MQ=255)                                                                                              gCTGGCCAGGGTGATCTGACGGGGGGTAAAGGTCATACAGCCTCGACGGCAGACAGCGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAGGCGCA‑GGTCATTTCTTCTTCAGGCAATCg                                                                            >  8:83531/1‑127 (MQ=255)                                                                                              gCTGGCCAGGGTGATCTGACGGGGGGTAAAGGTCATACAGCCTCGACGGCAGACAGCGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAGGCGCA‑GGTCATTTCTTCTTCAGGCAATCg                                                                            <  7:83531/127‑1 (MQ=255)                                                                                              gCTGGCCAGGGTGATCTGACGGGGGGTAAAGGTCATACAGCCTCGACGGCAGACAGCGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAGGCGCA‑GGTCATTTCTTCTTCAGGCAATCGCTCATG‑AAGGCCTTGCGCTCAt                                                     >  8:319725/1‑149 (MQ=255)                                                                                                  gCCAGGGTGATCTGACGGGGGGTAAAGGTCATACAGCCTCGACGGCAGACAGCGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAGGCGCA‑GGTCATTTCTTCTTCAGGCAATCGCTCATG‑AAGGCCTTGCGCTCATcccc                                                 >  5:171397/1‑149 (MQ=255)                                                                                                  gCCAGGGTGATCTGACGGGGGGTAAAGGTCATACAGCCTCGACGGCAGACAGCGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAGGCGCA‑GGTCATTTCTTCTTCAGGCAATCGCTCATG‑AAGGCCTTGCGCTCATcc                                                   <  2:60055/147‑1 (MQ=255)                                                                                                   ccAGGGTGATCTGACGGGGGGTAAAGGTCATACAGCCTCGACGGCAGACAGCGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAGGCGCA‑GGTCATTTCTTCTTCAGGCAATCGCTCATG‑AAGGCCTTGCGCTCATccccc                                                >  3:50429/1‑149 (MQ=255)                                                                                                    cAGGGTGATCTGACGGGGGGTAAAGGTCATACAGCCTCGACGGCAGACAGCGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAGGCGCA‑GGTCATTTCTTCTTCAGGCAATCGCTCATG‑AAGGCCTTGCGCTCATCCCCCt                                               >  8:141240/1‑149 (MQ=255)                                                                                                    cAGGGTGATCTGACGGGGGGTAAAGGTCATACAGCCTCGACGGCAGACAGCGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAGGCGCA‑GGTCATTTCTTCTTCAGGCAATCGCTCATG‑AAGGCCTTGCGCTCATCCCCCt                                               >  3:88545/1‑149 (MQ=255)                                                                                                     aGGGTGATCTGACGGGGGGTAAAGGTCATACAGCCTCGACGGCAGACAGCGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAGGCGCA‑GGTCATTTCTTCTTCAGGCAATCGCTCATG‑AAGGCCTTGCGCTCATCCCCCtt                                              <  4:338661/149‑1 (MQ=255)                                                                                                             cTGACGGGGGGTAAAGGTCATACAGCCTCGACGGCAGACAGCGGGAATGCCGCAGTCTGTCAGCCTTGGCGGCGCAGGCATAGGCGCA‑GGTCATTTCTTCTTCAGGCAATCGCTCATG‑AAGGCCTTGCGCTCATCCCCCTTCAGGGCt                                       <  7:298834/148‑1 (MQ=255)                                                                                                              tGACGGGGGGTAAAGGTCATACAGCCTCGACGGCAGACAGCGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAGGCGCA‑GGTCATTTCTTCTTCAGGCAATCGCTCATG‑AAGGCCTTGCGCTCATCCCCCTTCAGGGCttt                                     >  2:426218/1‑149 (MQ=255)                                                                                                               gACGGGGGGTAAAGGTCATACAGCCTCGACGGCAGACAGCGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAGGCGCA‑GGTCATTTCTTCTTCAGGCAATCGCTCATG‑AAGGCCTTGCGCTCATCCCCCTTCAGGGCTTtg                                    >  5:309513/1‑149 (MQ=255)                                                                                                               gACGGGGGGTAAAGGTCATACAGCCTCGACGGCAGACAGCGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAGGCGCA‑GGTCATTTCTTCTTCAGGCAATCGCTCATG‑AAGGCCTTGCGCTCATCCCCCTTCAGGGCTTtg                                    >  1:352048/1‑149 (MQ=255)                                                                                                                aCGGGGGGTAAAGGTCATACAGCCTCGACGGCAGACAGCGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAGGCGCA‑GGTCATTTCTTCTTCAGGCAATCGCTCATG‑AAGGCCTTGCGCTCATCCCCCTTCAGGGCTTtgg                                   <  8:27287/149‑1 (MQ=255)                                                                                                                aCGGGGGGTAAAGGTCATACAGCCTCGACGGCAGACAGCGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAGGCGCA‑GGTCATTTCTTCTTCAGGCAATCGCTCATG‑AAGGCCTTGCGCTCATCCCCCTTCAGGGCTTtg                                    <  1:217252/148‑1 (MQ=255)                                                                                                                aCGGGGGGTAAAGGGCATACAGCCTCGACGGCAGACAGCGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAGGCGCA‑GGTCATTTCTTCTTCAGGCAATCGCTCATG‑AAGGCCTTGCGCTCATCCCCCTTCAGGGCttt                                     <  3:85982/147‑1 (MQ=255)                                                                                                                  ggggggTAAAGGTCATACAGCCTCGACGGCAGACAGCGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAGGCGCA‑GGTCATTTCTTCTTCAGGCAATCGCTCATGAAAGGCCTTGCGCTCATCCCCCTTCAGGGCTTtggt                                  <  1:284688/149‑1 (MQ=255)                                                                                                                  ggggggTAAAGGTCATACAGCCTCGACGGCAGACAGCGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAGGCGCA‑GGTCATTTCTTCTTCAGGCAATCGCTCATG‑AAGGCCTTGCGCTCATCCCCCTTCAGGGCTTtggtg                                 >  4:198496/1‑149 (MQ=255)                                                                                                                    ggggTAAAGGTCATACAGCCTCGACGGCAGACAGCGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAGGCGCA‑GGTCATTTCTTCTTCAGGCAATCGCTCATG‑AAGGCCTTGCGCTCATCCCCCTTCAGGGCTTtggtggt                               <  6:309513/149‑1 (MQ=255)                                                                                                                      ggTAAAGGTCATACAGCCTCGACGGCAGACAGCGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAGGCGCA‑GGTCATTTCTTCTTCAGGCAATCGCTCATG‑AAGGCCTTGCGCTCATCCCCCTTCAGGGCTTtggtggtgg                             >  4:228603/1‑149 (MQ=255)                                                                                                                       gTAAAGGTCATACAGCCTCGACGGCAGACAGCGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAGGCGCA‑GGTCATTTCTTCTTCAGGCAATCGCTCATG‑AAGGCCTTGCGCTCATCCCCCTTCAGGGCTTtggtggt                               <  6:229372/146‑1 (MQ=255)                                                                                                                       gTAAAGGTCATACAGCCTCGACGGCAGACAGCGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAGGCGCA‑GGTCATTTCTTCTTCAGGCAATCGCTCATG‑AAGGCCTTGCGCTCATCCCCCTTCAGGGCTTtggtggt                               >  5:229372/1‑146 (MQ=255)                                                                                                                       gTAAAGGTCATACAGCCTCGACGGCAGACAGCGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAGGCGCA‑GGTCAGTTCTACTTCAGGAAATCGCTCATG‑AAGGCCTTGCGCTCATCCCCCTTCAGGGATTTGGTGGTGGc                            >  3:313158/1‑149 (MQ=255)                                                                                                                          aaGGTCATACAGCCTCGACGGCAGACAGCGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAGGCGCA‑GGTCATTTCTTCTTCAGGCAATCGCTCATG‑AAGGCCTTGCGCTCATCCCCCTTCAGGGCTTTGGTGGTGGCGTc                         <  4:50429/149‑1 (MQ=255)                                                                                                                            ggTCATACAGCCTCGACGGCAGACAGCGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAGGCGCA‑GGTCATTTCTTCTTCAGGCAATCGCTCATG‑AAGGCCTTGCGCTCATCCCCCTTCAGGGCTTTGGTGGTGGCGTCAg                       <  5:138048/149‑1 (MQ=255)                                                                                                                            ggTCATACAGCCTCGACGGCAGACAGCGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAGGCGCA‑GGTCATTTCTTCTTCAGGCAATCGCTCATG‑AAGGCCTTGCGCTCATCCCCCTTCAGGGCTTTGGTCGTGGCGTCa                        >  1:240407/1‑148 (MQ=255)                                                                                                                                aTACAGCCTCGACGGCAGACAGCGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAGGCGCA‑GGTCATTTCTTCTTCAGGCAATCGCTCATG‑AAGGCCTTgcgc                                                         >  7:183377/1‑111 (MQ=255)                                                                                                                                aTACAGCCTCGACGGCAGACAGCGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAGGCGCA‑GGTCATTTCTTCTTCAGGCAATCGCTCATG‑AAGGCCTTgcgc                                                         <  8:183377/111‑1 (MQ=255)                                                                                                                                 tACAGCCTCGACGGCAGACAGCGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAGGCGCA‑GGTCATTTCTTCTTCAGGCAATCGCTCATG‑AAGGCCTTGCGCTCATCCCCCTTCAGGGCTTTGGTGGTGGCGTCAGCAtt                   >  8:314105/1‑148 (MQ=255)                                                                                                                                 tACAGCCTCGACGGCAGACAGCGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAGGCGCA‑GGTCATTTCTTCTTCAGGCAATCGCTCATG‑AAGGCCTTGCGCTCATCCCCCTTCAGGGCTTTGGTGGTGGCGTCAGCATTg                  >  7:60432/1‑149 (MQ=255)                                                                                                                                  aCAGCCTCGACGGCAGACAGCGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAGGCGCA‑GGTCATTTCTTCTTCAGGCAATCGCTCATG‑AAGGCCTTGCGCTCATCCCCCTTCAGGGCTTTGGTGGTGGCGTCAGCATTg                  >  6:339158/1‑148 (MQ=255)                                                                                                                                    aGCCTCGACGGCAGACAGCGGGAATGCCGCAGTCTGCCAGCCTTGGAGGCGCAGGCATATGCGCA‑GGTCATTTCTTCTTCAGGCAATCGCTCATG‑AAGGCCATGCGCTCATCACCCTACAGGGCTTTGGTGGTGGCGTCAGCATTGCAg               >  6:281726/1‑149 (MQ=255)                                                                                                                                          gACGGCAGACAGCGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAGGCGCA‑GGTCATTTCTTCTTCAGGCAATCGCTCATG‑AAGGCCTTGCGCTCATCCCCCTTCAGGGCTTTGGTGGTGGCGTCAGCATTGCAGGTcgt          >  5:350753/1‑146 (MQ=255)                                                                                                                                            tggCAGACTGCGGCATTCCCGCTGTCTGCCAGCCTTGGCGGCGCAGGCATAGGCGCAGGGTCATTTCTTCTTCAGGCAATCGCTCATG‑AAGGCCTTGCGCTCATCCCCCTTCAGGGCTTTGGTGGTGGCGTCAGCATTGCAGGTCTTCa        >  1:343192/2‑149 (MQ=255)                                                                                                                                            tggCAGACTGCGGCATTCCCGCTGTCTGCCAGCCTTGGCGGCGCAGGCATAGGCGCA‑GGTCATTTCTTCTTCAGGCAATCGCTCATG‑AAGGCCTTGCGCTCATCCCCCTTCAGGGCTTTGGTGGTGGCGTCAGCATTGCAGGTCTTCAt       >  3:239184/2‑149 (MQ=255)                                                                                                                                              gCAGACAGCGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAGGCGCA‑GGTCATTTCTTCTTCAGGCAATCGCTCATG‑AAGGCCTTGCGCTCATCCCCCTTCAGGGCTTTGGTGGTGGCGTCAGCATTGCAGGTCTTCAttt     >  8:22375/1‑149 (MQ=255)                                                                                                                                               cAGACAGCGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGGAGGCATAGGCGCA‑GGTCATTTCTTCTTCAGGCAATCGCTCATG‑AAGGCCTTGCGCTCATCCCCCTTCAGGGCTTTGGTGGTGGCGTCAGCATTGCAGGTCTTCAtttt    <  7:112240/149‑1 (MQ=255)                                                                                                                                                 gACAGCGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAGGCGCA‑GGTCATTTCTTCTTCAGGCAATCGCTCATG‑AAGGCCTTGCGCTCATCCCCCTTCAGGGCTTTGGTGGTGGCGTCAGCATTGCAGGTCTTCAtttttt  <  2:352048/149‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                                    GTTGACCAGTTCGAACACCAGCAGACCGGCCAGCAGGCTGGGCAGCAACTTCAGGGGCAAAGCCAGCAGCAGCCCGGCCATGACGATGATCAGGCTGGCCAGGGTGATCTGACGGGGGGTAAAGGTCATACAGCCTCGACGGCAGACAGCGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAGGCGCA‑GGTCATTTCTTCTTCAGGCAATCGCTCATG‑AAGGCCTTGCGCTCATCCCCCTTCAGGGCTTTGGTGGTGGCGTCAGCATTGCAGGTCTTCATTTTTT  >  NC_002947/5518392‑5518684

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