Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,897,778	C→G	7.7%	intergenic (‑178/+56)	PP_5172 ← / ← PP_5173	hypothetical protein/RND family transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,897,778	0	C	G	7.7%	186.1 / 5.5	64	intergenic (‑178/+56)	PP_5172/PP_5173	hypothetical protein/RND family transporter
Reads supporting (aligned to +/- strand):  ref base C (27/31);  new base G (2/3);  total (30/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.24e-01

GGTTTGAAGGGGTTTTGCGCCTGGCTTGAGGGCAATAGTAGTGAGGCGGGATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTC‑GATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCA‑‑GGGGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAGGCCACGTACAGCGCCGGCAGGAACAGCAGCGTCAGCAAGGTCGCCACGATGATCCC  >  NC_002947/5897635‑5897926                                                                                                                                                |                                                                                                                                                        ggTTTGAAGGGGTTTTGCGCCTGGCTTGAGGGCAATAGTAGTGAGGCGGGATTTTTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTC‑GATg                                                                                                                                                    <  1:325109/148‑1 (MQ=255) ggTTTGAAGGGGTTTTGCGCCTGGCTTGAGGGCAATAGTAGTGAGGCGGGATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTC‑GATg                                                                                                                                                    <  5:61786/148‑1 (MQ=255)      gtaGGGGTTTTGCGCCTGGCTTTAGGGCAATACTAGTGAGGCGCGATTATTCTAAAAAATACTGGTTAAGAATTTTTATATAACCATTTAAGAGAACACCCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTC‑GATTCACCa                                                                                                                                               <  7:39823/146‑1 (MQ=255)      gAAGGGGTTTTGCGCCTGGCTTGAGGGCAATAGTAGTGAGGCGGGATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTC‑GATGCACCAt                                                                                                                                              >  4:254696/1‑149 (MQ=255)      gAAGGGGTTTTGCGCCTGGCTTGAGGGCAATAGTAGTGAGGCGGGATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTC‑GATGCACCAt                                                                                                                                              >  5:355400/1‑149 (MQ=255)       aaGGGGTTTTGCGCCTGGCTTGAGGGCAATAGTAGTGAGGCGGGATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTC‑GATGCACCATc                                                                                                                                             >  8:288909/1‑149 (MQ=255)        aGGGGTTTTGCGCCTGGCTTGAGGGCAATAGTAGTGAGGCGGGATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTC‑GATGCACCATc                                                                                                                                             <  7:63247/148‑1 (MQ=255)         ggggTTTTGCGCCTGGCTTGAGGGCAATAGTAGTGAGGCGGGATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTATGAGCTGGTTTACACGCGAATAATGTGGTA‑GATTCAACATCg                                                                                                                                            >  3:244014/1‑148 (MQ=255)         ggggTTTTGCGCCTGGCTTGAGGGCAATAGTAGTGAGGCGGGATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCATCTGGTAGGAGCGGGTTTACCAGCTAGTACGGTGGTC‑GATGACCCATCg                                                                                                                                            >  8:183337/1‑148 (MQ=255)                 gggcagggCTTGAGGGCAATAGTAGTGCGGCGGGATTATTATTAAAAATACTGTTTAAGAATTTTTAAATAACCATTTAAGAGACCACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATGCGATGGTG‑CATCGACCACCGTATTcgcg                                                                                                                                    <  6:48598/142‑1 (MQ=21)                 gcgcCTGGCTTGAGGGCAATAGTAGTGAGGCGGGATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATGCGATGGTG‑CATCGACCACCGTATTcgcg                                                                                                                                    <  2:7755/148‑1 (MQ=25)                 gcgcCTGGCTTGAGGGCAATAGTAGTGAGGCGGGATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTC‑GATGCACCATCGCATTCGCgg                                                                                                                                   >  1:234220/1‑149 (MQ=255)                  cgcCTGGCTTGAGGGCAATAGTAGTGAGGCGGGATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTC‑GATGCACCATCGCATTCGCggg                                                                                                                                  <  6:355400/149‑1 (MQ=255)                  cgcCTGGCTTGAGGGCAATAGTAGTGAGGCGGGATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTC‑GATGCACCATCGCATTCGCgg                                                                                                                                   <  6:230403/148‑1 (MQ=255)                  cgcCTGGCTTGAGGGCAATAGTAGTGAGGCGGGATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTC‑GATGCACCATCGCATTCGCgg                                                                                                                                   <  7:288909/148‑1 (MQ=255)                  cgcCTGGCTTGAGGGCAATAGTAGTGAGGCGGGATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTC‑GATGCACCATCGCATTCGCgg                                                                                                                                   <  3:254696/148‑1 (MQ=255)                    ccTGGCTTGAGGGCAATAGTAGTGAGGCGGGATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTC‑GATGCACCATCGCATTCGCGGGTa                                                                                                                                <  4:266872/149‑1 (MQ=255)                      tGGCTTGAGGGCAATAGTAGTGAGGCGGGATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTACGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTC‑GATGCACCATCGCATTCGCGGGTaaa                                                                                                                              <  7:86616/149‑1 (MQ=255)                          ttGAGGGCAATAGTAGTGAGGCGGGATTATTCTTAAAAATACTGTTTAAGAATTTTTATTTAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTTACCGCGAATACGGTGGTG‑GATGAACCATAGCATTAGCGGGTAAAccag                                                                                                                          >  5:249578/1‑147 (MQ=255)                          ttGAGGGCAATAGTAGTGAGGCGGGATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTC‑GATGCACCATCGCATTCGCGGGTAAACCCg                                                                                                                          >  6:271476/1‑149 (MQ=255)                          ttGAGGGCAATAGTAGTGAGGCGGGATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTC‑GATGCACCATCGCATTCGCGGGTAAACCCg                                                                                                                          >  1:19732/1‑149 (MQ=255)                            gAGGGCAATAGTAGTGAGGCGGGATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTACGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTC‑GATGCACCATCGCATTCGCGGGTAAACCCGc                                                                                                                         <  7:331827/148‑1 (MQ=255)                                     agtagtGAGGCGGGATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTC‑GATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCA‑‑gg                                                                                                             >  4:118480/1‑149 (MQ=255)                                            aGGCGGGATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCAGCGAATACGGTGGTC‑GATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCA‑‑GGGGCCata                                                                                                      >  5:282254/1‑149 (MQ=255)                                              gCGGGATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTC‑GATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCA‑‑GGGGCCATATg                                                                                                    >  6:385975/1‑149 (MQ=255)                                                gggATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTC‑GATGCACCATCGCATTCGCgg                                                                                                                                   >  2:166345/1‑118 (MQ=255)                                                gggATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTC‑GATGCACCATCGCATTCGCgg                                                                                                                                   <  1:166345/118‑1 (MQ=255)                                                          ttAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTC‑GATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCA‑‑GGGGCCATATGGGCCTTCAGGCg                                                                                        >  2:315040/1‑149 (MQ=255)                                                            aaaaaTAATGATTAAGAATTTTTATATAACCATTTAAGAGAACACGACTCTGGTAGTAGCGGGTTTACCCGCGAATACGGTGGGC‑GATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCA‑‑GGGGCCATATGGGCCTTCAGTCgt                                                                                       <  7:183337/148‑1 (MQ=38)                                                             aaaaTACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCTAATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCA‑‑GGGGCCATATGGGCCTTCAGGCgt                                                                                       <  1:116971/148‑1 (MQ=255)                                                                    tGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTAACCCGCGAATACGGTGGTC‑GATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCA‑‑GGGGCCATATGGGCCTTCAGGCGTGTTTCTGCg                                                                              <  3:307021/149‑1 (MQ=255)                                                                          aGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTC‑GATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCA‑‑GGGGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTcgc                                                                        <  2:234220/149‑1 (MQ=255)                                                                             aTTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTC‑GATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCA‑‑GGGGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGa                                                                     <  6:282254/149‑1 (MQ=255)                                                                              tttttATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCCCGAAAACGGTGGTC‑GATTCACCATAGCATTCGCGGGTAAACCCGCTCCTACCA‑‑GGGGCCATATGGGCCTTCAAGCGTGTTTCTTCGGAACGCGGa                                                                     >  4:83753/1‑148 (MQ=37)                                                                               ttttATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTAACCGCGAATACGGTGGTC‑GATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCA‑‑GGGGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATc                                                                   >  7:299961/1‑149 (MQ=255)                                                                                    tataACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTC‑GATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCA‑‑GGGGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTa                                                              <  5:271476/149‑1 (MQ=255)                                                                                      taACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTC‑GATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCA‑‑GGGGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTAc                                                             <  2:19732/148‑1 (MQ=255)                                                                                       aaCCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTC‑GATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCA‑‑GGGGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTAccc                                                           >  7:214824/1‑148 (MQ=255)                                                                                             ttAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTC‑GATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCA‑‑GGGGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAGGCCAc                                                     >  6:314507/1‑149 (MQ=255)                                                                                                     acacGCGTGTGGTTGGAGGGGGTTTACCGTCGAATACGTTGGTC‑GATGCACCATTGCATTCGGGGTTAAACCCGCTCCTACCA‑‑GGGGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAGGCCACGTACAgc                                              <  3:83753/148‑1 (MQ=255)                                                                                                     acacGCCTCTGGTTGGAGCGGGTTTACCCGCGAATACGGTGGTC‑GATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCATCGGGGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAGGCCACGTACAg                                               <  1:315040/149‑1 (MQ=255)                                                                                                      cacaCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTC‑GATGCACCATCGCATTCGCGGGTAAACACGCTCCTACCA‑‑GGGGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTAACAGGCCACGTACAgcgc                                            >  6:370972/1‑149 (MQ=255)                                                                                                       acGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTC‑GATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCA‑‑GGGGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAGGCCACGTACAGCGcc                                           <  5:385975/149‑1 (MQ=255)                                                                                                         gCCTCTGGTAGGAGCGGGTTTACCGGCGAATACGGTGGGC‑GATGCACCATCGCATTCGGGGGTAAACCCGCTCCTACCA‑‑GGGGCGATATGGGCCTTCAGGCGTGTTTCTGCg                                                                              <  5:358360/112‑1 (MQ=37)                                                                                                              ttgTAGGTGCGGGTTTACCCGCGAATACGGTGGTC‑GATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCA‑‑GGGGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAGGCCACGTACAGCGCCGGCAg                                      <  7:62790/145‑1 (MQ=255)                                                                                                                gTAGGAGCGGGTTTACCCGCGAATACGGTGGTC‑GATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCA‑‑GGGGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAGGCCACGTACAGCGCCGGCAGGaa                                   <  5:314507/148‑1 (MQ=255)                                                                                                                gTAGGAGCGGGTTTACCCGCGAATACGGTGGTC‑GATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCA‑‑GGGGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAGGCCACGTACAGCGCCGGCAGGaa                                   <  3:118480/148‑1 (MQ=255)                                                                                                                gTAGGAGCGGGTTTACCCGCGAATACGGTGGTC‑GATGCACCATCCCATTCGCGGGTAAACCCGCTCCTACCA‑‑GGGGCCATATGGGCCTTCACGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAGGCCACGTACAGCGCCGGCAGGaa                                   <  2:86498/148‑1 (MQ=255)                                                                                                                    gAGCGGGTTTACCCGCGAATACGGTGGTC‑GATGCAGCATCGCAATCGCGGGTAAACCCGCTCCAACCA‑‑GGGCCCAAATGGGCCTTCAGGCGTGGTTCTGCGGCTCCCCGATCATGAACCAGGCCACGGTCAGCGCCCGCAGGAAcag                                >  8:272411/1‑147 (MQ=37)                                                                                                                      gCGGGTTTACCCGCGAAAACGGTGGTC‑GATGCACCATCGCAATCGCGGGTAAACCCGCTCCAACCA‑‑GGGGCCATATGGGCCATCAGGCGTCTTTCTGGGGCTCGCGGAACTTGTACCAGGCCACGGACAGCGCCGGCAGGCACagccgc                            >  8:127639/1‑149 (MQ=255)                                                                                                                                ccGCGAATGCGATGGTG‑CATCGACCACCGTATTCGCGGGTAAACCCGCTCCTACCA‑‑GGGGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAg                                                          <  7:182204/109‑1 (MQ=18)                                                                                                                                ccGCGAATGCGATGGTG‑CATCGACCACCGTATTCGCGGGTAAACCCGCTCCTACCA‑‑GGGGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAg                                                          >  8:182204/1‑109 (MQ=18)                                                                                                                                ccGCGAATACGGTGGTC‑GATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCA‑‑GGGGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAGGCCACGTACAGCGcc                                           >  3:25454/1‑124 (MQ=255)                                                                                                                                ccGCGAATACGGTGGTC‑GATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCA‑‑GGGGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAGGCCACGTACAGCGcc                                           <  4:25454/124‑1 (MQ=255)                                                                                                                                ccGCGAATACGGTGGTC‑GATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCA‑‑GGGGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAGGCCACGTACAGCGCCGGCAGGAACAGCAGCGTCAGCAAgg                  >  1:110785/1‑149 (MQ=255)                                                                                                                                ccGCGAATACGGTGGTC‑GATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCA‑‑GGGGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAGGCCACGTACAGCGCCGGCAGGAACAGCAGCGTCAGCAAgg                  >  3:156446/1‑149 (MQ=255)                                                                                                                                ccGCGAATACGGTGGTC‑GATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCA‑‑GGGGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAGGCCACGTACAGCGCCGGCAGGAACAGCAGCGTCAGCAAgg                  >  6:163442/1‑149 (MQ=255)                                                                                                                                 cgcgAATACGGTGGTC‑GATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCA‑‑GGGGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAGGCCACGTACAGCGCCGGCAGGAACAGCAGCGTCAGCAAGGt                 <  3:88406/149‑1 (MQ=255)                                                                                                                                 cgcgAATACGGTGGTC‑GATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCA‑‑GGGGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAGGCCACGTACAGCGCCGGCAGGAACAGCAGCGTCAGCAAGGt                 >  4:88406/1‑149 (MQ=255)                                                                                                                                  gcgAATACGGTGGTC‑GATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCA‑‑GGGGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAGGCCACGTACAGCGCCGGCAGGAACAGCAGCGTCAGCAAGGTc                >  3:215382/1‑149 (MQ=255)                                                                                                                                    gAATACGGTGGTC‑GATGCACCATCGCATTCGCGGGTAAACCCGCTc                                                                                                                       >  3:79796/1‑46 (MQ=255)                                                                                                                                    gAATACGGTGGTC‑GATGCACCATCGCATTCGCGGGTAAACCCGCTc                                                                                                                       <  4:79796/46‑1 (MQ=255)                                                                                                                                    gAATACGGTGGTC‑GATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCA‑‑GGGGCCATATGGGCCTTCAGGCgt                                                                                       <  3:98511/76‑1 (MQ=255)                                                                                                                                    gAATACGGTGGTC‑GATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCA‑‑GGGGCCATATGGGCCTTCAGGCgt                                                                                       >  4:98511/1‑76 (MQ=255)                                                                                                                                                c‑GATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCA‑‑GGGGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAGGCCACGTACAGCGCCGGCAGGAACAGCAGCGTCAGCAAGGTCGCCACGATGATccc  <  2:110785/149‑1 (MQ=255)                                                                                                                                                |                                                                                                                                                        GGTTTGAAGGGGTTTTGCGCCTGGCTTGAGGGCAATAGTAGTGAGGCGGGATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTC‑GATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCA‑‑GGGGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAGGCCACGTACAGCGCCGGCAGGAACAGCAGCGTCAGCAAGGTCGCCACGATGATCCC  >  NC_002947/5897635‑5897926

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