Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,796,837	A→T	10.2%	intergenic (+21/‑111)	PP_t43 → / → PP_2449	tRNA‑Val/hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,796,837	0	A	T	10.2%	164.7 / 7.9	59	intergenic (+21/‑111)	PP_t43/PP_2449	tRNA‑Val/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (24/29);  new base T (3/3);  total (27/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.29e-01

AAAAGAGCTTTACAGGTACGCATCGGGCCACTAATATGCGCACCACACGACAGGCACATAGCTCAGTTGGTTAGAGCACCACCTTGACATGGTGGGGGTCGTTGGTTCGAGTCCAATTGTGCCTACCAAATCGAAAAAGGGCGATCAGAAATGATCGCCCTTTTTGC‑ATTTGCGCCCGAAGCAGCGTGTAACGTTTGGTCGCGGTCCGGGCATACCCCGGCGTGCTAGAATCCGCCCCTTCGAATCTGGAGGTACACGCGTGCGCAAACTGGCTGTGGTAATGGCAGTGCTCGCCCT  >  NC_002947/2796689‑2796985                                                                                                                                                     |                                                                                                                                                      aaaaGAGCGTTACAGCTACGCCACGGGCCACTAATATGCGCACCACACGACAGGCACATAGCTCAGTTGGTTAGAGCACCACCTTGACATGGTGGGGGTCGTTGGTTCGAGTCCAATTGTGCCTACCAAATCGAAAAAGGGCGATCAGa                                                                                                                                                       <  8:262641/149‑1 (MQ=255)        cTTTACAGGTACGCATCGGGCCACTAATATGCGCACCACACGACAGGCACATAGCTCAGTTGGTTAGAGCACCACCTTGACATGGTGGGGGTCGTTGGTTCGAGTCCAATTGTGCCTACCAAATCGAAAAAGGGCGATCAGAAATGATc                                                                                                                                                <  4:23684/149‑1 (MQ=255)         tttACAGGTACGCATCGGGCCACTAATATGCGCACCACACGACAGGCACATAGCTCAGTTGGTTAGAGCACCACCTTGACATGGTGGGGGTCGTTGGTTCGAGTCCAATTGTGCCTACCAAATCGAAAAAGGGCGATCAGAAATGATCg                                                                                                                                               >  7:251808/1‑149 (MQ=255)            aCAGGTACGCATCGGGCCACTAATATGCGCACCACACGACAGGCACATAGCTCAGTTGGTTAGAGCACCACCTTGACATGGTGGGGGTCGTTGGTTCGAGTCCAATTGTGCCTACCAAATCGAAAAAGGGCGATCATTTCTGATCGcc                                                                                                                                             <  6:133456/148‑1 (MQ=255)              aGGTACGCATCGGGCCACTAATATGCGCACCACACGACAGGCACATAGCTCAGTTGGTTAGAGCACCACCTTGACATGGTGGGGGTCGTTGGTTCGAGTCCAATTGTGCCTACCAAATCGAAAAAGGGCGATCAGAAATGATCGCCCtt                                                                                                                                          >  6:317364/1‑149 (MQ=255)                gTACGCATCGGGCCACTAATATGCGCCCCACACGACAGGCACATAGCTCAGTTGGATAGAGCACCACCTTGACATGGTGGGGGTCGTTGGTTCGAGTCCAATTGTGCCTACCAAATCGAAAAAGGTCGATCAGAAATGATCGCCCtttt                                                                                                                                        >  5:266476/1‑149 (MQ=255)                 tACGCATCGGGCCACTAATATGCGCACCACACGACAGGCACATAGCTCAGTTGGTTAGAGCACCACCTTGACATGGTGGGGGTCGTTGGTTCGAGTCCAATTGTGCCTACCAAATCGAAAAAGGGCGATCAGAAATGATCGCCCttttt                                                                                                                                       <  1:95531/149‑1 (MQ=255)                      aTCTGGCCACTAATATGCGCACAACACGACAGGCACATAGCTCAGATGGTTAGAGCACCACCATGAAATGGTGGGGGTCGTTGGTTCGAGTCAAATTGTGCCTACCAAATCGAAAAAAGGCGATCAGAAATTATAGAACTTTTTGA‑Att                                                                                                                                 >  8:114235/1‑149 (MQ=17)                         gggCCACTAATATGCTCACCACACGACAGGCACATAGCTCAGTTGGTTAGTGCACCACCTTGACATGGTGGGGGTCGTTGGTTCGAGTCCAATTGTGCCTACCAAATCGAAAAAGGGCGATCAGAAATGATCGCCCTTTTTGC‑ATTTgc                                                                                                                              <  8:296918/149‑1 (MQ=255)                             cACTAATATGCGCACCACACGACAGGCACATAGCTCAGTTGGTTAGAGCACCACCTTGACATGGTGGGGGTCGTTGGTTCGAGTCCAATTGTGCCTACCAAATGCAAAAAGGGCGATCATTTCTGATCGCCCTTTTTCG‑ATTTggtagg                                                                                                                          <  5:97888/149‑6 (MQ=14)                                 ccgatGCGTACCACACGACAGGCACATAGCTCAGTTGGTTAGAGCACCACCTTGACATGGTGGGGGTCGTTGGTTCGAGTCCAATTGTGCCTACCAAATCGAAAAAGGGCGATCAGAAATGATCGCCCTTTTTGC‑ATTTGCGCCCGAag                                                                                                                      <  5:324204/146‑1 (MQ=255)                                          accacACGACAGGCACATAGCTCAGTTGGTTAGAGCACCACCTTGACATGGTGGGGGTCGTTGGTTCGAGTCCAATTGTGCCTACCAAATCGAAAAAGGGCGATCAGAAATGATCGCCCTTTTTGC‑ATTTGCGCCCGAAGCAGCgtgt                                                                                                              >  5:129889/1‑148 (MQ=255)                                             acacGACAGGCACATAGCTCAGTTGGTTAGAGCACCACCTTGACATGGTGGGGGTCGTTGGTTCGAGTCCAATTGTGCCTACCAAATCGAAAAAGGGCGATCAGAAATGATCGCCCTTTTTGC‑ATTTGCGCCCGAAGCAGCGTGTAACg                                                                                                          <  1:161947/149‑1 (MQ=255)                                                 gACAGGCACATAGCTCAGTTGGTTAGAGCACCACCTTGACATGGTGGGGGTCGTTGGTTCGAGTCCAATTGTGCCTACCAAATCGAAAAAGGGCGATCAGAAATGATCGCCCTTTTTGC‑ATTTGCGCCCGAAGCAGCGTGTAACGTTTg                                                                                                      <  1:158517/149‑1 (MQ=255)                                                 gACAGGCACATAGCTCAGTTGGTTAGAGCACCACCTTGACATGGTGGGGGTCGTTGGTTCGAGTCCAATTGTGCCTACCAAATCGAAAAAGGGCGATCAGAAATGATCGCCCTTTTTGC‑ATTTGCGCCCGAAGCAGCGTGTAACGTTTg                                                                                                      >  2:250509/1‑149 (MQ=255)                                                  aCAGTCACATAGCTCAGTTGGTTAGAGCACCACCTTGACATGGTGGGGGTCGTTGGTTCGAGTCCAATTGTGCCTACCAAATCGAAAAAGGGCGATCAGAAATGATCGCCCTTTTTGC‑ATTTGCGCCCGAAGCAGCGTGTAACGTTTgg                                                                                                     <  4:213714/149‑1 (MQ=255)                                                  aCAGGCACATAGCTCAGTTGGTTAGAGCACCACCTTGACATGGTGGGGGTCGTTGGTTCGAGTCCAATTGTGCCTACCAAATCGAAAAAGGGCGATCAGAAATGATCGCCCTTTTTGC‑ATTTGCGCCCGAAGCAGCGTGTAACGTTTgg                                                                                                     <  3:72195/149‑1 (MQ=255)                                                          aTAGCTCAGTTGGTTAGAGCACCACCTTGACATGGTGGGGGTCGTTGGTTCGAGTCCAATTGTGCCTACCAAATCGAAAAAGGGCGATAAGAAATGATCGCCCTTTTTGC‑ATTTGCGCCCGAAGCAGCGTGTAACGTTTGGTCGCGGTc                                                                                             >  8:355219/1‑149 (MQ=255)                                                          aTAGCTCAGTTGGTTAGAGCACCACCTTGACATGGTGGGGGTCGTTGGCTCGAGTCCAATTGTGCCTACCAAATCGAAAAAGGGCGATCATTTCTGATCGCCCt                                                                                                                                           <  1:169600/104‑1 (MQ=14)                                                          aTAGCTCAGTTGGTTAGAGCACCACCTTGACATGGTGGGGGTCGTTGGCTCGAGTCCAATTGTGCCTACCAAATCGAAAAAGGGCGATCATTTCTGATCGCCCt                                                                                                                                           >  2:169600/1‑104 (MQ=14)                                                            aGCTCAGTTGGTTAGAGCACCACCTTGACATGGTGGGGGTCGTTGGTTCGAGTCCAATTGTGCCTACCAAATCGAAAAAGGGCGATCAGAAATGATCGCCCTTTTTGC‑ATTTGCGCCCg                                                                                                                         <  3:372338/119‑1 (MQ=255)                                                            aGCTCAGTTGGTTAGAGCACCACCTTGACATGGTGGGGGTCGTTGGTTCGAGTCCAATTGTGCCTACCAAATCGAAAAAGGGCGATCAGAAATGATCGCCCTTTTTGC‑ATTTGCGCCCg                                                                                                                         >  4:372338/1‑119 (MQ=255)                                                              cTCAGTTGGTTAGAGCACCACCTTGACATGGTGGGGGTCGTTGGTTCGAGTCCAATTGTGCCTACCAAATCGAAAAAGGGCGATCAGAAATGATCGCCCTTTTTGC‑ATTTGCGCCCGAAGCAGCGTGTAACGTTTGGTCGCGGTCCggg                                                                                         <  8:2419/149‑1 (MQ=255)                                                                 aGTTGGTTAGAGCACCACCTTGACATGGTGGGGGTCGTTGGTTCGAGTCCAATTGTGCCTACCAAATCGAAAAAGGGCGATCAGAAATGATCGCCCTTTTTGC‑ATTTGCGCCCGAAGCAGCGTGTAACGTTTGGTCGCGGTCCGGGc                                                                                        <  6:45436/147‑1 (MQ=255)                                                                  gTTGGTTAGAGCACCACCTTGACATGGTGGGGGTCGTTGGTTCGAGTCCAATTGTGCCTACCAAATCGAAAAAGGGCGCTCAGAAATGATCGCCCTTTTTGC‑ATTTGCGCCCGACGCGGCGTGTATCGTGTGGTCGCGGTCCGGGCaga                                                                                     >  4:160413/1‑147 (MQ=21)                                                                         agagCACCACCTTGACATGGTGGGGGTCGTTGGTTCGAGTCCAATTGTGCCTACCAAATCGAAAAAGGGCGATCAGAAATGATCGCCCTTTTTGC‑ATTTGCGCCCGAAGCAGCGTGTAACGTTTGGTCGCGGTCCGGGCATACCCCGGc                                                                              >  1:86916/1‑149 (MQ=255)                                                                                  ccTTGACATGGTGGGGGTCGTTGGTTCGAGTCCAATTGTGCCTACCAAATCGAAAAAGGGCGATCAGAAATGATCGCCCTTTTTGC‑ATTTGCGCCCGAAGCAGCGTGTAACGTTTGGTCGCGGTCCGGGCATACCCCGGCGTGCTAGaa                                                                     >  7:90731/1‑149 (MQ=255)                                                                                      gACATGGTGGGGGTCGTTGGTTCGAGTCCAATTGTGCCTACCAAATCGAAAAAGGGCGATCAGAAATGATCGCCCTTTTTGC‑ATTTGCGCCCGAAGCAGCGTGTAACGTTTGGTCGCGGTCCGGGCATACCCCGGCGTGCTAGAATc                                                                   <  1:250509/147‑1 (MQ=255)                                                                                       aCATGGTGGGGGTCGTTGGTTCGAGTCCAATTGTGCCTACCAAATCGAAAAAGGGCGATCAGAAATGATCGCCCTTTTTGC‑ATTTGCGCCCGAAGCAGCGTGTAACGTTTGGTCGCGGTCCGGGCATACCCCGGCGTGCTAGAATCCg                                                                 <  5:317364/148‑1 (MQ=255)                                                                                       aCATGGTGGGGGTCGTTGGTTCGAGTCCAATTGTGCCTACCAAATCGAAAAAGGGCGATCAGAAATGATCGCCCTTTTTGC‑ATTTGCGCCCGAAGCAGCGTGTAACGTTTGGTCGCGGTCCGGGCATACCCCGGCGTGCTAGAATCCGc                                                                <  6:224221/149‑1 (MQ=255)                                                                                        cATGGTGGGGGTCGTTGGTTCGAGTCCAATTGTGCCTACCAAATCGAAAAAGGGCGATCAGAAATGATCGCCCTTTTTGC‑ATTTGCGCCCGAAGCAGCGTGTAACGTTTGGTCGCGGTCCGGGCATACCCCGGCGTGCTAGAATCCGcc                                                               <  7:355219/149‑1 (MQ=255)                                                                                            gtggGGGTCGTTGGTTCGAGTCCAATTGTGCCTACCAAATCGAAAAAGGGCGATCAGAAATGATCGCCCTTTTTGC‑ATTTGCGCCCGAAGCAGCGTGTAACGTTTGGTCGCGGTCCGGGCATACCCCGGCGTGCTAGAATCCGCCCCtt                                                           >  5:371825/1‑149 (MQ=255)                                                                                               ggggTCGTTGGTTCGAGTCCAATTGTGCCTACCAAATCGAAAAAGGGCGATCAGAAATGATCGCCCTTTTTGC‑ATTTGCGCCCGAAGCAGCGTGTAACGTTTGGTCGCGGTCCGGGCATACCCCGGCGTGCTa                                                                        <  1:359167/133‑1 (MQ=255)                                                                                               ggggTCGTTGGTTCGAGTCCAATTGTGCCTACCAAATCGAAAAAGGGCGATCAGAAATGATCGCCCTTTTTGC‑ATTTGCGCCCGAAGCAGCGTGTAACGTTTGGTCGCGGTCCGGGCATACCCCGGCGTGCTa                                                                        >  2:359167/1‑133 (MQ=255)                                                                                                     gTTGGTTCGAGTCCAATTGTGCCTACCAAATCGAAAAAGGGCGATCAGAAATGATCGCCCTTTTTGC‑ATTTGCGCCCGAAGCAGCGTGTAACGTTTGGTCGCGGTCCGGGCATACCCCGGCGTGCTAGAATCCGCCCCTTCGAATCTg                                                   <  1:150661/148‑1 (MQ=255)                                                                                                        ggTTCGAGTCCAATTGTGCCTACCAAATCGAAAAAGGGCGATCAGAAATGATCGCCGTTTTTGC‑ATTTGCGCCCGAAGCAGCGTGTAACTTTTGGTCGCGGTCCGGGCATACCCCGGCGTGCTAGAATCCGCCCCTTCTAATCTGGAg                                                <  5:4126/148‑1 (MQ=255)                                                                                                         gTTCGAGTCCAATTGTGCCTACCAAATCGAAAAAGGGCGATCAGAAATGATCGccc                                                                                                                                            <  5:73206/56‑1 (MQ=255)                                                                                                         gTTCGAGTCCAATTGTGCCTACCAAATCGAAAAAGGGCGATCAGAAATGATCGccc                                                                                                                                            >  6:73206/1‑56 (MQ=255)                                                                                                                tCCAATTGTGCCTACCAAATCGAAAAAGGGCGATCAGAAATGATCGCCCTTTTTGC‑ATTTGCGCCCGAAGCAGCGTGTAACGTTTGGTCGCGGTCCGGGCATACCCCGGCGTGCTAGAATCCGCCCCTTCGAATCTGGAGGTACAcgcg                                       >  4:62/1‑149 (MQ=255)                                                                                                                  cAATTGTGCCTACCAAATCGAAAAAGGGCGATCAGAAATGATCGCCCTTTTTGC‑ATTTGCGCCCGAAGCAGCGTGTAACGTTTGGTCGCGGTCCGGGCATACCCCGGCGTGCTAGAATCCGCCCCTTCGAATCTGGAGGTACACGCGTg                                     >  1:63843/1‑149 (MQ=255)                                                                                                                     ttGTGCCTACCAAATCGAAAAAGGTCGATCAGAAATGCTCGCCCTTTTTGC‑ATTTGCGCCCGAAGCAGCGTGTAACGTTTGGTCGCGGTCCGGGCATACCCCGGCGTGCTAGAATCCGCCCCTTCGAATCTGGAGGTACACGCGTgcgc                                  <  6:266476/149‑1 (MQ=255)                                                                                                                           cTACCAAATCGAAAAAGGGCGATCAGAAATGATCGCCCTTTTTGC‑ATTTGCGCCCGAAGCAGCGTGTAACGTTTGGTCGCGGTCCGGGCATACCCCGGCGTGCTAGAATCCGCCCCTTCGAATCTGGAg                                                >  5:348857/1‑129 (MQ=255)                                                                                                                           cTACCAAATCGAAAAAGGGCGATCAGAAATGATCGCCCTTTTTGC‑ATTTGCGCCCGAAGCAGCGTGTAACGTTTGGTCGCGGTCCGGGCATACCCCGGCGTGCTAGAATCCGCCCCTTCGAATCTGGAg                                                <  6:348857/129‑1 (MQ=255)                                                                                                                            tACCAAATCGAAAAAGGGCGATCAGAAATGATCGCCCTTTTTGC‑ATTTGCGCCCGAAGCAGCGTGTAACGTTTGGTCGCGGTCCg                                                                                           <  8:137507/85‑1 (MQ=255)                                                                                                                            tACCAAATCGAAAAAGGGCGATCAGAAATGATCGCCCTTTTTGC‑ATTTGCGCCCGAAGCAGCGTGTAACGTTTGGTCGCGGTCCg                                                                                           >  7:137507/1‑85 (MQ=255)                                                                                                                            tACCAAATCGAAAAAGGGCGATCAGAAATGATCGCCCTTTTTGC‑ATTTGCGCCCGAAGCAGCGTGTAACGTTTGGTCGCGGTCCGGGCATACCCCGGCGTGCTAGAATCCGCCCCTTCGAATCTGGAGGTACACGCGTGCGCAAACTgg                           >  7:117034/1‑149 (MQ=255)                                                                                                                             aCCAAATCGAAAAAGGGCGATCAGAAATGATCGCCCTTTTTGC‑ATTTGCGCCCGAAGCAGCGTGTAACGTTTGGTcgcg                                                                                                <  5:51323/79‑1 (MQ=255)                                                                                                                             aCCAAATCGAAAAAGGGCGATCAGAAATGATCGCCCTTTTTGC‑ATTTGCGCCCGAAGCAGCGTGTAACGTTTGGTcgcg                                                                                                >  6:51323/1‑79 (MQ=255)                                                                                                                              ccAAATCGAAAAAGGGCGATCAGAAATGATCGCCCTTTTTGC‑ATTTGCGCCCGAAGCAGCGTGTAACGTTTGGTCGCGGTCCGGGCATACCCCGGCGTGCTAGAATCCGCCCCTTCGAATCTGGAGGTACACGCGTGCGCAAACTgg                           <  4:169906/147‑1 (MQ=255)                                                                                                                              ccAAATCGAAAAAGGGCGATCAGAAATGATCGCCCTTTTTGC‑ATTTGCGCCCGAAGCAGCGTGTAACGTTTGGTCGCGGTCCGGGCATACCCCGGCGTGCTAGAATCCGCCCCTTCGAATCTGGAGGTACACGCGTGCGCAAACTGGc                          <  2:86916/148‑1 (MQ=255)                                                                                                                                      aaaaaGGGCGATCAGAAATGATCGCCCTTTTTGC‑ATTTGCGCCCGAAGCAGCGTGTAACGTTTGGTCGCGGTCCGGGCATACCCCGGCGTGCTAGAATCCGCCCCTTCGAATCTGGAGGTACACGCGTGCGCAAACTGGCTGTGGTAAt                 <  6:371825/149‑1 (MQ=255)                                                                                                                                       aaaaGGGCGATCAGAAATGATCGCCCTTTTTGC‑ATTTGCGCCCGAAGCAGCGTGTAACGTTTGGTCGCTGTCCGGGCATACCCCGGCGTGCTAGAATCCGCCCCTTCGAATCTGGAGGTACACGCGTGCGCAAACTGGCTGTGGTAATg                <  8:251808/149‑1 (MQ=255)                                                                                                                                          aGGGCGATCATTTCTGATCGCCCTTTTTGC‑ATTTGCGCCCGAAGCAGCGTGTAACGTTTGGTCGCGGTCCGGGCATACCCCGGCGTGCTAGAATCCGCCCCTTCGAATCTGGAGGTACACGCGTGCGCAAACTGGCTGTGGTAATGGCa             >  4:206070/1‑149 (MQ=255)                                                                                                                                          aGGGCGATCAGAAATGATCGCCCTTTTTGC‑ATTTGCGCCCGAAGCAGCGTGTAACGTTTGGTCGCGGTCCGGGCATACCCCGGCGTGCTAGAATCCGCCCCTTCGAATCTGGAGGTACACGCGTGCGCAAACTGGCTGTGGTAATGGCa             <  6:129889/149‑1 (MQ=255)                                                                                                                                           gggCGATCATTTCTGATCGCCCTTTTTGCGATTTGCGCCCGAAGCAGCGTGTAACGTTTGGTCGCGGTCCGGGCATACCCCGGCGTGCTAGAATCCGCCCCTTCGAATCTGGAGGTACACGCGTGCGCAAACTGGCTGTGGTAATgg               >  1:28537/1‑147 (MQ=255)                                                                                                                                           gggCGATCAGAAATGATCGCCCTTTTTGC‑ATTTGCGCCCGAAGCAGCGTGTAACGTTTGGTCGCGGTCCGGGCATACCCCGGCGTGCTAGAATCCGCCCCTTCGAATCTGGAGGTACACGCGTGCGCAAACTGGCTGTGGTAATGGCAg            >  6:245362/1‑149 (MQ=255)                                                                                                                                              cGATCAGAAATGATCGCCCTTTTTGC‑ATTTGCGCCCGAAGCAGCGTGTAACGTTTGGTCGCGGTCCGGGCATACCCCGGCGTGCTAGAATCCGCCCCTTCGAATCTGGAGGTACACGCGTGCGCAAACTGGCTGTGGTAATGGCAGTGc         >  6:256182/1‑149 (MQ=255)                                                                                                                                                aTCAGAAATGATCGCCCTTTTTGC‑ATTTGCGCCCGAAGCAGCGTGTAACGTTTGGTCGCGGTCCGGGCATACCCCGGCGTGCTAGAATCCGCCCCTTCGAATCTGGAGGTACACGCGTGCGCAAACTGGCTGTGGTAATGGCa             <  3:344937/143‑1 (MQ=255)                                                                                                                                                aTCAGAAATGATCGCCCTTTTTGC‑ATTTGCGCCCGAAGCAGCGTGTAACGTTTGGTCGCGGTCCGGGCATACCCCGGCGTGCTAGAATCCGCCCCTTCGAATCTGGAGGTACACGCGTGCGCAAACTGGCTGTGGTAATGGCa             >  4:344937/1‑143 (MQ=255)                                                                                                                                                 tCAGAAATGATCGCCCTTTTTGC‑ATTTGCGCCCGAAGCAGCGTGTAACGTTTGGTCGCGGTCCGGGCATACCCCGGCGTGCTAGAATCCGCCCCTTCGAATCTGGAGGTACACGCGTGCGCAAACTGGCTGTGGTAATGGCAGTGCTCg      >  5:223883/1‑149 (MQ=255)                                                                                                                                                     aaaTGATCGCCCTTTTTGC‑ATTTGCGCCCGAAGCAGCGTGTAACGTTTGGTCGCGGTCCGGGCATACCCCGGCGTGCTAGAATCCGCCCCTTCGAATCTGGAGGTACACGCGTGCGCAGACTGGCTGTGGTAATGGCAGTGCTCGCCCt  >  5:203181/1‑149 (MQ=255)                                                                                                                                                     aaaTGATCGCCCTTTTTGC‑ATTTGCGCCCGAAGCAGCGTGTAACGTTTGGTCGCGGTCCGGGCATACCCCGGCGTGCTAGAATCCGCCCCTTCGAATCTGGAGGTACACGCGTGCGCAAACTGGCTGTGGTAATGGCAGTGCTCGCCCt  >  7:148488/1‑149 (MQ=255)                                                                                                                                                     aaaTGATCGCCCTTTTTGC‑ATTTGCGCCCGAAGCAGCGTGTAACGTTTGGTCGCGGTCCGGGCATACCCCGGCGTGCTAGAATCCGCCACTTCGAATCTGGAGGTACACGCGTGCGCAAACTGGCTGTGGTAATGGCAGTGCTCGCCCt  >  5:359294/1‑149 (MQ=255)                                                                                                                                                     |                                                                                                                                                      AAAAGAGCTTTACAGGTACGCATCGGGCCACTAATATGCGCACCACACGACAGGCACATAGCTCAGTTGGTTAGAGCACCACCTTGACATGGTGGGGGTCGTTGGTTCGAGTCCAATTGTGCCTACCAAATCGAAAAAGGGCGATCAGAAATGATCGCCCTTTTTGC‑ATTTGCGCCCGAAGCAGCGTGTAACGTTTGGTCGCGGTCCGGGCATACCCCGGCGTGCTAGAATCCGCCCCTTCGAATCTGGAGGTACACGCGTGCGCAAACTGGCTGTGGTAATGGCAGTGCTCGCCCT  >  NC_002947/2796689‑2796985

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