Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,969,993	T→G	10.3%	T84P (ACT→CCT)	fliC ←	flagellin

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,969,993	0	T	G	10.3%	164.8 / 8.9	57	T84P (ACT→CCT)	fliC	flagellin
Reads supporting (aligned to +/- strand):  ref base T (24/27);  new base G (3/3);  total (27/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.91e-01

CCGAAAGTGGTGGTGTTAGCAATACGGGTGAGCTCGCCCGACAGCGCGGTAAACTCCTGCTGCATGGAGGCGCGGTCGTCATCGCTGTTCGAGCCGTTCGCCGATTGCAGGGCCAGTTCACGCATACGCTGCATGATGTTGCCTGAAGTCGCCATGGCGCCTTCAGCAGTTTGCGCAATAGAGATACCATCGTTGGCGTTCTTGAC‑‑CGCAACGCTCAGACCTTTGATTTGAGTGGTCAGACGGTTGGAGATCTGCAGGCCTGCGGCGTCGTCTTTGGCACTGTTGATTTTCAGGCCA  >  NC_002947/4969845‑4970141                                                                                                                                                     |                                                                                                                                                       ccGAAAGTGGTGGTGTTAGCAATACGGGTGAGCTCGCCCGACAGCGCGGTAAACTCCTGCTGCATGGAGGCGCGGTCGTCATCGCTGTTCGAGCCGTTCGCCGATTGCAGGGCCAGTTCACGCATACGCTGCATGATGTTGCCTGAAGt                                                                                                                                                        >  5:215179/1‑149 (MQ=255)   gAAAGTGGTGGTGTTAGCAATACGGGTGAGCTCGCCCGACAGCGCGGTAAACTCCTGCTGCATGGAGGCGCGGTCGTCATCGCTGTTCGAGCCGTTCGCCGATTGCAGGGCCAGTTCACGCATACGCTGCATGATGTTGCCTGAAGTCg                                                                                                                                                      <  8:11566/149‑1 (MQ=255)          gtggtgTTAGCAATACGGGTGAGCTCGCCCGACAGCGCGGTAAACTCCTGCTGCATGGAGGCGCGGTCGTCATCGCTGTTCGAGCCGTTCGCCGATTGCAGGGCCAGTTCACGCATACGCTGCATGATGTTGCCTGAAGTCGCCATGgc                                                                                                                                               <  8:169313/149‑1 (MQ=255)             gtgtTAGCAATACGGGTGAGCTCGCCCGACAGCGCGGTAAACTCCTGCTGCATGGAGGCGCGGTCGTCATCGCTGTTCGAGCCGTTCGCCGATTGCAGGGCCAGTTCACGCATACGCTGCATGATGTTGCCTGAAGTCGCCATGgcgc                                                                                                                                             <  2:221486/148‑1 (MQ=255)             gtgtTAGCAATACGGGTGAGCTCGCCCGACAGCGCGGTAAACTCCTGCTGCATGGAGGCGCGGTCGTCATCGCTGTTCGAGCCGTTCGCCGATTGCAGGGCCAGTTCACGCATACGCTGCATGATGTTGCCTGAAGTCGCCATGGCGcc                                                                                                                                            <  5:263963/149‑1 (MQ=255)                  aGCAATACGGGTGAGCTCGCCCGACAGCGCGGTAAACTCCTGCTGCATGGAGGCGCGGTCGTCATCGCTGTTCGAGCCGTTCGCCGATTGCAGGGCCAGTTCACGCATACGCTGCATGATGTTGCCTGAAGTCGCCATGGCGCCTTcag                                                                                                                                       >  8:53181/1‑149 (MQ=255)                   gCAATACGGGTGAGCTCGCCCGACAGCGCGGTAAACTCCTGCTGCATGGAGGCGCGGTCGTCATCGCTGTTCGAGCCGTTCGCCGATTGCAGGGCCAGTTCACGCATACGCTGCATGATGTTGCCTGAAGGCGCCATGGCGACTTcag                                                                                                                                       <  7:349232/148‑1 (MQ=255)                   gCAATACGGGTGAGCTCGCCCGACAGCGCGGTAAACTCCTGCTGCATGGAGGCGCGGTCGTCATCGCTGTTCGAGCCGTTCGCCGATTGCAGGGCCAGTTCACGCATACGCTGCATGATGTTGCCTGAAGGCGCCATGGCGACTTcag                                                                                                                                       <  2:235048/148‑1 (MQ=255)                                gCTCGCCCGACAGCGCGGTAAACTCCTGCTGCATGGAGGCGCGGTCGTCATCGCTGTTCGAGCCGTTCGCCGATTGCAGGGCCAGTTCACGCATACGCTGCATGATGTTGCCTGAAGTCGCCATGGCGCCTTCAGCAGTTTGCGCAATa                                                                                                                         <  5:196565/149‑1 (MQ=255)                                    gCCCGACAGCGCGGTAAACTCCTGCTGCATGGAGGCGCGGTCGTCATCGCTGTTCGAGCCGTTCGCCGATTGCAGGGCCAGTTCACGCATACGCTGCATGATGTTGCCTGAAGTCGCCa                                                                                                                                                   >  3:163169/1‑119 (MQ=255)                                    gCCCGACAGCGCGGTAAACTCCTGCTGCATGGAGGCGCGGTCGTCATCGCTGTTCGAGCCGTTCGCCGATTGCAGGGCCAGTTCACGCATACGCTGCATGATGTTGCCTGAAGTCGCCa                                                                                                                                                   <  4:163169/119‑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                                                                      <  4:179228/146‑1 (MQ=255)                                                                                    gCTGTTCGAGCCGTTCGCCGATTGCAGGGCCAGTTCACGCATACGCTGCATGATGTTGCCTGAAGTCGCCATGGCGCCTTCAGCAGTTTGCGCAATAGAGATACCATCGTTGGCGTTCTTGAC‑‑CGCAACGCTCAGACCtttgattt                                                                      >  3:179228/1‑146 (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                                                                          <  4:44707/135‑1 (MQ=255)                                                                                           gAGCCGTTCGCCGATTGCAGGGCCAGTTCACGCATACGCTGCATGATGTTGCCTGAAGTCGCCATGGCGCCTTCAGCAGTTTGCGCAATAGAGATACCATCGTTGGCGTTCTTGAC‑‑CGCAACGCTCAGACCtttg                                                                          >  3:44707/1‑135 (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                                                               >  1:130577/1‑126 (MQ=255)                                                                                                               ggCCAGTTCACGCATACGCTGCATGATGTTGCCTGAAGTCGCCATGGCGCCTTCAGCAGTTTGCGCAATAGAGATACCATCGTTGGCGTTCTTGAC‑‑CGCAACGCTCAGACCTTTGATTTGAGTGGt                                                               <  2:130577/126‑1 (MQ=255)                                                                                                               ggCCAGTTCACGCATACGCTGCATGATGTTGCCTGAAGTCGCCATGGCGCCTTCAGCAGTTTGCGCAATAGAGATACCATCGTTGGCGTTCTTGAC‑‑CGCAACGCTCAGACCTTTGATTTGAGTGGTCAGACGGTTGga                                                   <  5:313104/138‑1 (MQ=255)                                                                                                               ggCCAGTTCACGCATACGCTGCATGATGTTGCCTGAAGTCGCCATGGCGCCTTCAGCAGTTTGCGCAATAGAGATACCATCGTTGGCGTTCTTGAC‑‑CGCAACGCTCAGACCTTTGATTTGAGTGGTCAGACGGTTGga                                                   >  6:313104/1‑138 (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                                                    >  6:311608/1‑120 (MQ=255)                                                                                                                                gCTGCATGATGTTGCCTGAAGTCGCCATGGCGCCTTCAGCAGTTTGCGCAATAGAGATACCATCGTTGGCGTTCTTGAC‑‑CGCAACGCTCAGACCTTTGATTTGAGTGGTCAGACGGTTgg                                                    <  5:311608/120‑1 (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)                                                                                                                                                     |                                                                                                                                                       CCGAAAGTGGTGGTGTTAGCAATACGGGTGAGCTCGCCCGACAGCGCGGTAAACTCCTGCTGCATGGAGGCGCGGTCGTCATCGCTGTTCGAGCCGTTCGCCGATTGCAGGGCCAGTTCACGCATACGCTGCATGATGTTGCCTGAAGTCGCCATGGCGCCTTCAGCAGTTTGCGCAATAGAGATACCATCGTTGGCGTTCTTGAC‑‑CGCAACGCTCAGACCTTTGATTTGAGTGGTCAGACGGTTGGAGATCTGCAGGCCTGCGGCGTCGTCTTTGGCACTGTTGATTTTCAGGCCA  >  NC_002947/4969845‑4970141

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