Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,518,546	A→T	6.0%	I12N (ATT→AAT)	PP_4850 ←	membrane protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,518,546	0	A	T	6.0%	258.2 / 4.9	83	I12N (ATT→AAT)	PP_4850	membrane protein
Reads supporting (aligned to +/- strand):  ref base A (44/34);  new base T (3/2);  total (47/36)
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.46e-01

GTTCGAACACCAGCAGACCGGCCAGCAGGCTGGGCAGCAACTTCAGGGGCAAAGCCAGCAGCAGCCCGGCCATGACGATGATCAGGCTGGCCAGGGTGATCTGACGGGGGGTAAAGGTCATACAGCCTCGACGGCAGACAGCGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAGGCGCA‑GGTCATTTCTTCTTCAGGCAATCGCTCATG‑AAGGCCTTGCGCTCATCCCCCTTCAGGGCTTTGGTGGTGGCGTCAGCATTGCAGGTCTTCATTTTTTCTTGCTGG  >  NC_002947/5518400‑5518692                                                                                                                                                   |                                                                                                                                                     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                                                                                                                                       >  4:425548/1‑139 (MQ=255)                        agcagGCTGGGCAGCAACTTCAGGGGCAAAGCCAGCAGCAGCCCGGCCATGACGATGATCAGGCTGGCCAGGGTGATCTGACGGGGGGTAAAGGTCATACAGCCTCGACGGCAGACAGCGGGAATGCCGCAGTCTGCCa                                                                                                                                       <  3:425548/139‑1 (MQ=255)                        agcagGCTGGGCAGCAACTTCAGGGGCAAAGCCAGCAGCAGCCCGGCAATGACGATGATCAGGCTGGCCAGGGTGATCTGACGGGGGGTAAAGGTCATACAGCCTCGACGGCAGACAGCGGGAATGCCGCAGTCTGCCAGCCTTggcgg                                                                                                                             >  7:27287/1‑149 (MQ=255)                             gCTGGGCAGCAACTTCAGGGGCAAAGCCAGCAGCAGCCCGGCCATGACGATGATCAGGCTGGCCAGGGTGATCTGACGGGGGGTAAAGGTCATACAGCCTCGACGGCAGACAGCGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAg                                                                                                                        >  1:169965/1‑149 (MQ=255)                             gCTGGGCAGCAACTTCAGGGGCAAAGCCAGCAGCAGCCCGGCCATGACGATGATCAGGCTGGCCAGGGTGATCTGACGGGGGGTAAAGGTCATACAGCCTCGACGGCAGACAGCGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAg                                                                                                                        <  7:232520/149‑1 (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                                                       >  3:88545/1‑149 (MQ=255)                                                                                            cAGGGTGATCTGACGGGGGGTAAAGGTCATACAGCCTCGACGGCAGACAGCGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAGGCGCA‑GGTCATTTCTTCTTCAGGCAATCGCTCATG‑AAGGCCTTGCGCTCATCCCCCt                                                       >  8:141240/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                                                                 >  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)                                                                                                                                              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)                                                                                                                                                   |                                                                                                                                                     GTTCGAACACCAGCAGACCGGCCAGCAGGCTGGGCAGCAACTTCAGGGGCAAAGCCAGCAGCAGCCCGGCCATGACGATGATCAGGCTGGCCAGGGTGATCTGACGGGGGGTAAAGGTCATACAGCCTCGACGGCAGACAGCGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAGGCGCA‑GGTCATTTCTTCTTCAGGCAATCGCTCATG‑AAGGCCTTGCGCTCATCCCCCTTCAGGGCTTTGGTGGTGGCGTCAGCATTGCAGGTCTTCATTTTTTCTTGCTGG  >  NC_002947/5518400‑5518692

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