Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,050,727	A→C	7.9%	N187H (AAC→CAC)	PP_2663 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,050,727	0	A	C	7.9%	182.3 / 4.1	61	N187H (AAC→CAC)	PP_2663	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (30/26);  new base C (2/3);  total (32/29)
Fisher's exact test for biased strand distribution p-value = 6.62e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCTGCTCGATGGCTTGTCCAGCCGCGAGGAGATGGTGCTTGC‑CGCACTCAGCGCAGCACTGGGGGATATCCCGCACTTCGGCGGCTCGGCCGGTGACGACAATTACCTGAACCGCACCCATGTGTATTACGGCGGTGCATTCCACAACGGCGCAGCCGTGGT‑GGTGCTGGTCAACACCTGGCTCGACTTCGAAGTGTTCACCACGCATCACATCCTGCCGCGCCAGGAAAAGCTGGTGGTAACCGGTGCCGACAGCGCCCAGCGCCGGGTTTTCGAACTGAATGCCGAGC  >  NC_002947/3050582‑3050871                                                                                                                                                   |                                                                                                                                                  gctgctCGATGGCTTGTCCAGCCGCGAGGAGATGGTGCTTGC‑CGCACTCAGCGCAGCACTGGGGGATATCCCGCACTTCGGCGGCTCGGCCGGTGACGACAATTACCTGAACCGCACCCATGTGTATTACGGCGGTGCATTCCACCACg                                                                                                                                                >  8:7296/1‑149 (MQ=255)     ctCGATGGCTTGTCCAGCCGCGAGGAGATGGTGCTTGC‑CGCACTCAGCGCAGCACTGGGGGATATCCCGCACTTCGGCGGCTCGGCCGGTGACGACAATTACCTGAACCGCACCCATGTGTATTACGGCGGTGCATTCCACAACGgcg                                                                                                                                             >  2:110102/1‑148 (MQ=255)       cGATGGCTTGTCCAGCCGCGAGGAGATGGTGCTTGC‑CGCACTCAGCGCAGCACTGGGGGATATCCCGCACTTCGGCGGCTCGGCCGGTGACGACAATTACCTGAACCGCACCCATGTGTATTACGGCGGTGCATTCCACAACGGCGCAg                                                                                                                                          >  3:248752/1‑149 (MQ=255)        gATGGCTTGTCCAGCCGCGAGGAGATGGTGCTTGC‑CGCACTCAGCGCAGCACTGGGGGATATCCCGCACTTCGGCGGCTCGGCCGGTGACGACAATTACCTGAACCGCACCCATGTGTATTACGGCGGTGCATTCCACAACGGCGCAGc                                                                                                                                         >  1:92351/1‑149 (MQ=255)         aTGGCTTGTCCAGCCGCGAGGAGATGGTGCTTGC‑CGCACTCAGCGCAGCACTGGGGGATATCCCGCACTTCGGCGGCTCGGCCGGTGACGACAATTACCTGAACCGCACCCATGTGTATTACGGCGGTGCATTCCACAACGGCGCAGcc                                                                                                                                        >  1:310478/1‑149 (MQ=255)             cTTGTCCAGCCGCGAGGAGATGGTGCTTGC‑CGCACTCAGCGCAGCACTGGGGGATATCCCGCACTTCGGCGGCTCGGCCGGTGACGACAATTACCTGAACCGCACCCATGTGTATTACGGCGGTGCATTCCACAACGGCGCAGCCgtgg                                                                                                                                    <  8:26408/149‑1 (MQ=255)              ttGTCCAGCCGCGAGGAGATGGTGCTTGC‑CGCACTCAGCGCAGCACTGGGGGATATCCCGCACTTCGGCGGCTCGGCCGGTGACGACAATTACCTGAACCGCACCCATGTGTATTACGGCGGTGCATTCCACAACGGCGCAGCCgtggt                                                                                                                                   >  7:137446/1‑149 (MQ=255)                gTCCAGCCGCGAGGAGATGGTGCTTGCGCGCACTCAGCGCAGCACTGGGGGATATCCCGCACTTCGGCGGCTCGGCCGGTGACGACAATTACCTGAACCGCACCCATGTGTATTACGGCGGTGCATTCCACAACGGCGCAGCCgtggt‑g                                                                                                                                 >  1:17116/1‑149 (MQ=255)                gTCCAGCCGCGAGGAGATGGTGCTTGC‑CGCACTCAGCGCAGCACTGGGGGATATCCCGCACTTCGGCGGCTCGGCCGGTGACGACAATTACCTGAACCGCACCCATGTGTATTACGGCGGTGCATTCCACCACGGCTGCGCCgttgt‑g                                                                                                                                 <  7:179857/148‑1 (MQ=255)                   cAGCCGCGAGGAGATGGTGCTTGC‑CGCACTCAGCGCAGCACTGGGGGATATCCCGCACTTCGGCGGCTCGGCCGGTGACGACAATTACCTGAACCGCACCCATGTGTATTACGGCGGTGCATTCCACCACGGCTGCGCCgttgt‑ggaa                                                                                                                              <  1:244006/148‑3 (MQ=255)                    aGCCGCGAGGAGATGGTGCTTGC‑CGCACTCAGCGCAGCACTGGGGGATATCCCGCACTTCGGCGGCTCGGCCGGTGACGACAATTACCTGAACCGCACCCATGTGTATTACGGCGGTGCATTCCACAACGGCGCAGCCGTGGT‑GGTGCt                                                                                                                            >  4:152586/1‑149 (MQ=255)                     gCCGCGAGGAGATGGTGCTTGC‑CGCACTCAGCGCAGCACTGGGGGATATCCCGCACTTCGGCGGCTCGGCCGGTGTCGACAATTACCTGAACCGCACCCATGTGTATTACGGCGGTGCATTCCACAACGGCGCAGCCGTGGT‑GGTGCTg                                                                                                                           <  4:156797/149‑1 (MQ=255)                     gCCGCGAGGAGATGGTGCTTGC‑CGCACTCAGCGCAGCACTGGGGGATATCCCGCACTTCGGCGGCTCGGCCGGTGACGACAATTACCTGAACCGCACCCATGTGTATTACGGCGGTGCATTCCACAACGGCGCAGCCGTGGT‑GGTGCTg                                                                                                                           <  1:110102/149‑1 (MQ=255)                      ccGCGAGGAGATGGTGCTTGC‑CGCACTCAGCGCAGCACTGGGGGATATCCCGCACTTCGGCGGCTCGGCCGGTGACGACAATTACCTGAACCGCACCCATGTGTATTACGGCGGTGCATTCCACAACGGCGCAGCCGTGGT‑GGTGCTgg                                                                                                                          <  5:373678/149‑1 (MQ=255)                        gcgAGGAGATGGTGCTTGC‑CGCGATCAGCGCAGCAGTGGGGGATAACCCGCAGTTCGGCGGCGCGGCCGGTGACGACAATTACCTGAACCGCACCCATGTGTATTACGGCGGTGCATTCCACAACGGCGCAGCCGTGGT‑GGAggtggt                                                                                                                         >  7:240447/1‑148 (MQ=255)                        gcgAGGAGATGGTGCTTGC‑CGCACTCAGCGCAGCACTGGGGGATATCCCGCACTTCGGCGGCTCGGCCGGTGACGACAATTACCTGAACCGCACCCATGTGTATTACGGCGGTGCATTCCACAACGGCGCAGCCGTGGT‑GGTGCTGGt                                                                                                                         >  5:230139/1‑148 (MQ=255)                        gcgAGGAGATGGTGCTTGC‑CGCACTCAGCGCAGCACTGGGGGATATCCCGCACTTCGGCGGCTCGGCCGGTGACGACAATTACCTGAACCGCACCCATGTGTATTACGGCGGTGCATTCCACAACGGCGCAGCCGTGGT‑GGTGCTGGTc                                                                                                                        <  6:289916/149‑1 (MQ=255)                        gcgAGGAGATGGTGCTTGC‑CGCACTCAGCGCAGCACTGGGGGATATCCCGCACTTCGGCGGCTCGGCCGGTGACGACAATTACCTGAACCGCACCCATGTGTATTACGGCGGTGCATTCCACAACGGCGCAGCCGTGGT‑GGTGCTGGTc                                                                                                                        <  4:213547/149‑1 (MQ=255)                          gaggagATGGTGCTTGC‑CGCACTCAGCGCAGCACTGGGGGATATCCCGCACTTCGGCGGCTCGGCCGGTGACGACAATTACCTGAACCGCACCCATGTGTATTACGGCGGTGCATTCCACAACGGCGCAGCCGTGGT‑GGTGCTGGTCaa                                                                                                                      <  6:230139/149‑1 (MQ=255)                                aTGGTGCTTGC‑CGCACTCAGCGCAGCACTGGGGGATATCCCGCACTTCGGCGGCTCGGCCGGTGACGACAATTACCTGAACCGCACCCATGTGTATTACGGCGGTGCATTCCACCACGGCTGCGCCgttgt‑ggaat                                                                                                                             >  2:304971/1‑133 (MQ=255)                                aTGGTGCTTGC‑CGCACTCAGCGCAGCACTGGGGGATATCCCGCACTTCGGCGGCTCGGCCGGTGACGACAATTACCTGAACCGCACCCATGTGTATTACCGCGGTGCATTCCACCACGGCTGCGCCgttgt‑ggaat                                                                                                                             <  1:304971/136‑4 (MQ=255)                                 tGGTGCTTGC‑CGCACTCAGCGCAGCACTGGGGGATATCCCGCACTTCGGCGGCTCGGCCGGTGACGACAATTACCTGAACCGCACCCATGTGTATTACGGCGGTGCATTCCACAACGGCGCAGCCGTGGT‑GGTGCTGGTCAACACCTgg                                                                                                               <  6:389620/149‑1 (MQ=255)                                 tGGTGCTTGC‑CGCACTCAGCGCAGCACTGGGGGATATCCCGCACTTCGGCGGCTCGGCCGGTGACGACAATTACCTGAACCGCACCCATGTGTATTACGGCGGTGCATTCCACAACGGCGCAGCCGTGGT‑GGTGCTGGTCAACACCTgg                                                                                                               >  6:31996/1‑149 (MQ=255)                                  ggTGCTTGC‑CGCACTCAGCGCAGCACTGGGGGATATCCCGCACTTCGGCGGCTCGGCCGGTGACGACAATTACCTGAACCGCACCCATGTGTATTACGGCGGTGCATTCCACAACGGCGCAGCCGTGGT‑GGTGCTGGTCAACACCTGGc                                                                                                              >  3:346778/1‑149 (MQ=255)                                   gTGCTTGC‑CGCACTCAGCGCAGCACTGGGGGATATCCCGCACTTCGGCGGCTCGGCCGGTGACGACAATTACCTGAACCGCACCCATGTGTATTACGGCGGTGCATTCCACAACGGCGCAGCCGTGGT‑GGTGCTGGTCAACACCTGGc                                                                                                              <  2:92351/148‑1 (MQ=255)                                    tGCTTGC‑CGCACTCAGCGCAGCACTGGGGGATATCCCGCACTTCGGCGGCTCGGCCGGTGACGACAATTACCTGAACCGCACCCATGTGTATTACGGCGGTGCATTCCACAACGGCGCAGCCGTGGT‑GGTGCTGGTCAACACCTGGCTc                                                                                                            <  8:137446/149‑1 (MQ=255)                                    tGCTTGC‑CGCACTCAGCGCAGCACTGGGGGATATCCCGCACTTCGGCGGCTCGGCCGGTGACGACAATTACCTGAACCGCACCCATGTGTATTACGGCGGTGCATTCCACAACGGCGCAGCCGTGGT‑GGTGCTGGTCAACACCTGGCTc                                                                                                            <  4:346778/149‑1 (MQ=255)                                      cTTGC‑CGCACTCAGCGCAGCACTGGGGGATATCCCGCACTTCGGCGGCTCGGCCGGTGACGACAATTACCTGAACCGCACCCATGTGTATTACGGCGGTGCATTCCACAACGGCGCAGCCGTGGT‑GGTGCTGGTCAACACCTGGCTc                                                                                                            >  2:378610/1‑147 (MQ=255)                                      cTTGC‑CGCACTCAGCGCAGCACTGGGGGATATCCCGCACTTCGGCGGCTCGGCCGGTGACGACAATTACCTGAACCGCACCCATGTGTATTACGGCGGTGCATTCCACAACGGCGCAGCCGTGGT‑GGTGCTGGTCAACACCTGGCTc                                                                                                            <  1:378610/147‑1 (MQ=255)                                              cACTCAGCGCAGCACTGGGGGATATCCCGCACTTCGGCGGCTCGGCCGGTGACGACAATTACCTGAACCGCACCCATGTGTATTACGGCGGTGCATTCCACAACGTCGCAGCCGTGGT‑GGTGCTGGTCAACACCTGGCTCGACTTCGaa                                                                                                   <  2:17116/149‑1 (MQ=255)                                              cACTCAGCGCAGCACTGGGGGATATCCCGCACTTCGGCGGCTCGGCCGGTGACGACAATTACCTGAACCGCACCCATGTGTATTACGGCGGTGCATTCCACAACGGCGCAGCCGTGGT‑GGTGCTGGTCAACACCTGGCTCGACTTCGaa                                                                                                   <  5:31996/149‑1 (MQ=255)                                              cACTCAGCGCAGCACTGGGGGATATCCCGCACTTCGGCGGCTCGGCCGGTGACGACAATTACCTGAACCGCACCCATGTGTATTACGGCGGTGCATTCCACAACGGCGCAGCCGTGGT‑GGTGCTGGTCAACACCTGGCTCGACTTCGaa                                                                                                   <  2:310478/149‑1 (MQ=255)                                                cTCAGCGCAGCACTGGGGGATATCCCGCACTTCGGCGGCTCGGCCGGTGACGACAATTACCTGAACCGCACCCATGTGTATTACGGCGGTGCATTCCACAACGGCGCAg                                                                                                                                          >  2:340377/1‑109 (MQ=255)                                                cTCAGCGCAGCACTGGGGGATATCCCGCACTTCGGCGGCTCGGCCGGTGACGACAATTACCTGAACCGCACCCATGTGTATTACGGCGGTGCATTCCACAACGGCGCAg                                                                                                                                          <  1:340377/109‑1 (MQ=255)                                                 tCAGCGCAGCACTGGGGGATATCCCGCACTTCGGCGGCTCGGCCGGTGACGACAATTACCTGAACCGCACCCATGTGTATTACGGCGGTGCATTCCACAACGGCGCAGCCGTGGT‑GGTGc                                                                                                                             <  6:210418/120‑1 (MQ=255)                                                 tCAGCGCAGCACTGGGGGATATCCCGCACTTCGGCGGCTCGGCCGGTGACGACAATTACCTGAACCGCACCCATGTGTATTACGGCGGTGCATTCCACAACGGCGCAGCCGTGGT‑GGTGc                                                                                                                             >  5:210418/1‑120 (MQ=255)                                                 tCAGCGCAGCACTGGGGGATATCCCGCACTTCGGCGGCTCGGCCGGTGACGACAATTACCTGAACCGCACCCATGTGTATTACGGCGGTGCATTCCACAACGGCGCAGCCGTGGT‑GGTGCTGGTCAACACCTGGCTCGACTTCGAAgtg                                                                                                >  8:61870/1‑149 (MQ=255)                                                                       cccGCACTTCGGCGGCTCGGCCGGTGACGACAATTACCTGAACCGCACCCATGTGTATTACGGCGGTGCATTCCACAACGGCGCAGCCGTGGT‑GGTGCTGGTCAACACCTGGCTCGACTTCGAAGTGTTCACCACGCATCACATCCTGc                                                                          >  7:361571/1‑149 (MQ=255)                                                                        ccGCACTTCGGCGGCTCGGCCGGTGACGACAATTACCTGAACCGCACCCATGTGTATTACGGCGGTGCATTCCACAACGGCGCAGCCGTGGT‑GGTGCTGGTCAACACCTGGCTCGCCTTCGAAGTGTTCACCACGCATCACATCCTGcc                                                                         >  6:258980/1‑149 (MQ=255)                                                                        ccGCACTTCGGCGGCTCGGCCGGTGACGACAATTACCTGAACCGCACCCATGTGTATTACGGCGGTGCATTCCACAACGGCGCAGCCGTGGT‑GGTGCTGGTCAACACCTGGCTCGACTTCGAAGTGTTCACCACGCATCACATCCTGcc                                                                         >  6:389909/1‑149 (MQ=255)                                                                        ccGCACTTCGGCGGCTCGGCCGGTGACGACAATTACCTGAACCGCACCCATGTGTATTACGGCGGTGCATTCCACAACGGCGCAGCCGTGGT‑GGTGCTGGTCAACACCTGGCTCGACTTCGAAGTGTTCACCACGCATCACATCCTGcc                                                                         >  2:106315/1‑149 (MQ=255)                                                                              ttCGGCGGCTCGGCCGGTGACGACAATTACCTGAACCGCACCCATGTGTATTACGGCGGTGCATTCCACAACGGCGCAGCCGTGGT‑GGTGCTGGTCAACACCTGGCTCGACTTCGAAGTGTTCACCACGCATCACATCCTGCCgcgc                                                                     >  2:200993/1‑147 (MQ=255)                                                                                cggcggCTCGGCCGGTGACGACAATTACCTGAACCGCACCCATGTGTATTACGGCGGTGCATTCCACAACGGCGCAGCCGTGGT‑GGTGCTGGTCAACACCTGGCTCGACTTCGAAGTGTTCACCACGCATCACATCCTGCCGCGCCAg                                                                  <  4:248752/148‑1 (MQ=255)                                                                                  gcggcTCGGCCGGTGACGACAATTACCTGAACCGCACCCATGTGTATTACGGCGGTGCATTCCACAACGGCGCAGCCGTGGT‑GGTGCTGGTCAACACCTGGCTCGACTTCGAAGTGTTCACCACGCATCACATc                                                                              >  4:82650/1‑134 (MQ=255)                                                                                  gcggcTCGGCCGGTGACGACAATTACCTGAACCGCACCCATGTGTATTACGGCGGTGCATTCCACAACGGCGCAGCCGTGGT‑GGTGCTGGTCAACACCTGGCTCGACTTCGAAGTGTTCACCACGCATCACATc                                                                              <  3:82650/134‑1 (MQ=255)                                                                                        cGGCCGGTGACGACAATTACCTGAACCGCACCCATGTGTATTACGGCGGTGCATTCCACAACGGCGCAGCCGTGGT‑GGTGCTGGTCAACACCTGGCTCGACTTCGAAGTGTTCACCACGCATCACATCCTGCCGCGCCAGGAAAAGCtg                                                         <  8:240447/149‑1 (MQ=255)                                                                                          gCCGGTGACGACAATTACCTGAACCGCACCCATGTGTATTACGGCGGTGCATTCCACAACGGCGCAGCCGTGGT‑GGTGCTGGTCAACACCTGGCTCGACTTCGAAGTGTTCACCACGCATCACATCCTGCCGCGCCAGGAAAAGCtggt                                                       >  7:300131/1‑149 (MQ=255)                                                                                                 acgacAATTACCTGAACCGCACCCATGTGTATTACGGCGGTGCATTCCACAACGGCGCAGCCGTGGT‑GGTGCTGGTCAACACCTGGCTCGACTTCGAAGTGTTCACCACGCATCACATCCTGCCGCGCCAGGAAAAGCTGGTGGTAAc                                                 <  8:361571/148‑1 (MQ=255)                                                                                                   gacAATTACCTGAACCGCACCCATGTGTATTACGGCGGTGCATTCCACAACGGCGCAGCCGTGGT‑GGTGCTGGTCAACACCTGGCTCGACTTCGAAGTGTTCACCACGCATCACATCCTGCCGCGCCAGGAAAAGCTGGTGGTAACCgg                                              >  7:210304/1‑149 (MQ=255)                                                                                                       aTTACCTGAACCGCACCCATGTGTATTACGGCGGTGCATTCCACAACGGCGCAGCCGTGGT‑GGTGCTGGTCAACACCTGGCTCGACTTCGAAGTGTTCACCACGCATCACATCCTGCCGCGCCAGGAAAAGCTGGTGGTAACCGGTGcc                                          <  3:152586/149‑1 (MQ=255)                                                                                                           ccTGAACCGCACCCATGTGTATTACGGCGGTGCATTCCACAACGGCGCAGCCGTGGT‑GGTGCTGGTCAACACCTGGCTCGACTTCGAAGTGTTCACCACGCATCACATCCTGCCGCGCCAGGAAAAGCTGGTGGTAACCGGTGCCGa                                        >  3:167953/1‑147 (MQ=255)                                                                                                                aCCGCACCCATGTGTATTACGGCGGTGCATTCCACAACGGCGCAGCCGTGGT‑GGTGCTGGTCAACACCTGGCTCGACTTCGAAGTGTTCACCACGCATCACATCCTGCCGCGCCAGGAAAAGCTGGTGGTAACCGGTGCCGAcagcgcc                                 <  8:300131/149‑1 (MQ=255)                                                                                                                       ccATGTGTATTACGGCGGTGCATTCCACAACGGCGCAGCCGTGGTGGGTGCTGGTCAACACCTGGCTCGACTTCGAAGTGTTCACCACGCATCACATCCTGCCGCGCCAGGAAAAGCTGGTGGTAACCGGTGCCGAcagcgcccagcgc                           >  1:78363/1‑149 (MQ=255)                                                                                                                         aTGTGTATTACGGCGGTGCATTCCACAACGGCGCAGCCGTGGT‑GGTGCTGGTCAACACCTGGCTCGACTTCGAAGTGTTCACCACGCATCACATCCTGCCGCGCCAGGAAAAGCTGGTGGTAACCGGTGCCGACAGCGCCCAGCGCCgg                        >  1:228286/1‑149 (MQ=255)                                                                                                                               aTTACGGCGGTGCATTCCACAACGGCGCAGCCGTGGT‑GGTGCTGGTCAACACCTGGCTCGACTTCGAAGTGTTCACCACGCATCACATCCTGCCGCGCCAGGAAAAGCTGGTGGTAACCGGTGCCGACAGCGCCCAGCGCCGGGtttt                   <  2:228286/148‑1 (MQ=255)                                                                                                                                ttACGGCGGTGCATTCCACAACGGCGCAGCCGTGGT‑GGTGCTGGTCAACACCTGGCTCGACTTCGAAGTGTTCACCACGCATCACATCCTGCCGCGCCAGGAAAAGCTGGTGGTAACCGGTGCCGACAGCGCCCAGCGCCGGGTTTTc                  >  5:188690/1‑148 (MQ=255)                                                                                                                                   cggcggTGCATTCCACAACGGCGCAGCCGTGGT‑GGTGCTGGTCAACACCTGGCTCGACTTCGAAGTGTTCACCACGCATCACATCCTGCCGCGCCAGGAAAAGCTGGTGGTAACCGGTGCCGACAGCGCCCAGCGCCGGGTTTTCGaa               <  6:379058/148‑1 (MQ=255)                                                                                                                                          gCATTCCACAACGGCGCAGCCGTGGT‑GGTGCTGGTCAACACCTGGCTCGACTTCGAAGTGTTCACCACGCATCACATCCTGCCGCGCCAGGAAAAGCTGGTGGTAACCGGTGCCGACAGCGCCCAGCGCCGGGTTTTCGAACTGAATGc       >  2:77013/1‑149 (MQ=255)                                                                                                                                             ttCCACAACGGCGCAGCCGTGGT‑GGTGCTGGTCAACACCTGGCTCGACTTCGAAGTGTTCACCACGCATCACATCCTGCCGCGCCAGGAAAAGCTGGTGGTAACCGGTGCCGACAGCGCCCAGCGCCGGGTTTTCGAACTGAATGcc      >  2:237155/1‑147 (MQ=255)                                                                                                                                             ttCCACAACGGCGCAGCCGTGGT‑GGTGCTGGTCAACACCTGGCTCGACTTCGAAGTGTTCACCACGCATCACATCCTGCCGCGCCAGGAAAAGCTGGTGGTAACCGGTGCCGACAGCGCCCAGCGCCGGGTTTTCGAACTGAATGCCGa    <  5:389909/149‑1 (MQ=255)                                                                                                                                               ccACAACGGCGCAGCCGTGGT‑GGTGCTGGTCAACACCTGGCTCGACTTCGAAGTGTTAACCACGCATCACATCCTGCCGCGCCAGGAAAAGCTGGTGGTAACCGCTGCCGACAGCGCCCAGCGCCGGGTTTTCGAACTGAATGCCGAGc  >  2:327003/1‑149 (MQ=255)                                                                                                                                                cacaACGGCGCAGCCGTGGT‑GGTGCTGGTCAACACCTGGCTCGACTTCGAAGTGTTCACCACGCATcac                                                                                 <  2:358639/69‑1 (MQ=255)                                                                                                                                                cacaACGGCGCAGCCGTGGT‑GGTGCTGGTCAACACCTGGCTCGACTTCGAAGTGTTCACCACGCATcac                                                                                 >  1:358639/1‑69 (MQ=255)                                                                                                                                                   |                                                                                                                                                  GCTGCTCGATGGCTTGTCCAGCCGCGAGGAGATGGTGCTTGC‑CGCACTCAGCGCAGCACTGGGGGATATCCCGCACTTCGGCGGCTCGGCCGGTGACGACAATTACCTGAACCGCACCCATGTGTATTACGGCGGTGCATTCCACAACGGCGCAGCCGTGGT‑GGTGCTGGTCAACACCTGGCTCGACTTCGAAGTGTTCACCACGCATCACATCCTGCCGCGCCAGGAAAAGCTGGTGGTAACCGGTGCCGACAGCGCCCAGCGCCGGGTTTTCGAACTGAATGCCGAGC  >  NC_002947/3050582‑3050871

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