Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,204,754	G→C	7.5%	A20P (GCA→CCA)	PP_2811 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,204,754	0	G	C	7.5%	192.1 / 4.7	64	A20P (GCA→CCA)	PP_2811	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (32/26);  new base C (2/3);  total (35/29)
Fisher's exact test for biased strand distribution p-value = 6.54e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.24e-01

GCAACCCTCGCCTGAAGCGCGACTTCTTCACCTCGGATAACCTGCGAGCCATCAGCGAGCGCTAAAACCGACGTCAGATCGAGAGCCAGCAT‑GAAAAAAATAATTACGAAAGCGTCAGGCCTTGTCATCGTGTGCGCAACCCCGTTGGCAATGGGGTTCACCTTTGAAACCGAAACCCTCAATGGTTCCTTCGACTCCACCATCACCGCAGGCATGGGCCTGCGTACGGAGTCGCGGGGCTGCAACCTGATCAACCAGGGGCCGACCGGGCACAACGTGCCCTCCGGCTGC  >  NC_002947/3204607‑3204897                                                                                                                                                     |                                                                                                                                                gCAACCCTCGCCTGAAGCGCGACTTCTTCACCTCGGATAACCTGCGAGCCATCAGCGAGCGCTAAAACCGACGTCAGATCGAGAGCCAGCAT‑GAAAAAAATAATTACGAAAGCGTCAGGCCTTGTCATCGTGTGCGCAACCCCGTTgg                                                                                                                                                 >  4:188601/1‑148 (MQ=255)  cAACCCTCGCCTGAAGCGCGACTTCTTCACCTCGGATAACCTGCGAGCCATCAGCGAGCGCTAAAACCGACGTCAGATCGAGAGCCAGCAT‑GAAAAAAATAATTACGAAAGCGTCAGGCCTTGTCATCGTGTGCGCAACCCCGTTGGc                                                                                                                                                >  3:44923/1‑148 (MQ=255)      ccTCGCCTGAAGCGCGACTTCTTCACCTCGGATAACCTGCGAGCCATCAGCGAGCGCTAAAACCGACGTCAGATCGAGAGCCAGCAT‑GAAAAAAATAATTACGAAAGCGTCAGGCCTTGTCATTGTGTGCGCAACCCCGTTGGCAATgg                                                                                                                                           <  2:82006/149‑1 (MQ=255)       cTCGCCTGAAGCGCGACTTCTTCACCTCGGATAACCTGCGAGCCATCAGCGAGCGCTAAAACCGACGTCAGATCGAGAGCCAGCAT‑GAAAAAAATAATTACGAAAGCGTCAGGCCTTGTCATCGTGTGCGCAACCCCGTTGGCAATggg                                                                                                                                          <  8:100983/149‑1 (MQ=255)       cTCGCCTGAAGCGCGACTTCTTCACCTCGGATAACCTGCGAGCCATCAGCGAGCGCTAAAACCGACGTCAGATCGAGAGCCAGCAT‑GAAAAAAATAATTACGAAAGCGTCAGGCCTTGTCATCGTGTGCGCAACCCCGTTGGCAATggg                                                                                                                                          >  7:242639/1‑149 (MQ=255)          gCCTGAAGCGCGACTTCTTCACCTCGGATAACCTGCGAGCCATCAGCGAGCGCTAAAACCGACGTCAGATCGAGAGCCAGCAT‑GAAAAAAATAATTACGAAAGCGTCAGGCCTTGTCATCGTGTGCGCAACCCCATTGCCAACGGGGtt                                                                                                                                       <  4:24864/149‑1 (MQ=255)          gCCTGAAGCGCGACTTCTTCACCTCGGATAACCTGCGAGCCATCAGCGAGCGCTAAAACCGACGTCAGATCGAGAGCCAGCAT‑GAAAAAAATAATTACGAAAGCGTCAGGCCTTGTCATCGTGTGCGCAACCCCATTGCCAACGGGGtt                                                                                                                                       <  6:186534/149‑1 (MQ=255)                 gcgcGACTTCTTCACCTCGGATAACCTGCGAGCCATCAGCGAGCGCTAAAACCGACGTCAGATCGAGAGCCAGCAT‑GAAAAAAATAATTACGAAAGCGTCAGGCCTTGTCATCGTGTGCGCAACCCCGTTGGCAATGGGGTTCACCttt                                                                                                                                <  5:281960/149‑1 (MQ=255)                 gcgcGACTTCTTCACCTCGGATAACCTGCCAGCCATCAGCGAGCGCTAAAACCGACGTCAGATCGAGAGCCAGCAT‑GAAAAAAATAATTACGAAAGCGTCAGGCCTTGTCATCGTGTGCGCAACCCCGTTGGCAATGGGGTTCACCttt                                                                                                                                <  7:262061/149‑1 (MQ=255)                   gcgACTTCTTCACCTCGGATAACCTGCGAGCCATCAGCGAGCGCTAAAACCGACGTCAGATCGAGAGCCAGCAT‑GAAAAAAATAATTACGAAAGCGTCAGGCCTTGTCATCGTGTGCGCAACCCCGTTGGCAATGGGGTTCACCTTTg                                                                                                                               >  6:337899/1‑148 (MQ=255)                     gACTTCTTCACCTCGGATAACCTGCGAGCCATCAGCGAGCGCTAAAACCGACGTCAGATCGAGAGCCAGCAT‑GAAAAAAATAATTACGAAAGCGTCAGGCCTTGTCATCGTGTGCGCAACCCCGTTGGCAATGGGGTTCACCTTTgaaa                                                                                                                            >  6:238583/1‑149 (MQ=255)                     gACTTCTTCACCTCGGATAACCTGCGAGCCATCAGCGAGCGCTAAAACCGACGTCAGATCGAGAGCCAGCAT‑GAAAAAAATAATTACGAAAGCGTCAGGCCTTGTCATCGTGTGCGCAACCCCGTTGGCAATGGGGTTCACCTTTgaaa                                                                                                                            <  4:44923/149‑1 (MQ=255)                       cttcttCACCTCGGATAACCTGCGAGCCATCAGCGAGCGCTAAAACCGACGTCAGATCGAGAGCCAGCAT‑GAAAAAAATAATTACGAAAGCGTCAGGCCTTGTCATCGTGTGCGCAACCCCGTTGGCAATGGGGTTCACCTTTgaaacc                                                                                                                          <  1:281239/149‑1 (MQ=255)                       cttcttCACCTCGGATAACCTGCGAGCCATCAGCGAGCGCTAAAACCGACGTCAGATCGAGAGCCAGCAT‑GAAAAAAATAATTACGAAAGCGTCAGGCCTTGTCATCGTGTGCGCAACCCCGTTGGCAATGGGGTTCACCTTTgaaa                                                                                                                            <  4:51992/147‑1 (MQ=255)                         tcttcACCTCGGATAACCTGCGAGCCATCAGCGAGCGCTAAAACCGACGTCAGATCGAGAGCCAGCAT‑GAAAAAAATAATTACGAAAGCGTCAGGCCTTGTCATCGTGTGCGCAACCCCGTTGGCAATGGGGTTCACCTTTgaaaccg                                                                                                                         >  3:89884/1‑148 (MQ=255)                             cACCTCGGATAACCTGCGAGCCATCAGCGAGCGCTAAAACCGACGTCAGATCGAGAGCCAGCAT‑GAAAAAAATAATTACGAAAGCGTCAGGCCTTGTCATCGTGTGCGCAACCCCGTTGGCAATGGGGTTCACCTTTgaaaccgaaac                                                                                                                     <  5:238583/148‑1 (MQ=255)                              aCCTCGGATAACCTGCGAGCCATCAGCGAGCGCTAAAACCGACGTCAGATCGAGAGCCAGCAT‑GAAAAAAATAATTACGAAAGCGTCAGGCCTTGTCATCGTGTGCGCAACCCCGTTGGCAATGGGGTTCACCTTTGAAACCGAAAccc                                                                                                                   >  5:295540/1‑149 (MQ=255)                                cTCGGATAACCTGCGAGCCATCAGCGAGCGCTAAAACCGACGTCAGATCGAGAGCCAGCAT‑GAAAAAAATAATTACGAAAGCGTCAGGCCTTGTCATCGTGTGCGCAACCCCGTTGGCAATGGGGTTCACCTTTGAAACCGAAACCCt                                                                                                                  >  2:383894/1‑148 (MQ=255)                                cTCGGATAACCTGCGAGCCATCAGCGAGCGCTAAAACCGACGTCAGATCGAGAGCCAGCAT‑GAAAAAAATAATTACGAAAGCGTCAGGCCTTGTCATCGTGTGCGCAACCCCGTTGGCAATGGGGTTCACCTTTGAAACCGAAACCCt                                                                                                                  >  2:188920/1‑148 (MQ=255)                                   ggATAACCTGCGAGCCATCAGCGAGCGCTAAAACCGACGTCAGATCGAGAGCCAGCAT‑G‑AAAAAATAATTACGAAAGCGTCAGGCCTTGTCATCGTGTGCGCAACCCCGTTGGCAATGGGGTTCACCTTTGAAACCGAAACCCTCAATg                                                                                                             >  7:168650/1‑149 (MQ=255)                                    gATAACCTGCGAGCCATCAGCGAGCGCTAAAACCGACGTCAGATCGAGAGCCAGCAT‑GAAAAAAATAATTACGAAAGCGTCAGGCCTTGTCATCGTGTGCGCAACCCCGTTGGCAATGGGGTTCACCTTTGAAACCGAAACCCTCAATg                                                                                                             <  8:242639/149‑1 (MQ=255)                                     aTAACCTGCGAGCCATCAGCGAGCGCTAAAACCGACGTCAGATCGAGAGCCAGCAT‑GAAAAAAATAATTACGAAAGCGTCAGGCCTTGTCATCGTGTGCGCAACCCCGTTGGCAATGGGGTTCACCTTTGAAACCGAAACCCTCAATgg                                                                                                            >  3:278038/1‑149 (MQ=255)                                        aCCTGCGAGCCATCAGCGAGCGCTAAAACCGACGTCAGATCGAGAGCCAGCAT‑GAAAAAAATAATTACGAAAGCGTCAGGCCTTGTCATCGTGTGCGCAACCCCGTTGGCAATGGGGTTCACCTTTGAAACCGAAACCCTCAATGGtt                                                                                                          <  8:171842/148‑1 (MQ=255)                                         ccTGCGAGCCATCAGCGAGCGCTAAAACCGACGTCAGATCGAGAGCCAGCAT‑GAAAAAAATAATTACGAAAGCGTCAGGCCTTGTCATCGTGTGCGCAACCCCGTTGGCAATGGGGTTCACCTTTGAAACCGAAACCCTCAATGGTTcc                                                                                                        >  5:49604/1‑149 (MQ=255)                                                  cATCAGCGAGCGCTAAAACCGACGTCAGATCGAGAGCCAGCAT‑GAAAAAAATAATTACGAAAGCGTCAGGCCTTGTCATCGTGTGCGCAACCCCGTTGGCAATGGGGTTCACCTTTGAAACCGAAACCCTCAATGGTTCCTTCGACTc                                                                                                >  5:150256/1‑148 (MQ=255)                                                    tCAGCGAGCGCTAAAACCGACGTCAGATCGAGAGCCAGCAT‑GAAAAAAATAATTACGAAAGCGTCAGGCCTTGTCATCGTGTGCGCAACCCCGTTGGCAATGGGGTTCACCTTTGAAACCGAAACCCTCAATGGTTCCTTCGACTccac                                                                                             <  7:229171/149‑1 (MQ=255)                                                       gcgagcgCTAAAACCGACGTCAGATCGAGAGCCAGCAT‑GAAAAAAATAATTACGAAAGCGTCAGGCCTTGTCATCGTGTGCGCAACCCCGTTGGCAATGGGGTTCACCTTTGAAACCGAAACCCTCAATGGTTCCTTCGACTccacca                                                                                           >  4:259899/1‑148 (MQ=255)                                                       gcgagcgCTAAAACCGACGTCAGATCGAGAGCCAGCAT‑GAAAAAAATAATTACGAAAGCGTCAGGCCTTGTCATCGTGTGCGCAACCCCGTTGGCAATGGGGTTCACCTTTGAAACCGAAACCCTCAATGGTTCCTTCGACTCCACCAt                                                                                          >  2:753/1‑149 (MQ=255)                                                         gagcgCTAAAACCGACGTCAGATCGAGAGCCAGCAT‑GAAAAAAATAATTACGAAAGCGTCAGGCCTTGTCATCGTGTGCGCAACCCCGATGTCAATGGGGTTAACCTTTGAAACCGAAACCCTCAATGGTTCCTTCGACTCCACCATCa                                                                                        >  1:152964/1‑149 (MQ=255)                                                               tAAAACCGACGTCAGATCGAGGGCCAGCAT‑GAAAAAAATAATTACGAAAGCGTCAGGCCTTGTCATCGTGTGCGCAACCCCGTTGGCAATGGGGTTCACCTTTGAAACCGAAACCCTCAATGGTTCCTTCGACTCCACCATCACCGCAg                                                                                  >  1:378275/1‑149 (MQ=255)                                                                 aaaCCGACGTCAGATCGAGAGCCAGCATGGAAAAAAATAATTACGAAAGCGTCAGGCCTTGTCATCGTGTGCGCAACCCCGTTGGCAATGGGGTTCACCTTTGAAACCGAAACCCTCAATGGTTCCTTCGACTCCACCATCACCGCAgg                                                                                 >  1:23579/1‑149 (MQ=255)                                                                 aaaCCGACGTCAGATCGAGAGCCAGCAT‑GAAAAAAATAATTACGAAAGCGTCAGGCCTTGTCATCGTGTGCGCAACCCCGTTGGCAATGGGGTTCACCTTTGAAACCGAAACCCTCAATGGTTCCTTCGACTCCACCATCACCGCAgg                                                                                 <  1:262250/148‑1 (MQ=255)                                                                   aCCGACGTCAGATCGAGAGCCAGCAT‑GAAAAAAATAATTACGAAAGCGTCAGGCCTTGTCATCGTGTGCGCAACCCCGTTGGCAATGGGGTTCACCTTTGAAACCGAAACCCTCAATGGTTCCTTCGACTCCACCATCACCGCAGGCAt                                                                              >  4:32777/1‑149 (MQ=255)                                                                      gACGTCAGATCGAGAGCCAGCAT‑GAAAAAAATAATTACGAAAGCGTCAGGCCTTGTCATCGTGTGCGCAACCCAGTTGGCAATGGGGTTCACCTTTGAAACCGAAACCCTCAATGGTTCCTTCGACTCCACCATCACCGCAGGCATggg                                                                           >  7:81107/1‑149 (MQ=255)                                                                          tCAGATCGAGAGCCAGCAT‑GAAAAAAATAATTACGAAAGCGTCAGGCCTTGTCATCGTGTGCGCAACCCCGTTGGCAATGGGGTTCACCTTTGAGACCGAAACCCTCAATGGTTGCTTCGACTCCACCATCACCGCAGGCATGGGCCTg                                                                       >  6:375567/1‑149 (MQ=255)                                                                          tCAGATCGAGAGCCAGCAT‑GAAAAAAATAATTACGAAAGCGTCAGGCCTTGTCATCGTGTGCGCAACCCCGTTGGCAATGGGGTTCACCTTTGAAACCGAAACCCTCAATGGTTCCTTCGACTCCACCATCACCGCAGGCATGGGCCTg                                                                       >  3:2128/1‑149 (MQ=255)                                                                          tCAGATCGAGAGCCAGCAT‑GAAAAAAATAATTACGAAAGCGACAGGCCTTCTCATCGTGTGCGCAACCCCGTTGCCAATGGCGTTCACCATTGAAACCCAAACCCTCAATGGTTCCTTCGACTCAACCATCACCGCAGGCATCGCCCTg                                                                       >  7:29222/1‑149 (MQ=255)                                                                              aTCGAGAGCCAGCAT‑GAAAAAAATAATTACGAAAGCGTCAGGCCTTGTCATCGTGTGCGCAACCCCGTTGGCAATGGGGTTCGCCTTTGAAACCGAAACCCTCAATGGTTCCTTCGACTCCACCATCACCGCAGGCATGGGCCTGCGTa                                                                   >  6:282475/1‑149 (MQ=255)                                                                               tCGAGAGCCAGCAT‑GAAAAAAATAATTACGAAAGCGTCAGGCCTTGTCATCGTGTGCGCAACCCCGTTGGCAATGGGGTTCACCTTTGAAACCGAAACCCTCAATGGTTCCTTCGACTCCACCATCACCGCAGGCATGGGCCTGCGTAc                                                                  >  8:63060/1‑149 (MQ=255)                                                                                  agagCCAGCAT‑GAAAAAAATAATTACGAAAGCGTCAGGCCTTGTCATCGTGTGCGCAACCCCGTTGGCAATGGGGTTCACCTTTGAAACCGAAACCCTCAATGGTTCCTTCGACTCCACCATCACCGCAGGCATGGGCCTGCGTACGGa                                                               <  1:753/149‑1 (MQ=255)                                                                                       cAGCAT‑GAAAAAAATAATTACGAAAGCGTCAGGCCTTGTCATCGTGTGCGCAACCCCGTTGGCAATGGGGTTCACCTTTGAAACCGAAACCCTCAATGGTTCCTTCGACTCCACCATCACCGCAGGCATGGGCCTGCGTACGGAGTcgc                                                          <  3:32777/149‑1 (MQ=255)                                                                                         gCAT‑GAAAAAAATAATTACGAAAGCGTCAGGCCTTGTCATCCTGTGCGCAACCCCGTTGGCAATGGGGTTCACCTTTGAAACCGAAACCCTCAATGGTTCCTTCGACTCCACCATCACCGCAGGCATGGGCCTGCGTACGGAGTCGCgg                                                        <  5:337899/149‑1 (MQ=255)                                                                                               aaaaaaaTAATTACGAAAGCGTCAGGCCTTGTCATCGTGTGCGCAACCCCGTTGGCAATGGGGTTCACCTTTGAAACCGAAACCCTCAATGGTTCCTTCGACTCCACCATCACCGCAGGCATGGGCCTGCGTACGGAGTCGCGGGGCt                                                    <  2:378275/148‑1 (MQ=255)                                                                                               aaaaaaaTAATTACGAAAGCGTCAGGCCTTGTCATCGTGTGCGCAACCCCGTTGGCAATGGGGTTCACCTTTGAAACCGAAACCCTCAATGGTTCCTTCGACTCCACCATCACCGCAGGCATGGGCCTGCGTACGGAGTCGCGGGGCt                                                    <  2:23579/148‑1 (MQ=255)                                                                                                aaaaaaTAATTACGAAAGCGTCAGGCCTTGTCATCGTGTGCGCAACCCCGTTGGCAATGGGGTTCACCTTTGAAACCGAAACCCTCAATGGTTCCTTCGACTCCACCATCACCGCAGGCATGGGCCTGCGTACGGAGTCGCGGGGCTGc                                                  <  1:188920/149‑1 (MQ=255)                                                                                                 aaaaaTAATTACGAAAGCGTCAGGCCTTGTCATCGTGTGCGCAACCCCGTTGGCAATGGGGTTCACCTTTGAAACCGAAACCCTCAATGGTTCCTTCGACTCCACCATCACCGCAGGCATGGGCCTGCGTACGGAGTCGCGGGGCTGc                                                  >  1:25070/1‑148 (MQ=255)                                                                                                  aaaaTAATTACGAAAGCGTCAGGCCTTGTCATCGTGTGCGCAACCCCGTTGGCAATGGGGTTCACCTTTGAAACCGAAACCCTCAATGGTTCCTTCGACTCCACCATCACCGCAGGCATGGGCCTGCGTACGGAGTCGCGGGGCTGCaa                                                <  8:168650/149‑1 (MQ=255)                                                                                                    aataatTACGAAAGCGTCAGGCCTTGTCATCGTGTGCGCAACCCCGTTGGCAATGGGGTTCACCTTTGAAACCGAAACCCTCAATGGTTCCTTCGACTCCACCATCACCGCAGGCATGGGCCTGCGTACGGAGTCGCGGGGCTGCAAcc                                              <  8:29222/149‑1 (MQ=255)                                                                                                                aGCGTCAGGCCTTGTCATCGTGTGCGCAACCCCGTTGGCAATGGGGTTCACCTTTGAAACCGAAACCCTCAATGGTTCCTTCGACTCCACCATCACCGCAGGCATGGGCCTGCGTACGGAGTCGCGGGGCTGCAACCTGATCAACCAgg                                  <  3:259899/149‑1 (MQ=255)                                                                                                                  cGTCAGGCCTTGTCATCGTGTGCGCAACCCCGTTGGCAATGGGGTTCACCTTTGAAACCGAAACCCTCAATGGTTCCTTCGACTCCACCATCACCGCAGGCATGGGCCTGCGTACGGAGTCGCGGGGCTGCAACCTGATCAACCAgggg                                <  6:295540/149‑1 (MQ=255)                                                                                                                   gTCAGGCCTTGTCATCGTGTGCGCAACCCCGTTGGCAATGGGGTTCACCTTTGAAACCGAAACCCTCAATGGTTCCTTCGACTCCACCATCACCGCAGGCATGGGCCTGCGTACGGAGTCGCgg                                                        <  5:395023/124‑1 (MQ=255)                                                                                                                   gTCAGGCCTTGTCATCGTGTGCGCAACCCCGTTGGCAATGGGGTTCACCTTTGAAACCGAAACCCTCAATGGTTCCTTCGACTCCACCATCACCGCAGGCATGGGCCTGCGTACGGAGTCGCgg                                                        >  6:395023/1‑124 (MQ=255)                                                                                                                        gCCTTGTCATCGTGTGCGCAACCCCGTTGGCAATGGGGTTCACCTTTgaaaccgaaac                                                                                                                     <  2:207633/58‑1 (MQ=255)                                                                                                                        gCCTTGTCATCGTGTGCGCAACCCCGTTGGCAATGGGGTTCACCTTTgaaaccgaaac                                                                                                                     >  1:207633/1‑58 (MQ=255)                                                                                                                           ttGTCATCGTGTGCGCAACCCCGTTGGCAATGGGGTTCACCTTTGAAACCGAAACCCTCAATGGTTCCTTCGACTCCACCATCACCGCAGGCATGGGCCTGCGTACGGAGTCGCGGGGCTGCAACCTGATCAACCAGGGGCCGACCggg                       >  1:214233/1‑149 (MQ=255)                                                                                                                             gTCATCGTGTGCGCAACCCCGTTGGCAATGGGGTTCACCTTTGAAACCGAAACCCTCAATGGTTCCTTCGACTCCACCATCACCGCAGGCATGGGCCTGCGTACGGAGTCGCGGGGCTGCAACCTGATCAACCAGGGGCCGACCggg                       <  1:383894/147‑1 (MQ=255)                                                                                                                             gTCATCGTGTGCGCAACCCCGTTGGCAATGGGGTTCACCTTTGAAACCGAAACCCTCAATGGTTCCTTCGACTCCACCATCACCGCAGGCATGGGCCTGCGTACGGAGTCGCGGGGCTGCAACCTGATCAACCAGGGGCCGACCGGGca                     >  7:244239/1‑149 (MQ=255)                                                                                                                              tCATCGTGTGCGCAACCGCGTTGGCAATGGGGTTCACCTTTGAAACCGAAACCCTCAATGGTTCCTTCGACTCCACAATCACCGCAGGCAGGGGCCTGCGTACGGAGTCGCGGGGCTGCAACCTGATCAACAAGGGGCCGACCGGGca                     >  8:280603/1‑148 (MQ=255)                                                                                                                                aTCGTGTGCGCAACCCCGTTGGCAATGGGGTTCACCTTTGAAACCGAAACCCTCAATGGTTCCTTCGACTCCACCATCACCGCAGGCATGGGCCTGCGTACGGAGTCGCGGGGCTGCAACCTGATCAACCAGGGGCCGACCGGGCACaa                  >  1:270397/1‑149 (MQ=255)                                                                                                                                    tgtgCGCAACCCCGTTGGCAATGGGGTTCACCTTTGAAACCGAAACCCTCAATGGTTCCTTCGACTCCACCATCACCGCAGGCATGGGCCTGCGTACGGAGTCGCGGGGCTGCAACCTGATCAACCAGGGGCCGACCGGGCACAACGTg              <  4:2128/149‑1 (MQ=255)                                                                                                                                      aggtgaaCCCCATTGCCAACGGGGTTCACCTTTGAAACCGAAACCCTCAATGGTTCCTTCGACTCCACCATCACCGCAGGCATGGGCCTGCGTACGGAGTCGCgg                                                        <  2:175973/100‑1 (MQ=255)                                                                                                                                      aggtgaaCCCCATTGCCAACGGGGTTCACCTTTGAAACCGAAACCCTCAATGGTTCCTTCGACTCCACCATCACCGCAGGCATGGGCCTGCGTACGGAGTCGCgg                                                        >  1:175973/6‑105 (MQ=255)                                                                                                                                           aaCCCCGTTGGCAATGGGGTTCACCTTTGAAACCGAAACCCTCAATGGTTCCTTCGACTCCACCATCACCGCAGGCATGGGCCTGCGTACGGAGTCGCGGGGCTGCAACCTGATCAACCAGGGGCCGACCGGGCACAACGTGCCCTCCg       >  2:322490/1‑149 (MQ=255)                                                                                                                                             ccccGTTGGCAATGGGGTTCACCTTTGAAACCGAAACCCTCAATGGTTCCTTCGACTCCACCATCACCGCAGGCATGGGCCTGCGTACGGAGTCGCGGGGCTGCAACCTGATCAACCAGGGGCCGACCGGGCACAACGTGCCCTCCGGc     >  6:297613/1‑149 (MQ=255)                                                                                                                                                 gTTGGCAATGGGGTTCACCTTTGAAACCGAAACCCTCAATGGTTCCTTCGACTCCACCATCACCGCAGGCATGGGCCTGCGTACGGAGTCGCGGTTCTTCAACCTGATCAACCAGGGGCCGACCGGGCACAACGTGCCCTCCGGcatc  >  8:312214/1‑145 (MQ=255)                                                                                                                                                     |                                                                                                                                                GCAACCCTCGCCTGAAGCGCGACTTCTTCACCTCGGATAACCTGCGAGCCATCAGCGAGCGCTAAAACCGACGTCAGATCGAGAGCCAGCAT‑GAAAAAAATAATTACGAAAGCGTCAGGCCTTGTCATCGTGTGCGCAACCCCGTTGGCAATGGGGTTCACCTTTGAAACCGAAACCCTCAATGGTTCCTTCGACTCCACCATCACCGCAGGCATGGGCCTGCGTACGGAGTCGCGGGGCTGCAACCTGATCAACCAGGGGCCGACCGGGCACAACGTGCCCTCCGGCTGC  >  NC_002947/3204607‑3204897

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