Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,209,223	A→G	5.9%	A392A (GCA→GCG)	PP_3692 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,209,223	0	A	G	5.9%	206.8 / 2.4	68	A392A (GCA→GCG)	PP_3692	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (40/24);  new base G (2/2);  total (42/26)
Fisher's exact test for biased strand distribution p-value = 6.33e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.97e-01

CCGCACTTGAGCTAGTCAACTTGAGTGAAAAAGCCGCGAGAACACGCAAAGCCAGTGAAGAGGCTGCGGCCGCTATTTTCGGTAAGGCAGAGCTGGAGGGCGTTGGCTCATCGACGTGGCAGCAGATGTGGGAGTACGCTCGCGCATA‑CTCCGCCGCCGTGGCGTATCCTGAGTCACCCTTCCCCAACGTGTCAGCTGAGTCACGATGCGTGCTTTGCCATCAGGATCTCGACGATCCCGCAAAAGCACGCCTGGCCGAATTTGAAAAGTTTGTCACTGCTGGTTTGGAGAC  >  NC_002947/4209078‑4209369                                                                                                                                                  |                                                                                                                                                    ccGCACTTGAGCTAGTCAACTTGAGTGAAAAAGCCGCGAGAACACGCAAAGCCAGTGAAGAGGCTGCGGCCGCTATTTTCGGTAAGGCAGAGCTGGAGGGCGTTGGCTCATCGACGTGGCAGCAGATGTGGGAGTACGCTCGCGCATA‑c                                                                                                                                                 >  8:115477/1‑149 (MQ=255) ccGCACTTGAGCTAGTCAACTTGAGTGAAAAAGCCGCGAGAACACGCAAAGCCAGTGAAGAGGCTGCGGCCGCTATTTTCGGTAAGGCAGAGCTGGAGGGCGTTGGCTCATCGACGTGGCAGCAGATGTGGGAGTACGCTCGCGCATA‑c                                                                                                                                                 >  2:269829/1‑149 (MQ=255)  cGCACTTGAGCTAGTCAACTTGAGTGAAAAAGCCGCGAGAACACGCAAAGCCAGTGAAGAGGCTGCGGCCGCTATTTTCGGTAAGGCAGAGCTGGAGGGCGTTGGCTCATCGACGTGGCAGCAGATGTGGGAGTACGCTCGCGCATA‑Ct                                                                                                                                                >  6:127214/1‑149 (MQ=255)         gAGCTAGTCAACTTGAGTGAAAAAGCCGCGAGAACACGCAAAGCCAGTGAAGAGGCTGCGGCCGCTATTTTCGGTAAGGCAGAGCTGGAGGGCGTTGGCTCATCGACGTGGCAGCAGATGTGGGAGTACGCTCGCGCATA‑CTccgccgc                                                                                                                                         >  1:222886/1‑149 (MQ=255)               gTCAACTTGAGTGAAAAAGCCGCGAGAACACGCAAAGCCAGTGAAGAGGCTGCGGCCGCTATTTTCGGTAAGGCAGAGCTGGAGGGCGTTGGCTCATCGACGTGGCAGCAGATGTGGGAGTACGCTCGCGCATA‑CTGCGCCGCCGTGGc                                                                                                                                   >  1:129072/1‑149 (MQ=255)                     ttGAGTGAAAAAGCCGCGAGAACACGCAAAGCCAGTGAAGAGGCTGCGGCCGCTATTTTCGGTAAGGCAGAGCTGCAGGGCGTTGGCTCATCGACGTGGCAGCAGATGTGGGAGTACGCTCGCGCATA‑CTCCGCCGCCATGGCGtctac                                                                                                                             >  6:284512/1‑145 (MQ=255)                          tGAAAAAGCCGCGAGAACACGCAAAGCCAGTGAAGAGGCTGCGGCCGCTATTTTCGGTAAGGCAGAGCTGGAGGGCGTTGGCTCATCGACGTGGCAGCAGATGTGGGAGTATGCGCGAGCGTA‑CTcccaca                                                                                                                                          <  1:228147/131‑5 (MQ=255)                          tGAAAAAGCCGCGAGAACACGCAAAGCCAGTGAAGAGGCTGCGGCCGCTATTTTCGGTAAGGCAGAGCTGGAGGGCGTTGGCTCATCGACGTGGCAGCAGATGTGGGAGTATGCGCGAGCGTA‑CTcccaca                                                                                                                                          >  2:228147/1‑127 (MQ=255)                          tGAAAAAGCCGCGAGAACACGCAAAGCCAGTGAAGAGGCTGCGGCCGCTATTTTCGGTAAGGCAGAGCTGGAGGGCGTTGGCTCATCGACGTGGCAGCAGATGTGGGAGTACGCTCGCGCATA‑CTCCGCCGCCGTGGCGTATCCTGAGt                                                                                                                        >  2:82401/1‑149 (MQ=255)                             aaaaGCCGCGAGAACACGCAAAGCCAGTGAAGAGGCTGCGGCCGCTATTTTCGGTAAGGCAGAGCTGGAGGGCGTTGGCTCATCGACGTGGCAGCAGATGTGGGAGTACGCTCGCGCATA‑CTCCGCCGCCGTGGCGTATCCTGAGTCAc                                                                                                                     <  1:269829/149‑1 (MQ=255)                               aaGCCGCGAGAACACGCAAAGCCAGTGAAGAGGCTGCGGCCGCTATTTTCGGTAAGGCAGAGCTGGAGGGCGTTGGCTCATCGACGTGGCAGCAGATGTGGGAGTACGCTCGCGCATA‑CTCCGCCGCCGTGGCGTATCCTGAGTCAccc                                                                                                                   >  5:280683/1‑149 (MQ=255)                                 gCCGCGAGAACACGCAAAGCCAGTGAAGAGGCTGCGGCCGCTATTTTCGGTAAGGCAGAGCTGGAGGGCGTTGGCTCATCGACGTGGCAGCAGATGTGGGAGTACGCTCGCGCATA‑CTCCGCCGCCGTGGCGTATCCTGAGTCACCCtt                                                                                                                 <  1:135375/149‑1 (MQ=255)                                             cGCAAAGCCAGTGAAGAGGCTGCGGCCGCTATTTTCGGTAAGGCAGAGCTGGAGGGCGTTGGCTCATAGACGTGGCAGCAGATGTGGGAGTACGCTCGCGCATA‑CTCCGCCGCCGTGGCGTATCCTGAGTCCCCCTTCCACATCgtgt                                                                                                      >  4:55017/1‑148 (MQ=255)                                              gCAAAGCCAGTGAAGAGGCTGCGGCCGCTATTTTCGGTAAGGCAGAGCTGGAGGGCGTTGGCTCATCGACGTGGCAGCAGATGTGGGAGTACGCTCGCGCATA‑CTCCGCCGCCGTGGCGTATCCTGAGTCACCCTTcc                                                                                                               <  8:88499/138‑1 (MQ=255)                                              gCAAAGCCAGTGAAGAGGCTGCGGCCGCTATTTTCGGTAAGGCAGAGCTGGAGGGCGTTGGCTCATCGACGTGGCAGCAGATGTGGGAGTACGCTCGCGCATA‑CTCCGCCGCCGTGGCGTATCCTGAGTCACCCTTcc                                                                                                               >  7:88499/1‑138 (MQ=255)                                               cAAAGCCAGTGAAGAGGCTGCGGCCGCTATTTTCGGTAAGGCATAGCTGGAGGGCGTTGGCTCATCGACGTGGCAGCAGATGTGGGAGTACGCTCGCGCATACCTCCGCCGCCGTGGCGTATCCTGAGTCACCCTTCCCCAACGTGTCa                                                                                                    <  1:82401/149‑1 (MQ=255)                                                aaaGCCAGTGAAGAGGCTGCGGCCGCTATTTTCGGTAAGGCAGAGCTGGAGGGCGTTGGCTCATCGACGTGGCAGCAGATGTGGGAGTACGCTCGCGCATA‑CTCCGCCGCCGTGGCGTATCCTGAGTCACCCTTCCCCAACGTGTCAGc                                                                                                  <  2:222886/149‑1 (MQ=255)                                                      aGTGAAGAGGCTGCGGCCGCTATTTTCGGTAAGGCAGAGCTGGAGGGCGTTGGCTCATCGACGTGGCAGCAGATGTGGGAGTACGCTCGCGCATA‑CTCCGCCGCCGTGGCGTATCCTGAGTGACCCTTCCCCAACGTGTCAGCTGAGTc                                                                                            <  6:280683/149‑1 (MQ=255)                                                      aGTGAAGAGGCTGCGGCCGCTATTTTCGGTAAGGCAGAGCTGGAGGGCGTTGGCTCATCGACGTGGCAGCAGATGTGGGAGTACGCTCGCGCATA‑CTCCGCCGCCGTGGCGTATCCTGAGTCACCCTTCCCCAACGTGTCAGCTGAGt                                                                                             >  3:269857/1‑148 (MQ=255)                                                                 tGCGGCCGCTATTTTCGGTAAGGCAGAGCTGGAGGGCGTTGGCTCATCGACGTGGCAGCAGATGTGGGAGTACGCTCGCGCATA‑CTCCGCCGCCGTGGCGTATCCTGAGTCACCCTTCCCCAACGTGTCAGCTGAGTCACGATGCGTGc                                                                                 <  3:55017/149‑1 (MQ=255)                                                                       cGCTATTTTCGGTAAGGCAGAGCTGGAGGGCGTTGGCTCATCGACGTGGCAGCAGATGTGGGAGTACGCTCGCGCATA‑CTCCGCCGCCGTGGCGTATCCTGAGTCACCCTTCCCCAACGTGTCAGCTGAGTCACGATGCGTGCTTTGc                                                                            >  4:302422/1‑148 (MQ=255)                                                                            ttttCGGTAAGGCAGAGCTGGAGGGCGTTGGCTCATCGACGTGGCAGCAGATGTGGGAGTACGCTCGCGCATA‑CTCCGCCGCCGTGGCGTATCCTGAGTCACCCTTCCCCAACGTGTCAGCTGAGTCACGATGCGTGCTTTGCCATCAg                                                                      <  1:165054/149‑1 (MQ=255)                                                                            ttttCGGTAAGGCAGAGCTGGAGGGCGTTGGCTCATCGACGTGGCAGCAGATGTGGGAGTACGCTCGCGCATA‑CTCCGCCGCCGTGGCGTATCCTGAGTCACCCTTCCCCAACGTGTCAGCTGAGTCACGATGCGTGCTTTGCCATCAg                                                                      >  3:21840/1‑149 (MQ=255)                                                                                 ggTAAGGCAGAGCTGGAGGGCGTTGGCTCATCGACGTGGCAGCAGATGTGGGAGTACGCTCGCGCATA‑CTCCGCCGCCGTGGCGTATCCTGAGTCACCCTTCCCCAACGTGTCAGCTGAGTCACGATGCGTGCTTTGCCATCAGGAtct                                                                 >  7:138395/1‑149 (MQ=255)                                                                                  gTAAGGCAGAGCTGGAGGGCGTTGGCTCATCGACGTGGCAGCAGATGTGGGAGTACGCTCGCGCATA‑CTCCGCCGCCGTGGCGTATCCTGCGTCACCCTTCCCCAACGTGTCAGCTGAGTCACGATGCGTGCTTTGCCATCAGGAtctc                                                                >  4:310985/1‑149 (MQ=255)                                                                                  gTAAGGCAGAGCTGGAGGGCGTTGGCTCATCGACGTGGCAGCAGATGTGGGAGTACGCTCGCGCATA‑CTCCGCCGCCGTGGCGTATCCTGAGTCACCCTTCCCCAACGTGTCAGCTGAGTCACGATGCGTGCTTTGCCATCAGGAtctc                                                                >  5:315190/1‑149 (MQ=255)                                                                                   tAAGGCAGAGCTGGAGGGCGTTGGCTCATCGACGTGGCAGCAGTGGTGGGAGTACGCTCGCGCATA‑CTCCGCCGCCATGGCGTATCCTGTGTCACCCTTCCCCAACGTGTCAGCTGAGTCACGATGCGTGCTTTGCCATCAGGATCTcg                                                               <  5:284512/149‑1 (MQ=255)                                                                                    aaGGCAGAGCTGGAGGGCGTTGGCTCATCGACGTGGCAGCAGATGTGGGAGTACGCTCGCGCATA‑CTCCGCCGCCGTGGCGTATCCTGAGTCACCCTTCCCCAACGTGTCAGCTGAGTCACGATGCGTGCTTTGCCATCAGGATCTcga                                                              >  6:116337/1‑149 (MQ=255)                                                                                       gCAGAGCTGGAGGGCGTTGGCTCATCGACGTGGCAGCAGATGTGGGAGTACGCTCGCGCATA‑CTCCGCCGCCGTGGCGTATCCTGAGTCACCCTTCCCCAACGTGTCAGCTGAGTCACGATGCGTGCTTTGCCATCAGGATCTcgacga                                                           >  5:94260/1‑149 (MQ=255)                                                                                        cAGAGCTGGGGGGCGTTGGCTCATCGACGTGGCAGCAGATGTGGGAGTACGCTCGCGCATA‑CTCCGCCGCCGTGGCGTATCCTGAGTCACCCTTCCCCAACGTGTCAGCTGAGTCACGATGCGTGCTTTGCCATCa                                                                       >  7:268655/1‑136 (MQ=255)                                                                                        cAGAGCTGGGGGGCGTTGGCTCATCGACGTGGCAGCAGATGTGGGAGTACGCTCGCGCATA‑CTCCGCCGCCGTGGCGTATCCTGAGTCACCCTTCCCCAACGTGTCAGCTGAGTCACGATGCGTGCTTTGCCATCa                                                                       <  8:268655/136‑1 (MQ=255)                                                                                             cTGGAGGGCGTTGGCTCATCGACGTGGCAGCAGATGTGGGAGTACGCTCGCGCATA‑CTCCGCCGCCGTGGCGTATCCTGAGTCACCCTTCCCCAACGTGTCAGCTGAGTCACGATGCGTGCTTTGCCATCAGGATCTCGACGATCCCGc                                                     <  7:115477/149‑1 (MQ=255)                                                                                                gAGGGCGTTGGCTCATCGACGTGGCAGCAGATGTGGGAGTACGCTCGCGCATA‑CTCCGCCGCCGTGGCGTATCCTGAGTCACCCTTCCCCAACGTGTCAGCTGAGTCACGATGCGTGCttt                                                                              >  4:225040/1‑121 (MQ=255)                                                                                                gAGGGCGTTGGCTCATCGACGTGGCAGCAGATGTGGGAGTACGCTCGCGCATA‑CTCCGCCGCCGTGGCGTATCCTGAGTCACCCTTCCCCAACGTGTCAGCTGAGTCACGATGCGTGCTTTGCCATCAGGAt                                                                   >  5:272064/1‑132 (MQ=255)                                                                                                gAGGGCGTTGGCTCATCGACGTGGCAGCAGATGTGGGAGTACGCTCGCGCATA‑CTCCGCCGCCGTGGCGTATCCTGAGTCACCCTTCCCCAACGTGTCAGCTGAGTCACGATGCGTGCTTTGCCATCAGGAt                                                                   <  6:272064/132‑1 (MQ=255)                                                                                                 aGGGCTTTGGCTCATCGACGTGGCAGCAGATGTGGGAGTACGCTCGCGCATA‑CTCCGCCGCCGTGGCGTATCCTGAGTCACCCTTCCCCAACGTGTCAGCTGAGTCACGATGCGTGCTTTGCCATCAGGATCTCGACGATCCCGCaaaa                                                 <  7:202836/149‑1 (MQ=255)                                                                                                      gTTGGCTCATCGACGTGGCAGCAGATGTGGGAGTACGCTCGCGCATA‑CTCCGCCGCCGTCGCGTATCCTGAGTCACCCTTCCCCAACGTGTCAGCTGAGTCACGATGCGTGCTTTGCCATCAGGATCTCGACGATCCCGCAAAAGCACg                                            >  6:11738/1‑149 (MQ=255)                                                                                                              aTCGACGTGGCAGCAGATGTGGGAGTACGCTCGCGCATA‑CTCCGCCGCCGTGGCGTATCCTGAGTCACCCTTCCCCAACGTGTCAGCTGAGTCACGATGCGTGCTTTGCCATCAGGATCTCGACGATCCCGCAAAAGCACGCCTGGcc                                     <  2:129072/148‑1 (MQ=255)                                                                                                              aTCGACGTGGCAGCAGATGTGGGAGTACGCTCGCGCATA‑CTCCGCCGCCGTGGCGTATCCTGAGTCACCCTTCCCCAACGTGTCAGCTGAGTCACGATGCGTGCTTTGCCATCAGGATCTCGACGATCCCGCAAAAGCACGCCTGGcc                                     >  6:217665/1‑148 (MQ=255)                                                                                                                 gACGTGGCAGCAGATGTGGGAGTACGCTCGCGCATA‑CTCCGCCGCCGTGGCGTATCCTGAGTCACCCTTCCCCAACGTGTCAGCTGAGTCACGATGCGTGCTTTGCCATCAGGATCTCGACGATCCCGCAAAAGCACGCCTGGCCGaa                                  <  5:217665/148‑1 (MQ=255)                                                                                                                 gACGTGGCAGCAGATGTGGGAGTACGCTCGCGCATA‑CTCCGCCGCCGTGGCGTATCCTGAGTCACCCTTCCCCAACGTGTCAGCTGAGTCACGATGCGTGCTTTGCCATCAGGATCTCGACGATCCCGCAAAAGCACGCCTGGCCGAAt                                 >  5:143122/1‑149 (MQ=255)                                                                                                                   cGTGGCAGCAGATGTGGGAGTACGCTCGCGCATA‑CTCCGCCGCCGTGGCGTATCCTGAGTCACCCTTCCCCAACGTGTCAGCTGAGTCACGATGCGTGCTTTGCCATCAGGATCTCGACGATCCCGCAAAAGCACGCCTGGCCGAAttt                               <  7:201087/149‑1 (MQ=255)                                                                                                                      ggCAGCAGATGTGGGAGTACGCTCGCGCATA‑CTCCGCCGCCGTGGCGTATCCTGAGTCACCCTTCCCCAACGTGTCAGCTGAGTCACGATGCGTGCTTTg                                                                             >  8:276373/1‑100 (MQ=255)                                                                                                                         agcagATGTGGGAGTACGCTCTCGCATA‑CTCCGCCGCCGTGGCGTATCCTGAGTCACCCTTCCCCAACGTGTCAGCTGAGTCACGATGCGTGCTTTGCCATCAGGATCTCGACGATCCCGCAAAAGCACGCCTGGCCGAATTTGAAAAg                         >  6:127703/1‑149 (MQ=255)                                                                                                                         agcagATGTGGGAGTACGCTCGCGCATT‑CTCCGCCGCCGTGGCGTATCCTGAGTCACCCTTCCCCAACGTGTCAGCTGAGTCACGATGCGTGCTTTGCCATCAGGATCTCGACGATCCCGCAAAAGCACGCCTGGCCGAATTTGAAAAg                         >  8:228339/1‑149 (MQ=255)                                                                                                                            agATGTGGGAGTACGCTCGCGCATA‑CTCCGCCGCCGTGGCGTATCCTGATTCACCCTTCCCCAACGTGTCAGCTGAGTCACGATGCGTGCTTTGCCATCAGGATCTCGACGATCCAGCAAAAGCACGCCTGGCCGAATTTGAAAAGttt                      <  5:127214/149‑1 (MQ=255)                                                                                                                              aTGTGGGAGTACGCTCGCGCATA‑CTCCGCCGCCGTGGCGTATCCTGAGTCACCCTTCCCCACCGTGTCAGCTGAGTCACGATGCGTGCTTTGCCATCAGGATCTCGACGATCCCGCAAAAGCACGCCTGGCCGAATTTGAAAAGTTTGt                    <  4:117394/149‑1 (MQ=255)                                                                                                                              aTGTGGGAGTACGCTCGCGCATA‑CTCCGCCGCCGTGGCGTATCCTGAGTCACCCTTCCCCAACGTTTCAGCTGAGTCACGATGCGTGCTTTTCCATCAGGATCTCGACGATCCCGCAAAAGCACGCCTGGCCGAATTTGAAAAGTTTg                     >  7:231076/1‑148 (MQ=255)                                                                                                                                gtgGGAGTACGCTCGCGCATA‑CTCCGCCGCCGTGGCGTATCCTGAGTCACCCTTCCCCAACGTGTCAGCTGAGTCACGATGCGTGCTTTGCCATCAGGATCTCGACGATCCCGCAAAAGCACGCCTGGCCGAATTTGAAAAGTTTGTc                   >  7:180950/1‑148 (MQ=255)                                                                                                                                  gggAGTACGCTCGCGCATA‑CTCCGCCGCCGTGGCGTATCCTGAGTCACCCTTCCCCAACGTGTCAGCTGAGTCACGATGCGTGCTTTGCCATCAGGATCTCGACGATCCCGc                                                     >  2:221138/1‑112 (MQ=255)                                                                                                                                  gggAGTACGCTCGCGCATA‑CTCCGCCGCCGTGGCGTATCCTGAGTCACCCTTCCCCAACGTGTCAGCTGAGTCACGATGCGTGCTTTGCCATCAGGATCTCGACGATCCCGc                                                     <  1:221138/112‑1 (MQ=255)                                                                                                                                  cggAGTATGCGCGAGCGTA‑CTCCGCCGCCGTGGCGTATCCTGAGTCACCCTTCCCCAACGTGTCAGCTGAGTCACGATGCGt                                                                                   <  6:28322/81‑1 (MQ=38)                                                                                                                                  cggAGTATGCGCGAGCGTA‑CTCCGCCGCCGTGGCGTATCCTGAGTCACCCTTCCCCAACGTGTCAGCTGAGTCACGATGCGt                                                                                   >  5:28322/2‑82 (MQ=38)                                                                                                                                    gAGTACGCTCGCGCATA‑CTCCGCCGCCGTGGCGTATCCTGAGTCACCCTTCCCCAACGTGTCAGCTGAGTCACGATGCGTGCTTTGCCATCAGGATCTCGACGATCCCGCAAAAGCACGCCTGGCCGAATTTGAAAAGTTTGTCActgc              <  4:269857/149‑1 (MQ=255)                                                                                                                                      gTACGCTCGCGCATA‑CTCCGCCGCCGTGGCGTATCCTGAGTCACCCTTCCCCAACGTGTCAGCTGAGTCACGATGCGTGCTTTGCCATCAGGATCTCGACGCTCCCGCAAAAGCACGCCTGGCCGAATTTGAAAAGTTTGTCActgctg            >  6:282167/1‑149 (MQ=255)                                                                                                                                      gTACGCTCGCGCATA‑CTCCGCCGCCGTGGCGTATCCTGAGTCACCCTTCCCCAACGTGTCAGCTGAGTCACGATGCGTGCTTTGCCATCAGGATCTCGACGATCCCGCAAAAGCACGCCTGGCCGAATTTGAAAAGTTTGTCActgctg            >  4:118847/1‑149 (MQ=255)                                                                                                                                          gCTCGCGCATA‑CTCCGCCGCCGTGGCGTATCCTGAGTCACCCTTCCCCAACGTGTCAGCTGAGTCACGATGCGTGCTTTGCCATCAGGATCTCGACGATCCCGCaa                                                   <  6:299089/106‑1 (MQ=255)                                                                                                                                          gCTCGCGCATA‑CTCCGCCGCCGTGGCGTATCCTGAGTCACCCTTCCCCAACGTGTCAGCTGAGTCACGATGCGTGCTTTGCCATCAGGATCTCGACGATCCCGCaa                                                   >  5:299089/1‑106 (MQ=255)                                                                                                                                           cTCGCGCATA‑CTCCGCCGCCGTGGCGTATCCTGAGTCACCCTTCCCCAACGTGTCAGCTGAGTCACGATGCGTGCTTTGCCATCAGGATCTCGACGATCCCGCAAAAGCACGCCTGGCCGAATTTGAAAAGTTTGTCACTGCTGGTTTg       <  6:143122/149‑1 (MQ=255)                                                                                                                                            tCGCGCATA‑CTCCGCCGCCGTGGCGTATCCTGAGTCACCCTTCCCCAACGTGTCAGCTGAGTCACGATGCGTGCTTTGCCATCAGGATCTCGACGATCCCGCAAAAGCACGCCTGGCCGAATTTGAAAAGTTTGTCACTGCTGGTTTg       >  8:213587/1‑148 (MQ=255)                                                                                                                                            tCGCGCATA‑CTCCGCCGCCGTGGCGTATCCTGAGTCACCCTTCCCCAACGTGTCAGCTGAGTCACGATGCGTGCTTTGCCATCAGGATCTCGACGATCCCGCAAAAGCACGCCTGGCCGAATTTGAAAAGTTTGTCACTGCTGGTTTg       >  3:249140/1‑148 (MQ=255)                                                                                                                                              gcgcATA‑CTCCGCCGCCGTGGCGTATCCTGAGTCACCCTTCCCCAACGTGTCAGCTGAGTCACGATGCGTGCTTTGCCATCAGGATCTCGACGATCCCGCAAAAGCACGCCTGGCCGAATTTGAAAAGTTTGTCACTGCTGGTTTg       >  5:34727/1‑146 (MQ=255)                                                                                                                                              gcgcATA‑CTCCGCCGCCGTGGCGTATCCTGAGTCACCCTTCCCCAACGTGTCAGCTGAGTCACGATGCGTGCTTTGCCATCAGGATCTCGACGATCCCGCAAAAGCACGCCTGGCCGAATTTGAAAAGTTTGTCACTGCTGGTTTg       >  5:137473/1‑146 (MQ=255)                                                                                                                                              gcgcATA‑CTCCGCCGCCGTGGCGTATCCTGAGTCACCCTTCCCCAACGTGTCAGCTGAGTCACGATGCGTGCTTTGCCATCAGGATCTCGACGATCCCGCAAAAGCACGCCTGGCCGAATTTGAAAAGTTTGTCACTGCTGGTTTGga     <  5:11738/148‑1 (MQ=255)                                                                                                                                                gcATA‑CTCCGCCGCCGTGGCGTATCCTGAGTCACCCTTCCCCAACGTGTCAGCTGAGTCACGATGCGTGCTTTGCCATCAGGATCTCGACGATCCCGCAAAAGCACGCCTGGCCGAATTTGAAAAGTTTGTCACTGCTGGTTTGgaga   >  7:212554/1‑148 (MQ=255)                                                                                                                                                gcATA‑CTCCGCCGCCGTGGCGTATCCTGAGTCACCCTTCCCCAACGTGTCAGCTGAGTCACGATGCGTGCTTTGCCATCAGGATCTCGACGATCCCGCAAAAGCACGCCTGGCCGAATTTGAAAAGTTTGTCACTGCTGGTTTGgaga   >  5:225538/1‑148 (MQ=255)                                                                                                                                                gcATA‑CTCCGCCGCCGTGGCGTATCCTGAGTCACCCTTCCCCAACGTGTCAGCTGAGTCACGATGCGTGCTTTGCCATCAGGATCTCGACGATCCCGCAAAAGCACGCCTGGCCGAATTTGAAAAGTTTGTCACTGCTGGTTTGGAGAc  <  5:127703/149‑1 (MQ=255)                                                                                                                                                 cATA‑CTCCGCCGCCGTGGCGTATCCTGAGTCACCCTTCCCCAACGTGTCAGCTGAGTCACGATGCGTGCTTTGCCATCAGGATCTCGACGATCCCGCAAAAGCACGCCTGGCCGAATTTGAAAAGTTTGTCACTGCTGGTTTGgaga   >  6:359466/1‑147 (MQ=255)                                                                                                                                                  |                                                                                                                                                    CCGCACTTGAGCTAGTCAACTTGAGTGAAAAAGCCGCGAGAACACGCAAAGCCAGTGAAGAGGCTGCGGCCGCTATTTTCGGTAAGGCAGAGCTGGAGGGCGTTGGCTCATCGACGTGGCAGCAGATGTGGGAGTACGCTCGCGCATA‑CTCCGCCGCCGTGGCGTATCCTGAGTCACCCTTCCCCAACGTGTCAGCTGAGTCACGATGCGTGCTTTGCCATCAGGATCTCGACGATCCCGCAAAAGCACGCCTGGCCGAATTTGAAAAGTTTGTCACTGCTGGTTTGGAGAC  >  NC_002947/4209078‑4209369

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