Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	396,848	Δ1 bp	6.4%	coding (98/924 nt)	PP_0329 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	396,848	0	C	.	6.4%	200.3 / 7.0	62	coding (98/924 nt)	PP_0329	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (30/28);  new base . (2/2);  total (32/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.65e-01

GTCGTTGTATCCAAAGGGGGTGTAGGTGTGATTTTTTTGCTTTCCGTCAGGCGTGCTCACGACAAGCACTGAGTCCATCTCATCTGTTGCGAAAAGCTGGGAACTTCCTGCTTCTTTGCCAGCGGACTCAGCAAGAGGGAGCGTTCCAATCCGCAGGATT‑GAACGCTTCTGATCGTTTCCAATGACAGTGGTCAGTTTGCTGGTCTGGTAGAAATATTTGGGGGTATCAGGACTGTTTTTCACCGACGTGCTCCCGTTGTCTGCAGTTGGCCCCGGAAATAGATGGCCTGGGC  >  NC_002947/396703‑396995                                                                                                                                                  |                                                                                                                                                     gTCGTTGTATCCAAAGGGGGTGTAGGTGTGATTTTTTTGCTTTCCGTCAGGCGTGCTCACGACAAGCACTGAGTCCATCTCATCTGTTGCGAAAAGCTGGGAACTTCCTGCTTCTTTGCCAGCGGACTCAGCAAGAGGGAGCGTTc                                                                                                                                                      >  4:207947/1‑146 (MQ=255) gTCGTTGTATCCAAAGGGGGTGTAGGTGTGATTTTTTTGCTTTCCGTCAGGCGTGCTCACGACAAGCACTGAGTCCATCTCATCTGTTGCGAAAAGCTGGGAACTTCCTGCTTCTTTGCCAGCGGACTCAGCAAGAGGGAGCGTTc                                                                                                                                                      <  3:207947/146‑1 (MQ=255)  tCGTTGTATCCAAAGGGGGTGGAGGTGTGATTTTTTTGCTTTCCGTCAGGCGTGCTCACGACAAGCACTGAGTCCATCTCATCTGTTGCGAAAAGCTGGGAACTTCCTGCTTCTTTGCCAGCGGACTCAGCAAGAGGGAGCGTTCCAAt                                                                                                                                                  <  2:17363/149‑1 (MQ=255)             aaaGGGGGTGTAGGTGTGATTTTTTTGCTTTCCGTCAGGCGTGCTCACGACAAGCACTGAGTCCATCTCATCTGTTGCGAAAAGCTGGGAACTTCCTGCTTCTTTGCCAGCGGACTCAGCAAGAGGGAGCGTTCCAATCCGCAGGAtt                                                                                                                                        >  2:215475/1‑148 (MQ=255)               aGGGGGTGTAGGTGTGATTTTTTTGCTTTCCGTCAGGCGTGCTCACGACAAGCACTGAGTCCATCTCATCTGTTGCGAAAAGCTGGGAACTTCCTGCTTCTTTGCCAGCGGACTCAGCAAGAGGGAGCGTTCCAATCCGCAGGATT‑Gaa                                                                                                                                    >  6:283884/1‑149 (MQ=255)                          gtgtGATTTTTTTGCTTTCCGTCAGGCGTGCTCACGACAAGCACTGAGTCCATCTCATCTGTTGCGAAAAGCTGGGAACTTCCTGCTTCTTTGCCAGCGGACTCAGCAAGAGGGAGCGTTCCAATCCGCAGGATT‑GAACGCTTCTGATc                                                                                                                         <  4:115861/149‑1 (MQ=255)                              gATTTTTTTGCTTTCCGTCAGGCGTGCTCACGACAAGCACTGAGTCCATCTCATCTGTTGCGAAAAGCTGGGAACTTCCTGCTTCTTTGCCAGCGGACTCAGCAAGAGGGAGCGTTCCAATCCGCAGGATT‑GAACGCTTCTGATCGtt                                                                                                                      >  4:199858/1‑148 (MQ=255)                               aTTTTTTTGCTTTCCGTCAGGCGTGCTCACGACAAGCACTGAGTCCATCTCATCTGTTGCGAAAAGCTGGGAACTTCCTGCTTCTTTGCCAGCGGACTCAGCAAGAGGGAGCGTTCCAATCCGCAGGATT‑GAACGCTTCTGATCGttt                                                                                                                     <  4:185646/148‑1 (MQ=255)                                            ccGTCAGGCGTGCTCACGACAAGCACTGAGTCCATCTCATCTGTTGCGAAAAGCTGGGAACTTCCTGCTTCTTTGCCAGCGGACTCAGCAAGAGGCAGCGTTCCAATCCGCAGGATT‑GAACGCTTCTGATCGTTTCCAATGACAGTgg                                                                                                        >  3:215516/1‑148 (MQ=255)                                             cGTCAGGCGTGCTCACGACAAGCACTGAGTCCATCTCATCTGTTGCGAAAAGCTGGGAACTTCCTGCTTCTTTGCCAGCGGACTCAGCAAGAGGGAGCGTTCCAATCCGCAGGATT‑GAACGCTTCTGATCGTTTCCAATGACAGTGGt                                                                                                       <  1:163652/148‑1 (MQ=255)                                              gTCAGGCGTGCTCACGACAAGCACTGAGTCCATCTCATCTGTTGCGAAAAGCTGGGAACTTCCTGCTTCTTTGCCAGCGGACTCAGCAAGAGGGAGCGTTCCAATCCGCAGGATT‑GAACGCTTCTGATCGTTTCCAATGACAGTGGt                                                                                                       >  2:221384/1‑147 (MQ=255)                                                    cGTGCTCACGACAAGCACTGAGTCCATCTCATCTGTTGCGAAAAGCTGGGAACTTCCTGCTTCTTTGCCAGCGGACTCAGCAAGAGGGAGCGTTCCAATCCGCAGGATT‑GAACGCTTCTGATCGTTTCCAATGACAGTGGTCAGTTTGc                                                                                               <  5:283884/149‑1 (MQ=255)                                                      tGCTCACGACAAGCACTGAGTCCATCTCATCTGTTGCGAAAAGCTGGGAACTTCCTGCTTCTTTGCCAGCGGACTCAGCAAGAGGGAGCGTTCCAATCCGCAGGATT‑GAACGCTTCTGATCGTTTCCAATGACAGTGGTCAGTTTGctg                                                                                             <  1:215475/149‑1 (MQ=255)                                                        cTCACGACAAGCACTGAGTCCATCTCATCTGTTGCGAAAAGCTGGGAACTTCCGGCTTCTTTGCCAGCGGACTCAGCAAGAGGGAGCGTTCCAATCCGCAGGATT‑GAACGCTTCTGATCGTTTCCAATGACAGTGGTCAGTTTGctggt                                                                                           >  1:245638/1‑149 (MQ=255)                                                           aCGACAAGCACTGAGTCCATCTCATCTGTTGCGAAAAGCTGGGAACTTCCTGCTTCTTTGCCAGCGGACTCAGCAAGAGGGAGCGTTCCAATCCGCAGGATT‑GAACGCTTCTGATCGTTTCCAATGACAGTGGTCa                                                                                                     >  2:2504/1‑136 (MQ=255)                                                             gACAAGCACTGAGTCCATCTCATCTGTTGCGAAAAGCTGGGAACTTCCTGCTTCTTTGCCAGCGGACTCAGCAAGAGGGAGCGTTCCAATCCGCAGGATT‑GAACGCTTCTGATCGtt                                                                                                                      >  2:290517/1‑117 (MQ=255)                                                                aaGCACTGAGTCCATCTCATCTGTTGCGAAAAGCTGGGAACTTCCTGCTTCTTTGCCAGCGGACTCAGCAAGAGGGAGCGTTCCAATCCGCAGGATT‑GAACGCTTCTGATCGTTTCCAATGACAGTGGTCAGTTTGCTGGTCTGGTAg                                                                                    >  8:2230/1‑148 (MQ=255)                                                                  gCACTGAGTCCATCTCATCTGTTGCGAAAAGCTGGGAACTTCCTGCTTCTTTGCCAGCGGACTCAGCAAGAGGGAGCGTTCCAATCCGCAGGATT‑GAACGCTTCTGATCGTTTCCAATGACAGTGGTCAGTTTGCTGGTCTGGTAGaaa                                                                                 >  4:287943/1‑149 (MQ=255)                                                                         gTCCATCTCATCTGTTGCGAAAAGCTGGGAACTTCCTGCTTCTTTGCCAGCGGACTCAGCAATAGTGAGCGTTCCAATCCGGAGGATT‑GAACGCTTCTGATCGTTTCCAATGACAGTGGTCAGTTTGCTGGTCTGGTAGAAATATTTgg                                                                          <  4:311496/149‑1 (MQ=255)                                                                            catctcatctGTTGCGAAAAGCTGGGAACTTCCTGCTTCTTTGCCAGCGGACTCAGCAAGAGGGAGCGTTCCAATCCGCAGGATT‑GAACGCTTCTGATCGTTTCCAATGACAGTGGTCAGTTTGCTGGTCTGGTAGAAATATTTgggg                                                                        <  5:149343/148‑1 (MQ=255)                                                                               ctcatctGTTGCGAAAAGCTGGGAACTTCCTGCTTCTTTGCCAGCGGACTCAGCAAGAGGGAGCGTTCCAATCCGCAGGATT‑GAACGCTTCTGATCGTTTCCAATGACAGTGGTCAGt                                                                                                   >  4:222419/1‑118 (MQ=255)                                                                               ctcatctGTTGCGAAAAGCTGGGAACTTCCTGCTTCTTTGCCAGCGGACTCAGCAAGAGGGAGCGTTCCAATCCGCAGGATT‑GAACGCTTCTGATCGTTTCCAATGACAGTGGTCAGt                                                                                                   <  3:222419/118‑1 (MQ=255)                                                                                    ctGTTGCGAAAAGCTGGGAACTTCCTGCTTCTTTGCCAGCGGACTCAGCAAGAGGGAGCGTTCCAATCCGCAGGATT‑GAACGCTTCTGATCGTTTCCAATGACAGTGGTCa                                                                                                     >  2:201283/1‑111 (MQ=255)                                                                                    ctGTTGCGAAAAGCTGGGAACTTCCTGCTTCTTTGCCAGCGGACTCAGCAAGAGGGAGCGTTCCAATCCGCAGGATT‑GAACGCTTCTGATCGTTTCCAATGACAGTGGTCa                                                                                                     <  1:201283/111‑1 (MQ=255)                                                                                       ttGCGAAAAGCTGGGAACTTCCTGCTTCTTTGCCAGCGGACTCAGCAAGAGGGAGCGTTCCAATCCGCAGGATT‑GAACGCTTCTGATCGTTTCCAATGACAGTGGTCAGTTTGCTGGTCTGGTAg                                                                                    >  5:215126/1‑125 (MQ=255)                                                                                       ttGCGAAAAGCTGGGAACTTCCTGCTTCTTTGCCAGCGGACTCAGCAAGAGGGAGCGTTCCAATCCGCAGGATT‑GAACGCTTCTGATCGTTTCCAATGACAGTGGTCAGTTTGCTGGTCTGGTAg                                                                                    <  6:215126/125‑1 (MQ=255)                                                                                       ttGCGAAAAGCTGGGAACTTCCTGCTTCTTTGCCAGCGGACTCAGCAAGAGGGAGCGTTCCAATCCGCAGGATT‑GAACGCTTCTGATCGTTTCCAATGACAGTGGTCAGTTTGCTGGTCTGGTAGAAATATTTGGGGGTATCAGGa                                                               <  3:199858/146‑1 (MQ=255)                                                                                       ttGCGAAAAGCTGGGAACTTCCTGCTTCTTTGCCAGCGGACTCAGCAAGAGGGAGCGTTCCAATCCGCAGGATT‑GAACGCTTCTGATCGTTTCCAATGACAGTGGTCAGTTTGCTGGTCTGGTAGAAATATTTGGGGGTATCAGGACTg                                                            >  8:251703/1‑149 (MQ=255)                                                                                       ttGCGAAAAGCTGGGAACTTCCTGCTTCTTTGCCAGCGGACTCAGCAAGAGGGAGCGTTCCAATCCGCAGGATT‑GAACGCTTCTGATCGTTTCCAATGACAGTGGTCAGTTTGCTGGTCTGGTAGAAATATTTGGGGGTATCAGGACTg                                                            <  7:2230/149‑1 (MQ=255)                                                                                             aaaGCTGGGAACTTCCTGCTTCTTTGCCAGCGGACTCAGCAAGAGGGAGCGTTCCAATCCGCAGGATT‑GAACGCTTCTGATCGTTTCCAATGACAGTGGTCAGTTTGCTGGTCTGGTAGAAATATTTGGGGGTATCAg                                                                 <  1:174440/138‑1 (MQ=255)                                                                                             aaaGCTGGGAACTTCCTGCTTCTTTGCCAGCGGACTCAGCAAGAGGGAGCGTTCCAATCCGCAGGATT‑GAACGCTTCTGATCGTTTCCAATGACAGTGGTCAGTTTGCTGGTCTGGTAGAAATATTTGGGGGTATCAg                                                                 >  2:174440/1‑138 (MQ=255)                                                                                              aaGCTGGGAACTTCCTGCTTCTTTGCCAGCGGACTCAGCAAGAGGGAGCGTTCCAATCCGCAGGATT‑GAACGCTTCTGATCGTTTCCAATGACAGTGGTCAGTTTGCTGGTCTGGTAGAAATATTTGGGGGTATCAGGACTGTTTTTCa                                                     <  3:287943/149‑1 (MQ=255)                                                                                                  tGGGAACTTCCTGCTTCTTTGCCAGCGGACTCAGCAAGAGGGAGCGTTCCAATCCGCAGGATT‑GAACGCTTCTGATCGTTTCCAATGACAGTGGTCAGTTTGctg                                                                                             >  8:97986/1‑105 (MQ=255)                                                                                                  tGGGAACTTCCTGCTTCTTTGCCAGCGGACTCAGCAAGAGGGAGCGTTCCAATCCGCAGGATT‑GAACGCTTCTGATCGTTTCCAATGACAGTGGTCAGTTTGctg                                                                                             <  7:97986/105‑1 (MQ=255)                                                                                                    ggAACTTCCTGCTTCTTTGCCAGCGGACTCAGCAAGAGGGAGCGTTCCAATCCGCAGGATT‑GAACGCTTCTGATCGTTTCCAATGACAGTGGTCAGTTTGCTGGTCTGGTAGAAATATTTGGGGGTATCAGGACTGTTTTTCACCGAc                                                <  1:221384/148‑1 (MQ=255)                                                                                                    ggAACTTCCTGCTTCTTTGCCAGCGGACTCAGCAAGAGGGAGCGTTCCAATCCGCAGGATT‑GAACGCTTCTGATCGTTTCCAATGACAGTGGTCAGTTTGCTGGTCTGGTAGAAATATTTGGGGGTATCAGGACTGTTTTTCACCGAc                                                <  7:251703/148‑1 (MQ=255)                                                                                                      aaCTTCCTGCTTCTTTGCCAGCGGACTCAGCAAGAGGGAGCGTTCCAATCCGCAGGATT‑GAACGCTTCTGATCGTTTCCAATGACAGTGGTCAGTTTGCTGGTCTGGTAGAAATATTTGGGGGTATCAGGACTGTTTTTCACCGACGTg                                             >  3:268545/1‑149 (MQ=255)                                                                                                         ttCCTGCTTCTTTGCCAGCGGACTCAGCAAGAGGGAGCGTTCCAATCCGCAGGATT‑GAACGCTTCTGATCGTTTCCAATGACAGTGGTCAGTTTGc                                                                                               >  8:204146/1‑96 (MQ=255)                                                                                                         ttCCTGCTTCTTTGCCAGCGGACTCAGCAAGAGGGAGCGTTCCAATCCGCAGGATT‑GAACGCTTCTGATCGTTTCCAATGACAGTGGTCAGTTTGc                                                                                               <  7:204146/96‑1 (MQ=255)                                                                                                           ccTGCTTCTTTGCCAGCGGACTCAGCAAGAGGGAGCGTTCCAATCCGCAGGATT‑GAACGCTTCTGATCGTTTCCAATGACAGTGGTCAGTTTGCTGGTCTGGTAGaa                                                                                  >  1:293238/1‑107 (MQ=255)                                                                                                           ccTGCTTCTTTGCCAGCGGACTCAGCAAGAGGGAGCGTTCCAATCCGCAGGATT‑GAACGCTTCTGATCGTTTCCAATGACAGTGGTCAGTTTGCTGGTCTGGTAGaa                                                                                  <  2:293238/107‑1 (MQ=255)                                                                                                           ccTGCTTCTTTGCCAGCGGACTCAGCAAGAGGGAGCGTTCCAATCAGCAGGATT‑GACCGCTTCTGATCGTTTCCAATGACAGTGGTCAGTTTGCTGGTCTGGTAGAAATATTTGGGGGTATCAGGACTGTTTTTCCCCGACGTGCTccc                                        >  6:261886/1‑149 (MQ=255)                                                                                                             tGCTTCTTTGCCAGCGGACTCAGCAAGAGGGAGCGTTCCAATCCGCAGGATT‑GAACGCTTCTGATCGTTTCCAATGACAGTGGTCAGTTTGCTGGTCTGGTAGAAATATTTGGGGGTATCAGGACTGTTTTTCACCGACGTGCTCCCGt                                      >  5:187968/1‑149 (MQ=255)                                                                                                             tGCTTCTTTGCCAGCGGACTCAGCAAGAGGGAGCGTTCCAATCCGCAGGATT‑GAACGCTTCTGATCGTTTCCAATGACAGTGGGCAGTTTGGTGGTCTGGTAGAAATATTTGGGGGTATCAGGAATGTTTTTCACCGCCGTGGTcaggt                                      >  4:20511/1‑145 (MQ=255)                                                                                                               cttcttTGCCATCGGACTCAGCAAGAGGGATCTTTCCAATCCGCAGGATT‑GAACGCTTCTGATCGTTTCCAATGACAGTGGTCAGTTTGCTGGTCTTGTAGAAATATTTGGGGGTATCAGGACTGTTTTTCACCGACGTGCTCCCGtt                                     <  6:187968/148‑1 (MQ=255)                                                                                                               cttcttTGCCAGCGGACTCAGCAAGAGGGAGCGTTCCATTCCGCAGGATT‑GAACGCTTCTGATCGTTTCCAATGACAGTGGTCAGTTTGCTGGTCTGGTAGAAATATTTGGGGGTATCAGGACTGTTTTTCACCGACGTGCTCCCGtt                                     <  6:81518/148‑1 (MQ=255)                                                                                                                     tGCCAGCGGACTCAGCAAGAGGGAGCGTTCCAATCCGCAGGATT‑GAACGCTTCTGATCGTTTCCAATGACAGTGGTCAGTTTGctggtctg                                                                                        <  6:108503/91‑1 (MQ=255)                                                                                                                     tGCCAGCGGACTCAGCAAGAGGGAGCGTTCCAATCCGCAGGATT‑GAACGCTTCTGATCGTTTCCAATGACAGTGGTCAGTTTGctggtctg                                                                                        >  5:108503/1‑91 (MQ=255)                                                                                                                       ccAGCGGACTCAGCAAGAGGGAGCGTTCCAATCCGCAGGATT‑GAACGCTTCTGATCGTTTCCAATGACAGTGGTCAGTTTGCTGGTCTGGTAGAAATATTTggg                                                                         >  2:309490/1‑104 (MQ=255)                                                                                                                       ccAGCGGACTCAGCAAGAGGGAGCGTTCCAATCCGCAGGATT‑GAACGCTTCTGATCGTTTCCAATGACAGTGGTCAGTTTGCTGGTCTGGTAGAAATATTTggg                                                                         <  1:309490/104‑1 (MQ=255)                                                                                                                        cAGCGGACTCAGCAAGAGGGAGCGTTCCAATCCGCAGGATT‑GAACGCTTCTGATCGTTTCCAATGACAGTGGTCAGTTTGCTGGTCTGGTAGAAATATTTGGGGGTATCAGGACTGTTTTTCACCGACGTGCTCCCGTTGTCTGCAGtt                           <  4:215516/149‑1 (MQ=255)                                                                                                                          gCGGACTCAGCAAGAGGGAGCGTTCCAATCCGCAGGATT‑GAACGCTTCTGATCGTTTCCAATGACAGTGGTCAGTTTGCTGGTCTGGTAGAAATATTTGGGGGTATCAGGACTGTTTTTCACCGACGTGCTCCCGTTGTCTGCAGTTgg                         >  3:77929/1‑149 (MQ=255)                                                                                                                                   gCAAGAGGGAGCGTTCCAATCCGCAGGATT‑GAACGCTTCTGATCGTTTCCAATGACAGTGGTCAGTTTGCTGGTCTGGTAGAAATATTTGGGGGTATCAGGACTGTTTTTCACCGACGTGGTCCCGTTGTCTGCAGTTGGCCCCGGaaa                >  1:280136/1‑149 (MQ=255)                                                                                                                                   gCAAGAGGGAGCGTTCCAATCCGCAGGATT‑GAACGCTTCTGATCGTTTCCAATGACAGTGGTCAGTTTGCTGGTCTGGTAGAAATATTTGGGGGTATCAGGACTGTTTTTCACCGACGTGCTCCCGTTGTCTGCAGTTGGCCCCGGaaa                <  5:68077/149‑1 (MQ=255)                                                                                                                                    cAAGATGGAGCGTTCAAATCCGCATGATT‑GAACGCTTCTGATCGTTTCCAATGACAGTGGTCAGTTTGCTGGTCTGGTAGAAATATTTGGGGGTATCAGGACTGTTTTTCACCGACGTGCTCCCGTTGTCTg                                <  7:125434/132‑1 (MQ=255)                                                                                                                                    cAAGAGGGAGCGTTCCAATCCGCAGGATT‑GAACGCTTCTGATCGTTTCCAATGACAGTGGTCAGTTTGCTGGTATGGTATAAATATTTGTGGGTATCAGTACTGTTTTTCACCGACGTTCTCCCGTTGTCTg                                >  8:125434/1‑132 (MQ=255)                                                                                                                                       gagGGATCGTTCCAATCCGCCGGATT‑GAACGCTTCTGATCGTTTCCAATGCCAGTGGTCAGTTTGCTGGTCTGGTTGAAATATTTGGGGGTATCAGGACTGTTTTTCACCGACGTGCTCCCGTTTTATGCAGTTGGCCCCGGAAATAGa            <  7:48810/149‑1 (MQ=255)                                                                                                                                        cagaaGCGTT‑CAATCCTGCGGATTGGAACGCTTCTGATCGTTTCCAATGACAGTGGTCAGTTTGCTGGTCTGGTAGAAATATTTGGGGGTATCAGGACTGTTTTTCa                                                     >  4:297252/5‑107 (MQ=37)                                                                                                                                        cagaaGCGTT‑CAATCCTGCGGATTGGAACGCTTCTGATCGTTTCCAATGACAGTGGTCAGTTTGCTGGTCTGGTAGAAATATTTGGGGGTATCAGGACTGTTTTTCa                                                     <  3:297252/103‑1 (MQ=37)                                                                                                                                          ggAGCGTTCCAATCCGCAGGATT‑GAACGCTTCTGATCGTTTCCAATGACAGTGGTCAGTTTGCTGGTCTGGTAGAAATATTTGGGGGTATCAGGACTGTTTTTCACCGACGTGCTCCCGTTGTCTGCAGTTGGCCCCGGAAATAGATgg         >  5:16322/1‑149 (MQ=255)                                                                                                                                          ggAGCGTTCCAATCCGCAGGATT‑GAACGCTTCTGATCGTTTCCAATGACAGTGGTCAGTTTGCTGGTCTGGTAGAAATATTTGGGGGTATCAGGACTGTTTTTCACCGACGTGCTCCCGTTGTCTGCAGTTGGCCCCGGAAATAGATgg         >  6:196248/1‑149 (MQ=255)                                                                                                                                           aaGCGTT‑CAATCCTGCGGATTGGAACGCTTCTGATCGTTTCCAATGACAGTGGTCAGTTTGCTGGTCTGGTAGAAATATTTGGGGGTATCAGGACTGTTTTTCACCGACGTGCTCCCGTTGTCTGCAGTTgg                         >  3:280657/2‑132 (MQ=255)                                                                                                                                           aaGCGTT‑CAATCCTGCGGATTGGAACGCTTCTGATCGTTTCCAATGACAGTGGTCAGTTTGCTGGTCTGGTAGAAATATTTGGGGGTATCAGGACTGTTTTTCACCGACGTGCTCCCGTTGTCTGCAGTTgg                         <  4:280657/131‑1 (MQ=255)                                                                                                                                               gTTCCAATCCGCAGGATT‑GAACGCTTCTGATCGTTTCCAATGACAGTGGTCAGTTTGCTGGTCTGGTAGAAATATTTGGGGGTATCAGGACTGTTTTTCACCGACGTGCTCCCGTTGTCAGCAGTTGGCCCCGCAAATAGATGGCCTgg    >  7:289732/1‑149 (MQ=255)                                                                                                                                                  ccAATCCGCAGGATT‑GAACGCTTCTGATCGTTTCCAATGACAGTGGTCAGTTTGCTGGTCTGGTAGAAATATTTGGGGGTATCAGGACTGTTTTTCACAGACGTGCTCCAGTTGTCTGCAGTTGGCCCCTGAAATAGATGGCATGGGc  >  2:216847/1‑148 (MQ=255)                                                                                                                                                  |                                                                                                                                                     GTCGTTGTATCCAAAGGGGGTGTAGGTGTGATTTTTTTGCTTTCCGTCAGGCGTGCTCACGACAAGCACTGAGTCCATCTCATCTGTTGCGAAAAGCTGGGAACTTCCTGCTTCTTTGCCAGCGGACTCAGCAAGAGGGAGCGTTCCAATCCGCAGGATT‑GAACGCTTCTGATCGTTTCCAATGACAGTGGTCAGTTTGCTGGTCTGGTAGAAATATTTGGGGGTATCAGGACTGTTTTTCACCGACGTGCTCCCGTTGTCTGCAGTTGGCCCCGGAAATAGATGGCCTGGGC  >  NC_002947/396703‑396995

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