Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,218,755	G→T	11.6%	A181E (GCG→GAG)	PP_3698 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,218,755	0	G	T	11.6%	162.8 / 3.5	62	A181E (GCG→GAG)	PP_3698	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (30/24);  new base T (6/2);  total (36/26)
Fisher's exact test for biased strand distribution p-value = 4.50e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CAGAACCTCCCCTCTCGTGCACTGATCACTGATAATCGATTGCCGAAGCAGGCGGTCAGC‑ATTACGCCATGAGTCCTCAACTCTTCGGAAATCAGCCCCCAGGCTCGTCGCCAGCGCGTAGATCTCGCGGATGCCTTCAATCGCTTTGGCTTCATCAAAGATCGGTACATGGCCGGCTTTGGCCCGGTCAAGCTCGATCTGCAGTTCACGGATCCTGCGCTTCAAGCTGAGTGCGCGGCTCTTGGGGTCTGAGTTGAGGTTCGTCTCGAGATTTTCGATCTCACGCTGCA  >  NC_002947/4218613‑4218902                                                                                                                                                |                                                                                                                                                    cAGAACCTCCCCTCTCGTGCACTGATCACTGATAATCGATTGCCGAAGCAGGCGGTCAGCAATAACGCCATTAGTCCTAAACTGTTCGGAAATAAGCCCCCAGGCTCGTCCCCAGCGCGTATATATCGCGGTTGCCTTAAATCGCGTTg                                                                                                                                                >  1:125029/1‑149 (MQ=255)  aGAACCTCCCCTCTCGTGCACTGATCACTGATAATCGATTGCCGAAGCAGGCGGTCAGC‑ATTACGCCATGAGTCCTCAACTCTTCGGAAATCAGCCCCCAGGCTCGTCGCCAGCGCGTAGATCTCGCGGATGCCTTCAATCGCTTTGGc                                                                                                                                              <  3:92286/149‑1 (MQ=255)  aGAACCTCCCCTCTCGTGCACTGATCACTGATAATCGATTGCCGAAGCAGGCGGTCAGC‑ATTACGCCATGAGTCCTCAACTCTTCGGAAATCAGCCCCCAGGCTCGTCGCCAGCGCGTAGATCTCGCGGATGCCTTCAATCGCTTTGGc                                                                                                                                              >  2:142066/1‑149 (MQ=255)     aCCTCCCCTCTCGTGCACTGATCACTGATAATCGATTGCCGAAGCAGGCGGTCAGC‑ATTACGCCATGAGTCCTCAACTCTTCGGAAATCAGCCCCCAGGCTCGTCGCCAGCGCGTAGATCTCGCGGATGCCTTCAATCGCTTTGGCTtc                                                                                                                                           <  4:387892/149‑1 (MQ=255)     aCCTCCCCTCTCGTGCACTGATCACTGATAATCGATTGCCGAAGCAGGCGGTCAGC‑ATTACGCCATGAGTCCTCAACTCTTCGGAAATCAGCCCCCAGGCTCGTCGCCAGCGCGTAGATCTCGCGGATGCCTTCAATCGCTTTGGCTtc                                                                                                                                           <  2:195800/149‑1 (MQ=255)     aCCTCCCCTCTCCTGCACTGATCACTGATCATCGATTGCCGAAGCAGGCGGGCAGC‑ATTACGCCATGAGTCCTCAACTCTTCGGAAATCAGCCCCCAGGCTCGTCGCCAGCGCGTAGATCTCGCGGATGCCTTCAATCGCTTTGGCTtc                                                                                                                                           <  6:445110/149‑1 (MQ=255)        tCCCCTCTCGTGCACTGATCACTGATAATCGATTGCCGAAGCAGGCGGTCAGC‑ATTACGCCATGAGTCCTCAACTCTTCGGAAATCAGCCCCCAGGCTCGTCGCCAGCGCGTAGATCTCGCGGATGCCTTCAATCGCTTTGGCTtcat                                                                                                                                         >  5:242611/1‑148 (MQ=255)        tCCCCTCTCGTGCACTGATCACTGATAATCGATTGCCGAAGCAGGCGGTCAGC‑ATTACGCCATGAGTCCTCAACTCTTCGGAAATCAGCCCCCAGGCTCGTCGCCAGCGCGTAGATCTCGCGGATGCCTTCAATCGCTTTGGCTtc                                                                                                                                           <  4:9228/146‑1 (MQ=255)        tCCCCTCTCGTGCACTGATCACTGATAATCGATTGCCGAAGCAGGCGGTCAGC‑ATTACGCCATGAGTCCTCAACTCTTCGGAAATCAGCCCCCAGGCTCGTCGCCAGCGCGTAGATCTCGCGGATGCCTTCAATCGCTTTGGCTtc                                                                                                                                           >  3:9228/1‑146 (MQ=255)             tctcGTGCACTGATCACTGATAATCGATTGCCGAAGCAGGCGGTCAGC‑ATTACGCCATGAGTCCTCAACTCTTCGGAAATCAGCCCCCAGGCTCGTCGCCAGCGCGTAGATCTCGCGGATGCCTTCAATCGCTTTGGCTtcatca                                                                                                                                       >  7:83602/1‑145 (MQ=255)                 gTGCACTGCACCCTGATATACGCTGGCCGAAGCAGGCTGGCAGC‑ATTACGCCCTGAGTCCTCATCTCTTCGGTAATCAGCCCCCACGCTCGTCGCCAGCGCGTAGAGCTCGCCGATGCCTTCTATCTCTTTGGCTTCATCAAAGATCg                                                                                                                                <  8:98430/148‑1 (MQ=255)                 gTGCACTGATCACTGATAATCGATTGCCGAAGCAGGCGGTCAGC‑ATTACGCCATGAGTCCTCAACTCTTCGGAAATCAGCCCCCAGGCTCGTCGCCAGCGCGTAGATCTCGCGGATGCCTTCAATCGCTTTGGCTTCATCAAAGATCgg                                                                                                                               >  7:376615/1‑149 (MQ=255)                 gTGCACTGATCACTGATAATCGATTGCCGAAGCAGGCGGTCAGC‑ATTACGCCATGAGTCCTCAACTCTTCGGAAATCAGCCCCCAGGCTCGTCGCCAGCGCGTAGATCTCGCGGATGCCTTCAATAGCTTTGGCTtcatca                                                                                                                                       >  7:98430/1‑141 (MQ=255)                  tGCACTGATCACTGATAATCGTTTGCCGAAGCAGGCGGTCAGC‑ATTACGCTATGAGTCCTCAACTCTTCGGAAATCAGCCCCCAGGCTCGTCGCCAGCGCGTAGATCTCGCGGATGCCTTCAATCGCTTTGGCTTCATCAAAGATCGGt                                                                                                                              <  6:99132/149‑1 (MQ=255)                   gCACTGATCACTGATAATCGATTGCCGAAGCAGGCGGTCAGC‑ATTACGCCATGAGTCCTCAACTCTTCGGAAATCAGCCCCCAGGCTCGTCGCCAGCGCGTAGATCTCGCGGATGCCTTCAATCGCTTTGGCTTCATCaaa                                                                                                                                     >  3:323235/1‑141 (MQ=255)                   gCACTGATCACTGATAATCGATTGCCGAAGCAGGCGGTCAGC‑ATTACGCCATGAGTCCTCAACTCTTCGGAAATCAGCCCCCAGGCTCGTCGCCAGCGCGTAGATCTCGCGGATGCCTTCAATCGCTTTGGCTTCATCaaa                                                                                                                                     <  4:323235/141‑1 (MQ=255)                      ctgatcactgatAATCGATTGCCGAAGCAGGCGGTCAGC‑ATTACGCCATGAGTCCTCAACTCTTCGGAAATCAGCCCCCAGGCTCGTCGCCAGCGCGTAGATCTCGCGGATGCCTTCAATCGCTTTGGCTTCATCAAAGATCGGTACAt                                                                                                                          >  2:104408/1‑149 (MQ=255)                          tcactgatAATCGATTGCCGAAGCAGGCGGTCAGC‑ATTACGCCATGAGTCCTCAACTCTTCGGAAATCAGCCCCCAGGCTCGTCGCCAGCGCGTAGATCTCGCGGATGCCTTCAATCGCTTTGGCTTCATCAAAGATCGGTACATGGc                                                                                                                       <  5:69686/148‑1 (MQ=255)                               gatAATCGATTGCCGAAGCAGGCGGTCAGC‑ATTACGCCATGAGTCCTCAACTCTTCGGAAATCAGCCCCCAGGCTCGTCGCCAGCGCGTAGATCTCGCGGATGCCTTCAATCGCTTTgg                                                                                                                                               <  2:298264/119‑1 (MQ=255)                               gatAATCGATTGCCGAAGCAGGCGGTCAGC‑ATTACGCCATGAGTCCTCAACTCTTCGGAAATCAGCCCCCAGGCTCGTCGCCAGCGCGTAGATCTCGCGGATGCCTTCAATCGCTTTgg                                                                                                                                               >  1:298264/1‑119 (MQ=255)                                  aaTCGATTGCCGAAGCAGGCGGTCAGC‑ATTACGCCATGAGTCCTCAACTCTTCGGAAATCAGCCCCCAGGCTCGTCGCCAGCGCGTAGATCTCGCGGATGCCTTCAATCGCTTTGGCTTCATCAAAGATCGGTACATGGCCGGCTTTgg                                                                                                              >  8:201517/1‑149 (MQ=255)                                   aTCGATTGCCGAAGCAGGCGGTCAGC‑ATTACGCCATGAGTCCTCAACTCTTCGGAAATCAGCCCCCAGGCTCGTCGCCAGCGCGTAGATCTCGCGGCTGCCTTCAATCGCTTTGGCTTCATCAAAGATCGGTACATGGCCGGCTTTGGc                                                                                                             >  8:84570/1‑149 (MQ=255)                                          gCCGAAGCAGGCGGTCAGC‑ATTACGCCATGAGTCCTCAACTCTTCGGAAATCAGCCCCCAGGCTCGTCGCCAGCGCGTAGATCTCGCGGATGCCTTCAATCGCTTTGGCTTCATCAAAGATCGGTACATGGcc                                                                                                                      <  5:400761/133‑1 (MQ=255)                                          gCCGAAGCAGGCGGTCAGC‑ATTACGCCATGAGTCCTCAACTCTTCGGAAATCAGCCCCCAGGCTCGTCGCCAGCGAGTAGATCTCGCTGATCCCATAAAAATATAAAAAAACATCAAAGATCGGTACATGGcc                                                                                                                      >  6:400761/1‑133 (MQ=255)                                           ccGAAGCAGGCGGTCAGC‑ATTACGCCATGAGTCCTCAACTCTTCGGAAATCAGCCCCCAGGCTCGTCGCCAGCGCGTAGATCTCGCGGATGCCTTCAATCGCTTTGGCTTCATCAAAGATCGGTACATGGCCGGCTTTGGCCCGGTCa                                                                                                      <  6:306966/148‑1 (MQ=255)                                             gAAGCAGGCGGTCAGC‑ATTACGCCATGAGTCCTCAACTCTTCGGAAATCAGCCCCCAGGCTCGTCGCCAGCGAGTAGATCTCGCGGATGCCTTCAATATATTTTACTTCATCAAAGATCGGTACATGGCCGGCTTTGGCCCGGTCAAGc                                                                                                   >  2:89976/1‑149 (MQ=255)                                                  aGGCGGTCAGC‑ATTACGCCATGAGTCCTCAACTCTTCGGAAATCAGCCCCCAGGCTCGTCGCCAGCGCGTAGATCTCGCGGATGCCTTCAATCGCTTTGGCTTCATCAAAGATCGGTACATGGCCGGCTTTGGCCCGGTCAAGCTCGAt                                                                                              <  4:276507/149‑1 (MQ=255)                                                       gTCAGC‑ATTACGCCATGAGTCCTCAACTCTTCGGAAATCAGCCCCCAGGCTCGTCGCCAGCGCGTAGATCTCGCGGATGCCTTCAATCGCTTTGGCTTCATCAAAGATCGGTACATGGCCGGCTTTGGCCCGGTCAAGCTCGATCTGCa                                                                                         <  2:125029/149‑1 (MQ=255)                                                           gC‑ATTACGCCATGAGTCCTCAACTCTTCGGAAATCAGCCCCCAGGCTCGTCGCCAGCGCGTAGATCTCGCGGATGCCTTCAATCGCTTTGGCTTCATCAAAGATCGGTACATGGCCGGCTTTGGCCCGCTCAAGCTCGATCTGCAGtt                                                                                      >  1:19011/1‑148 (MQ=255)                                                                tACGCCATGAGTCCTCAACTCTTCGGAAATCAGCCCCCAGGCTCGTCGCCAGCGAGTAGATCTCGCTGATCCCATAAAAATAAAAATAAAAATCAAAGATCGGTACATGGCCGGCTATGGCCCGGTCAAGCTCGATCTGCAGTTCACg                                                                                  >  8:194633/1‑148 (MQ=255)                                                                  cGCCATGAGTCCTCAACTCTTCGGAAATCAGCCCCCAGGCTCGTCGCCAGCGCGTAGATCTCGCGGATGCCTTCAATCGCTTTGGCTTCATCAAAGATCGGTACATGGCCGGCTTTGGCCCGGTCAAGCTCGATCTGCAGTTCACGGa                                                                                <  6:242611/148‑1 (MQ=255)                                                                     cATGAGTCCTCAACTCTTCGGAAATCAGCCCCCAGGCTCGTCGCCAGCGAGTAGATCTCGCGGATGCCATAAAACTATTTATCTACATCAAAGATCGGTACATGGCCGGCTTTGGCCCGGTCAAGCTCGATCTGCAGTTCACGGCTCCt                                                                            >  6:107002/1‑149 (MQ=255)                                                                      aTGAGTCCTCAACTCTTCGGAAATCAGCCCCCAGGCTCGTCGCCAGCGCGTAGATCTCGCGGATGCCTTCAATCGCTTTGGCTTAATAAAAGATCGGTACATGGCCGGCTTTGGCCCGGTCAAGCTCGATCTGCAGTTCACGGATCCTg                                                                           >  7:88211/1‑149 (MQ=255)                                                                       tGAGTCCTCAACTCTTCGGAAATCAGCCCCCAGGCTCGTCGCCAGCGCGTAGATCTCGCGGATGCCTTCAATCGCTTTGGCTTCATCAAAGATCGGTACATGGCCGGCTTTGGCCCGGTCAAGCTCGATCTGCAGTTCACGGATCCTgc                                                                          <  7:201517/149‑1 (MQ=255)                                                                        gAGTCCTCAACTCTTCGGAAATCAGCCCCCAGGCTCGTCGCCAGCGCGTAGATCTCGCGGATGCCTTCAATCGCTTTGGCTTCATCAAAGATCg                                                                                                                                <  1:377811/94‑1 (MQ=255)                                                                        gAGTCCTCAACTCTTCGGAAATCAGCCCCCAGGCTCGTCGCCAGCGCGTAGATCTCGCGGATGCCTTCAATCGATTTATCTTCATCAAAGATCg                                                                                                                                >  2:377811/1‑94 (MQ=255)                                                                         aGTCCTCAACTCTTCGGAAATCAGCCCCCAGGCTCGTCGCCAGCGCGTAGATCTCGCGGATGCCTTCAATCGCTTTGGCTTCATCAAAGATCGGTACATGGCCGGCTTTGGCCCGGTCAAGCTCGATCTGCAGTTCACGGATCCTgcg                                                                         <  8:376615/148‑1 (MQ=255)                                                                           tCCTCAACTCTTCGGAAATCAGCCCCCAGGCTCGTCGCCAGCGCGTAGATCTCGCGGATGCCTTCAATCGCTTTGGCTTCATCAAAGATCGGTACATGGCCGGCTTTGGCCCGGTCAAGCTCGATCTGCAGTTCACGGATCCTGCGCt                                                                       >  7:99334/1‑148 (MQ=255)                                                                           tCCTCAACTCTTCGGAAATCAGCCCCCAGGCTCGTCGCCAGCGCGTAGATCTCGCGGATGCCTTCAATCGCTTTGGCTTCATCAAAGATCGGTACATGGCCGGCTTTGGCCCGGTCAAGCTCGATCTGCAGTTCACGGATCCTGCGCt                                                                       >  4:265730/1‑148 (MQ=255)                                                                            ccTCAACTCTTCGGAAATCAGCCCCCAGGCTCGTCGCCAGCGAGTAGATCTCGCGGATGCCTTCAATCGCTTTGGCTTCATCAAAGATCGGTACATGGCCGGCTTTGGCCCGGTCAAGCTCGATCTGCAGTTCACGGATCCTGCGCtt                                                                      >  2:205576/1‑148 (MQ=255)                                                                              tCAACTCTTCGGAAATCAGCCCCCAGGCTCGTCGCCAGCGCGTAGATCTCGCGGATGCCTTCAATCGCTTTGGCTTCATCAAAGATCGGTACATGGCCGGCTTTGGCCCGGTCAAGCTCGATCTGCAGTTCACGGATCCTGCGCTTc                                                                     >  3:382125/1‑147 (MQ=255)                                                                                aaCTCTTCGGAAATCAGCCCCCAGGCTCGTCGCCAGCGCGTAGATCTCGCGGATGCCTTCAATCGCTTTGGCTTCATCAAAGATCGGTACATGGCCGGCTTTGGCCCGGTCAAGCTCGATCTGCAGTTCACGGATCCTGCGCTTCaa                                                                   <  4:382125/147‑1 (MQ=255)                                                                                   tctTCGGAAATCAGCCCCCAGGCTCGTCGCCAGCGAGTAGATCTCGCGGATCCCATCAAAAGATTTAAATAAATCAAAGATAGGTACATGGCCGTCTTTGGCCCGGTCAAGATCTAACAGCAGTACACGGATCATGCGCTTCAAGCt                                                                >  8:405990/1‑147 (MQ=255)                                                                                     ttCGGAAATCAGCCCCCAGGCTCGTCGCCAGCGCGTAGATCTCGCGGATGCCTTCAATCGCTTTGGCTTCATCAAAGATCGGTACATGGCCGGCTTTGGCCCGGTCAAGCTCGATCTGCAGTTCACGGATCCTGCGCTTCAAGCTGa                                                              <  1:239546/147‑1 (MQ=255)                                                                                       cGGAAATCAGCCCCCAGGCTCGTCGCCAGCGCGTAGATCTCGCGGATGCCTTCAATCGCTTTGGCTTCATCAAAGATCGGTACATGGCCGGCTTTGGCCCGGTCAAGCTCGATCTGCAGTTCACGGATCCTGCGCTTCAAGCTGAGTgc                                                          >  7:110560/1‑149 (MQ=255)                                                                                       cGGAAATCAGCCCCCAGGCTCGTCGCCAGCGAGTAGATCTCGCGGATGCCATAAATATATTTATCTTAATCAAATATCGGTACATGGCCGGATTTGGCCCGTTCAAGCTCGATCTGCAGTTCACGGATCCTGCGCTTCAAGCTCAGTgc                                                          >  8:336766/1‑149 (MQ=255)                                                                                         gAAATCAGCCCCCAGGCTCGTCGCCAGCGCGTAGATCTCGCGGATGCCTTCAATCGCTTTGGCTTCATCAAAGATCGGTACATGGCCGGCTTTGGCCCGGTCAAGCTCGATCTGCAGTTCACGGATCCTGCGCTTCAAGCTGAGTgcgc                                                        <  8:110560/149‑1 (MQ=255)                                                                                             tCAGCCCCCAGGCTCGTCGCCAGCGAGTAGATCTCGCTGATCCCATAAAAATAAATAAAATAATAAAAGATCGGTACATGGCCGGCTTTGGCCCGGTCAAGCTCGATCTGCAGTTCACGGATCCTGCGCTTCAAGCTGAGTGCGCgga                                                     >  4:343834/1‑147 (MQ=255)                                                                                                   cccAGGCTCGTCGCCAGCGCGTAGATCTCGCGGATGCCTTCAATCGCTTTGGCTTCATCAAAGATCGGTACATGGCCGGCTTTGGCCCGGTCAAGCTCGATCTGCAGTTCACGGATCCTGCGCTTCAAGCTGAGTGCGCGGCTCTTggg                                              >  3:122861/1‑149 (MQ=255)                                                                                                         cTCGTCGCCAGCGCGTAGATCTCGCGGATGCCTTCAATCGCTTTGGCTTCATCAAAGATCGGTACATGGCCGGCTTTGGCCCGGTCAAGCTCGATCTGCAGTTCACGGATCCTGCGCTTCAAGCTGAGTGCGCGGCTCTTGGGGTCTGa                                        <  3:265730/149‑1 (MQ=255)                                                                                                           cgtcgCCAGCGCGTAGATCTCGCGGATGCCTTCAATCGCTTTGGCTTCATCAAAGATCGGTACATGGCCGGCTTTGGCCCGGTCAAGCTCGATCTGCAGTTCACGGATCCTGCGCTTCAAGCTGAGTGCGCGGCTCTTGGGGTCTGAGt                                      >  3:378213/1‑149 (MQ=255)                                                                                                                        tAGATCTCGCGGATGCCTTCAATCGCTTTGGCTTCATCAAAGATCGGTACATGGCCGGCTTTGGCCCGGTCAAGCTCGATCTGCAGTTCACGGATCCTGCGCTTCAAGCTGAGTGCGCGGCTCTTGGGGTCTGAGTTGAGGTTCGtcgc                         >  7:380166/1‑147 (MQ=255)                                                                                                                          gATCTCGCGGATGCCTTCAATCGCTTTGGCTTCATCAAAGATCGGTACATGGCCGGCTTTGGCCCGGTCAAGCTCGATCTGCAGTTCACGGATCCTGCGCTTCAAGCTGAGTGCGCGGCTCTTGGGGTCTGAGTTGAGGTTCGTCTCga                       <  1:205576/149‑1 (MQ=255)                                                                                                                           aTCTCGCGGATGCCTTCAATCGCTTTGGCTTCATCAAAGATCGGTACATGGCCGGCTTTGGCCCGGTCAAGCTCGATCTGCAGTTCACGGATCCTGCGCTTCAAGCTGAGTGCGCGGCTCTTGGGGTCTGAGTTGAGGTTCGtctc                         >  6:295483/1‑146 (MQ=255)                                                                                                                           aTCTCGCGGATGCCTTCAATCGCTTTGGCTTCATCAAAGATCGGTACATGGCCGGCTTTGGCCCGGTCAAGCTCGATCTGCAGTTCACGGATCCTGCGCTTCAAGCTGAGTGCGCGGCTCTTGGGGTCTGAGTTGAGGTTCGTCACga                       >  1:374773/1‑148 (MQ=255)                                                                                                                              tcGCGGATGCCTTCAATCGCTTTGGCTTCATCAAAGATCGGTACATGGCCGGCTTTGGCCCGGTCAAGCTCGATCTGCAGTTCACGGATCCTGCGCTTCAAGCTGAGTGCGCGGCTCTTGGGGTCTGAGTTGAGGTTCGTCTCGAGAtt                   >  3:132758/1‑149 (MQ=255)                                                                                                                                gcgGATGCCTTCAATCGCTTTGGCTTCATCAAAGATCGGTACATGGCCGGCTTTGGCCCGGTCAAGCTCGATCTGCAGTTCACGGATCCTGCGCTTCAAGCTGAGTGCGCGGCTCTTGGGGTCTGAGTTGAGGTTCGTCTCGAGAtttt                 >  4:122248/1‑149 (MQ=255)                                                                                                                                 cgGATGCCTTCAATCTCTTTGGCTTCATCAAAGATCGGTACATGGCCGGCTTTGGCCCGGTCAAGCTCGATCTGCAGTTCACGGATCCTGCGCTTCAAGCTGAGTGCGCGGCTCTTGGGGTCTGAGTTGAGGTTCGTCTCGAGATTTTc                <  8:83602/149‑1 (MQ=255)                                                                                                                                          tCAATCGCTTTGGCTTCATCAAAGATCGGTACATGGCCGGCTTTGGCCCGGACAAGCTCGATCTGCAGTTCACGGATCCTGCGCTTCAAGCTGAGTGCGCGGCTCTTGGGGTCTGAGTTGAGGTTCGTCTCGAGATTTTCGATATCACg       >  3:58814/1‑149 (MQ=255)                                                                                                                                            aaTCGCTTTGGCTTCATCAAAGATCGGTACATGGCCGGCTTTGGCCCGGTCAAGCTCGATCTGCAGTTCACGGATCCTGCGCTTCAAGCTGAGTGCGCGGCTCTTGGGGTCTGAGTTGAGGTTCGTCTCGAGATTTTCGATCTCACGc      <  2:374773/148‑1 (MQ=255)                                                                                                                                               cGCTTTGGCTTCATCAAAGATCGGTACATGGCCGGCTTTGGCCCGGTCAAGCTCGATCTGCAGTTCACGGATCCTGCGCTTCAAGCTGAGTGCGCGGCTCTTGGGGTCTGAGTTGAGGTTCGTCTCGAGATTTTCg               <  6:293448/136‑1 (MQ=255)                                                                                                                                               cGCTTTGGCTTCATCAAAGATCGGTACATGGCCGGCTTTGGCCCGGTCAAGCTCGATCTGCAGTTCACGGATCCTGCGCTTCAAGCTGAGTGCGCGGCTCTTGGGGTCTGAGTTGAGGTTCGTCTCGAGATTTTCg               >  5:293448/1‑136 (MQ=255)                                                                                                                                                gCTTTGGCTTCATCAAAGATCGGTACATGGCCGGCTTTGGCCCGGTCAAGCTCGATCTGCAGTTCACGGATCCTGCGCTTCAAGCTGAGTGCGCGGCTCTTGGGGTCTGAGTTGAGGTTCGTCTCGAGATTTTCGATCTCACGCTGca  <  4:122861/148‑1 (MQ=255)                                                                                                                                                |                                                                                                                                                    CAGAACCTCCCCTCTCGTGCACTGATCACTGATAATCGATTGCCGAAGCAGGCGGTCAGC‑ATTACGCCATGAGTCCTCAACTCTTCGGAAATCAGCCCCCAGGCTCGTCGCCAGCGCGTAGATCTCGCGGATGCCTTCAATCGCTTTGGCTTCATCAAAGATCGGTACATGGCCGGCTTTGGCCCGGTCAAGCTCGATCTGCAGTTCACGGATCCTGCGCTTCAAGCTGAGTGCGCGGCTCTTGGGGTCTGAGTTGAGGTTCGTCTCGAGATTTTCGATCTCACGCTGCA  >  NC_002947/4218613‑4218902

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