Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,948,143	A→C	76.4%	F111V (TTC→GTC)	rho ←	transcription termination factor Rho

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,948,143	0	A	C	76.4%	130.7 / 37.2	67	F111V (TTC→GTC)	rho	transcription termination factor Rho
Reads supporting (aligned to +/- strand):  ref base A (7/9);  new base C (25/26);  total (32/35)
Fisher's exact test for biased strand distribution p-value = 7.80e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.84e-01

CCGGTCAAGTCCTCGGTGGAGCCGTTACCGGCCTCCATCTTCAGGCGCTTGTTCGGGAAC‑AGCGGCGTCAGGTTTTCGAACAGGATCTTGTTCTTCGCGTTTTCCGGGCGGTCGAAGTTGATGGTATCGACCTTCAGCAGGGCGAAGTAACGCTCCCCTTCCTTCGGTGGGCGGATCTTGCCGACGATGGTGTCGCCGGTACGCAGGTTGAAACGGCGGATCTGGCTGGGCGAGACATAGATATCGTCTGGGCCGGCCAGGTAGGACGCATCAGCCGAAC  >  NC_002947/5947998‑5948277                                                                                                                                                   |                                                                                                                                       ccGGTCAAGTCCTCGGTGGAGCCGTTACCGGCCTCCATCTTCAGGCGCTTGTTCGGGAAC‑AGCGGCGTCAGGTTTTCGAACAGGATCTTGTTCTTCGCGTTTTCCGGGCGGTCGAAGTTGCTGGTATCGACCTTCAGCAGGGCGACGt                                                                                                                                      >  6:209400/1‑148 (MQ=255)            cTCGGTGGAGCCGTTACCGGCCTCCATCTTCAGGCGCTTGTTCGGGAAC‑AGCGGCGTCAGGTTTTCGAACAGGATCTTGTTCTTCGCGTTTTCCGGGCGGTCGAAGTTGATGGTATCGACCTTCAGCAGGGCGAAGTAACGCTCCcctt                                                                                                                          >  6:324605/1‑149 (MQ=255)             tCGGTGGAGCCGTTACCGGCCTCCATCTTCAGGCGCTTGTTCGGGAAC‑AGCGGCGTCAGGTTTTCGAACAGGATCTTGTTCTTCGCGTTTTCCGGGCGGTCGAAGTTGATGGTATCGACCTTCAGCAGGGCGACGTAACGCTCCccttc                                                                                                                         >  1:342752/1‑149 (MQ=255)                       cGTTACCGGCCTCCATCTTCAGGCGCTTGTTCGGGAAC‑AGCGGCGTCAGGTTTTCGAACAGGATCTTGTTCTTCGCGTTTTCCGGGCGGTCGAAGTTGATGGTATCGACCTTCAGCAGGGCGACGTAACGCTCCCCTTCCTTCGGTggg                                                                                                               <  5:53319/149‑1 (MQ=255)                             cGGCCTCCATCTTCAGGCGCTTGTTCGGGAAC‑AGCGGCGTCAGGTTTTCGAACAGGATCTTGTTCTTCGCGTTTTCCGGGCGGTCGAAGTTGATGGTATCGACCTTCAGCAGGGCGAAGt                                                                                                                                      <  5:290727/120‑1 (MQ=255)                             cGGCCTCCATCTTCAGGCGCTTGTTCGGGAAC‑AGCGGCGTCAGGTTTTCGAACAGGATCTTGTTCTTCGCGTTTTCCGGGCGGTCGAAGTTGATGGTATCGACCTTCAGCAGGGCGAAGt                                                                                                                                      >  6:290727/1‑120 (MQ=255)                               gCCTCCATCTTCAGGCGCTTGTTCGGGAAC‑AGCGGCGTCAGGTTTTCGAACAGGATCTTGTTCTTCGCGTTTTCCGGGCGGTCGAAGTTGATGGTATCGACCTTCAGCAGGGCGACGTAACGCTCCCCTTCCTTCGGTGGGCGGATCtt                                                                                                       >  6:307616/1‑149 (MQ=255)                                ccTCTATCTTCAGGCGCTTGTTCGGGAAC‑AGCGGCGTCAGGTTTTCGTACAGGATCTTGTTCTTCGCGTTTTCCGGGCGGTCGAAGTTGATGGTATCGACCTTCAGCAGGGCGACGTAACTCTCCCCTTCCTTCGGTGGGCGGATCTTg                                                                                                      <  8:318888/149‑1 (MQ=255)                                ccTCCATCTTCAGGCGCTTGTTCGGGAAC‑AGCGGCGTCAGGTTTTCGAACAGGATCTTGTTCTTCGCGTTTTCCGGGCGGTCGAAGTTGATGGTATCGACCTTCAGCAGGGCGACGTAACGCTCCCCTTCCTTCGGTGGGCGGATCTTg                                                                                                      >  5:172537/1‑149 (MQ=255)                                 cTCCATCTTCAGGCGCTTGTTCGGGAAC‑AGCGGCGTCAGGTTTTCGAACAGGATCTTGTTCTTCGCGTTTTCCGGGCGGTCGAAGTTGATGGTATCGACCTTCAGCAGGGCGAAGTAACGCTCCCCTTCCTTCGGTGGGCGGATCTTGc                                                                                                     <  4:96065/149‑1 (MQ=255)                                  tCCATCTTCAGGCGCTTGTTCGGGAAC‑AGCGGCGTCAGGTTTTCGAACAGGATCTTGTTCTTCGCGTTTTCCGGGCGGTCGAAGTTGATGGTATCGACCTTCAGCAGGGCGACGTAACGCTCCCCTTCCTTCGGTGGGCGGATCTTGcc                                                                                                    >  7:81566/1‑149 (MQ=255)                                    cATCTTCAGGCGCTTGTTCGGGAAC‑AGCGGCGTCAGGTTTTCGAACAGGATCTTGTTCTTCGCGTTTTCCGGGCGGTCGAAGTTGATGGTATCGACCTTCAGCAGGGCGACGTAACGc                                                                                                                                 <  7:182888/118‑1 (MQ=255)                                    cATCTTCAGGCGCTTGTTCGGGAAC‑AGCGGCGTCAGGTTTTCGAACAGGATCTTGTTCTTCGCGTTTTCCGGGCGGTCGAAGTTGATGGTATCGACCTTCAGCAGGGCGACGTAACGc                                                                                                                                 >  8:182888/1‑118 (MQ=255)                                        ttCAGGCGCTTGTTCGGGAAC‑AGCGGCGTCAGGTTTTCGAGCAGGATCTTGTTCTTCGCGTTTTCCGGGCGGTCGAAGTTGATGGTATCGACCTTCAGCAGGGCGACGTAACGCTCCCCTTCCTTCGGTGGGCGGATCTTGCCGACGAt                                                                                              >  1:254650/1‑149 (MQ=255)                                          cAGGCGCTTGTTCGGGAAC‑AGCGGCGTCAGGTTTTCGAGCAGGATCTTGTTCTTCGCGTTTTCCGGGCGGTCGAAGTTGATGGTATCGACCTTCAGCAGGGCGACGTAACGCTCCCCTTCCTTCGGTGGGCGGATCTTGCCGACGATg                                                                                             <  2:254650/148‑1 (MQ=255)                                          cAGGCGCTTGTTCGGGAAC‑AGCGGCGTCAGGTTTTCGAACAGGATCTTGTTCTTCGCGTTTTCCGGGCGGTCGAAGTTGATGGTATCGACCTTCAGCAGGGCGAAGTAACGCTCCCCTTCCTTCGGTGGGCGGATCTTGCCGACGATg                                                                                             <  5:324605/148‑1 (MQ=255)                                          cAGGCGCTTGTTCGGGAAC‑AGCGGCGTCAGGTTTTCGAACAGGATCTTGTTCTTCGCGTTTTCCGGGCGGTCGAAGTTGATGGTATCGACCTTCAGCAGGACGAAGTAACGCTCCCCTTCCTTCGGTGGGCGGATCTTGCCGACGATg                                                                                             <  6:301664/148‑1 (MQ=255)                                           acgccgctGTTCCCGAACAAGCGGCGTCAGGTTTTCGAACAGGATCTTGTTCTTCGCGTTTTCCGGGCGGTCGAAGTTGATGGTATCGACCTTCAGCAGGGCGACGTAACGCTCCCCTTCCTTCGGTGGGCGGATCTTGCCGACGATgg                                                                                            >  1:30490/8‑149 (MQ=255)                                              cgcTTGTTCGGGAAC‑AGCGGCGTCAGGTTTTCGAACAGGATCTTGTTCTTCGCGTTTTCCGGGCGGTCGAAGTTGATGGTATCGACCTTCAGCAGGGCGACGTAACGCTCCCCTTCCTTCGGTGGGCGGATCTTGCCGACGATGGTGTc                                                                                        >  6:279999/1‑149 (MQ=255)                                               gcTTGTTCGGGAAC‑AGCGGCGTCAGGTTTTCGAACAGGATCTTGTTCTTCGCGTTTTCCGGGCGGTCGAAGTTGATGGTATCGACCTTCAGCAGGGCGACGTAACGCTCCCCTTCCTTCGGTGGGCGGATCTTGCCGACGATGGTGTCg                                                                                       >  6:82662/1‑149 (MQ=255)                                               gcTTGTTCGGGAAC‑AGCGGCGTCAGGTTTTCGAACAGGATCTTGTTCTTCGCGTTTTCCGGGCGGTCGAAGTTGATGGTATCGACCTTCAGCAGGGCGACGTAACGCTCCCCTTCCTTCGGTGGGCGGATCTTGCCGACGATGGTGTCg                                                                                       >  1:320295/1‑149 (MQ=255)                                                cTTGTTCGGGAAC‑AGCGGCGTCAGGTTTTCGAACAGGATCTTGTTCTTCGCGTTTTCCGGGCGGTCGAAGTTGATGGTATCGACCTTCAGCAGGGCGACGTAACGCTCCCCTTCCTTCGGTGGGCGGATCTTGCCGACGATGGTGTCGc                                                                                      <  2:320295/149‑1 (MQ=255)                                                   gTTCGGGAAC‑AGCGGCGTCAGGTTTTCGAACAGGATCTTGTTCTTCGCGTTTTCCGGGCGGTCGAAGTTGATGGTATCGACCTTCAGCAGGGCGACGTAACGCTCCCCTTCCTTCGGTGGGCGGATCTTGcc                                                                                                    >  3:168140/1‑132 (MQ=255)                                                   gTTCGGGAAC‑AGCGGCGTCAGGTTTTCGAACAGGATCTTGTTCTTCGCGTTTTCCGGGCGGTCGAAGTTGATGGTATCGACCTTCAGCAGGGCGACGTAACGCTCCCCTTCCTTCGGTGGGCGGATCTTGcc                                                                                                    <  4:168140/132‑1 (MQ=255)                                                    ttCGGGAAC‑AGCGGCGTCAGGTTTTCGAACAGGATCTTGTTCTTCGCGTTTTCCGGGCGGTCGAAGTTGATGGTATCGACCTTCAGCAGGGCGACGTAACGCTCCCCTTCCTTCGGTGGGCGGATCTTGCCGACGATGGTGTCGCCgg                                                                                   <  2:342752/148‑1 (MQ=255)                                                    ttCGGGAAC‑AGCGGCGTCAGGTTTTCGAACAGGATCTTGTTCTTCGCGTTTTCCGGGCGGTCGAAGTTGATGGTATCGACCTTCAGCAGGGCGACGTAACGCTCCCCTTCCTTCGGTGGGCGGATCTTGCCGACGATGGTGTCGCCgg                                                                                   <  5:279999/148‑1 (MQ=255)                                                    ttCGGGAAC‑AGCGGCGTCAGGTTTTCGAACAGGATCTTGTGCTTCGCGTTTTCCGGGCGGTCGAAGTTGATGGTATCGACCTTCAGCAGGGCGACGTAACGCTCCCCTTCCTTCGGTGGGCGGATCTTGCCGACGATGGTGTCGCCgg                                                                                   <  4:271338/148‑1 (MQ=255)                                                      cGGGAAC‑AGCGGCGTCAGGTTTTCGAACAGGATCTTGTTCTTCGCGTTTTCCGGGCGGTCGAAGTTGATGGTATCGACCTTCAGCAGGGCGACGTAACGCTCCCCTTCCTTCGGTGGGCGGATCTTGCCGACGATGGTGTCGCCGGTAc                                                                                >  1:235190/1‑149 (MQ=255)                                                       gggAAC‑AGCGGCGTCAGGTTTTCGAACAGGATCTTGTTCTTCGCGTTTTCCGGGCGGTCGAAGTTGATGGTATCGACCTTCAGCAGGGCGACGTAACGCTCCCCTTCCTTCGGTGGGCGGATCTTGCCGACGATGGTGTCGCCGGTAc                                                                                >  1:283734/1‑148 (MQ=255)                                                            caaGCGGCGTCAGGTTTTCGAACAGGATCTTGTTCTTCGCGTTTTCCGGGCGGTCGAAGTTGATGGTATCGACCTTCAGCAGGGCGACGTAACGCTCCCCTTCCTTCGGTGGGCGGATCTTGCCGACGATGGTGTCGCCGGTACGCAgg                                                                           <  2:30490/147‑1 (MQ=255)                                                              aGCGGCGTCAGGTTTTCGAACAGGATCTTGTTCTTCGCGTTTTCCGGGCGGTCGAAGTTGATGGTATCGACCTTCAGCAGGGCGACGTAACGCTCCCCTTCCTTCGGTGGGCGGATCTTGCCGACGATGGTGTCGCCGGTACGCAGGtt                                                                         <  5:82662/149‑1 (MQ=255)                                                                cggcgTCAGGTTTTCGAACAGGATCTTGTTCTTCGCGTTTTCCGGGCGGTCGAAGTTGATGGTATCGACCTTCAGCAGGGCGAAGTAACGCTCCCCTTCCTTCGGTGGGCGGATCTTGCCGACGATGGTGTCGCCGGTACGCAGGTTg                                                                        >  1:67261/1‑148 (MQ=255)                                                                   cgTCAGGTTTTCGAACAGGATCTTGTTCTTCGCGTTTTCCGGGCGGTCGAAGTTGATGGTATCGACCTTCAGCAGGGCGAAGTAACGCTCCCCTTCCTTCGGTGGGCGGATCTTGCCGACGATGGTGTCGCCGGTACGCAGGt                                                                          >  8:202030/1‑143 (MQ=255)                                                                   cgTCAGGTTTTCGAACAGGATCTTGTTCTTCGCGTTTTCCGGGCGGTCGAAGTTGATGGTATCGACCTTCAGCAGGGCGAAGTAACGCTCCCCTTCCTTCGGTGGGCGGATCTTGCCGACGATGGTGTCGCCGGTACGCAGGt                                                                          <  7:202030/143‑1 (MQ=255)                                                                       aGGTTTTCGAACAGGATCTTGTTCTTCGCGTTTTCCGGGCGGTCGAAGTTGATGGTATCGACCTTCAGCAGGGCGACGTAACGCTCCCCTTCCTTCGGTGGGCGGATCTTGCCGACGATGGTGTCGCCGGTACGCAGGTTGAAAcggcg                                                                <  8:81566/149‑1 (MQ=255)                                                                         gTTTTAGAACAGGATATTGTTCTTCGCTTTTTCCGTGCGGTATAAGTTGATGGTATCGAACTTAAGCAGGGCGAAGTAACGCTCCCCTTCCTTAGGTGGGAGGATCTTGCCGACTATGTTGTAGCCGGTACTCAGCGTGaaaagtagga                                                              >  8:330747/1‑142 (MQ=255)                                                                                  cAGGATCTTGTTCTTCGCGTTTTCCGGGCGGTCGAAGTTGAAGGTATCGACCTTCAGCAGGGCGACGTAACGCTCCCCTTCCTTCGGTGGGCGGATCTTGCCGACGATGGTGTCGCCGGTACGCAGGTTGAAACGGCGGATctggctgg                                                     <  5:307616/149‑1 (MQ=255)                                                                                    ggATCTTGTTCTTCGCGTTTTCCGGGCGGTCGAAGTTGATGGTATCGACCTTCAGCAGGGCGAAGTAACGCTCCCCTTCCTTCGGTGGGCGGATCTTGCCGACGATGGTGTCGCCGGTACGCAGGTTGAAACGGCGGATCTGGCTGGGc                                                   <  2:67261/149‑1 (MQ=255)                                                                                           gTTCTTCGCGTTTTCCGGGCGGTCGAAGTTGATGGTATCGACCTTCAGCAGGGCGACGTAACGCTCCCCTTCCTTCGGTGGGCGGATCTTGCCGACGATGGTGTCGCCGGTACGCAGGTTGAAACGGCGGATCTGGCTGGGCgag                                                <  2:235190/145‑1 (MQ=255)                                                                                           gTTCTTCGCGTTTTCCGGGCGGTCGAAGTTGATGGTATCGACCTTCAGCAGGGCGACGTAACGCTCCCCTTCCTTCGGTGGGCGGATCTTGCCGACGATGGTGTCGCCGGTACGCAGGTTGAAACGGCGGATCTGGCTGGGCga                                                 <  5:209400/144‑1 (MQ=255)                                                                                           gTTCTTCGCGTTTTCCGGGCGGTCGAAGTTGATGGTATCGACCTTCAGCAGGGCGACGTAACGCTCCCCTTCCTTCGGTGGGCGGATCTTGCCGACGATGGTGTCGCCGGTACGCAGGTTGAAACGGCGGATCTGGCTGGGCGAGACa                                             >  5:92665/1‑148 (MQ=255)                                                                                               ttcGCGTTTTCCGGGCGGTCGAAGTTGATGGTATCGACCTTCAGCAGGGCGAAGTAACGCTCCCCTTCCTTCGGTGGGCGGATCTTGCCGACGATGGTGTCGCCGGTACGCAGGTTGAAACGGCGGATCTGGCTGGGCGAGACATAGa                                         >  7:3111/1‑148 (MQ=255)                                                                                               ttcGCGTTGTCCGGGCGGTCGATTTTGATGGTATCGACCTTCAGCAGGGCGACGTAACGCTCCCCTTCCTTCGGTGGGGGGTTCTTGCcgacga                                                                                               <  6:266903/94‑1 (MQ=255)                                                                                                  gcgTTTTCCGGGCGGTCGAAGTTGATGGTATCGACCTTCAGCAGGGCGACGTAACGCTCCCCTTCCTTCGGTGGGCGGATCTTGCCGACGATGGTGTCGCCGGTACGCAGGTTGAAACGGCGGATCTGGCTGGGCGAGACATa                                           >  5:44788/1‑143 (MQ=255)                                                                                                  gcgTTTTCCGGGCGGTCGAAGTTGATGGTATCGACCTTCAGCAGGGCGACGTAACGCTCCCCTTCCTTCGGTGGGCGGATCTTGCCGACGATGGTGTCGCCGGTACGCAGGTTGAAACGGCGGATCTGGCTGGGCGAGACATa                                           <  6:44788/143‑1 (MQ=255)                                                                                                   cgTTTTCCGGGCGGTCGAAGTTGATGGTATCGACCTTCAGCAGGGCGACGTAACGCTCCCCTTCCTTCGGTGGGCGGATCTTGCCGACGATGGTGTCGCCGGTACGCAGGTTGAAACGGCGGATCTGGCTGGGCGAGACa                                             >  1:121964/1‑140 (MQ=255)                                                                                                   cgTTTTCCGGGCGGTCGAAGTTGATGGTATCGACCTTCAGCAGGGCGACGTAACGCTCCCCTTCCTTCGGTGGGCGGATCTTGCCGACGATGGTGTCGCCGGTACGCAGGTTGAAACGGCGGATCTGGCTGGGCGAGACa                                             <  2:121964/140‑1 (MQ=255)                                                                                                     ttttCCGGGCGGTCGAAGTTGATGGTATCGACCTTCAGCAGGGCGAAGTAACGCTCCCCTTCCTTCGGTGGGCGGATCTTGCCGACGATGGTGTCGCCGGTACGCAGGTTGAAACGGCGGATCTGGCTGGGCGAGACATAGATATCGTc                                  <  8:3111/149‑1 (MQ=255)                                                                                                       tgccGGGCGGGCGATGTTGATGGTATCGACCTTCAGCAGGGCGACGTAACGCTCCCCTTCCTTCGGTGGGCGGATCTTGCCGACGATGGTGTCGCCGGTACGCAGGTTGAAACGGCGGATCTGGCTGGGCGAGACATAGATATCGTCt                                 <  6:92665/146‑1 (MQ=255)                                                                                                         ccGGGCGGTCGAAGTTGATGGTATCGACCTTCAGCAGGGCGACGTAACGCTCCCCTTCCTTCGGTTGGCGGATCTTGCCGACGATGGTGTCGCCGGTACGCAGGTTGAAACGGCGGATCTGGCTGGGCGAGACATAGATATCGTCTgg                               <  2:283734/148‑1 (MQ=255)                                                                                                         ccGGGCGGTCGAAGTTGATGGTATCGACCTTCAGCAGGGCGACGTAACGCTCCCCTTCCTTCGGTGGGCGGATCTTGCCGACGATGGTGTCGCCGGTACGCAGGTTGAAACGGCGGATCTGGCTGGGCGAGACATAGATATCGTCTgg                               <  7:330747/148‑1 (MQ=255)                                                                                                                   gAAGTTGATGGTATCGACCTTCAGCAGGGCGACGTAACGCTCCCCTTCCTTCGGTGGGCGGATCTTGCCGACGATGGTGTCGCCGGTACGCAGGTTGAAACGGCGGATCTGGCTGGGCgag                                                >  7:244522/1‑121 (MQ=255)                                                                                                                   gAAGTTGATGGTATCGACCTTCAGCAGGGCGACGTAACGCTCCCCTTCCTTCGGTGGGCGGATCTTGCCGACGATGGTGTCGCCGGTACGCAGGTTGAAACGGCGGATCTGGCTGGGCgag                                                <  8:244522/121‑1 (MQ=255)                                                                                                                   gAAGTTGATGGTATCGACCTTCAGCAGGGCGAAGTAACGCTCCCCTTCCTTCGGTGGGCGGATCTTGCCGACGATCGTGTCGCCGGTACGCAGGTTGAAACGGCGGATCTGGCTGGGCGAGACATAGATATCGTCTGGGCCGGCCAGGt                    <  1:318729/149‑1 (MQ=255)                                                                                                                         gATGGTATCGACCTTCAGCAGGGCGACGTAACGCTCCCCTTCCTTCGGTGGGCGGATCTTGCCGACGATGGTGTCGCCGGTACGCAGGTTGAAACGGCGGATCTGGCTGGGCGAGACATAGATATCGTCTGGGCCGGCCAg                      <  2:318131/141‑1 (MQ=255)                                                                                                                         gATGGTATCGACCTTCAGCAGGGCGACGTAACGCTCCCCTTCCTTCGGTGGGCGGATCTTGCCGACGATGGTGTCGCCGGTACGCAGGTTGAAACGGCGGATCTGGCTGGGCGAGACATAGATATCGTCTGGGCCGGCCAg                      >  1:318131/1‑141 (MQ=255)                                                                                                                          aTGGTATCGACCTTCAGCAGGGCGACGTAACGCTCCCCTTCCTTCGGTGGGCGGATCTTGCCGACGATGGTGTCGCCGGTACGCAGGTTGAAACGGCGGATCTGGCTGGGCGAGACATAGATATCGTCTGGGCCGGCCAGGTAGGACg              >  6:151693/1‑148 (MQ=255)                                                                                                                           tGGTATCGACCTTCAGCAGGGCGACGTAACGCTCCCCTTCCTTCGGTGGGCGGATCTTGCCGACGATGGTGTCGCCGGTACGCAGGTTGAAACGGTGGATCTGGCTGGGCGAGACATAGATATCGTCTGGGCCGGCCAGGTAGGACGCa            >  2:304087/1‑149 (MQ=255)                                                                                                                            ggTATCGACCTTCAGCAGGGCGAAGTAACGCTCCCCTTCCTTCGGTGGGCGGATCTTGCCGACGATGGTGTCGCCGGTACGCAGGTTGAAACGGCGGATCTGGCTGGGCGAGACATAGATATc                                     <  5:56146/123‑1 (MQ=255)                                                                                                                            ggTATCGACCTTCAGCAGGGCGAAGTAACGCTCCCCTTCCTTCGGTGGGCGGATCTTGCCGACGATGGTGTCGCCGGTACGCAGGTTGAAACGGCGGATCTGGCTGGGCGAGACATAGATATc                                     >  6:56146/1‑123 (MQ=255)                                                                                                                               aTCGACCTTCAGCAGGGCGACGTAACGCTCCCCTTCCTTCGGTGGGCGGATCTTGCCGACGATGGTGTCGCCGGTACGCAGGTTGAAACGGCGGATCTGGCTGGGCGAGACATAGATATCGTCTGGGCCGGCCAGGTAGGACGCATCAg        >  7:298523/1‑149 (MQ=255)                                                                                                                                  gACCTTCAGCAGGGCGACGTAACGCTCCCCTTCCTTCGGTGGGCGGATCTTGCCGACGATGGTGTCGCCGGTACGCAGGTTGAAACGGCGGATCTGGCTGGGCGAGACATAGATATCGTCTGGGCCGGCCAGGTAGGACGCATCAGCCg     >  4:260953/1‑149 (MQ=255)                                                                                                                                   aCCTTCAGCAGGGCGACGTAACGCTCCCCTTCCTTCGGTGGGCGGATCTTGCCGACGATGGTGTCGCCGGTACGCAGGTTGAAACGGCGGATCTGGCTGGGCgaga                                               >  8:188913/1‑106 (MQ=255)                                                                                                                                   aCCTTCAGCAGGGCGACGTAACGCTCCCCTTCCTTCGGTGGGCGGATCTTGCCGACGATGGTGTCGCCGGTACGCAGGTTGAAACGGCGGATCTGGCTGGGCGAGACATAGATATCGTCTGGGCCGGCCAGGTAGGACGCATCAGCCGa    <  5:151693/149‑1 (MQ=255)                                                                                                                                     cTTCAGCAGGGCGACGTAACGCTCCCCTTCCTTCGGTGGGCGGATCTTGCCGACGATGGTGTCGCCGGTACGCAGGTTGAAACGGTGGATCTGGCTGGGCGAGACATAGATATCGTCTGGGCCGGCCAGGTAGGACGCATCAGCCGAAc  <  1:304087/149‑1 (MQ=255)                                                                                                                                              ggCGACGTAACGCTCCCCTTCCTTCGGTGGGCGGATCTTGCCGACGATGGTGTCGCCGGTACGCAGGTTGAAACGGCGGATCTGGCTGGGCGAGACATAGATATCGTCTGGGCCGGCCAGGTAGGACGCATCAGc       <  6:294362/135‑1 (MQ=255)                                                                                                                                              ggCGACGTAACGCTCCCCTTCCTTCGGTGGGCGGATCTTGCCGACGATGGTGTCGCCGGTACGCAGGTTGAAACGGCGGATCTGGCTGGGCGAGACATAGATATCGTCTGGGCCGGCCAGGTAGGACGCATCAGc       >  5:294362/1‑135 (MQ=255)                                                                                                                                                   |                                                                                                                                       CCGGTCAAGTCCTCGGTGGAGCCGTTACCGGCCTCCATCTTCAGGCGCTTGTTCGGGAAC‑AGCGGCGTCAGGTTTTCGAACAGGATCTTGTTCTTCGCGTTTTCCGGGCGGTCGAAGTTGATGGTATCGACCTTCAGCAGGGCGAAGTAACGCTCCCCTTCCTTCGGTGGGCGGATCTTGCCGACGATGGTGTCGCCGGTACGCAGGTTGAAACGGCGGATCTGGCTGGGCGAGACATAGATATCGTCTGGGCCGGCCAGGTAGGACGCATCAGCCGAAC  >  NC_002947/5947998‑5948277

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