Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,763,486	G→A	6.6%	A142V (GCC→GTC)	PP_2417 ←	iron ABC transporter permease

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,763,486	0	G	A	6.6%	180.6 / 2.1	60	A142V (GCC→GTC)	PP_2417	iron ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base G (31/25);  new base A (2/2);  total (33/27)
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.98e-01

GCCAGAACAGGATGCCGCGCGCCTGTTCGGCGGTGGCCGACTTGGTGATCAGGAAGGCGGTCAGTGC‑ATTGAACAGCTGCGAGCCGGCGATACCGGCAAGGATGACCTGTGCATTGTTGCTGCTGGGGCCGGCCGCGCGGGCCAGCACCAGCACCAGGGCAAATGCGGCCCCGGCGCCGATGAAGGCCCCAGCAGACATGCTTAATGCCAGGCTACCCAGCCCCAGCAGGCCGACCAGTACCGCGCCGGTCGAGGCACCGGCAGACAGCCCCAGCAGGTAGGGCTCG  >  NC_002947/2763346‑2763632                                                                                                                                              |                                                                                                                                                   gcCAGAACAGGATGCCGCGCGCCTGTTCGGCGGTGGCCGACTTGGTGATCAGGAAGGCGGTCAGTGC‑ATTGAACAGCTGCGAGCCGGCGATACCGGCAAGGATGACCTGTGCATTGTTGCTGCTGGGGCCGGCCGCGCGGAcc                                                                                                                                                  <  1:363378/143‑1 (MQ=255) gcCAGAACAGGATGCCGCGCGCCTGTTCGGCGGTGGCCGACTTGGTGATCAGGAAGGCGGTCAGTGC‑ATTGAACAGCTGCGAGCCGGCGATACCGGCAAGGATGACCTGTGCATTGTTGCTGCTGGGGCCGGCCGCGCGGAcc                                                                                                                                                  >  2:363378/1‑143 (MQ=255)  ccAGAACAGGATGCCGCGCGCCTGTTCGGCGGTGGCCGACTTGGTGATCAGGAAGGCGGTCAGTGC‑ATTGAACAGCTGCGAGCCGGCGATACCGGCAAGGATGACCTGTGCATTGTTGCTGCTGGGGCCGGCCGCGCGGGccagcac                                                                                                                                             >  8:215947/1‑147 (MQ=255)    aGAACAGGATGCCGCGCGCCTGTTCGGCGGTGGCCGACTTGGTGATCAGGAAGGCGGTCAGTGC‑ATTGAACAGCTGCGAGCCGGCGATACCGGCAAGGATGACCTGTGCATTGTTGCTGCTGGGGCCGGCCGCGCGGAccagcaccag                                                                                                                                          >  6:397890/1‑148 (MQ=255)      aaCAGGATGCCGCGCGCCTGTTCGGCGGTGGCCGACTTGGTGATCAGGAAGGCGGTCAGTGC‑ATTGAACAGCTGCGAGCCGGCGATACCGGCAAGGATGACCTGTGCATTGTTGCTGCTGGGGCCGGCCGCGCGGGCcagcaccagcac                                                                                                                                       >  8:40498/1‑149 (MQ=255)       aCAGGATGCCGCGCGCCTGTTCGGCGGTGGCCGACGGGGTGATCAGGAAGGCGGGCAGTGC‑ATTGAACAGCTGCGAGCCGGCGATACCGGCAGGGATGCGCTGTGCTTTGTTGCTGCTGGGGCCGGCCGCGCGGGCCagcgcccgcccc                                                                                                                                      >  2:112316/1‑149 (MQ=255)              gCCGCGCGCCTGTTCGGCGGTGGCCGACTTGGTGATCAGGAAGGCGGTCAGTGC‑ATTGAACAGCTGCGAGCCGGCGATACCGGCAAGGATGACCTGTGCATTGTTGCTGCTGGGGCCGGCCGCGCGGGCCAGCACCAGCACCAGGGCa                                                                                                                                >  6:191892/1‑148 (MQ=255)                cgcgcgCCTGTTCGGCGGTGGCCGACTTGGTGATCAGGAAGGCGGTCAGTGC‑ATTGAACAGCTGCGAGCCGGCGATACCGGCAAGGATGACCTGTGCATTGTTGCTGCTTGGGACGGCCGCGCGGGACAGCACCAGCACCAGCGCAAAt                                                                                                                             >  8:269482/1‑149 (MQ=255)                 gcgcgcCTGTTCGGCGGTGGCCGACTTGGTGATCAGGAAGGCGGTCAGTGC‑ATTGAACAGCTGCGAGCCGGCGATACCGGCAAGGATGACCTGTGCATTGTTGCTGCTGGGGCCGGCCGCGCGGGCCAGCACCAGCACCAGGGCAAATg                                                                                                                            >  6:199631/1‑149 (MQ=255)                  cgcgcCTGTTCGGCGGTGGCCGACTTGGTGATCAGGAAGGCGGTCAGTGC‑ATTGAACAGCTGCGAGCCGGCGATACCGGCAAGGATGACCTGTGCATTGTTGCTGCTGGGGCCGGCCGCGCGGGCCAGCACCAGCACCAGGGCAAATGc                                                                                                                           <  3:384566/149‑1 (MQ=255)                  cgcgcCTGTTCGGCGGTGGCCGACTTGGTGATCAGGAAGGCGGTCAGTGC‑ATTGAACAGCTGCGAGCCGGCGATACCGGCAAGGATGACCTGTGCATTGTTGCTGCTGGGGCCGGCCGCGCGGGCCAGCACCAGCACCAGGGCAAATGc                                                                                                                           <  6:387020/149‑1 (MQ=255)                       cTGTTCGGCGGTGGCCGACTTGGTGATCAGGAAGGCGGTCAGTGCGATTGAACAGCTGCGAGCCGGCGATACCGGCAAGGATGACCTGTGCATTGTTGCTGCTGGGGCCGGCCGCGCGGGCCAGCACCAGCACCAGGGCAAATGCGGcc                                                                                                                       >  1:80610/1‑149 (MQ=255)                       cTGTTCGGCGGTGGCCGACTTGGTGATCAGGAAGGCGGTCAGTGC‑ATTGAACAGCTGCGAGCCGGCGATACCGGCAAGGATGACCTGTGCATTGTTGCTGCTGGGGCCGGCCGCGCGGGCCAGCACCAGCACCAGGGCAAATGCGGccc                                                                                                                      >  6:177932/1‑149 (MQ=255)                                 gTGGCCGACTTGGTGATCAGGAAGGCGGTCAGTGC‑ATTGAACAGCTGCGAGCCGGCGATACCGGCAAGGATGACCTGTGCATTGTTGCTGCTGGGGCCGGCCGCGCGGACCAGCACCAGCACCAGGGCAAATGCGGCCCCGGCGCCGAt                                                                                                            <  5:397890/149‑1 (MQ=255)                                    gCCGACTTGGTGATCAGGAAGGCGGTCAGTGC‑ATTGAACAGCTGCGAGCCGGCGATACCGGCAAGGATGACCTGTGCATTGTTGCTGCTGGGGCCGGCCGCGCGGGCCAGCACCAGCACCAGGGCAAATGCGGCCCCGGCGCCGATGa                                                                                                          <  5:199631/148‑1 (MQ=255)                                      cGACTTGGTGATCAGGAAGGCGGTCAGTGC‑ATTGAACAGCTGCGAGCCGGCGATACCGGCAAGGATGACCTGTGCATTGTTGCTGCTGGGGCCGGCCGCGCGGGCCAGCACCAGCACCAGGGCAAATGCGGCCCCGGCGCCGATGAAgg                                                                                                       >  8:186992/1‑149 (MQ=255)                                          ttGGTGATCAGGAAGGCGGTCAGTGC‑ATTGAACAGCTGCGAGCCGGCGATACCGGCAAGGATGACCTGTGCATTGTTGCTGCTGGGGCCGGCCGCGCGGGCCAGCACCAGCACCAGGGCAAATGCGGCCCCGGCGCCGATGAAGGccc                                                                                                    >  8:209633/1‑148 (MQ=255)                                             gTGATCAGGAAGGCGGTCAGTGC‑ATTGAACAGCTGCGAGCCGGCGATACCGGCAAGGATGACCTGTGCATTGTTGCTGCTGGGGCCGGCCGCGCGGGCCAGCACCAGCACCAGGGCAAATGCGGCCCCGGCGCCGATGAAGGCCCcag                                                                                                 >  1:382134/1‑148 (MQ=255)                                               gATCAGGAAGGCGGTCAGTGC‑ATTGAACAGCTGCGAGCCGGCGATACCGGCAAGGATGACCTGTGCATTGTTGCTGCTGGGGCCGGCCGCGCGGGCCAGCACCAGCACCAGGGCAAATGCGGCCCCGGCGCCGATGAAGGCCCcagcag                                                                                              >  8:301183/1‑149 (MQ=255)                                                aTCAGGAAGGCGGTCAGTGC‑ATTGAACAGCTGCGAGCCGGCGATACCGGCAAGGATGACCTGTGCATTGTTGCTGCTGGGGCCGGCCGCGCGGGCCAGCACCAGCACCAGGGCAAATGCGGCCCCGGCGCCGATGAAGGCCCcagcag                                                                                              <  2:14099/148‑1 (MQ=255)                                                 tCAGGAAGGCGGTCAGTGC‑ATTGAACAGCTGCGAGCCGGCGATACCGGCAAGGATGACCTGTGCATTGTTGCTGCTGGGGCCGGCCGCGCGGGCCAGCACCAGCACCAGGGCAAATGCGGCCCCGGCGCCGATGAAGGCCCCAGCAGAc                                                                                            <  8:77224/149‑1 (MQ=255)                                                   aGGAAGGCGGTCAGTGC‑ATTGAACAGCTGCGAGCCGGCGATACCGGCAAGGATGACCTGTGCATTGTTGCTGCTGGGGCCGGCCGCGCGGGCCAGCACCAGCACCAGGGCAAATGCGGCCCCGGCGCCGATGAAGGCCCCAGCAGACa                                                                                           <  3:86449/148‑1 (MQ=255)                                                    ggAAGGCGGTCAGTGC‑ATTGAACAGCTGCGAGCCGGCGATACCGGCAAGGATGACCTGTGCATTGTTGCTGCTGGGGCCGGCCGCGCGGGCCAGCACCAGCACCAGGGCAAATGCGGCCCCGGCGCCGATGAAGGCCCCAGCAGACATg                                                                                         >  2:261639/1‑149 (MQ=255)                                                     gAAGGCGGTCAGTGC‑ATTGACAAGCTGCGATCAGTCGATACCGGAAAGGATGACCTGTGCATTGTTGCTTCTGGGGCCGGCCGCGCTGGCCAGAACCAGCACCAGGGCAAATGCGGCCCCGGCGCCGATGAAGGCCCCAGCAGACATGc                                                                                        <  5:61100/149‑1 (MQ=255)                                                     gAAGGCGGTCAGTGC‑ATTGAACAGCTGCGAGCCGGCGATACCGGCAAGGATGACCTGTGCATTGTTGCTGCTGGGGCCGGCCGCGCGGGCCAGCACCAGCACCAGGGCAAATGCGGCCCCGGCGCCGATGAAGGCCCCAGCAGACATGc                                                                                        <  8:354042/149‑1 (MQ=255)                                                      aaGGCGGTCAGTGC‑ATTGAACAGCTGCGAGCCGGCGATACCGGCAAGGATGACCTGTGCATTGTTGCTGCTGGGGCCGGCCGCGCGGGCCAGCACCAGCACCAGGGCAAATGCGGCCCCGGCGCCGATGAAGGCCCCAGCAGACATGCt                                                                                       <  7:269482/149‑1 (MQ=255)                                                          cGGTCAGTGC‑ATTGAACAGCTGCGAGCCGGCGATACCGGCAAGGATGACCTGTGCATTGTTGCTGCTGGGGCCGGCCGCGCGGGCCAGCACCAGCACCAGGGCAAATGCGGCCCCGGCGCCGATGAAGGCCCCAGCAGACATGCTTAAt                                                                                   <  7:209633/149‑1 (MQ=255)                                                          cGGGCAGTGC‑ATTGAACAGCTGCGAGCCGGCGATACCGGCAAGGATCACGTGTGCATTGTTGCTGCTGGGGCCGGCCGCGCGGGCGAGCACCAGCACCAGGGCAAATGCGGCCCCGGCGCCGATGAAGGCCCCAGCAGACATGCTTAAt                                                                                   <  7:216974/149‑1 (MQ=255)                                                               aGTGC‑ATTGAACAGCTGCGAGCCGGCGATACCGGCAAGGATGACCTGTGCATTGTTGCTGCTGGGGCCGGCCGCGCGGGCCAGCACCAGCACCAGGGCAAATGCGGCCCCGGCGCCGATGAAGGCCCCAGCAGACATGCTTAATGCCa                                                                               <  6:222074/148‑1 (MQ=255)                                                                gTGC‑ATTGAACAGCTGCGAGCCGGCGATACCGGCAAGGATGACCTGTGCATTGTTGCTGCTGGGGCCGGCCGCGCGGGCCAGCACCAGCACCAGGGCAAATGCGGCCCCGGCGCCGATGAAGGCCCCAGCAGACATGCTTAATGCCAgg                                                                             >  5:379670/1‑149 (MQ=255)                                                                     aTTGAACAGCTGCGAGCCGGCGATACCGGCAAGGATGACCTGTGCATTGTTGCTGCTGGGGCCGGCCGCGCGGGCCAGCACCAGCACCAGGGCAAATGCGGCCCCGGCGCCGATGAAGGCCCCAGCAGACATGCTTAATGCCAGGCTAc                                                                         >  3:138258/1‑149 (MQ=255)                                                                      ttGAACAGCTGCGAGCCGGCGATACCGGCAAGGATGACCTGTGCATTGTTGCTGCTGGGGCCGGCCGCGCGGGCCAGCACCAGCACCAGGGCAAATGCGGCCCCGGCGCCGATGAAGGCCCCAGCAGACATGCTTAATGCCAGGCTAcc                                                                        >  5:122887/1‑149 (MQ=255)                                                                      ttGAACAGCTGCGAGCCGGCGATACCGGCAAGGATGACCTGTGCATTGTTGCTGCTGGGGCCGGCCGCGCGGGCCAGCACCAGCACCAGGGCAAATGCGGCCCCGGCGCCGATGAAGGCCCCAGCAGACATGCTTAATGCCAGGCTAcc                                                                        >  4:286960/1‑149 (MQ=255)                                                                         aaCAGCTGCGAGCCGGCGATACCGGCAAGGATGACCTGTGCATTGTTGCTGCTGGGGCCGGCCGCGCGGGCCAGCACCAGCACCAGGGCAAATGCGGCCCCGGCGCCGAGGAAGGCCCCAGCAGACATCCATaa                                                                                    >  4:219192/1‑134 (MQ=255)                                                                         aaCAGCTGCGAGCCGGCGATACCGGCAAGGATGACCTGTGCATTGTTGCTGCTGGGGCCGGCCGCGAGGGCCAGCACCAGCACCAGGGCAAATGCGGCCCCGGCGCCGATGAAGGCCCCAGCAGACATGCTTaa                                                                                    <  3:219192/134‑1 (MQ=255)                                                                           cAGCTGCGAGCCGGCGATACCGGCAAGGATGACCTGTGCATTGTTGCTGCTGGGGCCGGCCGCGCGGGCCAGCACCAGCACCAGGGCAAATGCGGCCCCGGCGCCGATGAAGGCCCCAGCAGACATGCTTAATGCCAGGCTAcccagcc                                                                   >  3:389776/1‑149 (MQ=255)                                                                            aGCTGCGAGCCGGCGATACCGGCAAGGATGACCTGTGCATTGTTGCTGCTGGGGCCGGCCGCGCGGGCCAGCACCAGCACCAGGGCAAATGCGGCCCCGGCGCCGATGAAGGCCCCAGCAGACATGCTTACTGCCAGGCTAcccagccc                                                                  >  6:273924/1‑149 (MQ=255)                                                                            aGCTGCGAGCCGGCGATACCGGCAAGGATGACCTGTGCATTGTTGCTGCTGGGGCCGGCCGCGCGGGCCAGCACCAGCACCAGGGCAAATGCGGCCCCGGCGCCGATGAAGGCCCCAGCAGACATGCTTAATGCCAGGCTAcccagccc                                                                  <  5:191892/149‑1 (MQ=255)                                                                            aGCTGCGAGCCGGCGATACCGGCAAGGATGACCTGTGCATTGTTGCTGCTGGGGCCGGCCGCGCGGGCCAGCACCAGCACCAGGGCAAATGCGGCCCCGGCGCCGATGAAGGCCCCAGCAGACATGCTTAATGCCAGGCTAcccagccc                                                                  <  7:186992/149‑1 (MQ=255)                                                                                     ccGGCGATACCGGCAAGGATGACCTGTGCATTGTTGCTGCTGGGGCCGGCCGCGCGGGCCAGCACCAGCACCAGGGCAAATGCGGCCCCGGCGCCGATGAAGGCCCCAGCAGACATGCTTAATGCCAGGCTACCCAGCCCCAGCAGGcc                                                         >  8:343803/1‑149 (MQ=255)                                                                                        gCGATACCGGCAAGGATGACCTGTGCATTGTTGCTGCTGGGGCCGGCCGCGCGGGCCAGCACCAGCACCAGGGCAAATGCGGCCCCGGCGCCGATGAAGGCCCCAGCAGACATGCTTAATGCCAGGCTACCCAGCCCCAGCAGGCCGAc                                                      >  2:205144/1‑149 (MQ=255)                                                                                        gCGATACCGGCAAGGATGACCTGTGCATTGTTGCTGCTGGGGCCGGCCGCGCGGGCCAGCACCAGCACCAGGGCAAATGCGGCCCCGGCGCCGATGAAGGCCCCAGCAGACATGCTTAATGCCAGGCTACCCAGCCCCAGCAGGCCGAc                                                      <  6:379670/149‑1 (MQ=255)                                                                                              ccGGCAAGGATGACCTGTGCATTGTTGCTGCTGGGGCCGGCCGCGCGGGCCAGCACCAGCACCAGGGCAAATGCGGCCCCGGCGCCGATGAAGGCCCCAGCAGACATGCTTAATGCCAGGCTACCCAGCCCCAGCAGGCCGACCAGTa                                                 <  7:301183/148‑1 (MQ=255)                                                                                               cGGCAAGGATGACCTGTGCATTGTTGCTGCTGGGGCCGGCCGCGCGGGCCAGCACCAGCACCAGGGCAAATGCGGCCCCGGCGCCGATGAAGGCCCCAGCAGACATGCTTAATGCCAGGCTACCCAGCCCCAGCAGGCCGACCAGTAcc                                               >  2:305823/1‑149 (MQ=255)                                                                                                ggCAAGGATGACCTGTGCATTGTTGCTGCTGGGGCCGGCCGCGCGGGCCAGCACCAGCACCAGGGCAAATGCGGCCCCGGCGCCGATGAAGGCCCCAGCAGACATGCTTAATGCCAGGCTACCCAGCCCCAGCAGGCCGACCAGTACcg                                              <  1:261639/149‑1 (MQ=255)                                                                                                           ccTGTGCATTGTTGCTGCTGGGGCCGGCCGCGCGGGCCAGCACCAGCACCAGGGCAAATGCGGCCCCGGCGCCGATGAAGGCCCCAGCAGACATGCTTAATGCCAGGCTACCCAGCCCCAGCAGGCCGACAAGTACCGCGCCGGTCGAg                                   <  7:343803/149‑1 (MQ=255)                                                                                                            cTGTGCATTGTTGCTGCTGGGGCCGGCCGCGCGGGCCAGCACCAGCACCAGGGCAAATGCGGCCCCGGCGCCGATGAAGGCCCCAGCAGGCATGCTTAATGCCAGGCTACCCAGCCCCAGCAGGCCGACCAGTACcg                                              <  6:86525/137‑1 (MQ=255)                                                                                                            cTGTGCATTGTTGCTGCTGGGGCCGGCCGCGCGGGCCAGCACCAGCACCAGGGCAAATGCGGCCCCGGCGCCGATGAAGGCCCCAGCAGGCATGCTTAATGCCAGGCTACCCAGCCCCAGCAGGCCGACCAGTACcg                                              >  5:86525/1‑137 (MQ=255)                                                                                                                 cATTGTTGCTGCTGGGGCCGGCCGCGCGGGCCAGCACCAGCACCAGGGCAAATGCGGCCCCGGCGCCGATGAAGGCCCCAGCAGACATGCTTAATGCCAGGCTACCCAGCCCCAGCAGGCCGACCAGTACCGCGCCGGTCGAGGCACCg                             >  6:164432/1‑149 (MQ=255)                                                                                                                  aTTGTTGCTGCTGGGGCCGGCCGCGCGGGCCAGCACCAGCACCAGGGCAAATGCGGCCCCGGCGCCGATGAAGGCCCCAGCAGACATGCTTAATGCCAGGCTACCCAGCCCCAGCAGGCCGACCAGTACCGCGCCGGTCGAGGCACCgg                            <  5:177932/149‑1 (MQ=255)                                                                                                                   ttgttgCTGCTGGGGCCGGCCGCGCGGGCCAGCACCAGCCCCAGGGCAAATGCGGCCCCGGCGCCGATGAAGGCCCCAGCAGACATGCTTAATGCCAGGCTACCCAGCCCCCGCAGGCCGACCAGTACCGCGCCGGTCGAGGCACCGGc                           >  3:14942/1‑149 (MQ=255)                                                                                                                      ttgCTGCTGGGGCCGGCCGCGCGGGCCAGCACCAGCACCAGGGCAAATGCGGCCCCGGCGCCGATGAAGGCCCCAGCAGACATGCTTAATGCCAGGCTACCCAGCCCCAGCAGGCCGACCAGTACCGCGCCGGTCGAGGCACCGGCAGa                        >  1:379268/1‑149 (MQ=255)                                                                                                                               gggCCGGCCGCGCGGGCCAGCACCAGCACCAGGGCAAATGCGGCCCCGGCGCCGATGAAGGCCCCAGCAGACATGCTTAATGCCAGGCTACCCAGCCCCAGCAGGCCGACCAGTACCGCGCCGGTCGAGGCACCGGCAGACAGCCCcag               <  5:273924/149‑1 (MQ=255)                                                                                                                               gggCCGGCCGCGCGGGCCAGCACCAGCACCAGGGCAAATGCGGCCCCGGCGCCGATGAAGGCCCCAGCAGACATGCTTAATGCCAGGCTACCAAGCCCCAGCAGGCCGACAAGTACCGCGCCGGTCGAGTCACCGACAGACAGCCCcag               >  6:321179/1‑149 (MQ=255)                                                                                                                                ggccggccGCGCGGGCCAGCACCAGCACCAGGGCAAATGCGGCCCCGGCGCCGATGAAGGCCCCAGCAGACATGCTTAATGCCAGGCTACCCAGCCCCAGCAGGCCGCCCAGTACCGCGCCGGTCGAGGCACCGGCAGACAGCCCcag               >  1:246502/1‑148 (MQ=255)                                                                                                                                 gccggccgCGCGGGCCAGCACCAGCACCAGGGCAAATGCGGCCCCGGCGCCGATGAAGGCCCCAGCAGACATGCTTAATGCCAGGCTACCCAGCCCCAGCAGGCCGACCAGTACCGCGCCGGTCGAGGCACCGGCAGACAGcccccg               >  2:379606/1‑145 (MQ=255)                                                                                                                                 gccggccgCGCGGGCCAGCACCAGCACCAGGGCAAATGCGGCCCCGGCGCCGATGAAGGCCCCAGCAGACATGCTTAATGCCAGGCTACCCAGCCCCAGCAGGCCGACCAGTACCGCGCCGGTCGAGGCACCGGCAGACAGCCCcagca             <  5:164432/149‑1 (MQ=255)                                                                                                                                     gccgCGCGGGCCGGCACCAGCACCAGGGCAAATGCGGCCCCGGCGCCGATGAAGGCCCCAGCAGACATGCTTAATGCCAGGCTACCCAGCCCCAGCAGGCCGACCAGTACCGCGCCGGTCGAGGCACCGGCAGACAGCCCCAGCAGGTa         <  5:194021/149‑1 (MQ=255)                                                                                                                                     gccgCGCGGGCCAGCACCAGCACCAGGGCAAATGCGGCCCCGGCGCCGATGAAGGCCCCAGCAGACATGCTTAATGCCAGGCTACCCAGCCCCAGCAGGCCGACCAGTACCGCGCCGGTCGAGGCACCGGCAGACAGCCCCAGCAGGTa         <  2:80610/149‑1 (MQ=255)                                                                                                                                      ccgCGCGGGCCAGCACCAGCACCAGGGCAAATGCGGCCCCGGCGCCGATGAAGGCCCCAGCAGACATGCTTAATGCCAGGCTACCCAGCCCCAGCAGGCCGACCAGTACCGCGCCGGTCGAGGCACCGGCAGACAGCCCCAGCAGGTAg        >  4:352002/1‑149 (MQ=255)                                                                                                                                            gggCCAGCACCAGCACCAGGGCAAATGCGGCCCCGGCGCCGATGAAGGCCCCAGCAGACATGCTTAATGCCAGGCTACCCAGCCCCAGCAGGCCGACCAGTACCGCGCCGGTCGAGGCACCGGCAGACAGCCCCAGCAGGTAGGGCTCg  <  3:286960/149‑1 (MQ=255)                                                                                                                                            gggCCAGCACCAGCACCAGGGCAAATGCGGCCCCGGCGCCGATGAAGGCCCCAGCAGACATGCTTAATGCCAGGCTACCCAGCCCCAGCAGGCCGACCAGTACCGCGCCGGTCGAGGCACCGGCAGACAGCCCCAGCAGGTAGGGCTCg  <  4:389776/149‑1 (MQ=255)                                                                                                                                              |                                                                                                                                                   GCCAGAACAGGATGCCGCGCGCCTGTTCGGCGGTGGCCGACTTGGTGATCAGGAAGGCGGTCAGTGC‑ATTGAACAGCTGCGAGCCGGCGATACCGGCAAGGATGACCTGTGCATTGTTGCTGCTGGGGCCGGCCGCGCGGGCCAGCACCAGCACCAGGGCAAATGCGGCCCCGGCGCCGATGAAGGCCCCAGCAGACATGCTTAATGCCAGGCTACCCAGCCCCAGCAGGCCGACCAGTACCGCGCCGGTCGAGGCACCGGCAGACAGCCCCAGCAGGTAGGGCTCG  >  NC_002947/2763346‑2763632

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