Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	6,129,176	A→G	9.9%	intergenic (+38/+29)	PP_5376 → / ← PP_5377	hypothetical protein/hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	6,129,176	0	A	G	9.9%	165.4 / 8.1	62	intergenic (+38/+29)	PP_5376/PP_5377	hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (27/28);  new base G (3/3);  total (30/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.48e-01

GGAAGCTGTGACGTGATGCCTGCGCAAATGACCTATGAAGACTCGACAGGCAAGCTCAACATTCTCGAATACCGCGCAGCTGGCA‑TCAATTGCCGCGGCCAAAATTGAGCTCTGAATGCTGTCCAAACGAAAACGGCGCCTTAAAAGCGCCGTTTGCGTTGGGAAAACGGTGGGTCAGTGGCCTGGATACTCTGCAACCACTTCATCCCGTCCGCTCTTGTTGAGACCGATGACCTGGTATGCATGAGAGACACCATCCACTTCCATACCTGGAGAGCCTGCCGGCATA  >  NC_002947/6129031‑6129319                                                                                                                                                   |                                                                                                                                                ggAAGCTGTGACGTGATGCCTGCGCAAATGACCTATGAAGACTCGACAGGCAAGCTCAACATTCTCGAATACCGCGCAGCTGGCA‑TCAATTGCCGCGGCCAAAATTGAGCTCTGAATGCTGTCCAAACGAAAACGGCGCCTTAAAAgcg                                                                                                                                              >  2:52044/1‑149 (MQ=255)       tgtgACGTGATGCCTGCGCAAATGACCTATGAAGACTCGACAGGCAAGCTCAACATTCTCGAATACCGCGCAGCTGGCA‑TCAATTGCCGCGGCCAAAATTGAGCTCTGAATGCTGTCCAAACGCAAACGGCGCTTTTAAGGCGCCGttt                                                                                                                                        >  8:119123/1‑149 (MQ=255)           aCGTGATGCCTGCGCAAATGACCTATGAAGACTCGACAGGCAAGCTCAACATTCTCGAATACCGCGCAGCTGGCA‑TCAATTGCCGCGGCCAAAATTGAGCTCTGAATGCTGTCCAAACGCAAACGGCGCTTTTAAGGCGCCGTTTTCGt                                                                                                                                    <  1:86634/149‑1 (MQ=39)            cGTGATGCCTGCGCAAATGACCTATGAAGACTCGACAGGCAAGCTCAACATTCTCGAATACCGCGCAGCTGGCA‑TCAATTGCCGCGGCCAAAATTGAGCTCTGAATGCTGTCCAAACGCAAACGGCGCTTTTAAGGCGCCGTTTTCGtt                                                                                                                                   <  7:119123/149‑1 (MQ=39)              tGATGCCTGCGCAAATGACCTATGAAGACTCGACAGGCAAGCTCAACATTCTCGAATACCGCGCAGCTGGCA‑TCAATTGCCGCGGCCAAAATTGAGCTCTGAATGCTGTCCAAACGAAAACGGCGCCTTAAAAGCGCCGTTTGCGTTgg                                                                                                                                 <  8:281342/149‑1 (MQ=255)                  gCCTGCGCAAATGACCTATGAAGACTCGACAGGCAAGCTCAACATTCTCGAATACCGCGCAGCTGGCA‑TCAATTGCCGCGGCCAAAATTGAGCTCTGAATGCTGTCCAAACGAAAACGGCGCCTTAAAAGCGCCGttt                                                                                                                                        <  4:306467/138‑1 (MQ=255)                  gCCTGCGCAAATGACCTATGAAGACTCGACAGGCAAGCTCAACATTCTCGAATACCGCGCAGCTGGCA‑TCAATTGCCGCGGCCAAAATTGAGCTCTGAATGCTGTCCAAACGAAAACGGCGCCTTAAAAGCGCCGttt                                                                                                                                        >  3:306467/1‑138 (MQ=255)                     tGCGCAAATGACCTATGAAGACTCGACAGGCAAGCTCAACATTCTCGAATACCGCGCAGCTGGCA‑TCAATTGCCGCGGCCAAAATTGAGCTCTGAATGCTGTCCAAACGAAAACGGCGCCTTAAAAGCGCCGTTTGCGTTGGGAAAACg                                                                                                                          >  8:104899/1‑149 (MQ=255)                        gcAAATGACCTATGAAGACTCGACAGGCAAGCTCAACATTCTCGAATACCGCGCAGCTGGCA‑TCAATTGCCGCGGCCAAAATTGAGCTCTGAATGCTGTCCAAACGAAAACGGCGCCTTAAAAGCGCCGTTTGCGTTGGGAAAACGGTg                                                                                                                       <  7:249733/149‑1 (MQ=255)                          aaaTGACCTATGAAGACTCGACAGGCAAGCTCAACATTCTCGAATACCGCGCAGCTGGCA‑TCAATTGCCGCGGCCAAAATTGAGCTCTGAATGCTGTCCAAACGAAAACGGCGCCTTAAAAGCGCCGTTTGCGTTGGGAAAACGGTggt                                                                                                                     >  5:154040/1‑148 (MQ=255)                                    tGAAGACTCGACAGGCAAGCTCAACATTCTCGAATACCGCGCAGCTGGCA‑TCAATTGCCGCGGCCAAAATTGAGCTCTGAATGCTGTCCAAACGAAAACGGCGCCTTAAAAGCGCCGTTTGCGTTGGGAAAACGGTGGGTCAGTGGCCt                                                                                                           <  2:77346/149‑1 (MQ=255)                                    tGAAGACTCGACAGGCAAGCTCAACATTCTCGAATACCGCGCAGCTGGCA‑TCAATTGCCGCGGCCAAAATTGAGCTCTGAATGCTGTCCAAACGAAAACGGCGCCTTAAAAGCGCCGTTTGCGTTGGGAAAACGGTGGGTCAGTGGCCt                                                                                                           <  7:305238/149‑1 (MQ=255)                                     gAAGACTCGACAGGCAAGCTCAACATTCTCGAATTCCGCGCAGCTGGCA‑TCAATTGCCGCGGCCAAAATTGAGCTCTGAATGCTGTCCAAACGAAAACGGCGCCTTAAAAGCGCCGTTTGCGTTGGGAAAACGGTGGGTCAGTGGCCTg                                                                                                          >  1:60068/1‑149 (MQ=255)                                     gAAGACTCGACAGGCAAGCTCAACATTCTCGAATACCGCGCAGCTGGCATTCAATTGCCGCGGCCAAAATTGAGCTCTGAATGCTGTCCAAACGAAAACGGCGCCTTAAAAGCGCCGTTTGCGTTGGGAAAACGGTTGGTCAGTGGCCt                                                                                                           >  1:304585/1‑149 (MQ=255)                                      aaGACTCGACAGGCAAGCTCAACATTCTCGAATACCGCGCAGCTGGCA‑TCAATTGCCGCGGCCAAAATTGAGCTCTGAATGCTGTCCAAACGAAAACGGCGCCTTAAAAGCGCCGTTTGCGTTGGGAAAACGGt                                                                                                                        <  1:399116/134‑1 (MQ=255)                                      aaGACTCGACAGGCAAGCTCAACATTCTCGAATACCGCGCAGCTGGCA‑TCAATTGCCGCGGCCAAAATTGAGCTCTGAATGCTGTCCAAACGAAAACGGCGCCTTAAAAGCGCCGTTTGCGTTGGGAAAACGGt                                                                                                                        >  2:399116/1‑134 (MQ=255)                                       aGACTCGACAGGCAAGCTCAACATTCTCGAATACCGCGCAGCTGGCA‑TCAATTGCCGCGGCCAAAATTGAGCTCTGAATGCTGTCCAAACGAAAACGGCGCCTTAAAAGCGCCGTTTGCGTTGGGAAAACGGTGGGTCAGTGGCCTGGa                                                                                                        <  2:235748/149‑1 (MQ=255)                                         aCTCGACAGGCAAGCTCAACATTCTCGAATACCGCGCAGCTGGCA‑TCAATTGCCGCGGCCAAAATTGAGCTCTGAATGCTGTCCAAACGAAAACGGCGCCTTAAAAGCGCCGTTTGCGTTGGGAAAACGGTGGGTCAGTGGCCTGGATa                                                                                                      >  1:281595/1‑149 (MQ=255)                                             gACAGGCAAGCTCAACATTCTCGAATACCGCGCAGCTGGCA‑TCAATTGCCGCGGCCAAAATTGAGCTCTGAATGCTGTCCAAACGAAAACGGCGCCTTAAAAGCGCCGTTTGCGTTGGGAAAACGGTGGGTCAGTGGCCTGGATActct                                                                                                  >  8:153763/1‑149 (MQ=255)                                                aGGCAAGCTCAACATTCTCGAATACCGCGCAGCTGGCA‑TCAATTGCCGCGGCCAAAATTGAGCTCTGAATGCTGTCCAAACGAAAACGGCGCCTTAAAAGCGCCGTTTGCGTTGGGAAAACGGTGGGTCAGTGGCCTGGATACTCTGc                                                                                                >  5:51669/1‑148 (MQ=255)                                                  gCAAGCTCAACATTCTCGAATACCGCGCAGCTGGCA‑TCAATTGCCGCGGCCAAAATTGAGCTCTGAATGCTGTCCAAACGCAAACGGCGCTTTTAAGGCGcc                                                                                                                                            <  5:415385/102‑1 (MQ=39)                                                  gCAAGCTCAACATTCTCGAATACCGCGCAGCTGGCA‑TCAATTGCCGCGGCCAAAATTGAGCTCTGAATGCTGTCCAAACGCAAACGGCGCTTTTAAGGCGcc                                                                                                                                            >  6:415385/1‑102 (MQ=39)                                                  gCAAGCTCAACATTCTCGAATACCGCGCAGCTGGCA‑TCAATTGCCGCGGCCAAAATTGAGCTCTGAATGCTGTCCAAACGAAAACGGCGCCTTAAAAGCGCCGTTTGCGTTGGGAAAACGGTGGGTCAGTGGCCTGGATACTCTGCAAc                                                                                             <  3:63374/149‑1 (MQ=255)                                                  gCAAGCTCAACATTCTCGAATACCGCGCAGCTGGCA‑TCAATTGCCGCGGCCAAAATTGAGCTCTGAATGCTGTCCAAACGAAAACGGCGCCTTAAAAGCGCCGTTTGCGTTGGGAAAACGGTGGGTCAGTGGCCTGGATACTCTGCAAc                                                                                             >  4:292064/1‑149 (MQ=255)                                                          aaCATTCTCGAATACCGCGCAGCTGGCA‑TCAATTGCCGCGGCCAAAATTGAGCTCTGAATGCTGTCCAAACGAAAACGGCGCCTTAAAAGCGCCGTTTGCGTTGGGAAAACGGTGGGTCAGTGGCCTGGATACTCTGCAACCACTTc                                                                                       <  2:206140/147‑1 (MQ=255)                                                            cATTCTCGAATACCGCGCAGCTGGCA‑TCAATTGCCGCGGCCAAAATTGAGCTCTGAATGCTGTCCAAACGAAAACGGCGCCTTAAAAGCGCCGTTTGCGTTGGGAAAACGGTGGGTCAGTGGCCTGGATACTCTGCaa                                                                                              >  5:393741/1‑138 (MQ=255)                                                            cATTCTCGAATACCGCGCAGCTGGCA‑TCAATTGCCGCGGCCAAAATTGAGCTCTGAATGCTGTCCAAACGAAAACGGCGCCTTAAAAGCGCCGTTTGCGTTGGGAAAACGGTGGGTCAGTGGCCTGGATACTCTGCaa                                                                                              <  6:393741/138‑1 (MQ=255)                                                               tctcGAATACCGCGCAGCTGGCA‑TCAATTGCCGCGGCCAAAATTGAGCTCTGAATGCTGTCCAAACGAAAACGGCGCCTTAAAAGCGCCGTTTGCGTTGGGAAAACGGTGGGTCAGTGGCCTGGATACTCTGCAACCACTTCATCCCGt                                                                                >  8:262691/1‑149 (MQ=255)                                                                 tcGAATACCGCGCAGCTGGCA‑TCAATTGCCGCGGCCAAAATTGAGCTCTGAATGCTGTCCAAACGAAAACGGCGCCTTAAAAGCGCCGTTTGCGTTGGGAAAACGGTGGGTCAGTGGCCTGGATACTCTGCAACCACTTCATCCCGt                                                                                >  6:263171/1‑147 (MQ=255)                                                                 tcGAATACCGCGCAGCTGGCA‑TCAATTGCCGCGGCCAAAATTGAGCTCTGAATGCTGTCCAAACGAAAACGGCGCCTTAAAAGCGCCGTTTGCGTTGGGAAAACGGTGGGTCAGTGGCCTGGATACTCTGCAACCACTTCATCCCGTcc                                                                              >  6:267838/1‑149 (MQ=255)                                                                        ccGCGCAGCTGGCA‑TCAATTGCCGCGGCCAAAATTGAGCTCTGAATGCTGTCCAAACGAAAACGGCGCCTTAAAAGCGCCGTTTGCGTTGGGAAAACGGTGGGTCAGTGGCCTGGATACTCTGCAACCACTTCATCCCGTCCGCTCtt                                                                        <  5:46914/148‑1 (MQ=255)                                                                         cgcgCAGCTGGCA‑TCAATTGCCGCGGCCAAAATTGAGCTCTGAATGCTGTCCAAACGAAAACGGCGCCTTAAAAGCGCCGTTTGCGTTGGGAAAACGGTGGGTCAGTGGCCTGGATACTCTGCAACCACTTCATCCCGTCCGCTCttgt                                                                      >  1:390529/1‑149 (MQ=255)                                                                          gcgcAGCTGGCA‑TCAATTGCCGCGGCCAAAATTGAGCTCTGAATGCTGTCCAAACGAAAACGGCGCCTTAAAAGCGCCGTTTGCGTTGGGAAAACGGTGGGTCAGTGGCCTGGATACTCTGCAACCACTTCATCCCGTCCGCTCttgtt                                                                     <  7:153763/149‑1 (MQ=255)                                                                          gcgcAGCTGGCA‑TCAATTGCCGCGGCCAAAATTGAGCTCTGAATGCTGTCCAAACGAAAACGGCGCCTTAAAAGCGCCGTTTGCGTTGGGAAAACGGTGGGTCAGTGGCCTGGATACTCTGCAACCACTTCATCCCGTCCGCTCttgtt                                                                     <  1:52044/149‑1 (MQ=255)                                                                            gcAGCTGGCA‑TCAATTGCCGCGGCCAAAATTGAGCTCTGAATGCTGTCCAAACGAAAACGGCGCCTTAAAAGCGCCGTTTGCGTTGGGAAAACGGTGGGTCAGTGGCCTGGATACTCTGCAACCACTTCATCCCGTCCGCTCTTGTTga                                                                   <  3:104126/149‑1 (MQ=255)                                                                                cTGGCA‑TCAATTGCCGCGGCCAAAATTGAGCTCTGAATGCTGTCCAAACGAAAACGGCGCCTTAAAAGCGCCGTTTGCGTTGGGAAAACGGTGGGTCAGTGGCCTGGATACTCTGCAACCACTTCATCCCGTCCGCTCTTGTTGAGAcc                                                               <  5:263171/149‑1 (MQ=255)                                                                                   gCA‑TCAATTGCCGCGGCAAAAATTGAGCTCTGAATGCTGTCCAAACGAAAACGGCGCCTTAAAAGCGCCGTTTGCGTTGGGAAAACGGTGGGTCAGTGGCCTGGATACTCTGCAACCACTTCATCCCGTCCGCTCTTGTTGAGACCGAt                                                            <  2:345089/149‑1 (MQ=255)                                                                                       tCAATTGCCGCGGCCAAAATTGAGCTCTGAATGCTGTCCAAACGAAAACGGCGCCTTAAAAGCGCCGTTTGCGTTGGGAAAACGGTGGGTCAGTGGCCTGGATACTCTGCAACCACTTCATCCCGTCCGCTCTTGTTGAGACCGATGAc                                                         >  5:168824/1‑149 (MQ=255)                                                                                         aaTTGCCGCGGCCAAAATTGAGCTCTGAATGCTGTCCAAACGAAAACGGCGCCTTAAAAGCGCCGTTTGCGTTGGGAAAACGGTGGGTCAGTGGCCTGGATACTCTGCAACCACTTCATCCCGTCCGCTCTTGTTGAGAACGATGACCt                                                       >  3:109441/1‑149 (MQ=255)                                                                                              ccGCGGCCAAAATTGAGCTCTGAATGCTGTCCAAACGAAAACGGCGCCTTAAAAGCGCCGTTTGCGTTGGGAAAACGGTGGGTCAGTGGCCTGGATACTCTGCAACCACTTCATCCCGTCCGCTCTTGTTGAGACCGATGACCTGGTAt                                                  <  7:104899/149‑1 (MQ=255)                                                                                                    ccAAAATTGAGCTCTGAATGCTGTCCAAACGAAAACGGCGCCTTAAAAGCGCCGTTTGCGTTGGGAAAACGGTGGGTCAGTGGCCTGGATACTCTGCAACCACTTCATCCCGTCCGCTCTTGTTGAGACCGATGACCTGGTATGCATga                                            <  2:304585/149‑1 (MQ=255)                                                                                                      aaaaTTGAGCTCTGAATGCTGTCCAAACGAAAACGGCGCCTTAAAAGCGCCGTTTGCGTTGGGAAAACGGTGGGTCAGTGGCCTGGATACTCTGCAACCACTTCATCCCGTCCGCTCTTGTTGAGACCGATGACCTGGTATGCATgaga                                          <  2:281595/149‑1 (MQ=255)                                                                                                         aTTGAGCTCTGAATGCTGTCCAAACGAAAACGGCGCCTTAAAAGCGCCGTTTGCGTTGGGAAAACGGTGGGTCAGTGGCCTGGATACTCTGCAACCACTTCATCCCGTCCGCTCTTGTTGAGACCGATGACCTGGTATGCATGAGAGac                                       <  5:267838/149‑1 (MQ=255)                                                                                                                tctGAATGCTGTCCAAACGAAAACGGCGCCTTAAAAGCGCCGTTTGCGTTGGGAAAACGGTGGGTCAGTGGCCTGGATACTCTGCAACCACTTCATCCCGTCCGCTCTTGTTGAGACCGATGAc                                                         >  4:101951/1‑124 (MQ=255)                                                                                                                tctGAATGCTGTCCAAACGAAAACGGCGCCTTAAAAGCGCCGTTTGCGTTGGGAAAACGGTGGGTCAGTGGCCTGGATACTCTGCAACCACTTCATCCCGTCCGCTCTTGTTGAGACCGATGACCTGGTATGCATGAGAGACACCATc                                 >  8:57544/1‑148 (MQ=255)                                                                                                                              aaaCGAAAACGGCGCCTTAAAAGCGCCGTTTGCGTTGGGAAAACGGTGGGTCAGTGGCCTGGATACTCTGCAACCACTTCATCCCGTCCGCTCTTGTTGAGACCGATGACCTGGTATGCATgaga                                          >  1:249486/1‑125 (MQ=255)                                                                                                                              aaaCGAAAACGGCGCCTTAAAAGCGCCGTTTGCGTTGGGAAAACGGTGGGTCAGTGGCCTGGATACTCTGCAACCACTTCATCCCGTCCGCTCTTGTTGAGACCGATGACCTGGTATGCATgaga                                          <  2:249486/125‑1 (MQ=255)                                                                                                                              aaaCGAAAACGGCGCCTTAAAAGCGCCGTTTGCGTTGGGAAAACGGTGGGTCAGTGGCCTGGATACTCTGCAACCACTTCATCCCGTCCGCTCTTGTTGAGACCGATGACCTGGTATGCATGAGAGACAcc                                    <  1:390756/131‑1 (MQ=255)                                                                                                                              aaaCGAAAACGGCGCCTTAAAAGCGCCGTTTGCGTTGGGAAAACGGTGGGTCAGTGGCCTGGATACTCTGCAACCACTTCATCCCGTCCGCTCTTGTTGAGACCGATGACCTGGTATGCATGAGAGACAcc                                    >  2:390756/1‑131 (MQ=255)                                                                                                                               aaCGAAAACGGCGCCTTAAAAGCGCCGTTTGCGTTGGGAAAACGGTGGGTCAGTGGCCTGGATACTCTGCAACCACTTCATCCCGTCCGCTCTTGTTGAGACCGATGACCTGGTATGCATGAGAGACACCATCCACTTCCATACCTGga                 <  3:292064/149‑1 (MQ=255)                                                                                                                                aCGAAAACGGCGCCTTAAAAGCGCCGTTTGCGTTGGGAAAACGGTGGGTCAGTGGCCTGGATACTCTGCAACCACTTCATCCCGTCCGCTCTTGTTGAGACCGATGACCTGGTATGCATGAGAGACACCATCCACTTCCATAc                      <  1:146446/143‑1 (MQ=255)                                                                                                                                aCGAAAACGGCGCCTTAAAAGCGCCGTTTGCGTTGGGAAAACGGTGGGTCAGTGGCCTGGATACTCTGCAACCACTTCATCCCGTCCGCTCTTGTTGAGACCGATGACCTGGTATGCATGAGAGACACCATCCACTTCCATAc                      >  2:146446/1‑143 (MQ=255)                                                                                                                                 cGAAAACGGCGCCTTAAAAGCGCCGTTTGCGTTGGGAAAACGGTGGGTCAGTGGCCTGGATACTCTGCAACCACTTCATCCCGTCCGCTCTTGTTGAGACCGATGACCTGGTATGCATGAGAGACACCATCCACTTCCATACCTGgaga               <  2:60068/149‑1 (MQ=255)                                                                                                                                  gAAAACGGCGCCTTAAATGCGCAGTTTGCGTTGGGAAAACGGTTGGTCAGTGGCCTGGATACTCTGCAACCACTTCATCCCGTCCGCTCttgttg                                                                    <  6:272594/95‑1 (MQ=255)                                                                                                                                  gAAAACGGCGCCTTAAAAGCGCCGTTTGCGTTGGGAAAACGGTGGGTCAGTGGCCTGGATACTCTGCAACCACTTCATCCCGTCCGCTCTTGTTGAGACCGATGACCTGGTATGCATGAGAGACACCATCCACTTCCATACCTGGagag              >  6:273739/1‑149 (MQ=255)                                                                                                                                  gAAAACGGCGCCTTAAAAGCGCCGTTTGCGTTGGGAAAACGGTGGGTCAGTGGCCTGGATACTCTGCAACCACTTCATCACGTCCGCTCttgttg                                                                    >  5:272594/1‑95 (MQ=255)                                                                                                                                   caaaCGGCGCTTTTAAGGCGCCGTTTTCGTTTGGAAAACGGTGGGTCAGTGGCCTGGATACTCTGCAACCACTTCATCCCGTCCGCTCTTGTTGAGACCGATGACCTGGTATGCATGAGAGACACCATCCACTTCCATACCTGGAGAGc             >  7:219813/2‑149 (MQ=39)                                                                                                                                     aaCGGCGCCTTAAAAGCGCCGTTTGCGTTGGGAAAACGGTGGGTCAGTGGCCTGGATACTCTGCAACCACTTCATCCCGTCCGCTCTTGTTGAGACCGATGACCTGGTATGCATGAGAGACACCATCCACTTCCATACCTGGAGAGCCt           <  7:262691/149‑1 (MQ=255)                                                                                                                                         gcgcCTTAAAAGCGCCGTTTGCGTTGGGAAAACGGTGGGTCAGTGGCCTGGATACTCTGCAACCACTTCATCCCGTCCGCTCTTGTTGAGACCGATGACCTGGTATGCATGAGAGACACCATCCACTTCCATACCTGGAGAGCCTGCCg       <  4:109441/149‑1 (MQ=255)                                                                                                                                              ttAAAAGCGCCGTTTGCGTTGGGAAAACGGTGGGTCAGTGGCCTGGATACTCTGCAACCACTTCATCCCGTCCGCTCTTGTTGAGACCGATGACCTGGTATGCATGAGAGACACCATCCACTTCCATACCTGGAGAGCCTGCCGGCATa  >  4:135303/1‑149 (MQ=255)                                                                                                                                               tAAAAGCGCCGTTTGCGTTGGGAAAACGGTGGGTCAGTGGCCTGGATACTCTGCAACCACTTCATCCCGTCCGCTCTTGTTGAGACCGATGACCTGGTATGCATgaga                                          >  5:132618/1‑108 (MQ=255)                                                                                                                                               tAAAAGCGCCGTTTGCGTTGGGAAAACGGTGGGTCAGTGGCCTGGATACTCTGCAACCACTTCATCCCGTCCGCTCTTGTTGAGACCGATGACCTGGTATGCATgaga                                          <  6:132618/108‑1 (MQ=255)                                                                                                                                                   |                                                                                                                                                GGAAGCTGTGACGTGATGCCTGCGCAAATGACCTATGAAGACTCGACAGGCAAGCTCAACATTCTCGAATACCGCGCAGCTGGCA‑TCAATTGCCGCGGCCAAAATTGAGCTCTGAATGCTGTCCAAACGAAAACGGCGCCTTAAAAGCGCCGTTTGCGTTGGGAAAACGGTGGGTCAGTGGCCTGGATACTCTGCAACCACTTCATCCCGTCCGCTCTTGTTGAGACCGATGACCTGGTATGCATGAGAGACACCATCCACTTCCATACCTGGAGAGCCTGCCGGCATA  >  NC_002947/6129031‑6129319

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