Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,780,952	T→G	5.5%	S84A (TCT→GCT)	PP_2434 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,780,952	0	T	G	5.5%	220.5 / 3.9	73	S84A (TCT→GCT)	PP_2434	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (36/33);  new base G (2/2);  total (38/35)
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.32e-01

CGCCACCAGCCCTCCTGGGGCGATGCTGGAAGTGTTGGTGCACCTGGAGATCGACCCCGAGGATTTCCCTACCACCATGCGCTTGCTTCGCATCGAACTGCCCGACACGGTTTCCCAGGCGCAGTTGCCCGCTTTGCAGCCCGGCTGGTCTGCTCAGCCTGAGTTGACCCGAACGCTGGGCAATC‑GCTTTCTGGATGACTGTTCGGCATTGCTGCTGCCGGTGCCAAGCGCGATCATGCCCAGCACCACCAATTACCTGTTCAACCCACGGCACCCGCAGGCGCAGAGCGCGAAA  >  NC_002947/2780804‑2781098                                                                                                                                                     |                                                                                                                                                    cgccACCAGCCCTCCTGGGGCGATGCTGGAAGTGTTGGTGCACCTGGAGATCGACCCCGAGGATTTCCCTACCACCATGCGCTTGCTTCGCATCGAACTGCCCGACACGGTTTCCCAGGCGCAGTTGCCCGCTTTGTAGCCCGGCTGGt                                                                                                                                                     <  3:148491/149‑1 (MQ=255)   ccaccaGCCCTCCTGGGGCGATGCTGGAAGTGTTGGTGCACCTGGAGATCGACCCCGAGGATTTCCCTACCACCATGCGCTTGCTTCGCATCGAACTGCCCGACACGGTTTCCCAGGCGCAGTTGCCCGCTTTGCAGCCCGGCTGGTCt                                                                                                                                                   >  5:258928/1‑149 (MQ=255)      ccaGCCCTCCTGGGGCGATGCTGGAAGTGTTGGTGCACCTGGAGATCGACCCCGAGGATTTCCCTACCACCATGCGCTTGCTTCGCATCGAACTGCCCGACACGGTTTCCCAGGCGCAGTTGCCCGCTTTGCAGCCCGGCTGGTCTGCt                                                                                                                                                <  5:245110/149‑1 (MQ=255)        aGCCCTCCTGGGGCGATGCTGGAAGTGTTGGTGCACCTGGAGATCGACCCCGAGGATTTCCCTACCACCATGCGCTTGCTTCGCATCGAACTGCCCGACACGGTTTCCCAGGCGCAGTTGCCCGCTTTGCAGCCCGGCTGGTCTGCTCa                                                                                                                                              <  6:85685/149‑1 (MQ=255)                  gggCGATGCTGGAAGTGTTGGTGCACCTGGAGATCGACCCCGAGGATTTCCCTACCACCATGCGCTTGCTTCGCATCGAACTGCCCGACACGGTTTCCCAGGCGCAGTTGCCCGCTTTGCAGCCCGGCTGGTCTGCTCAGCCTGAGtt                                                                                                                                     >  4:259258/1‑148 (MQ=255)                    gCGATGCTGGAAGTGTTGGTGCACCTGGAGATCGACCCCGAGGATTTCCCTACCACCATGCGCTTGCTTCGCATCGAACTGCCCGACACGGTTTCCCAGGCGCAGTTGCCCGCTTTGCAGCCCGGCTGGTCTGCTCAGCCTGAGTTGAc                                                                                                                                  <  2:93713/149‑1 (MQ=255)                           tGGAAGTGTTGGTGCACCTGGAGATCGACCCCGAGGATTTCCCTACCACCATGCGCTTGCTTCGCATCGAACTGCCCGACACGGTTTCCCAGGCGCAGTTGCCCGCTTTGCAGCCCGGCTGGTCTGCTCAGCCTGAGTTGACCCGAACg                                                                                                                           >  6:191606/1‑149 (MQ=255)                            ggAAGTGTTGGTGCACCTGGAGATCGACCCCGAGGATTTCCCTACCACCATGCGCTTGCTTCGCATCGAACTGCCCGACACGGTTTCCCAGGCGCAGTTGCCCGCTTTGCAGCCCGGCTGGTCTGCTCAGCCTGAGTTGACCCGAACGc                                                                                                                          >  1:179446/1‑149 (MQ=255)                             gtattgtTGGTGCACCTGGAGATCGACCCCGAGGATTTCCCTACCACCATGCGCTTGCTTCGCATCGAACTCGCCGACACGGTTTCCCAGGCGCAGTTGCCCGCTTTGCAGCCCGGCTGGTCTGCTCAGCCTGAGTTGACCCGAACGCt                                                                                                                         <  8:157270/145‑1 (MQ=255)                                gtgtTGGTGCACCTGGAGATCGACCCCGAGGATTTCCCTACCACCATGCGCTTGCTTCGCATCGAACTGCCCGACACGGTTTCCCAGGCGCAGTTGCCCGCTTTGCAGCCCGGCTGGTCTGCTCAGCCTGAGTTGACCCGAACGCTggg                                                                                                                      <  1:264743/149‑1 (MQ=255)                                   ttGGTGCACCTGGAGATCGACCCCGAGGATTTCCCTACCACCATGCGCTTGCTTCGCATCGAACTGCCCGACACGGTTTCCCAGGCGCAGTTGCCCGCTTTGCAGCCCGGCTGGTCTGCTCAGCCTGAGTTGACCCGAACGCTGGGCaa                                                                                                                   >  8:174368/1‑149 (MQ=255)                                   ttGGTGCACCTGGAGATCGACCCCGAGGATTTCCCTACCACCATGCGCTTGCTTCGCATCGAACTGCCCGACACGGTTTCCCAGGCGCAGTTGCCCGCTTTGCAGCCCGGCTGGTCTGCTCAGCCTGAGTTGACCCGAACGCTGGGCaa                                                                                                                   >  2:373773/1‑149 (MQ=255)                                      gTGCACCTGGAGATCGACCCCGAGGATTTCCCTACCACCATGCGCTTGCTTCGCATCGAACTGCCCGACACGGTTTCCCAGGCGCAGTTGCCCGCTTTGCAGCCCGGCTGGTCTGCTCAGCCTGAGTTGACCCGAACGCTGGGCAATC‑g                                                                                                               >  5:348602/1‑149 (MQ=255)                                          aCCTGGAGATCGACCCCGAGGATTTCCCTACCACCATGCGCTTGCTTCGCATCGAACTGCCCGACACGGTTTCCCAGGCGCAGTTGCCCGCTTTGCAGCCCGGCTGGTCTGCTCAGCCTGAGTTGACCCGAACGCTGGGCAATC‑GCttt                                                                                                           >  8:276071/1‑149 (MQ=255)                                           ccTGGAGATCGACCCCGAGGATTTCCCTACCACCATGCGCTTGCTTCGCATCGAACTGCCCGACACGGTTTCCCAGGCGCAGTTGCCCGCTTTGCAGCCCGGCTGGTCTGCTCAGCCTGAGTTGACCCGAACGCTGGGCAATC‑GCTTTc                                                                                                          >  5:290821/1‑149 (MQ=255)                                            cTGGAGATCGACCCCGAGGATTTCCCTACCACCATGCGCTTGCTTCGCATCGAACTGCCCGACACGGTTTCCCAGGCGCAGTTGCCCGCTTTGCAGCCCGGCTGGTCTGCTCAGCCTGAGTTGACCCGAACGCTGGGCAATC‑GCTTTCt                                                                                                         >  2:222875/1‑149 (MQ=255)                                              ggAGATCGACCCCGAGGATTTCCCTACCACCATGCGCTTGCTTCGCATCGAACTGCCCGACACGGTTTCCCAGGCGCAGTTGCCCGCTTTGCAGCCCGGCTGGTCTGCTCAGCCTGAGTTGACCCGAACGCTGGGCAATC‑Gc                                                                                                              >  5:60084/1‑142 (MQ=255)                                              ggAGATCGACCCCGAGGATTTCCCTACCACCATGCGCTTGCTTCGCATCGAACTGCCCGACACGGTTTCCCAGGCGCAGTTGCCCGCTTTGCAGCCCGGCTGGTCTGCTCAGCCTGAGTTGACCCGAACGCTGGGCAATC‑Gc                                                                                                              <  6:60084/142‑1 (MQ=255)                                               gagaTCGACCCCGAGGATTTCCCTACCACCATGCGCTTGCTTCGCATCGAACTGCCCGACACGGTTTCCCAGGCGCAGTTGCCCGCTTTGCAGCCCGGCTGGTCTGCTCAGCCTGAGTTGACCCGAACGCTGGGCAATC‑GCTTTCTgg                                                                                                       <  6:258928/148‑1 (MQ=255)                                               gagaTCGACCCCGAGGATTTCCCTACCACCATGCGCTTGCTTCGCATCGAACTGCCCGACACGGTTTCCCAGGCGCAGTTGCCCGCTTTGCAGCCCGGCTGGTCTGCTCAGCCTGAGTTGACCCGAACGCTGGGCAATC‑GCTTTCTGGa                                                                                                      <  8:49708/149‑1 (MQ=255)                                                agaTCGACCCCGAGGATTTCCCTACCACCATGCGCTTGCTTCGCATCGAACTGCCCGACACGGTTTCCCAGGCGCAGTTGCCCGCTTTGCAGCCCGGCTGGTCTGCTCAGCCTGAGTTGACCCGAACGCTGGGCAATC‑GCTTTCTGGAt                                                                                                     <  3:259258/149‑1 (MQ=255)                                                   tCGACCCCGAGGATTTCCCTACCACCATGCGCTTGCTTCGCATCGAACTGCCCGACACGGTTTCCCAGGCGCAGTTGCCCGCTTTGCAGCCCGGCTGGTCTGCTCAGCCTGAGTTGACCCGAACGCTGGGCAATC‑GCTTTCTGGATGAc                                                                                                  >  7:69598/1‑149 (MQ=255)                                                       ccccGAGGATTTCCCTACCACCATGCGCTTGCTTCGCATCGAACTGCCCGACACGGTTTCCCAGGCGCAGTTGCCCGCTTTGCAGCCCGGCTGGTCTGCTCAGCCTGAGTTGACCCGAACGCTGGGCAATC‑GCTTTCTGGATGACTGtt                                                                                              <  8:44867/149‑1 (MQ=255)                                                        cccGAGGATTTCCCTACCACCATGCGCTTGCTTCGCATCGAACTGCCCGACACGGTTTCCCAGGCGCAGTTGCCCGCTTTGCAGCCCGGCTGGTCTGCTCAGCCTGAGTTGACCCGAACGCTGGGCAATC‑GCTTTCTGGATGACTGTTc                                                                                             >  5:100783/1‑149 (MQ=255)                                                          cGAGGATTTCCCTACCACCATGCGCTTGCTTCGCATCGAACTGCCCGACACGGTTTCCCAGGCGCAGTTGCCCGCTTTGCAGCCCGGCTGGTCTGCTCAGCCTGAGTTGACCCGAACGCTGGGCAATC‑GCTTTCTGGATGACTGTTCgg                                                                                           >  6:380628/1‑149 (MQ=255)                                                           gAGGATTTCCCTACCACCATGCGCTTGCTTCGCATCGCACTGCCCGACACGGTTTCCCAGGCGCAGTTGCCCGCTTTGCAGCCCGGCTGGTGTGCTCAGCCTGAGTTGACCCGAACGCTGGGCAATC‑GCGTTCTGGATGCCTGTTCgg                                                                                           >  6:315739/1‑148 (MQ=255)                                                                tttCCCTACCACCATGCGCTTGCTTCGCATCGAACTGCCCGACACGGTTTCCCAGGCGCAGTTGCCCGCTTTGCAGCCCGGCTGGTCTGCTCAGCCTGAGTTGACCCGAACGCTGGGCAATC‑GCTTTCTGGATGACTGTTCGGCATtgc                                                                                     >  4:45158/1‑149 (MQ=255)                                                                        ccaccaTGCGCTTGCGTCGCATCGAACTGCCCGACACGGTTTCCCAGGCGCAGTTGCCCGCTTTGCAGCCCGGCTGGTCTGCTCAGCCTGCGTTGCCCCGAACGCTGGGCAATC‑GCTTTCTGGATGACTGTTCGGCATTGCTGCTGCCg                                                                             >  8:370510/1‑149 (MQ=255)                                                                         caccaTGCGCTTGCTTCGCATCGAACTGCCCGACACGGTTTCCCAGGCGCAGTTGCCCGATTTGCAGCCCGGCTGGTCTGCTCAGCCTGAGTTGACCCGAACGCTGGGCAATC‑GCTTTCTGGATGACTGTTCGGCATTGCTGCTGCCgg                                                                            >  2:317280/1‑149 (MQ=255)                                                                             aTGCGCTTGCTTCGCATCGAACTGCCCGACACGGTTTCCCAGGCGCAGTTGCCCGCTTTGCAGCCCGGCTGGTCTGCTCAGCCTGAGTTGACCCGAACGCTGGGCAATC‑GCTTTCTGGATGACTGTTCGGCATTGCTGCTGCCGGTGcc                                                                        <  6:100783/149‑1 (MQ=255)                                                                                 gcttgcttCGCATCGAACTCCCCGACACGGTTTCCCAGGCGCAGTTGCCCGCTATGCAGCCCGGCTGGTATGCTCAGCCTGACTTGACCCGAACGCTGGGCAATC‑GCTTTCTGGATGACTGTTCTGCATTGCTGCTGCCGGTGCCAAgc                                                                    >  3:273301/1‑149 (MQ=255)                                                                                      cttCGCATCGAACTGCCCGACACGGTTTCCCAGGCGCAGTTGCCCGCTTTGCAGCCCGGCTGGTCTGCTCAGCCTGAGTTGACCCGAACGCTGGGCAATC‑GCTTTCTGGATGACTGTTCGGCATTGCTGCTGCCGGTGCCAAGCGCGAt                                                               <  3:45158/149‑1 (MQ=255)                                                                                          gCATCGAACTGCCCGACACGGTTTCCCAGGCGCAGTTGCCCGCTTTGCAGCCCGGCTGGTCTGCTCAGCCTGAGTTGACCCGAACGCTGGGCAATC‑GCTTTCTGGATGACTGTTCGGCATTGCTGCTGCCGGTGCCAAGCGCGATCAt                                                            <  7:276071/148‑1 (MQ=255)                                                                                          gCATCGAACTGCCCGACACGGTTTCCCAGGCGCAGTTGCCCGCTTTGCAGCCCGGCTGGTCTGCTCAGCCTGAGTTGACCCGAACGCTGGGCAATC‑GCTTTCTGGATGACTGTTCGGCATTGCTGCTGCCGGTGCCAAGCGCGATCAt                                                            <  7:174368/148‑1 (MQ=255)                                                                                          gCATCGAACTGCCCGACACGGTTTCCCAGGCGCAGTTGCCCGCTTTGCAGCCCGGCTGGTCTGCTCAGCCTGAGTTGACCCGAACGCTGGGCAATC‑GCTTTCTGGATGACTGTTCGGCATTGCTGCTGCCGGTGCCAAGCGCGATCATg                                                           <  1:222875/149‑1 (MQ=255)                                                                                           cATCGAACTGCCCGACACGGTTTCCCAGGCGCAGTTGCCCGCTTTGCAGCCCGGCTGGTCTGCTCAGCCTGAGTTGACCCGAACGCTGGGCAATC‑GCTTTCTGGATGACTGTTCGGCATTGCTGCTGCCGGTGCCAAGCGCGATCATGc                                                          <  1:373773/149‑1 (MQ=255)                                                                                             tCGAACTGCCCGACACGGTTTCCCAGGCGCAGTTGCCCGCTTTGCAGCCCGGCTGGTCTGCTCAGCCTGAGTTGACCCGAACGCTGGGCAATC‑GCTTTCTGGATGACTGTTCGGCATTTCTGCTGCCGGTGCCAAGCGCGATCATGccc                                                        <  6:290821/149‑1 (MQ=255)                                                                                             tCGAACTGCCCGACACGGTTTCCCAGGCGCAGTTGCCCGCTTTGCAGCCCGGCTGGTCTGCTCAGCCTGAGTTGACCCGAACGCTGGGCAATC‑GCTTTCTGGATGACTGTTCGGCATTGCTGCTGCCGGTGCCAAGCGCGATCATGccc                                                        <  8:69598/149‑1 (MQ=255)                                                                                               gAACTGCCCGACACGGTTTCCCAGGCGCAGTTGCCCGCTTTGCAGCCCGGCTGGTCTGCTCAGCCTGAGTTGACCCGAACGCTGGGCAATC‑GCTTTCTGGATGACTGTTCGGCATTGCTGCTGCCGGTGCCAAGCGCGATCATGCCCa                                                       <  2:179446/148‑1 (MQ=255)                                                                                               gAACTGCCCGACACGGTTTCCCAGGCGCAGTTGCCCGCTTTGCAGCCCGGCTGGTCTGCTCAGCCTGAGTTGACCCGAACGCTGGGCAATC‑GCTTTCTGGATGACTGTTCGGCATTGCTGCTGCCGGTGCCAAGCGCGATCATGCCCAg                                                      <  8:331241/149‑1 (MQ=255)                                                                                                  cTGCACGACACGGTTTCCCAGGCGCAGTTGCTCGCTTTGCAGCTCGGCTTGTCTGCTCAGCCTGAGTTGACCCGAACGCTTGGCAATC‑GCTTTCTGGATGACTGTTCGGCATTGCTGCTGCTGGTGCCAAGCGCGATCATGCCCAGcac                                                   <  4:281345/149‑1 (MQ=255)                                                                                                     cccGACACGGTTTCCCAGGCGCAGTTGCCCGCTTTGCAGCCCGGCTGGTCTGCTCAGCCTGAGTTGACCCGAACGCTGGGCAATC‑GCTTTCTGGATGACTGTTCGGCATTGCTGCTGCCGGTGCCAAGCGCGATCATGCCCAGcaccac                                                >  3:347040/1‑149 (MQ=255)                                                                                                      ccGACACGGTTTCCCAGGCGCAGTTGCCCGCTTTGCAGCCCGGCTGGTCTGCTCAGCCTGAGTTGACCCGAACGCTGGGCAATC‑GCTTTCTGGATGACTGTTCGGCATTGCTGCTGCCGGTGCCAAGCGCGATCATGCCCAGcaccac                                                >  6:290548/1‑148 (MQ=255)                                                                                                             ggTTTCCCAGGCGCAGTTGCCCGCTTTGCAGCCCGGCTGGTCTGCTCAGCCTGAGTTGACCCGAACGCTGGGCAATC‑GCTTTCTGGATGACTGTTCGGCATTGCTGCTGCCGGTGCCAAGCGCGATCATGCCCAGCACCACCAATTAcc                                        <  5:191606/149‑1 (MQ=255)                                                                                                             ggTTTCCCAGGCGCAGTTGCCCGCTTTGCAGCCCGGCTGGTCTGCTCAGCCTGAGTTGACCCGAACGCTGGGCAATC‑GCTTTCTGGATGACTGTTAAACATTGCTTCAGCCGGTGCCAAGCGCGATAATTCCCAGCACCACCAATTAcc                                        >  1:323860/1‑149 (MQ=255)                                                                                                              gTTTCCCAGGCGCAGTTGCCCGCTTTGCAGCCCGGCTGGTCTGCTCAGCCTGAGTTGACCCGAACGCTGGGCAATC‑GCTTTCTGGATGACTGTTCGGCATTGCTGCTGCCGGTGCCAAGCGCGATCATGCCCAGCACCACCAATTACCt                                       <  1:162701/149‑1 (MQ=255)                                                                                                                   ccAGGCGCAGTTGCCCGCTTTGCAGCCCGGCTGGTCTGCTCAGCCTGAGTTGACCCGAACGCTGGGCAATC‑GCTTTCTGGATGACTGTTCGGCATTGCTGCTGCCGGTGCCAAGCGCGATCATGCCCAGCACCACCAATTACCTGTTc                                   >  2:196133/1‑148 (MQ=255)                                                                                                                   ccAGGCGCAGTTGCCCGCTTTGCAGCCCGGCTGGTCTGCTCAGCCTGAGTTGACCCGAACGCTGGGCAATC‑GCTTTCTGGATGACTGTTCGGCATTGCTGCTGCCGGTGCCAAGCGCGATCATGCCCAGCACCACCAATTACCTGTTc                                   >  6:37502/1‑148 (MQ=255)                                                                                                                    cAGGCGCAGTTGCCCGCTTTGCAGCCCGGCTGGTCTGCTCAGCCTGAGTTGACCCGAACGCTGGGCAATC‑GCTTTCTGGATGACTGTTCGGCATTGCTGCTGCCGGTGCCAAGCGCGATCATGCCCAGCACCACCAAtt                                           <  8:284940/139‑1 (MQ=255)                                                                                                                    cAGGCGCAGTTGCCCGCTTTGCAGCCCGGCTGGTCTGCTCAGCCTGAGTTGACCCGAACGCTGGGCAATC‑GCTTTCTGGATGACTGTTCGGCATTGCTGCTGCCGGTGCCAAGCGCGATCATGCCCAGCACCACCAAtt                                           >  7:284940/1‑139 (MQ=255)                                                                                                                           aGTTGCCCGCTTTGCAGCCCGGCTGGTCTGCTCAGCCTGAGTTGACCCGAACGCTGGGCAATC‑GCTTTCTGGATGACTGTTCGGCATTGCTGCTGCCGGTGCCAAGCGCGATCATGCCCAGCACCACCAATTACCTGTTCAACCCACgg                          <  4:347040/149‑1 (MQ=255)                                                                                                                             ttGCCCGCTTTGCAGCCCGGCTGGTCTGCTCAGCCTGAGTTGACCCGAACGCTGGGCAATCTTCTTTCTGGATGACTGTTCGGCATTGCTGCTGCCGGTGCCAAGCGCGATCATGCCCAGCACCACCAATTACCTGTTCAACCCACGGc                         >  1:309389/1‑149 (MQ=255)                                                                                                                               gCCCGCTTTGCAGCCCGGCTGGTCTGCTCAGCCTGAGTTGACCCGAACGCTGGGCAATC‑GCTTTCTGGATGACTGTTCGGCATTGCTGCTGCCGGTGCCAAGCGCGATCATGCCCAGCACCACCAATTACCTGTTCAACCCACGGCAcc                      >  5:267289/1‑149 (MQ=255)                                                                                                                                cccGCTTTGCAGCCCGGCTGGTCTGCTCAGCCTGAGTTGACCCGAACGCTGGGCAATC‑GCTTTCTGGATGACTGTTCGGCATTGCTGCTGCCGGTGCCAAGCGCGATCATGCCCAGCACCACCAATTACCTGTTCACCCCACGGCAccc                     >  6:293457/1‑149 (MQ=255)                                                                                                                                cccGCTTTGCAGCCCGGCTGGTCTGCTCAGCCTGAGTTGACCCGAACGCTGGGCAATC‑GCTTTCTGGATGACTGTTCGGCATTGCTGCTGCCGGTGCCAAGCGCGATCATGCCCAGCACCACCAATTACCTGTTCAACCCACGGCAcc                      <  5:37502/148‑1 (MQ=255)                                                                                                                                cccGCTTTGCAGCCCGGCTGGTCTGCTCAGCCTGAGTTGACCCGAACGCTGGGCAATC‑GCTTTCTGGATGACTGTTCGGCATTGCTGCTGCCGGTGCCAAGCGCGATCATGCCCAGCACCACCAATTACCTGTTCAACCCACGGCAcc                      <  1:374106/148‑1 (MQ=255)                                                                                                                                 ccGCTTTGCAGCCCGGCTGGTCTGCTCAGCCTGAGTTGACCCGAACGCTGGGCAATC‑GCTTTCTGGATGACTGTTCGGCATTGCTGCTGCCGGTGCCAAGCGCGATCATGCCCAGCACCAACAATTACCTGTTCAACCAACGGCAcc                      >  2:374106/1‑147 (MQ=255)                                                                                                                                  cGCTTTGCAGCCCGGCTGGTCTGCTCAGCCTGAGTTGACCCGAACGCTGGGCAATC‑GCTTTCTGGATGACTGTTCGGCATTGCTGCTGCCGGTGCCAAGCGCGATCATGCCCAGCACCACCAATTACCt                                       <  5:146063/129‑1 (MQ=255)                                                                                                                                  cGCTTTGCAGCCCGGCTGGTCTGCTCAGCCTGAGTTGACCCGAACGCTGGGCAATC‑GCTTTCTGGATGACTGTTCGGCATTGCTGCTGCCGGTGCCAAGCGCGATCATGCCCAGCACCACCAATTACCt                                       >  6:146063/1‑129 (MQ=255)                                                                                                                                   gCTTTGCAGCCCGGCTGGTCTGCTCAGCCTGAGTTGACCCGAACGCTGGGCAATC‑GCTTTCTGGATGACTGTTCGGCATTGCTGCTGCCGGTGCCAAGCGCGATCATGCCCAGCACCACCAATTACCTGTTCAACCCACGGCACCCg                    >  2:103345/1‑147 (MQ=255)                                                                                                                                   gCTTTGCAGCCCGGCTGGTCTGCTCAGCCTGAGTTGACCCGAACGCTGGGCAATC‑GCTTTCTGGATGACTGTTCGGCATTGCTGCTGCCGGTGCCAAGCGCGATCATGCCCAGCACCACCAATTACCTGTTCAACCCACGGCACCCg                    >  8:337623/1‑147 (MQ=255)                                                                                                                                     tgagCAGACCAGCCGGGCTGCTCAGCCTGAGTTGACCCGAACGCTGGGCAATC‑GCTTTCTGGATGACTGTTCGGCATTGCTGCTGCCAGTGCCAAGCGCGATCATGCCCAGca                                                    <  6:227004/110‑1 (MQ=255)                                                                                                                                     tgagCAGACCAGCCGGGCTGCTCAGCCTGAGTTGACCCGAACGCTGGGCAATC‑GCTTTCTGGATGACTGTTCGGCATTGCTGCTGCCAGTGCCAAGCGCGATCATGCCCAGca                                                    <  6:232367/110‑1 (MQ=255)                                                                                                                                     tgagCAGACCAGCCGGGCTGCTCAGCCTGAGTTGACCCGAACGCTGGGCAATC‑GCTTTCTGGATGACTGTTCGGCATTGCTGCTGCCAGTGCCAAGCGCGATCATGCCCAGca                                                    >  5:227004/4‑113 (MQ=255)                                                                                                                                     tgagCAGACCAGCCGGGCTGCTCAGCCTGAGTTGACCCGAACGCTGGGCAATC‑GCTTTCTGGATGACTGTTCGGCATTGCTGCTGCCAGTGCCAAGCGCGATCATGCCCAGca                                                    >  5:232367/4‑113 (MQ=255)                                                                                                                                        gCAGCCCGGCTGGTCTGCTCAGCCTGAGTTGACCCGAACGCTGGGCAATC‑GCTTTCTGGATGACTGTTCGGCATTGCTGCTGCCGGTGCCAAGCGCGATCATGCCCAGCACCACCAATTACCTGTTCAACCCACGGCACCCGCAGgcg              >  1:111752/1‑148 (MQ=255)                                                                                                                                         cAGCCCGGCTGGTCTGCTCAGCCTGAGTTGACCCGAACGCTGGGCAATC‑GCTTTCTGGATGACTGTTCGGCATTGCTGCTGCCGGTGCCAAGCGCGATCATGCCCAGCACCACCAATTACCTGTTCAACCCACGGCACCCGCAGGCGCa            <  3:89375/149‑1 (MQ=255)                                                                                                                                           gCCCGGCTGGTCTGCTCAGCCTGAGTTGACCCGAACGCTGGGCAATC‑GCTTTCTGGATGACTGTTCGGCATTGCTGCTGCCGGTGCCAAGCGCGATCATGCCCAGCACCACCAATTACCTGTTCAACCCACGGCACCCGCAGGCGCaga          <  5:290548/149‑1 (MQ=255)                                                                                                                                            cccGGCTGGTCTGCTCAGCCTGAGTTGACCCGAACGCTGGGCAATC‑GCTTTCTGGATGACTGTTCGGCATTGCTGCTGCCGGTGCCAAGCGCGATCATGCCCAGcaccac                                                >  1:370033/1‑110 (MQ=255)                                                                                                                                            cccGGCTGGTCTGCTCAGCCTGAGTTGACCCGAACGCTGGGCAATC‑GCTTTCTGGATGACTGTTCGGCATTGCTGCTGCCGGTGCCAAGCGCGATCATGCCCAGcaccac                                                <  2:370033/110‑1 (MQ=255)                                                                                                                                              cGGCTGGTCTGCTCAGCCTGAGTTGACCCGAACGCTGGGCAATC‑GCTTTCTGGATGACTGTTCGGCATTGCTGCTGCCGGTGCCAAGCGCGATCATGCCCAGCACCACCAATTACCTGTTCAACCCACGGCACCCGCAGGCGCAGAgcg       <  1:196133/149‑1 (MQ=255)                                                                                                                                                   ggTCTGCTCAGCCTGAGTTGACCCGAACGATGGGCAATC‑GCTTTCTGGATGACAGTTCCGCATTGCTGCTGCCGGTGCCAAGCGAGATCATGCCCAGCACCAACAATAACCAGTTCAACCCACTGCACCCGCAGGCGCAGAGCGCGaa   >  6:291071/1‑148 (MQ=255)                                                                                                                                                    gTCTGCTCAGCCTGAGTTGACCCGAACGCTGGGCAATC‑GCTTTCTGGATGACTGTTCGGCATTGCTGCTGCCGGTGCCAAGCGCGATCATGCCCAGCACCACCAATTACCTGTTCAACCCACGGCACCCGCAGGCGCAGAGCGCGaaa  <  2:309389/148‑1 (MQ=255)                                                                                                                                                    gTCTGCTCAGCCTGAGTTGACCCGAACGCTGGGCAATC‑GCTTTCTGGATGACTGTTCGGCATTGCTGCTGCCGGTGCCAAGCGCGATCATGCCCAGCACCACCAATTACCTGTTCAACCCACGGCACCCGCAGGCGCAGAGCGCGaaa  >  4:121275/1‑148 (MQ=255)                                                                                                                                                     |                                                                                                                                                    CGCCACCAGCCCTCCTGGGGCGATGCTGGAAGTGTTGGTGCACCTGGAGATCGACCCCGAGGATTTCCCTACCACCATGCGCTTGCTTCGCATCGAACTGCCCGACACGGTTTCCCAGGCGCAGTTGCCCGCTTTGCAGCCCGGCTGGTCTGCTCAGCCTGAGTTGACCCGAACGCTGGGCAATC‑GCTTTCTGGATGACTGTTCGGCATTGCTGCTGCCGGTGCCAAGCGCGATCATGCCCAGCACCACCAATTACCTGTTCAACCCACGGCACCCGCAGGCGCAGAGCGCGAAA  >  NC_002947/2780804‑2781098

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