Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,970,004	C→A	8.5%	G80V (GGC→GTC)	fliC ←	flagellin

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,970,004	0	C	A	8.5%	182.7 / 6.7	59	G80V (GGC→GTC)	fliC	flagellin
Reads supporting (aligned to +/- strand):  ref base C (28/26);  new base A (3/2);  total (31/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.51e-01

GTGTTAGCAATACGGGTGAGCTCGCCCGACAGCGCGGTAAACTCCTGCTGCATGGAGGCGCGGTCGTCATCGCTGTTCGAGCCGTTCGCCGATTGCAGGGCCAGTTCACGCATACGCTGCATGATGTTGCCTGAAGTCGCCATGGCGCCTTCAGCAGTTTGCGCAATAGAGATACCATCGTTGGCGTTCTTGAC‑‑CGCAACGCTCAGACCTTTGATTTGAGTGGTCAGACGGTTGGAGATCTGCAGGCCTGCGGCGTCGTCTTTGGCACTGTTGATTTTCAGGCCAGAGGACAGAC  >  NC_002947/4969857‑4970151                                                                                                                                                    |                                                                                                                                                      gtgtTAGCAATACGGGTGAGCTCGCCCGACAGCGCGGTAAACTCCTGCTGCATGGAGGCGCGGTCGTCATCGCTGTTCGAGCCGTTCGCCGATTGCAGGGCCAGTTCACGCATACGCTGCATGATGTTGCCTGAAGTCGCCATGgcgc                                                                                                                                                       <  2:221486/148‑1 (MQ=255) gtgtTAGCAATACGGGTGAGCTCGCCCGACAGCGCGGTAAACTCCTGCTGCATGGAGGCGCGGTCGTCATCGCTGTTCGAGCCGTTCGCCGATTGCAGGGCCAGTTCACGCATACGCTGCATGATGTTGCCTGAAGTCGCCATGGCGcc                                                                                                                                                      <  5:263963/149‑1 (MQ=255)      aGCAATACGGGTGAGCTCGCCCGACAGCGCGGTAAACTCCTGCTGCATGGAGGCGCGGTCGTCATCGCTGTTCGAGCCGTTCGCCGATTGCAGGGCCAGTTCACGCATACGCTGCATGATGTTGCCTGAAGTCGCCATGGCGCCTTcag                                                                                                                                                 >  8:53181/1‑149 (MQ=255)       gCAATACGGGTGAGCTCGCCCGACAGCGCGGTAAACTCCTGCTGCATGGAGGCGCGGTCGTCATCGCTGTTCGAGCCGTTCGCCGATTGCAGGGCCAGTTCACGCATACGCTGCATGATGTTGCCTGAAGGCGCCATGGCGACTTcag                                                                                                                                                 <  2:235048/148‑1 (MQ=255)       gCAATACGGGTGAGCTCGCCCGACAGCGCGGTAAACTCCTGCTGCATGGAGGCGCGGTCGTCATCGCTGTTCGAGCCGTTCGCCGATTGCAGGGCCAGTTCACGCATACGCTGCATGATGTTGCCTGAAGGCGCCATGGCGACTTcag                                                                                                                                                 <  7:349232/148‑1 (MQ=255)                    gCTCGCCCGACAGCGCGGTAAACTCCTGCTGCATGGAGGCGCGGTCGTCATCGCTGTTCGAGCCGTTCGCCGATTGCAGGGCCAGTTCACGCATACGCTGCATGATGTTGCCTGAAGTCGCCATGGCGCCTTCAGCAGTTTGCGCAATa                                                                                                                                   <  5:196565/149‑1 (MQ=255)                        gCCCGACAGCGCGGTAAACTCCTGCTGCATGGAGGCGCGGTCGTCATCGCTGTTCGAGCCGTTCGCCGATTGCAGGGCCAGTTCACGCATACGCTGCATGATGTTGCCTGAAGTCGCCATGGCGCCTTCAGCAGTTTGCGCAATag                                                                                                                                  <  6:215179/146‑1 (MQ=255)                        gCCCGACAGCGCGGTAAACTCCTGCTGCATGGAGGCGCGGTCGTCATCGCTGTTCGAGCCGTTCGCCGATTGCAGGGCCAGTTCACGCATACGCTGCATGATGTTGCCTGAAGTCGCCATGGCGCCTTCAGCAGTTTGCGCAATAgaga                                                                                                                               <  7:53181/149‑1 (MQ=255)                          ccGACAGCGCGGTAAACTCCTGCTGCATGGAGGCGCGGTCGTCATCGCTGTTCGAGCCGTTCGCCGATTGCAGGGCCAGTTCACGCATACGCTGCATGATGTTGCCTGAAGTCGCCATGGCGCCTTCAGCAGTTTGCGCAATAGAGATa                                                                                                                             >  2:56462/1‑149 (MQ=255)                                     gTAAACTCCTGCTGCATGGAGGCGCGGTCGTCATCGCTGTTCGAGCCGTTCGCCGATTGCAGGGCCAGTTCACACATACGCTGCATGATGTTGCCTGAAGTCGCCATGGCGCCTTCAGCAGTTTGCGCAATAGAGATACCATCGTTGGc                                                                                                                  <  5:261850/149‑1 (MQ=255)                                           tCCTGCTGCATGGAGGCGCGGTCGTCATCGCTGTTCGAGCCGTTCGCCGATTGCAGGGCCAGTTCACGCATACGCTGCATGATGTTGCCTGAAGTCGCCATGGCGCCTTCAGCAGTTTGCGCAATAGAGATACCATCGTTGGCGTTCtt                                                                                                            >  4:40981/1‑149 (MQ=255)                                             ctgctgCATGGAGGCGCGGTCGTCATCGCTGTTCGAGCCGTTCGCCGATTGCAGGGCCAGTTCACGCATACGCTGCATGATGTTGCCTGAAGTCGCCATGGCGCCTTCAGCAGTTTGCGCAATAGAGATACCATCGTTGGCGTTCTTGa                                                                                                          <  3:40981/149‑1 (MQ=255)                                                   cATGGAGGCGCGGTCGTCATCGCTGTTCGAGCCGTTCGCCGATTGCAGGGCCAGTTCACGCATACGCTGCATGATGTTGCCTGAAGTCGCCATGGCGCCTTCAGCAGTTTGCGCAATAGAGATACCATCGTTGGCGTTCTTGAC‑‑CGCaa                                                                                                  <  6:350252/149‑1 (MQ=255)                                                         ggCGCGGTCGTCATCGCTGTTCGAGCCGTTCGCCGATTGCAGGGCCAGTTCACGCATACGCTGCATGATGTTGCCTGAAGTCGCCATGGCGCCTTCAGCAGTTTGCGCAATAGAGATACCATCGTTGGCGTTCTTGAC‑‑CGCAACGCTCa                                                                                            <  1:56462/149‑1 (MQ=255)                                                         ggCGCGGTCGTCATCGCTGTTCGAGCCGTTCGCCGATTGCAGGGCCAGTTCACGCATACGCTGCATGATGTTGCCTGAAGGCGCCATGGCGCCTTCAGCAGTTTGCGCAATAGAGATACCATCGTTGGCGTTCTTGAC‑‑CGCAACGCTCa                                                                                            <  5:84146/149‑1 (MQ=255)                                                               gtcgtcATCGCTGTTCGAGCCGTTCGCCGATTGCAGGGCCAGTTCACGCATACGCTGCATGATGTTGCCTGAAGTCGCCATGGCGCCTTCAGCAGTTTGCGCAATAGAGATACCATCGTTGGCGTTCTTGAC‑‑CGCAACGCTCAGACCtt                                                                                      >  7:301779/1‑149 (MQ=255)                                                                   tcATCGCTGTTCGAGCCGTTCGCCGATTGCAGGGCCAGTTCACGCATACGCTGCATGATGTTGCCTGAAGTCGCCATGGCGCCTTCAGCAGTTTGCGCAATAGAGATACCATCGTTGGCGTTCTTGAC‑‑CGCAACGCTCAGACCtttgat                                                                                  >  3:248646/1‑149 (MQ=255)                                                                     aTCGCTGTTCGAGCCGTTCGCCGATTGCAGGGCCAGTTCACGCATACGCTGCATGATGTTGCCTGAAGTCGCCATGGCGCCTTCAGCAGTTTGCGCAATAGAGATACCATCGTTGGCGTTCTTGAC‑‑CGCAACGCTCAGACCtttgattt                                                                                >  4:236214/1‑149 (MQ=255)                                                                     aTCGCTGTTCGAGCCGTTCGCCGATTGCAGGGCCAGTTCACGCATACACTGCATGATGTTGCCTGAAGTCGCCATGGCGCCTTCAGCAGTTTGCGCAATAGAGATACCATCGTTGGCGTTCTTGAC‑‑CGCAACGCTCAGACCtttgattt                                                                                <  1:266862/149‑1 (MQ=255)                                                                      tCGCTGTTCGAGCCGTTCGCCGATTGCAGGGCCAGTTCACGCATACGCTGCATGATGTTGCCTGAAGTCGCCATGGCGCCTTCAGCAGTTTGCGCAATAGAGATACCATCGTTGGCGTTCTTGAC‑‑CGCAACGCTCAGACCtttgattt                                                                                >  1:205586/1‑148 (MQ=255)                                                                        gCTGTTCGAGCCGTTCGCCGATTGCAGGGCCAGTTCACGCATACGCTGCATGATGTTGCCTGAAGTCGCCATGGCGCCTTCAGCAGTTTGCGCAATAGAGATACCATCGTTGGCGTTCTTGAC‑‑CGCAACGCTCAGACCtttgattt                                                                                >  3:179228/1‑146 (MQ=255)                                                                        gCTGTTCGAGCCGTTCGCCGATTGCAGGGCCAGTTCACGCATACGCTGCATGATGTTGCCTGAAGTCGCCATGGCGCCTTCAGCAGTTTGCGCAATAGAGATACCATCGTTGGCGTTCTTGAC‑‑CGCAACGCTCAGACCtttgattt                                                                                <  4:179228/146‑1 (MQ=255)                                                                         cTGTTCGAGCCGTTCGCCGATTGCAGGGCCAGTTCACGCATACGCTGCATGATGTTGCCTGAAGTCGCCATGGCGCCTTCAGCAGTTTGCGCAATAGAGATACCATCGTTGGCGTTCTTGAC‑‑CGCAACGCTCAGACCTTTGATTTGAGt                                                                            <  4:248646/149‑1 (MQ=255)                                                                          tGTTCGAGCCGTTCGCCGATTGCAGGGCCAGTTCACGCATACGCTGCATGATGTTGCCTGAAGTCGCCATGGCGCCTTCAGCAGTTTGCGCAATAGAGATACCATCGTTGGCGTTCTTGAC‑‑CGCAACGCTCAGACCTTTGATTTGAGTg                                                                           >  3:112445/1‑149 (MQ=255)                                                                              cGAGCCGTTCGCCGATTGCAGGGCCAGTTCACGCATACGCTGCATGATGTTGCCTGAAGTCGCCATGGCGCCTTCAGCAGTTTGCGCAATAGAGATACCATCGTTGGCGTTCTTGAC‑‑CGCAACGCTCAGACCTTTGATTTGAGTGGt                                                                         >  6:303840/1‑147 (MQ=255)                                                                              cGAGCCGTTCGCCGATTGCAGGGCCAGTTCACGCATACGCTGCATGATGTTGCCTGAAGTCGCCATGGCGCCTTCAGCAGTTTGCGCAATAGAGATACCATCGTTGGCGTTCTTGAC‑‑CGCAACGCTCAGACCTTTGATTTGAGTGGTc                                                                        <  2:205586/148‑1 (MQ=255)                                                                               gAGCCGTTCGCCGATTGCAGGGCCAGTTCACGCATACGCTGCATGATGTTGCCTGAAGTCGCCATGGCGCCTTCAGCAGTTTGCGCAATAGAGATACCATCGTTGGCGTTCTTGAC‑‑CGCAACGCTCAGACCtttg                                                                                    >  3:44707/1‑135 (MQ=255)                                                                               gAGCCGTTCGCCGATTGCAGGGCCAGTTCACGCATACGCTGCATGATGTTGCCTGAAGTCGCCATGGCGCCTTCAGCAGTTTGCGCAATAGAGATACCATCGTTGGCGTTCTTGAC‑‑CGCAACGCTCAGACCtttg                                                                                    <  4:44707/135‑1 (MQ=255)                                                                                  ccGTTCGCCGATTGCAGGGCCAGATCACGCATACGCTGCATGATGTTGCCTGAAGTCGCCATGGCGCCTTCAGCAGTTTGCGCAATAGAGATACCATCGTTGGCGTTCTTGAC‑‑CGCAACGCTCAg                                                                                           >  5:205664/1‑125 (MQ=255)                                                                                  ccGTTCGCCGATTGCAGGGCCAGATCACGCATACGCTGCATGATGTTGCCTGAAGTCGCCATGGCGCCTTCAGCAGTTTGCGCAATAGAGATACCATCGTTGGCGTTCTTGAC‑‑CGCAACGCTCAg                                                                                           <  6:205664/125‑1 (MQ=255)                                                                                           gATTGCAGGGCCAGTCCACGCATACGCTGCATGATGTTGCCTGAAGTCGCCATGGCGCCTTCAGCAGTTTGCGCAATAGAGATACCATCGTTGGCGTTCTTGAC‑‑CGCAACGCTCAGACCtttgatttga                                                                              >  8:106878/1‑129 (MQ=255)                                                                                                cAGGGCCAGTTCACGCATACGCTGCATGATGTTGCCTGAAGTCGCCATGGCGCCTTCAGCAGTTTGCGCAATAGAGATACCATCGTTGGCGTTCTTGAC‑‑CGCAACGCTCAGACCTTTGATTTGAGTGGTCAGACGGTTGGAGATCTGCa                                                     <  4:175254/149‑1 (MQ=255)                                                                                                 aGGGCCAGTTCACGCATACGCTGCATGATGTTGCCTGAAGTCGCCATGGCGCCTTCAGCAGTTTGCGCAATAGAGATACCATCGTTGGCGTTCTTGAC‑‑CGCAACGCTCAGACCTTTGATTTGAGTGGTCAGACGGTTGGAGATCTGCAg                                                    >  6:42405/1‑149 (MQ=255)                                                                                                   ggCCAGTTCACGCATACGCTGCATGATGTTGCCTGAAGTCGCCATGGCGCCTTCAGCAGTTTGCGCAATAGAGATACCATCGTTGGCGTTCTTGAC‑‑CGCAACGCTCAGACCTTTGATTTGAGTGGt                                                                         <  2:130577/126‑1 (MQ=255)                                                                                                   ggCCAGTTCACGCATACGCTGCATGATGTTGCCTGAAGTCGCCATGGCGCCTTCAGCAGTTTGCGCAATAGAGATACCATCGTTGGCGTTCTTGAC‑‑CGCAACGCTCAGACCTTTGATTTGAGTGGt                                                                         >  1:130577/1‑126 (MQ=255)                                                                                                   ggCCAGTTCACGCATACGCTGCATGATGTTGCCTGAAGTCGCCATGGCGCCTTCAGCAGTTTGCGCAATAGAGATACCATCGTTGGCGTTCTTGAC‑‑CGCAACGCTCAGACCTTTGATTTGAGTGGTCAGACGGTTGga                                                             >  6:313104/1‑138 (MQ=255)                                                                                                   ggCCAGTTCACGCATACGCTGCATGATGTTGCCTGAAGTCGCCATGGCGCCTTCAGCAGTTTGCGCAATAGAGATACCATCGTTGGCGTTCTTGAC‑‑CGCAACGCTCAGACCTTTGATTTGAGTGGTCAGACGGTTGga                                                             <  5:313104/138‑1 (MQ=255)                                                                                                        gTTCACGCATACGCTGCATGATGTTGCCTGAAGTCGCCATGGCGCCTTCAGCAGTTTGCGCAATAGAGATACCATCGTTGGCGTTCTTGAC‑‑CGCAACGCTCAGACCTTTGATTTGAGTGGTCAGACGGTTGGAGATCTGCAGGCCTgcg                                             >  6:119448/1‑149 (MQ=255)                                                                                                             cGCATACGCTGCATGATGTTGCCTGAAGTCGCCATGGCGCCTTCAGCAGTTTGCGCAATAGAGATACCATCGTTGGCGTTCTTGAC‑‑CGCAACGCTCAGACCTTTGATTTGAGTGGTCAGACGGTTGGAGATCTGCAGGCCTGCGGcgtc                                        <  5:303840/149‑1 (MQ=255)                                                                                                               cATACGCTGCATGATGTTTCCTGAAGTCGCCATGGCGCCTTCAGCAGTTTGCGGAATAGAGATACCATCGTTGGCGTTCTTGCCCTCTCAACGCTCAGACCTTTGATTTGAGTGGTCAGACGGTTGGAGATCTGCAGGCCTGCGGcgtc                                        <  1:296422/149‑1 (MQ=255)                                                                                                                  aCGCTGCATGATGTTGCCTGAAGTCGCCATGGCGCCTTCAGCAGTTTGCGCAATAGAGATACCATCGTTGGCGTTCTTGAC‑‑CTGAACGCTCAGACCtttgatttga                                                                              <  7:236862/106‑1 (MQ=255)                                                                                                                  aCGCTGCATGATGTTGCCTGAAGTCGCCATGGCGCCTTCAGCAGTTTGCGCAATAGAGATACCATCGTTGGCGTTCTTGAC‑‑CGCAACGCTCAGACCtttgatttga                                                                              >  8:236862/1‑106 (MQ=255)                                                                                                                  aCGCTGCATGATGTTGCCTGAAGTCGCCATGGCGCCTTCAGCAGTTTGCGCAATAGAGATACCATCGTTGGCGTTCTTGAC‑‑CGCAACGCTCAGACCTTTGATTTGAGTGGTCAGACGGTTGGAGATCTGCAGGCCTGCGGCGTCGTCtt                                   >  6:59768/1‑149 (MQ=255)                                                                                                                    gCTGCATGATGTTGCCTGAAGTCGCCATGGCGCCTTCAGCAGTTTGCGCAATAGAGATACCATCGTTGGCGTTCTTGAC‑‑CGCAACGCTCAGACCTTTGATTTGAGTGGTCAGACGGTTgg                                                              <  5:311608/120‑1 (MQ=255)                                                                                                                    gCTGCATGATGTTGCCTGAAGTCGCCATGGCGCCTTCAGCAGTTTGCGCAATAGAGATACCATCGTTGGCGTTCTTGAC‑‑CGCAACGCTCAGACCTTTGATTTGAGTGGTCAGACGGTTgg                                                              >  6:311608/1‑120 (MQ=255)                                                                                                                    gCTGCATGATGTTGCCTGAAGTCGCCATGGCGCCTTCAGCAGTTTGCGCAATAGAGATACCATCGTTGGCGTTCTTGAC‑‑CGCAACGCTCAGACCTTTGATTTGAGTGGTCAGACGGTTGGAGATCTGCAGGCCTGCGGCGTCGTCTTTg                                 <  5:42405/149‑1 (MQ=255)                                                                                                                      tGCATGATGTTGCCTGAAGTCGCCATGGCGCCTTCAGCAGTTTGCGCAATAGAGATACCATCGTTGGCGTTCTTGAC‑‑CGCAACGCTCAGACCtt                                                                                      >  3:281444/1‑94 (MQ=255)                                                                                                                      tGCATGATGTTGCCTGAAGTCGCCATGGCGCCTTCAGCAGTTTGCGCAATAGAGATACCATCGTTGGCGTTCTTGAC‑‑CGCAACGCTCAGACCtt                                                                                      <  4:281444/94‑1 (MQ=255)                                                                                                                               ttGCCTGAAGTCGCCATGGCGCCTTCAGCAGTTTGCGCAATAGAGATACCATCGTTGGCGTTCTTGAC‑‑CGCAACGCTCAGACCTTTGATTTGAGTGGTCAGACGGTTGGAGATCTGCAGGCCTGCGGCGTCGTCTTTGGCACTGTTGAt                      >  5:31268/1‑149 (MQ=255)                                                                                                                                 tgcTGAAGGCGCCATGGCGACTTCAGCAGTTTGCGCAATAGAGATACCATCGTTGGCGTTCTTGAC‑‑CGCAACGCTCAGACCTTTGATTTGAGTGGTCAGACGGTTGGAGATCTGCAGGCCTGCGGCGTCGTCTTTGGCACTGTTGAttt                    >  5:244233/3‑149 (MQ=255)                                                                                                                                 tgcTGAAGGCGCCATGGCGACTTCAGCAGTTTGCGCAATAGAGATACCATCGTTGGCGTTCTTGAC‑‑CGCAACGCTCAGACCTTTGATTTGAGTGGTCAGACGGTTGGAGATCTGCAGGCCTGCGGCGTCGTCTTTGGCACTGTTGAttt                    >  7:94839/3‑149 (MQ=255)                                                                                                                                  ccTGAAGTCGCCATGGCGCCTTCAGCAGTTTGCGCAATAGAGATACCATCGTTGGCGTTCTTGAC‑‑CGCAACGCTCAGACCTTTGATTTGAGTGGTCAGACGGTTGGAGATCTGCAGGCCTGCGGCGTCGTCTTTGGCACTGTTGAttt                    >  4:151464/1‑148 (MQ=255)                                                                                                                                     gAAGGCGCCATGGCGACTTCAGCAGTTTGCGCAATAGAGATACCATCGTTGGCGTTCTTGAC‑‑CGCAACGCTCAGACCTTTGATTTGAGTGGTCAGACGGTTGGAGATCTGCAGGCCTGCGGCGTCGTCTTTGGCACTGTTGATTTTCAg                >  3:54330/1‑149 (MQ=255)                                                                                                                                         tCGCCATGGCGCCTTCAGCAGTTTGCGCAATAGAGATACCATCGTTGGCGTTCTTGAC‑‑CGCAACGCTCAGACCTTTGATTTGAGTGGTCAGACGGTTGGAGATCTGCAGGCCTGCGGCGTCGTCTTTGGCACTGTTGATTTTCAGGc              >  1:98834/1‑147 (MQ=255)                                                                                                                                         tCGCCATGGCGCCTTCAGCAGTTTGCGCAATAGAGATACCATCGTTGGCGTTCTTGAC‑‑CGCAACGCTCAGACCTTTGATTTGAGTGGTCAGACGGTTGGAGATCTGCAGGCCTGCGGCGTCGTCTTTGGCACTGTTGATTTTCAGGCCa            <  4:112445/149‑1 (MQ=255)                                                                                                                                            ccATGGCGCCTTCAGCAGTTTGCGCAATAGAGATACCATCGTTGGCGTTCTTGAC‑‑CGCAACGCTCAGACCttt                                                                                     >  3:262113/1‑73 (MQ=255)                                                                                                                                            ccATGGCGCCTTCAGCAGTTTGCGCAATAGAGATACCATCGTTGGCGTTCTTGAC‑‑CGCAACGCTCAGACCttt                                                                                     <  4:262113/73‑1 (MQ=255)                                                                                                                                            ccATGGCGCCTTCAGCAGTTTGCGCAATAGAGATACCATCGTTGGCGTTCTTGAC‑‑CGCAACGCTCAGACCTTTGATTTGAGTGGTCAGACGGTTGGAGATCTGCAGGCCTGCGGCGTCGTCTTTgg                                <  4:26964/126‑1 (MQ=255)                                                                                                                                            ccATGGCGCCTTCAGCAGTTTGCGCAATAGAGATACCATCGTTGGCGTTCTTGAC‑‑CGCAACGCTCAGACCTTTGATTTGAGTGGTCAGACGGTTGGAGATCTGCAGGCCTGCGGCGTCGTCTTTgg                                >  3:26964/1‑126 (MQ=255)                                                                                                                                               tGGCGCCTTCAGCAGTTTGCGCAATAGAGATACCATCGTTGGCGTTCTTGAC‑‑CGCAACGCTCAGACCTTTGATTTGAGTGGTCAGACGGTTGGAGATCTGCAGGCCTGCGGCGTCGTCTTTGGCACTGTTGATTTTCAGGCCAGAGGa       >  2:351518/1‑148 (MQ=255)                                                                                                                                               tGGCGCCTTCAGCAGTTTGCGCAATAGAGATACCATCGTTGGCGTTCTTGAC‑‑CGCAACGCTCAGACCTTTGATTTGAGTGGTCAGACGGTTGGAGATCTGCAGGCCTGCGGCGTCGTCTTTGGCACTGTTGATTTTCAGGCCAGAGGa       >  8:205071/1‑148 (MQ=255)                                                                                                                                                   gcCTTCAGCAGTTTGCGCAATAGAGATACCATCGTTGGCGTTCTTGAC‑‑CGCAACGCTCAGACCTTTGATTTGAGTGGTCAGACGGTTGGAGATCTGCAGGCCTGCGGCGTCGTCTTTGGCACTGTTGATTTTCAGGCCAGAGGACAgac  >  5:195338/1‑149 (MQ=255)                                                                                                                                                    |                                                                                                                                                      GTGTTAGCAATACGGGTGAGCTCGCCCGACAGCGCGGTAAACTCCTGCTGCATGGAGGCGCGGTCGTCATCGCTGTTCGAGCCGTTCGCCGATTGCAGGGCCAGTTCACGCATACGCTGCATGATGTTGCCTGAAGTCGCCATGGCGCCTTCAGCAGTTTGCGCAATAGAGATACCATCGTTGGCGTTCTTGAC‑‑CGCAACGCTCAGACCTTTGATTTGAGTGGTCAGACGGTTGGAGATCTGCAGGCCTGCGGCGTCGTCTTTGGCACTGTTGATTTTCAGGCCAGAGGACAGAC  >  NC_002947/4969857‑4970151

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