Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,178,185	Δ1 bp	6.7%	coding (421/1437 nt)	PP_1931 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,178,185	0	G	.	6.7%	181.0 / 7.2	60	coding (421/1437 nt)	PP_1931	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (29/26);  new base . (2/2);  total (32/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 = 1.00e+00

GAGTTGATGTGGCTCACAGGACGTCTGGCTCCGGATTTCAAAACCATCGCCGACTTTCGCAAGGACAACGGTAAAGCCATTCGTAGCGTCTGCCGTCAGTTCGTGATCCTCTGCCGTAACCTCAACCTTTTCTCTGAGTCGATCATTGCCATCGACGGCAGCAAGTTCA‑‑‑AAGCAGTGAACAATCGCGACCGTAACTTCACCCAGGCCAAGGTGAAGGCCCGCATGCAGCAAATCGAGCAGAGCATTGAGCGCTACCTTACAGCGATGGATTCGGCAGACCGGGCAATGCCG  >  NC_002947/2178038‑2178328                                                                                                                                                    |                                                                                                                                                   gAGTTGTTGTGGCTCACTGGACGTCTGTCTCCGGATTTCAAAACCATCGCCGACTTTCTCAAGGACAACGGTAAATCCATTCGTAGCGTCTGCCGTCAGTTCGTGATCCTCTGCCGTAACCTCAACCTTTTCTCTGAGTCGATCATTGc                                                                                                                                                   <  6:117980/149‑1 (MQ=255)   gTTGATGTGGCTCACAGGACGTCTGGCTCCGGATTTCAAAACCATCGCCGACTTTCGCAAGGACAACGGTAAAGCCATTCGTAGCGTCTGCCGTCAGTTCGTGATCCTCTGCCGTAACCTCAACCTTTTCTCTGAGTCGATCATTGCCa                                                                                                                                                 <  7:161699/149‑1 (MQ=255)   gTTGATGTGGCTCACAGGACGTCTGGCTCCGGATTTCAAAACCATCGCCGACTTTCGCAAGGACAACGGTAAAGCCATTCGTAGCGTCAGCCGTCAGTTCGTGATCCTCTGCCGTAACCTCAACCTTTTCTCTGAGTCGATCATTGCCa                                                                                                                                                 <  3:137969/149‑1 (MQ=255)         gtgGCTCACAGGACGTCTGGCTCCGGATTTCAAAACCATCGCCGACTTTCGCAAGGACAACGGTAAAGCCATTCGTAGCGTCTGCCGTCAGTTCGTGATCCTCTGCCGTAACCTCAACCTTTTCTCTGAGTCGATCATTGCCATCGACg                                                                                                                                           >  5:50453/1‑149 (MQ=255)            gCTCACAGGACGTCTGGCTCCGGATTTCAAAACCATCGCCGACTTTCGCAAGGACAACGGTAAAGCCATTCGTAGCGTCTGCCGTCAGTTCGTGATCCTCTGCCGTAACCTCAACCTTTTCTCTGAGTCGATCATTGCCATCGACGgca                                                                                                                                        <  7:343801/149‑1 (MQ=255)              tCACAGGACGTCTGGCTCCGGATTTCAAAACCATCGCCGACTTTCGCAAGGACAACGGTAAAGCCATTCGTAGCGTCTGCCGTCAGTTCGTGATCCTCTGCCGTAACCTCAACCTTTTCTCTGAGTCGATCATTGCCATCGACGgcag                                                                                                                                       >  1:261927/1‑148 (MQ=255)                acaGGACGTCTGGCTCCGGATTTCAAAACCATCGCCGACTTTCGCAAGGACAACGGTAAAGCCATTCGTAGCGTCTGCCGTCAGTTCGTGATCCTCTGCCGTAACCTCAACCTTTTCTCTGAGTCGATCATTGCCATCGACGgcagc                                                                                                                                      <  7:282110/147‑1 (MQ=255)                    gACGTCTGGCTCCGGATTTCAAAACCATCGCCGACTTTCGCAAGGACAACGGTAAAGCCATTCGTAGCGTCTGCCGTCAGTTCGTGATCCTCTGCCGTAACCTCAACCTTTTCTCTGAGTCGATCATTGCCATCGACGGCAGCAAGtt                                                                                                                                 >  1:317713/1‑148 (MQ=255)                    gACGTCTGGCTCCGGATTTCAAAACCATCGCCGACTTTCGCAAGGACAACGGTAAAGCCATTCGTAGCGTCTGCCGTCAGTTCGTGATCCTCTGCCGTAACCTCAACCTTTTCTCTGAGTCGATCATTGCCATCGACGGCAGCAAGtt                                                                                                                                 >  6:34990/1‑148 (MQ=255)                      cGTGTGGCTCCGGATTTCAAAACCATCGCCGACTTTCGCAAGGACAACGGTAAAGCCATTCGTAGCGTCTGCCGTCAGTTCGTGATCCTCTGCCGTAACCTCAACCTTTTCTCTGAGTCGATCATTGCCATCGACGGCAGCAAGTTCa‑‑‑a                                                                                                                           <  4:326815/149‑1 (MQ=255)                       gTCTGGCTCCGGATTTCAAAACCATCGCCGACTTTCGCAAGGACAACGGTAAAGCCATTCGTAGCGTCTGCCGTCAGTTCGTGATCCTCTGCCGTAACATCAACCTTTTCTCTGAGTCGATCATTGCCATCGACGGCAGCAAGTTCa‑‑‑aa                                                                                                                          >  7:88941/1‑149 (MQ=255)                          tGGCTCCGGATTTCAAAACCATCGCCGACTTTCGCAAGGACAACGGTAAAGCCATTCGTAGCGTCTGCCGTCAGTTCGTGATCCTCTGCCGTAACCTCAACCTTTTCTCTGAGTCGATCATTGCCATCGACGGCAGCAAGTTCA‑‑‑AAGCa                                                                                                                       >  3:324612/1‑149 (MQ=255)                            gCTCCGGATTTCAAAACCATCGCCGACTTTCGCAAGGACAACGGTAAAGCCATTCGTAGCGTCTGCCGTCAGTTCGTGATCCTCTGCCGTAACCTCAACCTTTTCTCTGAGTCGATCATTGCCATCGACGGCAGCAAGTTCA‑‑‑AAGCa                                                                                                                       <  5:187952/147‑1 (MQ=255)                              tCCGGATTTCAAAACCATCGCCGACTTTCTCAAGGACAACGGTAAAGCCATTCGTAGCGACAGCCGTCAGTTCGTGATCCTCTGCCGTAACCTCAACCTTTTATCTGAGTCGATCATTGCCATCGACGGAAGCAAGTACA‑‑‑AAGCAGTg                                                                                                                    >  8:326440/1‑148 (MQ=255)                              tCCGGATTTCAAAACCATCGCCGACTTTCGCAAGGACAACGGTAAAGCCATTCGTAGCGTCTGCCGTCAGTTCGTGATCCTCTGCCGTAACCTCAACCTTTTCTCTGAGTCGATCATTGCCATCGACGGCAGCAAGTTCA‑‑‑AAGCAGTg                                                                                                                    >  4:368733/1‑148 (MQ=255)                                      tCAAAACCATCGCCGACTTTCGCAAGGACAACGGTAAAGCCATTCGTAGCGTCTGCCGTCAGTTCGTGATCCTCTGCCGTAACCTCAACCTTTTCTCTGAGTCGATCATTGCCATCGACGGCAGCAAGTTCA‑‑‑AAGCAGTGAACAATcgc                                                                                                           >  8:31036/1‑149 (MQ=255)                                        aaaaCCATCGCCGACTTTCGCAAGGACAACGGTAAAGCCATTCGTAGCGTCTGCCGTCAGTTCGTGATCCTCTGCCGTAACCTCAACCTTTTCTCTGAGTCGATCATTGCCATCGACGGCAGCAAGTTCA‑‑‑AAGCAGTGAACAATcgcg                                                                                                          <  8:88941/148‑1 (MQ=255)                                        aaaaCCATCGCCGACTTTCGCAAGGACAACGGTAAAGCCATTCGTAGCGTCTGCCGTCAGTTCGTGATCCTCTGCCGTAACCTCAACCTTTTCTCTGAGTCGATCATTGCCATCGACGGCAGCAAGTTCA‑‑‑AAGCAGTGAACAATCGCGa                                                                                                         <  6:50453/149‑1 (MQ=255)                                        aaaaCCATCGCCGACTTTCGCAAGGACAACGGTAAAGCCATTCGTAGCGTCTGCCGTCAGTTCGTGATCCTCTGCCGTAACCTCAACCTTTTCTCTGAGTCGATCATTGCCATCGACGGCAGCAAGTTCA‑‑‑AAGCAGTGAACAATCGCGa                                                                                                         <  2:59920/149‑1 (MQ=255)                                        aaaaCCATCGCCGACTTTCGCAAGGACAACGGTAAAGCCATTCGTAGCGTCTGCCGTCAGTTCGTGATCCTCTGCCGTAACCTCAACCTTTTCTCTGAGTCGATCATTGCCATCGACGGCAGCAAGTTCA‑‑‑AAGCAGTGAACAATCGCGa                                                                                                         <  4:324612/149‑1 (MQ=255)                                         aaaCCATCGCCGACTTTCGCAAGGACAACGGTAAAGCCATTCGTAGCGACTGCCGTCAGTTCGTGATCCTCTGCCGTAACCTCAACCTTTTATATGAGTCGATAATTGCCATCCACGGCAGCAAGTTCA‑‑‑AAGCAGTGAACAATCGCGAc                                                                                                        >  1:147427/1‑149 (MQ=255)                                               tCGCCGACTTTCGCAAGGACAACGGTAAAGCCATTCGTAGCGTCTGCCGTCAGTTCGTGATCCTCTGCCGTAACCTCAACCTTTTCTCTGAGTCGATCATTGCCATCGACGGCAGCAAGTTCA‑‑‑AAGCAGTGAACAATCGCGACCGTAAc                                                                                                  >  1:7161/1‑149 (MQ=255)                                                      cTTTCGCAAGGACAACGGTAAAGCCATTCGTAGCGTCTGCCGTCAGTTCGTGATCCTCTGCCGTAACCTCACCCTTTTCTCTGAGTCGATCATT‑CCATCGACGGCAGCAAGTTCA‑‑‑AAGCAGTGAACAATCGCGACCGTAACTTCAcc                                                                                            >  3:54040/1‑147 (MQ=255)                                                      cTTTCGCAAGGACAACGGTAAAGCCATTCGTAGCGTCTGCCGTCAGTTCGTGATCCTCTGCCGTAACCTCAACCTTTTCTCTGAGTCGATCATT‑ACATCGACGGCAGCAAGTTCA‑‑‑AAGCAGTGAACAATAGCGACCGTAACTTCAcc                                                                                            >  3:49100/1‑147 (MQ=255)                                                       tttCGCAAGGACAACGGTAAAGCCATTCGTAGCGTCTGCCGTCAGTTCGTGATCCTCTGCCGTAACCTCAACCTTTTCACTGAGTCGATCATTGCCATCGAGGGCAGCAAGTTCA‑‑‑AAGCAGTGAACAAAAGCGACCGTAACTTCAcc                                                                                            >  8:188099/1‑147 (MQ=255)                                                            cAAGGACAACGGTAAAGCCATTGGTAGCGTCTGCCGTCAGTTCGTGATCCTCTGCCGTAACCTCAACCTTTTCTCTGAGTCGATCATTGCCATCGACGGCAGCAAGTTCA‑‑‑AAGCAGTGAACAATCGCGACCGTAACTTCACCCAGGcc                                                                                      >  7:236866/1‑148 (MQ=255)                                                             aaGGACAACGGTAAAGCCATTCGTAGCGTCTGCCGTCAGTTCGTGATCCTCTGCCGTAACCTCAACCTTTTCTCTGAGTCGATCATTGCCATCGACGGCAGCAAGTTCA‑‑‑AAGCAGTGAACAATCGCGACCGTAACTTCACCCAGGcc                                                                                      >  6:307517/1‑147 (MQ=255)                                                              agtaCAACGGTAAAGCCATTCGTAGCGTCTGCCGTCAGTTCGTGATCCTCTGCCGTAACCTCAACCTTTTCTCTGAGTCGATCATTGCCATCGACGGCAGCAAGTTCA‑‑‑AAGCAGTGAACAATCGCGACCGTAACTTCACCCAGGCCAAg                                                                                   <  8:334999/146‑1 (MQ=255)                                                                   aaCGGTAAAGCCATTCGTAGCGTCTGCCGTCAGTTCGTGATCCTCTGCCGTAACCTCAACCTTTTCTCTGAGTCGATCATTGCCATCGACGGCAGCAAGTTCA‑‑‑AAGCAGTGAACAATCGCGACCGTAACTTCACCCAGGCCAAGGTGaa                                                                              <  2:7161/149‑1 (MQ=255)                                                                        tAAAGCCATTCGTAGCGTCTGCCGTCAGTTCGTGATCCTCTGCCGTAACCTCAACCTTTTCTCTGAGTCGATCATTGCCATCGACGGCAGCAAGTTCA‑‑‑AAGCAGTGAACAATCGCGACCGTAACTTCACCTTGGCCTGGGTGAAgttac                                                                         <  3:368733/149‑5 (MQ=255)                                                                        tAAAGCCATTCGTAGCGTCTGCCGTCAGTTCGTGATCCTCTGCCGTAACCTCAACCTTTTCTCTGAGTCGATCATTGCCATCGACGGCAGCAAGTTCA‑‑‑AAGCAGTGAACAATCGCGACCGTAACTTCACCCAGGCCAAGGTGAAGGccc                                                                         <  4:362095/149‑1 (MQ=255)                                                                           aGCCATTCGTAGCGTCTGCCGTCAGTTCGTGATCCTCTGCCGTAACCTCAACCTTTTCTCTGAGTCGATCATTGCCATCGACGGCAGCAAGTTCA‑‑‑AAGCAGTGAACAATCGCGACCGTAACTTCACCCAGGCCAAGGTGAAGGCCCGCa                                                                      >  1:207010/1‑149 (MQ=255)                                                                            gCCATTCGTAGCGTCTGCCGTCAGTTCGTGATCCTCTGCCGTAACCTCAACCTTTTCTCTGAGTCGATCATTGCCATCGACGGCAGCAAGTTCA‑‑‑AAGCAGTGAACAATCGCGACCGTAACTTCACCCAGGCCAAGATAAAGGCCCGCAt                                                                     >  1:129006/1‑149 (MQ=255)                                                                               aTTCGTAGCGTCTGCCGTCAGTTCGTGATCCTCTGCCGTAACCTTAACCTTTTCTCTGAGTCGATCATTGCCATCGACGGCAGCAAGTTCA‑‑‑AAGCAGTGAACAATCGCGACCGTAACTTCACCCAGGCCAAGGTGAAGGCCCGCATg                                                                    >  4:60787/1‑147 (MQ=255)                                                                                ttCGTATCGTCTGCCGTCAGTTCGTGATCCTCTTCCGTAACCTCAACCTTTTCTCTGAGTCGATCATT‑CCATCGACGGCAGCAAGTTCA‑‑‑AAGCAGTGAACAATCGCGACCGTAACTTCACCCAGGCCAAGGTGAAGGCCCGCATgcagc                                                                <  4:54040/149‑1 (MQ=255)                                                                                ttCGTAGCGTCTGCCGTCAGTTCGTGATCCTCTGCCGTAACCTCAACCTTTTCTCTGAGTCGATCATTGCCATCGACGGCAGCAAGTTCA‑‑‑AAGCAGTGAACAATCGCGACCGTAACTTCACCCAGGCCAAGGTGAAGGCCCGCATgcag                                                                 <  5:307517/149‑1 (MQ=255)                                                                                ttCGTAGCGTCTGCCGTCAGTTCGTGATCCTCTGCCGTAACCTCAACCTTTTCTCTGAGTCGATCATTGCCATCGACGGCAGCAAGTTCA‑‑‑AAGCAGTGAACAATCGCGACCGTAACTTCACCCAGGCCAAGGGGAAGGCCCGCATgcag                                                                 <  2:261927/149‑1 (MQ=255)                                                                                 tCGTAGCGTCTGCCGTCATTTCGTGATCCTCTCCCTTAACATCAACATTTTCTCTGATTCGATCATT‑CCATCGACGGCAGCAAGTTCA‑‑‑AAGCAGTGAACAATCGCGACCGTAACTTCACCCAGGCCAAGGTGAAGGCCCGCATgcagc                                                                <  4:49100/148‑1 (MQ=255)                                                                                             ccGTCAGTTCGTGATCCTCTGCCGTAACCTCAACCTTTTCTCTGAGTCGATCATTGCCATCGACGGCAGCAAGTTCA‑‑‑AAGCAGTGAACAATCGCGACCGTAACTTCACCCAGGCCAAGGTGAAGGCCCGCATGCAGCAAATCgagcaga                                                    <  2:207010/149‑1 (MQ=255)                                                                                                  aGTTCGTGATCCTCTGCCGTATCCTCAACCTTTTCTCTGAGTCGATCATTGCCATCGACGTCAGCAAGTTCA‑‑‑AAGCAGTGAACAATCGCTACCGTATCTTCACCCAGGCCAAGGTGAAGGCCCGCATGCAGCAAATCGAGCAGAGCAtt                                               <  7:326440/149‑1 (MQ=255)                                                                                                   gTTCGTGATCCTCTGCCGTAACCTCAACCTTTTCTCTGAGTCGATCATTGCCATCGACGGCAGCAAGTTCA‑‑‑AAGCAGTGAACAATCGCGACCGTAACTTCACCCAGGCCAAGGTGAAGGCCCGCATGCAGCAAATCGAGCAGAGCAtt                                               >  5:348520/1‑148 (MQ=255)                                                                                                    ttCGTGATCCTCTGCCGTAACCTCAACCTTTTCTCTGAGTCGATCATTGCCATCGACGGCAGCAAGTTCA‑‑‑AAGCAGTGAACAATCGCGACCGTAACTTCACCCAGGCCAAGGTGAAGGCCCGCATGCAGCAAATCGAGCAGAGCAtt                                               >  4:100352/1‑147 (MQ=255)                                                                                                     tCGTGATCCTCTGCCGTAACCTTAACCTTTTCTCTGAGTCGATCATTGCGATCGACGGCAGCAAGTTCA‑‑‑AAGCAGTGAACAATCGCGACCGTAACTTCACCCAGGCCAAGGTGAAGGCCCGCATGCAGCAAATCGAGCAGAGCATTGa                                             <  3:60787/148‑1 (MQ=255)                                                                                                       gTGATCCTCTGCCGTAACCTCAACCTTTTCTCTGAGTCGATCATTGCCATCGACGGCAGCAAGTTCAGAGAAGCAGTGAACAATCGCGACCGTAACTTCACCCAGGCCAAGGTGAAGGCCCGCATGCAGCAAATCGAGCAGAGCAtt                                               >  1:116400/1‑147 (MQ=255)                                                                                                       gTGATCCACTGCCGTAACCTCAACCTTTTCTCTGAGTCGATCATTGCCAGCGACGGCAGCAAGTTCA‑‑‑AAGCATTGAACAATCGCGACCGTAACTTCACCCAGGCCAAGGTGAAGGCCCGCATGCAGCAAATCGAGCAGAGCATTGAgc                                           <  6:20523/148‑1 (MQ=255)                                                                                                        tGATCCTCTGCCGTAACCTCAACCTTTTCTCTGAGTCGATCATTGCCATCGACGGCAGCAAGTTCA‑‑‑AAGCAGTGAACAATCGCGACCGTAACTTCACCCAGGCCAAGGTGAAGGCCCGCATGCAGCAAAt                                                            <  4:310526/130‑1 (MQ=255)                                                                                                        tGATCCTCTGCCGTAACCTCAACCTTTTCTCTGAGTCGATCATTGCCATCGACGGCAGCAAGTTCA‑‑‑AAGCAGTGAACAATCGCGACCGTAACTTCACCCAGGCCAAGGTGAAGGCCCGCATGCAGCAAAt                                                            >  3:310526/1‑130 (MQ=255)                                                                                                                  ccGTAACCTCAACCTTTTCTCTGAGTCGATCATTGCCATCGACGGCAGCAAGTTCA‑‑‑AAGCAGTGAACAATCGCGACCGTAACTTCACCCAGGCCAAGGTGAAGGCCCGCATGCAGCAAATCGAGCAGAGCATTGAGCGCTACCTTACa                                <  2:317713/148‑1 (MQ=255)                                                                                                                         cTCAACCTTTTCTCTGAGTCGATCATTGCCATCGACGGCAGCAAGTTCA‑‑‑AAGCAGTGAACAATCGAGACCGTAACTTCACCCAGGCCAAGGTGAAGGCCCGCATGCAGCAAATCGAGCAGAGCATTGAGCGCTACCTTACAGCGATGGa                        >  5:208371/1‑149 (MQ=255)                                                                                                                              ccTTTTCTCTGAGTCGATCATTGCCATCGACGGCAGCAAGTTCA‑‑‑AAGCAGTGAACAATCGCGACCGTAACTTCACCCAGGCCAAGGTGAAGGCCCGCATGCAGCAAATCGAGCAGAGCATTGAGCGCTACCTTACAGCGATGGATTCgg                   >  3:164529/1‑149 (MQ=255)                                                                                                                               cTTTTCTCTGAGTCGATCATTGCCATCGACGGCAGCAAGTTCA‑‑‑AAGCAGTGAACAATCGCGACCGTAACTTCACCCAGGCCAAGGTGAAGGCCCGCATGCAGCAAATCGAGCAGAGCATTGAGCGCTACCTTACAGCGATGGATTCgg                   >  6:80051/1‑148 (MQ=255)                                                                                                                                ttttCTCTGAGTCGATCATTGCCATCGACGGCAGCAAGTTCA‑‑‑AAGCAGTGAACAATCGCGACCGTAACTTCACCCAGGCCAAGGTGAAGGCCCGCATGCAGCAAATCGAGCAGAGCATTGAGCGCTACCTTa                                  >  6:33223/1‑132 (MQ=255)                                                                                                                                  ttCTCTGAGTCGATCATTGCCATCGACGGCAGCAAGTTCA‑‑‑AAGCAGTGAACAATCGCGACCGTAACTTCACCCAGGCCAAGGTGAAGGCCCGCATGCAGCAAATCGAGCAGAGCATTGAGCGCTACCTTACAGCGATGGATTCGGCa                 <  3:100352/147‑1 (MQ=255)                                                                                                                                     tctGAGTCGATCATTGCCATCGACGGCAGCAAGTTCA‑‑‑AAGCAGTGAACAATCGCGACCGTAACTTCACCCAGGCCAAGGTGAAGGCCCGCATGCAGCAAATCGAGCAGAGCATTGAGCGCTACCTTACAGCGATGGATTCGGCAGACCg            >  3:77018/1‑149 (MQ=255)                                                                                                                                     tctGAGTCGATCATTCCCATCGACTGCAGCAAGTTAA‑‑‑AAGCAATGAAAAATCGCGACCGTAACTTCACCCAGGTCAAGGTGAAGGCCCGCATCCAGCAAATCGACCAGAGCAATGAGCGCAACCTTACAGCGATTGATTCGGCAGACCg            >  2:113793/1‑149 (MQ=255)                                                                                                                                       tGAGTCGATCATTGGCATCGACGGCAGCAAGTTCA‑‑‑AAGCAGTGAACAATCGCGACCGTAACTTCACCCAGGTCATGGTGAAGGCCCGCATGCAGGAAATCGTGCAGAGCATTGAGCGCTACCTTACAGCGATGGATTCGGCAGACCggg          <  1:113793/149‑1 (MQ=255)                                                                                                                                       tGAGTCGATCATTGCCATCGACGGCAGCAAGTTCA‑‑‑AAGCAGTGAACAATCGCGACCGTAACTTCACCCAGGCCAAGGTGAAGGCCCGCATGCAGCAAATCGAGCAGAGCATTGAGCGCTACCTTACAGCGATGGATTCGGCAGACCggg          <  5:34990/149‑1 (MQ=255)                                                                                                                                       tGAGTCGATCATTGCCATCGACGGCAGCAAGTTCA‑‑‑AAGCAGTGAACAATCGCGACCGTAACTTCACCAAGGCCAAGGTGAAGGCCCGAATGCAGCAAATATACCAGACCATTGAGCGCTAACTTACAGCGATGGATTCGGCAGaaagg           >  3:208843/1‑144 (MQ=255)                                                                                                                                        gAGTCGATCATTGCCATCGACGGCAGCAAGTTCA‑‑‑AAGCAGTGAACAATCGCGACCGTAACTTCACCCAGGCCAAGGTGAAGGCCCGCATGCAGCAAATCGAGCAGAGCATTGAGCGCTACCTTACAGCGATGGATTCGGCAGACCggg          >  2:367692/1‑148 (MQ=255)                                                                                                                                                 aTTGCCATCGACGGCAGCAAGTTCA‑‑‑AAGCAGTGAACAATCGCGACCGTAACTTCACCCAGACCAAGGTGAAGGCCCGCATGCAGCAAATCGAGCAGAGCATTGAGCGCTACCTTACAGCGATGGATTCGGCAGACCGGGCAATGCCg  >  1:140177/1‑147 (MQ=255)                                                                                                                                                    |                                                                                                                                                   GAGTTGATGTGGCTCACAGGACGTCTGGCTCCGGATTTCAAAACCATCGCCGACTTTCGCAAGGACAACGGTAAAGCCATTCGTAGCGTCTGCCGTCAGTTCGTGATCCTCTGCCGTAACCTCAACCTTTTCTCTGAGTCGATCATTGCCATCGACGGCAGCAAGTTCA‑‑‑AAGCAGTGAACAATCGCGACCGTAACTTCACCCAGGCCAAGGTGAAGGCCCGCATGCAGCAAATCGAGCAGAGCATTGAGCGCTACCTTACAGCGATGGATTCGGCAGACCGGGCAATGCCG  >  NC_002947/2178038‑2178328

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