Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,388,037	G→T	7.6%	G151G (GGC→GGA)	PP_1209 ←	cold‑shock protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,388,037	0	G	T	7.6%	201.9 / 3.3	64	G151G (GGC→GGA)	PP_1209	cold‑shock protein
Reads supporting (aligned to +/- strand):  ref base G (32/27);  new base T (2/3);  total (34/30)
Fisher's exact test for biased strand distribution p-value = 6.59e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGTTGGTCAACGGCGGTTTACTGCAACCACGTCTTCGGCCTGCAGGCCCTTGTCGCGGTGCATCACCGAAAACTCCACGCGCTGGCCTTCGACCAGAATGCGGTGGCCCTCACCGCGGATGGCCCGAAAGTGCACGAAGATGTCATCGCCCGAGTCTCGGGAGATGAAACCAAAGC‑CTTTCGACGTGTTGAACCATTTGACGGTGCCGGTATCACGGTTGCCGGCATCTTGTGCCACGTGTTGAGTGCCGCGTGCCTTGCGCGTGCTGCGGGCAAAGCCCACGGCCAGATGCA  >  NC_002947/1387890‑1388182                                                                                                                                                    |                                                                                                                                                   cGTTGGTCAACGGCGGTTTACTGCAACCACGTCTTCGGCCTGCAGGCCCTTGTCGCGGTGCATCACCGAAAACTCCACGCGCTGGCCTTCGACCAGAATGCGGTGGCCCTCACCGCGGATGGCCCGAAAGTGCACGAAGATGTCATCGc                                                                                                                                                   >  4:258154/1‑149 (MQ=255)       tCAACGGCGGTTTACTGCAACCACGTCTTCGGCCTGCAGGCCCTTGTCGCGGTGCATCACCGAAAACTCCACGCGCTGGCCTTCGACCAGAATGCGGTGGCCCTCACCGCGGATGGCCCGAAAGTGCACGAAGATGTCATCGCCCGAGt                                                                                                                                             >  4:192442/1‑149 (MQ=255)       tCAACGGCGGTTTACTGCAACCACGTCTTCGGCCTGCAGGCCCTTGTCGCGGTGCATCACCGAAAACTCCACGCGCTGGCCTTCGACCAGAATGCGGTGGCCCTCACCGCGGATGGCCCGAAAGTGCACGAAGATGTCATCGCCCGAGt                                                                                                                                             >  2:346525/1‑149 (MQ=255)          aCGGCGGTTTACTGCAACCACGTCTTCGGCCTGCAGGCCCTTGTCGCGGTGCATCACCGAAAACTCCACGCGCTGGCCTTCGACCAGAATGCGGTGGCCCTCACCGCGGATGGCCCGAAAGTGCACGAAGATGTCATCTCCCGAGactc                                                                                                                                          >  8:52228/1‑149 (MQ=255)           cggcggTTTACTGCAACCACGTCTTCGGCCTGCAGGCCCTTGTCGCGGTGCATCACCGAAAACTCCGCGCGCTGGCCTTCGACCAGAATGCGGTGGCCCTCACCGCGGATGGCCCGAAAGTGCACGAAGATGTCATCGCCCGAGTCTCg                                                                                                                                         >  5:135942/1‑149 (MQ=255)            ggcggTTTACTGCAACCACGTCTTCGGCCTGCAGGCCCTTGTCGCGGTGCATCACCGAAAACTCCACGCGCTGGCCTTCGACCAGAATGCGGTGGCCCTCACCGCGGATGGCCCGAAAGTGCACGAAGATGTCATCGCCCGAGTCTCgg                                                                                                                                        >  7:168996/1‑149 (MQ=255)                 tttACTGCAACCACGTCTTCGGCCTGCAGGCCCTTGTCGCGGTGCATCACCGAAAACTCCACGCGCTGGCCTTCGACCAGAATGCGGTGGCCCTCACCGCGGATGGCCCGAAAGTGCACGAAGATGTCATCGCCCGAGTCTCGGGAGAt                                                                                                                                   <  7:251723/149‑1 (MQ=255)                     cTGCAACCACGTCTTCGGCCTGCAGGCCCTTGTCGCGGTGCATCACCGAAAACTCCACGCGCTGGCCTTCGACCAGAATGCGGTGGCCCTCACCGCGGATGGCCCGAAAGTGCACGAAGATGTCATCGCCCGAGTCTCGGGAGATGaa                                                                                                                                <  1:346525/148‑1 (MQ=255)                     cTGCAACCACGTCTTCGGCCTGCAGGCCCTTGTCGCGGTGCATCACCGAAAACTCCACGCGCTGGCCTTCGACCAGAATGCGGTGGCCCTCACCGCGGATGGCCCGAAAGTGCACGAAGATGTCATCGCCCGAGTCTCGGGAGATGaa                                                                                                                                <  3:192442/148‑1 (MQ=255)                      tGCAACCACGTCTTCGGCCTGCAGGCCCTTGTCGCGGTGCATCACCGAAAACTCCACGCGCTGGCCTTCGACCAGAATGCGGTGGCCCTCACCGCGGATGGCCCGAAAGTGCACGAAGATGTCATCGCCCGAGTCTCGGGAGATGAAAc                                                                                                                              >  3:141959/1‑149 (MQ=255)                          aCCACGTCTTCGGCCTGCAGGCCCTTGTCGCGGTGCATCACCGAAAACTCCACGCGCTGGCCTTCGACCAGAATGCGGTGGCCCTCACCGCGGATGGCCCGAAAGTGCACGACGATGTCATCGCCCGAGTCTCGGGAGATGAAACCaaa                                                                                                                          <  3:331034/149‑1 (MQ=255)                          aCCACGTCTTCGGCCTGCAGGCCCTTGTCGCGGTGCATCACCGAAAACTCCACGCGCTGGCCTTCGACCAGAATGCGGTGGCCCTCACCGCGGATGGCCCGAAAGTGCACGAAGATGTCATCTCCCGAGACTCGGGCGATGacatcttc                                                                                                                          <  7:52228/149‑8 (MQ=255)                          aCCACGTCTTCGGCCTGCAGGCCCTTGTCGCGGTGCATCACCGAAAACTCCACGCGCTGGCCTTCGACCAGAATGCGGTGGCCCTCACCGCGGATGGCCCGAAAGTGCACGAAGATGTCATCGCCCGAGTCTCGGGAGATGAAACCaaa                                                                                                                          <  3:258154/149‑1 (MQ=255)                          aCCACGTCTTCGGCCTGCAGGCCCTTGTCGCGGTGCATCACCGAAAACTCCACGCGCTGGCCTTCGACCAGAATGCGGTGGCCCTCACCGCGGATGGCCCGAAAGTGCACGAAGATGTCATCGCCCGAGTCTCGGGAGATGAAACCa                                                                                                                            >  5:143972/1‑147 (MQ=255)                                tCTTCGGCCTGCAGGCCCTTGTCGCGGTGCATCACCGAAAACTCCACGCGCTGGCCTTCGACCAGAATGCGGTGGCCCTCACCGCGGATGGCCCGAAAGTGCACGAAGATGTCATCGCCCGAGTCTCGGGAGATGAAACCAAAGC‑Cttt                                                                                                                   >  5:300670/1‑149 (MQ=255)                                 cTTCGGCCTGCAGGCCCTTGTCGCGGTGAATCACCGAAAACTCCACGCGCTGGCCTTCGACCAGAATGCGGTGGCCCTCACCGCAGATGGCCCGAAAGTGCACGAAGATGTCATCGCCCGAGTCTCGGGAGATGAAACCAAAGC‑CTTTc                                                                                                                  >  1:16118/1‑149 (MQ=255)                                        cTGCAGGCCCTTGTCGCGGTGCATCACCGAAAACTCCACGCGCTGGCCTTCGACCAGAATGCGGTGGCCCTCACCGCGGATGGCCCGAAAGTGCACGAAGATGTCATCGCCCGAGTCTCGGGAGATGAAACCAAAGC‑CTTTCGACgtgt                                                                                                           >  2:152855/1‑149 (MQ=255)                                         tGCAGGCCCTTGTCGCGGTGCATCACCGAAAACTCCACGCGCTGGCCTTCGACCAGAATGCGGTGGCCCTCACCGCGGATGGCCCGAAAGTGCACGAAGATGTCATCGCCCGAGTCTCGGGAGATGAAACCAAAGC‑CTTTCGACGTGtt                                                                                                          >  7:57956/1‑149 (MQ=255)                                           cAGGCCCTTGTCGCGGTGCATCACCGAAAACTCCACGCGCTGGCCTTCGACCAGAATGCGGTGGCCCTCACCGCGGATGGCCCGAAAGTGCACGAAGATGTCATCGCCCGAGTCTCGGGAGAt                                                                                                                                   <  7:3283/123‑1 (MQ=255)                                           cAGGCCCTTGTCGCGGTGCATCACCGAAAACTCCACGCGCTGGCCTTCGACCAGAATGCGGTGGCCCTCACCGCGGATGGCCCGAAAGTGCACGAAGATGTCATCGCCCGAGTCTCGGGAGAt                                                                                                                                   >  8:3283/1‑123 (MQ=255)                                                 cTTGTCGCGGTGCATCACCGAAAACTCCACGCGCTGGCCTTCGACCAGAATGCGGTGGCCCTCACCGCGGATGGCCCGAAAGTGCACGAAGATGTCATCGCCCGAGTCTCGGGAGATGAAACCAAAGC‑CTTTCGACGTGTTGAACCAtt                                                                                                  <  8:168996/149‑1 (MQ=255)                                                    gTCGCGGTGCATCACCGAAAACTCCACGCGCTGGCCTTCGACCAGAATGCGGTGGCCCTCACCGCGGATGGCCCGAAAGTGCACGAAGATGTCATCGCCCGAGTCTCGGGAGATGAAACCAAAGC‑CTTTCGACGTGTTGAACCAttt                                                                                                 <  7:262429/147‑1 (MQ=255)                                                    gTCGCGGTGCATCACCGAAAACTCCACGCGCTGGCCTTCGACCAGAATGCGGTGGCCCTCACCGCGGATGGCCCGAAAGTGCACGAAGATGTCATCGCCCGAGTCTCGGGAGATGAAACCAAAGC‑CTTTCGACGTGTTGAACCATTTGa                                                                                               >  7:158602/1‑149 (MQ=255)                                                     tCGCGGTGCATCACCGAAAACTCCACGCGCTGGCCTTCGACCAGAATGCGGTGGCCCTCACCGCGGATGGCCCGAAAGTGCACGAAGATGTCATCGCCCGAGTCTCGGGAGATGAAACCAAAGC‑CTTTCGACGTGTTGAACCATTTGa                                                                                               <  8:57956/148‑1 (MQ=255)                                                               tCACCGAAAACTCCACGCGCTGGCCTTCGACCAGAATGCGGTGGCCCTCACCGCGGATGGCCCGAAAGTGCACGAAGATGTCATCGCCCGAGTCTCGGGAGATGAAACCAAAGC‑CTTTCGAc                                                                                                               >  6:338024/1‑122 (MQ=255)                                                               tCACCGAAAACTCCACGCGCTGGCCTTCGACCAGAATGCGGTGGCCCTCACCGCGGATGGCCCGAAAGTGCACGAAGATGTCATCGCCCGAGTCTCGGGAGATGAAACCAAAGC‑CTTTCGAc                                                                                                               <  5:338024/122‑1 (MQ=255)                                                                  ccGAAAACTCCACGCGCTGGCCTTCGACCAGAATGCGGTGGCCCTCACCGCGGATGGCCCGAAAGTGCACGAAGATGTCATCGCCCGAGTCTCGGGAGATGAAACCAAAGC‑CTTTCGACGTGTTGaa                                                                                                       >  6:34711/1‑127 (MQ=255)                                                                  ccGAAAACTCCACGCGCTGGCCTTCGACCAGAATGCGGTGGCCCTCACCGCGGATGGCCCGAAAGTGCACGAAGATGTCATCGCCCGAGTCTCGGGAGATGAAACCAAAGC‑CTTTCGACGTGTTGAACCATTTGACGGTGCCGGTATCa                                                                                 >  3:42475/1‑149 (MQ=255)                                                                   cGAAAACTCCACGCGCTGGCCTTCGACCAGAATGCGGTGGCCCTCACCGCGGATGGCCCGAAAGTGCACGAAGATGTCATCGCCCGAGTCTCGGGAGATGAAACCAAAGC‑CTTTCGACGTGTTGAACCATTTGACGGTGCCGGTATCAc                                                                                >  8:156737/1‑149 (MQ=255)                                                                    gAAAACTCCACGCGCTGGCCTTCGACCAGAATGCGGTGGCCCTCACCGCGGATGGCCCGAAAGTGCACGAAGATGTCATCGCCCGAGTCTCGGGAGATGAAACCAAAGC‑CTTTCGACGTGTTGAACGATTTGACGGTGCCGGTATCACg                                                                               <  6:300670/149‑1 (MQ=255)                                                                            cccgcgCTGGCCTTCGACCAGAATGCGGTGGCCCTCACCGCGGATGGCCCGAAAGTGCACGAAGATGTCATCGCCCGAGTCTCGGGAGATGAAACCAAAGC‑CTTTCGACGTGTTGAACCATTTGACGGTGCCGGTATCACGGTTGcc                                                                         <  6:90995/145‑1 (MQ=255)                                                                            cACGCGCTGGCCTTCGACCAGAATGCGGTGGCCCTCACCGCGGATGGCCCGAAAGTGCACGAAGATGTCATCGCCCGAGTCTCGGGAGATGAAACCAAAGC‑CTTTCGACGTGTTGAACCATTTGACGGTGCCGGTATCACGGTTGCCg                                                                        <  8:158602/148‑1 (MQ=255)                                                                                       cTTCGACCAGAATGCGGTGGCCCTCACCGCGGATGGCCCGAAAGTGCACGAAGATGTCATCGCCCGAGTCTCGGGAGATGAAACCAAAGC‑CTTTCGACATGTTGAACCATTTGACGGTGCCGGTATCACGGTTGCCGGCATCTtgt                                                               <  6:10289/146‑1 (MQ=255)                                                                                       cTTCGACCAGAATGCGGTGGCCCTCACCGCGGATGGCCCGAAAGTGCACGAAGATGTCATCGCCCGAGTCTCGGGAGATGAAACCAAAGC‑CTTTCGACATGTTGAACCATTTGACGGTGCCGGTATCACGGTTGCCGGCATCTtgt                                                               >  5:10289/1‑146 (MQ=255)                                                                                        ttCGACCAGAATGCGGTGGCCCTCACCGCGGATGGCCCGAAAGTGCACGAAGATGTCATCTCCCGAGACTCGGGCGATGa                                                                                                                                 <  8:66379/80‑1 (MQ=38)                                                                                           gACCAGAATGCGGTGGCCCTCACCGCGGATGGCCCGAAAGTGCACGAAGATGTCATCGCCCGAGTCTCGGGAGATGAAACCAAAGC‑CTTTCGACGTGTTGAACCATTTGACGGTGCCGGTATCACGGTTGCCGGCATCTTGTGCCACgt                                                        <  4:141959/149‑1 (MQ=255)                                                                                            aCCAGAATGCGGTGGCCCTCACCGCGGATGGCCCGAAAGTGCACGAAGATGTCATCGCCCGAGTCTCGGGAGATGAAACCAAAGC‑CTTTCGACGTGTTGAACCATTTGACGGTGCCGGTATCACGGTTGCCGGCATc                                                                   >  4:7805/1‑137 (MQ=255)                                                                                            aCCAGAATGCGGTGGCCCTCACCGCGGATGGCCCGAAAGTGCACGAAGATGTCATCGCCCGAGTCTCGGGAGATGAAACCAAAGC‑CTTTCGACGTGTTGAACCATTTGACGGTGCCGGTATCACGGTTGCCGGCATc                                                                   <  3:2981/137‑1 (MQ=255)                                                                                            aCCAGAATGCGGTGGCCCTCACCGCGGATGGCCCGAAAGTGCACGAAGATGTCATCGCCCGAGTCTCGGGAGATGAAACCAAAGC‑CTTTCGACGTGTTGAACCATTTGACGGTGCCGGTATCACGGTTGCCGGCATc                                                                   >  4:2981/1‑137 (MQ=255)                                                                                            aCCAGAATGCGGTGGCCCTCACCGCGGATGGCCCGAAAGTGCACGAAGATGTCATCGCCCGAGTCTCGGGAGATGAAACCAAAGC‑CTTTCGACGTGTTGAACCATTTGACGGTGCCGGTATCACGGTTGCCGGCATc                                                                   <  3:7805/137‑1 (MQ=255)                                                                                            aCCAGAATGCGGTGGCCCTCACCGCGGATGGCCCGAAAGTGCACGAAGATGTCATCGCCCGAGTCTCGGGAGATGAAACCAAAGC‑CTTTCGACGTGTTGAACCATTTGACGGTGCCGGTATCACGGTTGCCGGCATCTTGTGCCACgtg                                                       >  1:161735/1‑149 (MQ=255)                                                                                                gAATGCGGTGGCCCTCACCGCAGATGGCCCGAAAGTGCACGAAGATGTCATCGCCCGAGTCTCGGGAGATGAAACCAAAGC‑CTTTCGACGTGTTGAACCATTTGACGGTGCCGGTATCACGGTTGCCGGCATCTTGTGCCACGTGTTGa                                                   <  2:16118/149‑1 (MQ=255)                                                                                                        tGGCCCTCACCGCGGATGGCCCGAAAGTGCACGAAGATGTCATCGCCCGAGTCTCGGGAGATGAAACCAAAGC‑CTTTCGACGTGTTGAACCATTTGACGGTGCCGGTATCACGGTTGCCGGCATCTTGTGCCACGTGTTGAGTGCcgcg                                           <  7:156737/149‑1 (MQ=255)                                                                                                         ggCCCTCACCGCGGATGGCCCGAAAGTGCACGAAGATGTCATCGCCCGAGTCTCGGGAGATGAAACCAAAGC‑CTTTCGACGTGTTGAACCATTTGACGGTGCCGGTATCACGGTTGCCGGCATCTTGTGCCACGTGTTGAGTGCCGCGt                                          <  1:152855/149‑1 (MQ=255)                                                                                                           cccTCACCGCGGATGGCCCGAAAGTGCACGAAGATGTCATCGCCCGAGTCTCGGGAGATGAAACCAAAGC‑CTTTCGACGTGTTGAACCATTTGACGGTGCCGGTATCACGGTTGCCGGCATCTtgt                                                               >  1:175360/1‑126 (MQ=255)                                                                                                           cccTCACCGCGGATGGCCCGAAAGTGCACGAAGATGTCATCGCCCGAGTCTCGGGAGATGAAACCAAAGC‑CTTTCGACGTGTTGAACCATTTGACGGTGCCGGTATCACGGTTGCCGGCATCTtgt                                                               <  2:175360/126‑1 (MQ=255)                                                                                                            ccTCACCGCGGATGGCCCGAAAGTGCACGAAGATGTCATCGCCCGAGTCTCGGGAGATGAAACCAAAGC‑CTTTCGACGTGTTGAACCATTTGACGGTGCCGGTATCACGGTTGCCGGCATCTTGTGCCACGTGTTGAGt                                                 >  3:330554/1‑139 (MQ=255)                                                                                                            ccTCACCGCGGATGGCCCGAAAGTGCACGAAGATGTCATCGCCCGAGTCTCGGGAGATGAAACCAAAGC‑CTTTCGACGTGTTGAACCATTTGACGGTGCCGGTATCACGGTTGCCGGCATCTTGTGCCACGTGTTGAGt                                                 <  4:330554/139‑1 (MQ=255)                                                                                                                aCCGCGGATGGCCCGAAAGTGCACGAAGATGTCATCGCCCGAGTCTCGGGAGATGAAACCAAAGC‑CTTTCGACGTCTTGAACCATTTGACGGTGCCGGTATCACGGTTGCCGGCATCTTGTGCCACGTGTTGAGTGCCGCGTGCCTTgc                                   <  5:352307/149‑1 (MQ=255)                                                                                                                 ccGCGGATGGCCCGAAAGTGCACGAAGATGTCATCGCCCGAGTCTCGGGAGATGAAACCAAAGC‑CTTTCGACGTGTTGAACCATTTGACGGTGCCg                                                                                       >  1:199490/1‑96 (MQ=255)                                                                                                                 ccGCGGATGGCCCGAAAGTGCACGAAGATGTCATCGCCCGAGTCTCGGGAGATGAAACCAAAGC‑CTTTCGACGTGTTGAACCATTTGACGGTGCCg                                                                                       <  2:199490/96‑1 (MQ=255)                                                                                                                      gATGGCCCGAAAGTGCACGAAGATGTCATCGCCCGAGTCTCGGGAGATGAAACCAAAGC‑CTTTCGACGTGTTGAACCATTTGACGGTGCCGGTATCACGGTTGCCGGCATCTTGTGCCACGTGTTGAGTGCCGCGTGCCTTGCGCGtgc                             <  4:42475/149‑1 (MQ=255)                                                                                                                           cccGAAAGTGCACGAAGATGTCATCGCCCGAGTCTCGGGAGATGAAACCAAAGC‑CTTTCGACGTGTTGAACCATTTGAGGGTGCCGGTATCACGGTTGCCGGCATCTTGTGCCACGTGTTGAGt                                                 <  1:245754/124‑1 (MQ=255)                                                                                                                           cccGAAAGTGCACGAAGATGTCATCGCCCGAGTCTCGGGAGATGAAACCAAAGC‑CTTTCGACGTGTTGAACCATTTGACGGTGCCGGTATCACGGTTGCCGGCATCTTGTGCCACGTGTTGAGt                                                 >  2:245754/1‑124 (MQ=255)                                                                                                                           cccGAAAGTGCACGAAGATGTCATCGCCCGAGTCTCGGGAGATGAAACCAAAGC‑CTTTCGACGTGTTGAACCATTTGACGGTGCCGGTATCACGGTTGCCGGCATCTTGTGCCACGTGTTGAGTGCCGCGTGCCTTgcgc                                 >  2:246142/1‑140 (MQ=255)                                                                                                                           cccGAAAGTGCACGAAGATGTCATCGCCCGAGTCTCGGGAGATGAAACCAAAGC‑CTTTCGACGTGTTGAACCATTTGACGGTGCCGGTATCACGGTTGCCGGCATCTTGTGCCACGTGTTGAGTGCCGCGTGCCTTGCGCGTGCTGCgg                        >  8:192701/1‑149 (MQ=255)                                                                                                                           cccGAAAGTGCACGAAGATGTCATCGCCCGAGTCTCGGGAGATGAAACCAAAGC‑CTTTCGACGTGTTGAACCATTTGACGGTGCCGGTATCACGGTTGCCGGCATCTTGTGCCACGTGTTGAGTGCCGCGTGCCTTGCGCGTGCTGCgg                        >  6:316055/1‑149 (MQ=255)                                                                                                                           cccGAAAGTGCACGAAGATGTCATCGCCCGAGTCTCGGGAGATGAAACCAAAGC‑CTTTCGACGTGTTGAACCATTTGACGGTGCCGGTATCACGCTTGCCGGCATCTTGTGCCACGTGTTGAGTGCCGCGTGCCTTgcgc                                 <  1:246142/140‑1 (MQ=255)                                                                                                                                     cACGAAGATGTCATCGCCCGAGTCTCGGGAGATGAAACCAAAGCGCTTTCGACGTGTTGAACCATTTGACGGTGCCGGTATCACGGTTGCCGGCATCTTGTGCCACGTGTTGAGTGCCGCGTGCCTTGCGCGTGCTGCGGGCAAAGcc                >  1:73744/1‑148 (MQ=255)                                                                                                                                            gtttCATCTCCCGAGACTCGGGCGATGAAACCAAAGC‑CTTTCGACGTGTTGAACCATTTGACGGTGCCGGTATCACGGTTGCCGGCATCTTGTGc                                                             >  1:249963/4‑95 (MQ=38)                                                                                                                                            gtttCATCTCCCGAGACTCGGGCGATGAAACCAAAGC‑CTTTCGACGTGTTGAACCATTTGACGGTGCCGGTATCACGGTTGCCGGCATCTTGTGc                                                             <  2:249963/92‑1 (MQ=38)                                                                                                                                            aTGTCATCGCCCGAGTCTCGGGAGATGAAACCAAAGC‑CTTTCGACGTGTTGAACCATTTGACGGTGCCGGTATCACGGTTGCCGGCATCTTGTGCCACGTGTTGAGTGCCGCGTGCCTTGCGCGTGCTGCGGGCAAAGCCCACGGc          >  3:159284/1‑146 (MQ=255)                                                                                                                                               tCATCGCCCGAGTCTCGGGAGATGAAACCAAAGC‑CTTTCGACGTGTTGAACCATTTGACGGTGCCGGTATCACGGTTGCCGGCATCTTGTGCCACGTGTTGAGTGCCGCGTGCCTTGCGCGTGCTGCGGGCAAAGCCCACGGCCAGATg    >  4:177524/1‑149 (MQ=255)                                                                                                                                                 aTCGCCCGAGTCTCGGGAGATGAAACCAAAGC‑CTTTCGACGTGTTGAACCATTTGACGGTGCCGGTATCACGGTTGCCGGCATCTTGTGCCACGTGTTGAGTGCCGCGTGCCTTGCGCGTGCTGCGGGCAAAGCCCACGGCCAGATGca  >  8:332948/1‑149 (MQ=255)                                                                                                                                                    |                                                                                                                                                   CGTTGGTCAACGGCGGTTTACTGCAACCACGTCTTCGGCCTGCAGGCCCTTGTCGCGGTGCATCACCGAAAACTCCACGCGCTGGCCTTCGACCAGAATGCGGTGGCCCTCACCGCGGATGGCCCGAAAGTGCACGAAGATGTCATCGCCCGAGTCTCGGGAGATGAAACCAAAGC‑CTTTCGACGTGTTGAACCATTTGACGGTGCCGGTATCACGGTTGCCGGCATCTTGTGCCACGTGTTGAGTGCCGCGTGCCTTGCGCGTGCTGCGGGCAAAGCCCACGGCCAGATGCA  >  NC_002947/1387890‑1388182

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