Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,518,553	A→T	5.7%	C10S (TGC→AGC)	PP_4850 ←	membrane protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,518,553	0	A	T	5.7%	271.3 / 3.3	86	C10S (TGC→AGC)	PP_4850	membrane protein
Reads supporting (aligned to +/- strand):  ref base A (43/38);  new base T (3/2);  total (46/40)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AACACCAGCAGACCGGCCAGCAGGCTGGGCAGCAACTTCAGGGGCAAAGCCAGCAGCAGCCCGGCCATGACGATGATCAGGCTGGCCAGGGTGATCTGACGGGGGGTAAAGGTCATACAGCCTCGACGGCAGACAGCGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAGGCGCA‑GGTCATTTCTTCTTCAGGCAATCGCTCATG‑AAGGCCTTGCGCTCATCCCCCTTCAGGGCTTTGGTGGTGGCGTCAGCATTGCAGGTCTTCATTTTTTCTTGCTGGGTTTGGGGC  >  NC_002947/5518405‑5518701                                                                                                                                                     |                                                                                                                                                       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)         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                                                                                                                   <  8:151708/136‑1 (MQ=255)                                                  ccAGCAGCAGCCCGGCCATGACGATGATCAGGCTGGCCAGGGTGATCTGACGGGGGGTAAAGGTCATACAGCCTCGACGGCAGACAGCGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAGGCGCA‑gg                                                                                                                   >  7:151708/1‑136 (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                                                                                             <  7:83531/127‑1 (MQ=255)                                                                                 gCTGGCCAGGGTGATCTGACGGGGGGTAAAGGTCATACAGCCTCGACGGCAGACAGCGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAGGCGCA‑GGTCATTTCTTCTTCAGGCAATCg                                                                                             >  8:83531/1‑127 (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                                                >  5:229372/1‑146 (MQ=255)                                                                                                          gTAAAGGTCATACAGCCTCGACGGCAGACAGCGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAGGCGCA‑GGTCATTTCTTCTTCAGGCAATCGCTCATG‑AAGGCCTTGCGCTCATCCCCCTTCAGGGCTTtggtggt                                                <  6:229372/146‑1 (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                                                                          <  8:183377/111‑1 (MQ=255)                                                                                                                   aTACAGCCTCGACGGCAGACAGCGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAGGCGCA‑GGTCATTTCTTCTTCAGGCAATCGCTCATG‑AAGGCCTTgcgc                                                                          >  7:183377/1‑111 (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)                                                                                                                                         cGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAGGCGCA‑GGTCATTTCTTCTTCAGGCAATCGCTCATG‑AAGGCCTTGCGCTCATCCCCCTTCAGGGCTTTGGTGGTGGCGTCAGCATTGCAGGTATTCATGTTTTCTTGc              >  6:217919/1‑149 (MQ=255)                                                                                                                                           ggAATGCCGCAGTCTGCCAGCCTTGCCGGCGCAGGCATAGGCGCA‑GGTCATTTCTTCTTCAGGCAATCGCTCATG‑AAGGCCTTGCGCTCATCCCCCTTCATGGCTTTGGTGGTGGCGTCAGCATTGCAGGTCTTCATTTTTTCTTGCTg            <  4:62705/149‑1 (MQ=255)                                                                                                                                            gAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAGGCGCA‑GGTCATTTCTTCTTCAGGCAATCGCTCATG‑AAGGCCTTGCGCTCATCCCCCTTCAGGGCTTTGGTGGTGGCGTCAGCATTGCAGGTCTTCATTTTTTCTTGCTgg           >  5:96161/1‑149 (MQ=255)                                                                                                                                                 ccGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAGGCGCA‑GGTCATTTCTTCTTCAGGCAATCGCTCATG‑AAGGCCTTGCGCTCATCCCCCTTCAGGGCTTTGGTGGTGGCGTCAGCATTGCAGGTCTTCATTTTTTCTTGCTGGGTTTg      <  7:20599/149‑1 (MQ=255)                                                                                                                                                  cGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAGGCGCA‑GGTCATTTCTTCTTCAGGCAATCGCTCATG‑AAGGCCTTGCGCTCATCCCCCTTCAGGGCTTTGGTGGTGGCGTCAGCATTGCAGGTCTTCATTTTTTCTTGCTGGGTTTgg     <  6:350753/149‑1 (MQ=255)                                                                                                                                                   gCAGTCTGCCAGCCTTGGCGGCGCAGGCATAGGCGCA‑GGTCATTTCTTCTTCAGGCAATCGCTCATG‑AAGGCCTTGCGCTCATCCCCCTTCAGGGCTTTGGTGGTGGCGTCAGCATTGCAGGTCTTCAttt                      >  6:330832/1‑131 (MQ=255)                                                                                                                                                   gCAGTCTGCCAGCCTTGGCGGCGCAGGCATAGGCGCA‑GGTCATTTCTTCTTCAGGCAATCGCTCATG‑AAGGCCTTGCGCTCATCCCCCTTCAGGGCTTTGGTGGTGGCGTCAGCATTGCAGGTCTTCAttt                      <  5:330832/131‑1 (MQ=255)                                                                                                                                                   gCAGTCTGCCAGCCTTGGCGGCGCAGGCATAGGCGCA‑GGTCATTTCTTCTTCAGGCAATCGCTCATG‑AAGGCCTTGCGCTCATCCCCCTTCAGGGCTTTGGTGGTGGCGTCAGCATTGCAGGTCTTCATTTTTTCTTGCTGGGTTTggg    <  4:88545/149‑1 (MQ=255)                                                                                                                                                   gCAGTCTGCCAGCCTTGGCGGCGCAGGCATAGGCGCA‑GGTCATTTCTTCTTCAGGCAATCGCTCATG‑AAGGCCTTGCGCTCATCCCCCTTCAGGGCTTTGGTGGTGGCGTCAGCATTGCAGGTCTTCATTTTTTCTTGCTGGGTTTggg    <  8:60432/149‑1 (MQ=255)                                                                                                                                                    cAGTCTGCCAGCCTTGGCGGCGCAGGCATAGGCGCA‑GGTCATTTCTTCTTCAGGCAATCGCTCATG‑AAGGCCTTGCGCTCATCCCCCTTCAGGGCTTTGGTGGTGGCGTCAGCATTGCAGGTCTTCATTTTTTCTTGCTGGGTTTgggg   <  3:198496/149‑1 (MQ=255)                                                                                                                                                    cAGTCTGCCAGCCTTGGCGGCGCAGGCATAGGCGCA‑GGTCATTTCTTCTTCAGGCAATCGCTCATG‑AAGGCCTTGCGCTCATCCCCCTTCAGGGCTTTGGTGGTGGCGTCAGCATTGAAGGTCTTCATTTTTTCTTGCTGGGTTTgggg   >  6:187878/1‑149 (MQ=255)                                                                                                                                                     aGTCTGCCAGCCTTGGCGGCTGAGGGATAGGCGCA‑GGTCATTTCTTCTTCAGGCAATCGCTCATG‑AAGGCCTTTCGCTCATCCCCCTTCAGGGCTTTGGTGGTGGCGTCATCATTGCAGGTCTTCATTTTTTCTTGCTGGGTTTGGGGc  <  4:183299/149‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                       AACACCAGCAGACCGGCCAGCAGGCTGGGCAGCAACTTCAGGGGCAAAGCCAGCAGCAGCCCGGCCATGACGATGATCAGGCTGGCCAGGGTGATCTGACGGGGGGTAAAGGTCATACAGCCTCGACGGCAGACAGCGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAGGCGCA‑GGTCATTTCTTCTTCAGGCAATCGCTCATG‑AAGGCCTTGCGCTCATCCCCCTTCAGGGCTTTGGTGGTGGCGTCAGCATTGCAGGTCTTCATTTTTTCTTGCTGGGTTTGGGGC  >  NC_002947/5518405‑5518701

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: ‑