Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,457,918	C→A	7.6%	H93N (CAC→AAC)	PP_1276 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,457,918	0	C	A	7.6%	181.9 / 4.2	62	H93N (CAC→AAC)	PP_1276	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (28/27);  new base A (3/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 = 8.97e-01

CTACGGGAGCGCCGAAAAATCAGCGTGGTTGTACGCGTACAAGGGGCAGGCCGAAGATGCGGACGTCTACATTGACTTCGATCTGGAACTACAGGAGGGCGTGCGGGATGTCGTCATCGGGGGTGAGGCGAACCGGGCGCTTGTTCACAAGCGCGGCACAACCTATGGCGGCTATGCGAAATCCGGGCGTATCCGGAAACTCGAAATGACTGCCACAAGCATCAGGGCTGAATCCTTCGAATTTGAAGGGGAGGATGATTTGCAGCGACCTTTCAGGGTGGTGGGTGGCCCAT  >  NC_002947/1457773‑1458065                                                                                                                                                  |                                                                                                                                                    cTACGGGAGCGCCGAAAAATCAGCGTGGTTGTACGCGTACAAGGGGCAGGCCGAAGATGCGGACGTCTACATTGACTTCGATCTGGAACTACAGGAGGGCGTGCGGGATGTCGTCATCGGGGGTGAGGCGAACCGGGCGCTTGTTcac                                                                                                                                                   >  7:348659/1‑148 (MQ=255)       gAGCGCCGAAAAATCAGCGTGGTTTTACGCGTACATGGGGCAGGCCGAAGATGCGGACGTCTACATTGACTTCGATCTGGAACTACAGGAGGGCGTGCGGGATGTCGTCATCGGGGGTGAGGCGAACCGGGCGCTTGTGAACAAgcgcc                                                                                                                                            <  7:49118/149‑2 (MQ=255)       gAGCGCCGAAAAATCAGCGTGGTTGTACGCGTACAAGGGGCAGGCCGAAGATGCGGACGTCTACATTGACTTCGATCTGGAACTACAGGAGGGCGTGCGGGATGTCGTCATCGGGGGTGAGGCGAACCGGGCGCTTGTTCACAAGcgcg                                                                                                                                            <  2:227344/149‑1 (MQ=255)        aGCGCCGAAAAATCAGCGTGGTTGTACGCGTACAAGGGGCAGGCCGAAGATGCGGACGTCTACATTGACTTCGATCTGGAACTACAGGAGGGCGTGCGGGATGTCGTCATCGGGGGTGAGGCGAACCGGGCGCTTGTTCACAAGCGCgg                                                                                                                                           <  7:14039/149‑1 (MQ=255)           gcCGAAAAATCAGCGTGGTTGTACGCGTACAAGGGGCAGGCCGAAGATGCGGACGTCTACATTGACTTCGATCTGGAACTACAGGAGGGCGTGCGGGATGTCGTCATCGGGGGTGAGGCGAACCGGGCGCTTGTTCACAAGCGCGGca                                                                                                                                         >  1:168055/1‑148 (MQ=255)                aaaaTCAGCGTGGTTGTACGCGTACAAGGGGCAGGCCGAAGATGCGGACGTCTACATTGACTTCGATCTGGAACTACAGGAGGGCGTGCGGGATGTCGTCATCGGGGGTGAGGCGAACCGGGCGCTTGTTCACAAGCGCGGCACAAcc                                                                                                                                    >  8:73393/1‑148 (MQ=255)                         gTGGTTGTACGCGTACAAGGGGCAGGCCGAAGATGCGGACGTCTACATTGACTTCGATCTGGAACTACAGGAGGGCGTGCGGGATGTCGTCATCGGGGGTGAGGCGAACCGGCCGGTTGTACACAGGCGCGGACCAACCTATgacgaca                                                                                                                          >  6:21877/1‑146 (MQ=255)                                tACGCGTACAAGGGGCAGGCCGAAGATGCGGACGTCTACATTGACTTCGATCTGGAACTACAGGAGGGCGTGCGGGATGTCGTCATCGGGGGTGAGGCGAACCGGGCGCTTGTTCACAAGCGCGGCACAACCTATGGCGGCTATGCGa                                                                                                                    >  4:72110/1‑148 (MQ=255)                                   gcctACAAGGGGCACGCCGAAGATGCGCAAGACTACATAGACATCCATCTGGAACTACAGCAGGGCGTGCGGGATGTCGTCATCGGGGGTGAGGCGAACCGGGCGCTTGATCACAAGCGCAGAACAACCTATGTCGGCAATGCGAAAta                                                                                                                >  4:323812/4‑148 (MQ=255)                                          aGGGGCAGGCCGAAGATGCGGACGTCTACATTGACTTCGATCTGGAACTACAGGAGGGCGTGCGGGATGTCGTCATCGGGGGTGAGGCGAACCGGGCGCTTGTTCACAAGCGCGGCACAACCTATGGCGGCTATGCGAAATCCGGGc                                                                                                           >  6:349275/1‑147 (MQ=255)                                          aGGGGCAGGCCGAAGATGCGGACGTCTACATTGACTTCGATCTGGAACTACAGGAGGGCGTGCGGGATGTCGTCATCGGGGGTGAGGCGAACCGGGCGCTTGTTCACAAGCGCGGCACAACCTATGGCGGCTATGCGAAATCCGGGc                                                                                                           <  5:349275/147‑1 (MQ=255)                                          aGGGGCAGGCCGAAGATGCGGACGTCTACATTGACTTCGATCTGGAACTACAGGAGGGCGTGCGGGATGTCGTCATCGGGGGTGAGGCGAACCGGGCGCTTGTTCACAAGCGCGGCACAACCTATGGCGGCTATGCGAAATCCGGGCGt                                                                                                         <  8:348659/149‑1 (MQ=255)                                                 ggCCGAAGATGCGGACGTCTACATTGACTTCGATCTGGAACTACAGGAGGGCGTGCGGGATGTCGTCATCGGGGGTGAGGCGAACCGGGCGCTTGTTCACAAGCGCGGCACAACCTATGGCGGCTATGCGAAATCCGGGCGTATCCGGa                                                                                                  <  1:347474/149‑1 (MQ=255)                                                   ccGAAGATGCGGACGTCTACATTGACTTCGATCTGGAACTACAGGAGGGCGTGCGGGATGTCGTCATCGGGGGTGAGGCGAACCGGGCGCTTGTTCACAAGCGCGGCACAACCTATGGCGGCTATGCGAAATCCGGGCGTATCCGGaa                                                                                                 >  6:88617/1‑148 (MQ=255)                                                    cGAAGATGCGGACGTCTACATTGACTTCGATCTGGAACTACAGGAGGGCGTGCGGGATGTCGTCATCGGGGGTGAGGCGAACCGGGCGCTTGTGAACAAgcgccc                                                                                                                                           <  4:374861/105‑3 (MQ=255)                                                    cGAAGATGCGGACGTCTACATTGACTTCGATCTGGAACTACAGGAGGGCGTGCGGGATGTCGTCATCGGGGGTGAGGCGAACCGGGCGCTTGTGAACAAgcgccc                                                                                                                                           >  3:374861/1‑103 (MQ=255)                                                      aaGATGCGGACGTCTACATTGACTTCGATCTGGAACTACAGGAGGGCGTGCGGGATGTCGTCATCGGGGGTGAGGCGAACCGGGCGCTTGTTCACAAGCGCGGCACAACCTATGGCGGCTATGCGAAATCCGGGCGTATCCGGAAACTc                                                                                             >  6:208327/1‑149 (MQ=255)                                                            cGGACGTCTACATTGACTTCGATCTGGAACTACAGGAGGGCGTGCGGGATGTCGTCATCGGGGGTGAGGCGAACCGGGCGCTTGTTCACAAGCGCGGCACAACCTATGGCGGCTATGCGAAATCCGGGCGTATCCGGAAACTCGAAAt                                                                                        <  2:241627/148‑1 (MQ=255)                                                              gACGTCTACATTGACTTCGATCTGGAACTACAGGAGGGCGTGCGGGATGTCGTCATCGGGGGTGAGGCGAACCGGGCGCTTGTTCACAAGCGCGGCACAACCTATggc                                                                                                                              <  8:73700/108‑1 (MQ=255)                                                              gACGTCTACATTGACTTCGATCTGGAACTACAGGAGGGCGTGCGGGATGTCGTCATCGGGGGTGAGGCGAACCGGGCGCTTGTTCACAAGCGCGGCACAACCTATggc                                                                                                                              >  7:73700/1‑108 (MQ=255)                                                              gACGTCTACATTGACTTCGATCTGGAACTACAGGAGGGCGTGCGGGATGTCGTCATCGGGGGTGAGGCGAACCGGGCGCTTGTTCACAAGCGCGGCACAACCTATGGCGGCTATGCGAAATCCGGGCGTATCCGGAAACTCGAAATGAc                                                                                     <  2:168055/149‑1 (MQ=255)                                                               aCGTCTACATTGACTTCGATCTGGAACTACAGGAGGGCGTGCGGGATGTCGTCATCGGGGGTGAGGCGAACCGGGCGCTTGTTCACAAGCGCGGCACAACCTATGGCGGCTATGCGAAATCCGGGCGTATCCGGAAACTCGAAATGACt                                                                                    <  1:377136/149‑1 (MQ=255)                                                                    tACATTGACTTCGATCTGGAACTACAGGAGGGCGTGCGGGATGTCGTCATCGGGGGTGAGGCGAACCGGGCGCTTGTTCACAAGCGCGGCACAACCTATGGCGGCTATGCGaa                                                                                                                   >  3:100854/1‑113 (MQ=255)                                                                    tACATTGACTTCGATCTGGAACTACAGGAGGGCGTGCGGGATGTCGTCATCGGGGGTGAGGCGAACCGGGCGCTTGTTCACAAGCGCGGCACAACCTATGGCGGCTATGCGaa                                                                                                                   <  4:100854/113‑1 (MQ=255)                                                                       aTTGACTTCGATCTGGAACTACAGGAGGGCGTGCGGGATGTCGTCATCGGGGGTGAGGCGAACCGGGCGCTTGTTCACAAGCGCGGCACAACCTATGGCGGCTATGCGAAATCCGGGCGTATCCGGAAACTCGAAATGACTGCCACAAg                                                                            <  3:238713/149‑1 (MQ=255)                                                                          gACTTCGATCTGGAACTACAGGAGGGCGTGCGGGATGTCGTCATCGGGGGTGAGGCGAACCGGGCGCTTGTTCACAAGCGCGGCACAACCAATGGCGGCTATGCGAAATCCGGGCGTATCCGGAAACTCGAAATGACTGCCACAAGCAt                                                                         >  3:375890/1‑149 (MQ=255)                                                                           aCTTCGATCTGGAACTACAGGAGGGCGTGCGGGATGTCGTCATCGGGGGTGAGGCGAACCGGGCGGGTGGTGACAACCGCGGCACAACCGATGGGGGCGTAGCGAAAACCGGGCCTATACGGAAACTCGAAATGACTGCCACAACCAt                                                                         >  2:386598/1‑148 (MQ=255)                                                                           aCTTCGATCTGGAACTACAGGAGGGCGTGCGGGATGTCGTCATCGGGGGTGAGGCGAACCGGGCGCTTGTTCACAAGCGCGGCACAACCTATGGCGGCTATGCGAAATCCGGGCGTATCCGGAAACTCGAAATGACTGCCACAAGCAt                                                                         <  1:54028/148‑1 (MQ=255)                                                                            cTTCGATCTGGAACTACAGGAGGGCGTGCGGGATGTCGTCATCGGGGGTGAGGCGAACCGGGCGCTTGTTCACAAGCGCGGCACAACCTATGGCGGCTATGCGAAATCCGGGCGTATCCGGAAACTCGAAATGACTGCCACAAGCATCa                                                                       <  3:247640/149‑1 (MQ=255)                                                                             ttCGATCTGGAACTACAGGAGGGCGTGCGGGATGTCGTCATCGGGGGTGAGGCGAACCGGGCGCTTGTTCACAAGCGCGGCACAACCTATGGCGGCTATGCGAAATCCGGGCGTATCCGGa                                                                                                  <  6:8421/121‑1 (MQ=255)                                                                             ttCGATCTGGAACTACAGGAGGGCGTGCGGGATGTCGTCATCGGGGGTGAGGCGAACCGGGCGCTTGTTCACAAGCGCGGCACAACCTATGGCGGCTATGCGAAATCCGGGCGTATCCGGa                                                                                                  >  5:8421/1‑121 (MQ=255)                                                                             ttCGATCTGGAACTACAGGAGGGCGTGCGGGATGTCGTCATCGGGGGTGAGGCGAACCGGGCGCTTGTTCACAAGCGCGGCACAACCTATGGCGGCTATGCGAAATCCGGGCGTATCCGGAAACTCGa                                                                                           <  2:357684/128‑1 (MQ=255)                                                                             ttCGATCTGGAACTACAGGAGGGCGTGCGGGATGTCGTCATCGGGGGTGAGGCGAACCGGGCGCTTGTTCACAAGCGCGGCACAACCTATGGCGGCTATGCGAAATCCGGGCGTATCCGGAAACTCGa                                                                                           >  1:357684/1‑128 (MQ=255)                                                                                  tCTGGAACTACAGGAGGGCGTGCGGGATGTCGTCATCGGGGGTGAGGCGAACCGGGCGCTTGTTCACAAGCGCGGCACAACCTATGGCGGCTATGCGAAATCCGGGCGTATCCGGAAACTCGAAATGACTGCCACAAGCATCAGGGCt                                                                  >  7:267315/1‑148 (MQ=255)                                                                                      gAACTACAGGAGGGCGTGCGGGATGTCGTCATCGGGGGTGAGGCGAACCGGGCGCTTGTTCACAAGCGCGGCACAACCTATGGCGGCTATGCGAAATCCGGGCGTATCCGGAAACTCGAAATGACTGCCACAAGCATCAGGGCTGAATc                                                             <  5:88617/149‑1 (MQ=255)                                                                                         cTACAGGAGGGCGTGCGGGATGGCGTCATCGGGGGTGAGGCGAACCGGGCGCTTGTTCACAAGCGCGGCACAACCTATGGCGGCTATGCGAAATCCGGGCGTATCCGGAAACTCGAAATGACTGCCACAAGCATCAGGGCTGAATCCtt                                                          <  5:359639/149‑1 (MQ=255)                                                                                          tACAGGAGGGCGTGCGGGATGTCGTCATCGGGGGTGAGGCGAACCGGGCGCTTGTTCACAAGCGCGGCACAACCTATGGCGGCTATGCGAAATCCGGGCGTATCCGGAAACTCGAAATGACTGCCACAAGCATCAGGGCTGAATCCTTc                                                         >  2:9464/1‑149 (MQ=255)                                                                                          tACAGGAGGGCGTGCGGGATGTCGTCATCGGGGGTGAGGCGAACCGGGCGCTTGTTCACAAGCGCGGCACAACCTATGGCGGCTATGCGAAATCCGGGCGTATCCGGAAACTCGAAATGACTGCCACAAGCATCAGGGCTGAATCCTTc                                                         >  4:170809/1‑149 (MQ=255)                                                                                           aCAGGAGGGCGTGCGGGATGTCGTCATCGGGGGTGAGGCGAACCGGGCGCTTGTTCACAAGCGCGGCACAACCTATGGCGGCTATGCGAAATCCGGGCGTATCCGGAAACTCGAAATGACTGCCACAAGCATCAGGGCTGAATCCTTc                                                         >  5:3667/1‑148 (MQ=255)                                                                                            cAGGAGGGCGTGCGGGATGTCGTCATCGGGGGTGAGGCGAACCGGGCGCTTGTTCACAAGCGCGGCACAACCTATGGCGGCTATGCGAAATCCGGGCGTATCCGGAAACTCGAAATGACTGCCACAAGCATCAGGGCTGAATCCTTc                                                         <  2:131835/147‑1 (MQ=255)                                                                                            cAGGAGGGCGTGCGGGATGTCGTCATCGGGGGTGAGGCGAACCGGGCGCTTGTTCACAAGCGCGGCACAACCTATGGCGGCTATGCGAAATCCGGGCGTATCCGGAAACTCGAAATGACTGCCACAAGCATCAGGGCTGAATCCTTc                                                         <  5:153190/147‑1 (MQ=255)                                                                                              ggaggGCGTGCGGGATGTCGTCATCGGGGGTGAGGCGAACCGGGCGCTTGTTCACAAGCGCGGCACAACCTATGGCGGCTATGCGAAATCCGGGCGTATCCGGAAACTCGAAATGACTGCCACAAGCATCAGGGCTGAATCCTTCGAAt                                                     <  1:9464/149‑1 (MQ=255)                                                                                                tgggCGTGCGGGATGTCGTCTTCGGGGGTGAGGCGAACCGGGCGCTTGTTCACATGCGCGGCACAACCTATGGCGGCTATGCGAAATCCGGGCGTATCCGGAAACTCGAAATGACTGCCACAAGCATCAGGGCTGAATCCTTCGAAttt                                                   <  7:277929/148‑1 (MQ=255)                                                                                                      tGCGGGATGTCGTCATCGGGGGTGGGGCGAACCGGGCGCTTGTTCACAAGCGCGGCACAACCTATGGCGGCTATGCGAAATCCGGGCGTATCCGGAAACTCGAAATGACTGCCACAAGCATCAGGGCTGAATCCTTCGAATTTGAAggg                                             >  6:263688/1‑149 (MQ=255)                                                                                                                 gtcATCGGGGGTGAGGCGAACCGGGCGCTTGTTCACAAGCGCGGCACAACCTATGGCGGCTATGCGAAATCCGGGCGTATCCGGAAACTCGAAATGACTGCCACAAGCATCAGGGCTGAATCCTTCGAATTTGAAGGGGAGGATGAttt                                  >  8:300745/1‑149 (MQ=255)                                                                                                                    aTCGGGGGTGAGGCGAACCGGGCGCTTGTTCACAAGCGCGGCACAACCTATGGCGGCTATGCGAAATCCGGGCGTAGCCGGAAACTCGAAATGACTGCCACAAGCATCAGGGCTGAATCCTTCGAATTTGAAGGGGAGGATGATTTGCa                               <  3:72110/149‑1 (MQ=255)                                                                                                                         gggTGAGGCGAACCGGGCGCTTGTTCACAAGCGCGGCACAACCTATGGCGGCTATGCGAAATCCGGGCGTATCCGGAAACTCGAAATGACTGCCACAAGCATCAGGGCTGAATCCTTCGAATTTGAAGGGGAGGATGAGTTGCAGCGAc                          >  6:357532/1‑149 (MQ=255)                                                                                                                              aGGCGAACCGGGCGCTTGTTCACAAGCGCGGCACAACCTATGGCGGCTATGCGAAATCCGGGCGTATCCGGAAACTCGAAATGACTGCCACAAGCATCAGGGCTGAATCCTTCGAATTTGAAGGGGAGGATGATTTGCAGCGACCttt                      <  5:208327/148‑1 (MQ=255)                                                                                                                               ggCGAACCGGGCGCTTGTTCACAAGCGCGGCACAACCTATGGCGGCTATGCGAAATCCGGGCGTATCCGGAAACTCGAAATGACTGCCACAAGCATCAGGGCTGAATCCTTCGAATTTGAAGGGGAGGATGATTTGCAGCGACCttt                      >  8:104824/1‑147 (MQ=255)                                                                                                                               ggCGAACCGGGCGCTTGTTCACAAGCGCGGCACAACCTATGGCGGCTATGCGAAATCCGGGCGTATCCGGAAACTCGAAATGACTGCCACAAGCATCAGGGCTGAATCCTTCGAATTTGAAGGGGAGGATGATTTGCAGCGACCTTTCa                    <  7:73393/149‑1 (MQ=255)                                                                                                                                 cGAACCGGGCGCTTGTTCACAAGCGCGGCACAACCTATGGCGGCTATGCGAAATCCGGGCGTATCCGGAAACTCGAAATGACTGCCACAAGCATCAGGGCTGAATCCTTCGAATTTGAAGGGGAGGATGATTTGCAGCGACCTTTCAgg                  <  6:3667/149‑1 (MQ=255)                                                                                                                                 cGAACCGGGCGCTTGTTCACAAGCGCGGCACAACCTATGGCGGCTATGCGAAATCCGGGCGTATCCGGAAACTCGAAATGACTGCCACAAGCATCAGGGCTGAATCCTTCGAATTTGAAGGGGAGGATGATTTGCAGCGACCTTTCAgg                  <  3:170809/149‑1 (MQ=255)                                                                                                                                  ttgtgccgcgcTTGTGAACAAGCGCGGCACAACCTATGGCGGCTATGCGAAATCCGGGCGTATCCGGAAACTCGAAATGACTGCCACAAGCATCAGGGCTGAATCCTTCGAATTTGAAGGGGAGGATGATTTGCAGCGACCTTTCAggg                 >  8:302623/8‑149 (MQ=255)                                                                                                                                     tgccgcgcTTGTGAACAAGCGCGGCACAACCTATGGCGGCTATGCGAAATCCGGGCGTATCCGGAAACTCGAAATGACTGCCACAAGCATCAGGGCTGAATCCTTCGAATTTGAAGGGGAGGATGATTTGCAGCGACCTTTCAGggtgg              >  4:189352/5‑149 (MQ=255)                                                                                                                                       gggCGCTTGTTCACAAGCGCGGCACAACCTATGGCGGCTATGCGAAATCCGGGCGTATCCGGAAACTCGAAATGACTGCCACAAGCATCAGGGCTGAATCCTTCGAATTTGAAGGGGAGGATGATTTGCAg                              >  6:59109/1‑131 (MQ=255)                                                                                                                                       gggCGCTTGTTCACAAGCGCGGCACAACCTATGGCGGCTATGCGAAATCCGGGCGTATCCGGAAACTCGAAATGACTGCCACAAGCATCAGGGCTGAATCCTTCGAATTTGAAGGGGAGGATGATTTGCAg                              <  5:59109/131‑1 (MQ=255)                                                                                                                                        ggCGCTTTTTGACAAGCGCGGCACAACCTATTGGGGCTATGCGAAATCCGGGCGTATCCGGAAAATCGAAATGACTGCCACAAGCATCAGGGCTGAATCCTTCGAATTTGaa                                                <  1:302754/112‑1 (MQ=255)                                                                                                                                        ggCGCTTGTTCACAAGCGCGGCACAACCTATGGCGGCTATGCGAAATCCGGGCGTATCCGGAAACTCGAAATGACTGCCACAAGCATCAGGGCTGAATCCTTCGAATTTGAAGGGGAGGATGATTTGCAGCGACCTTTCAGGGtggtgg           >  8:8520/1‑149 (MQ=255)                                                                                                                                             ttGTTCACAAGCGCGGCACAACCTATGGCGGCTATGCGAAATCCGGGCGTATCCGGAAACTCGAAATGACTGCCACAAGCATCAGGGCTGAATCCTTCGAATTTGAAGGGGAGGATGATTTGCAGCGACCTTTCAGGGTGGTGGGTGGc      >  4:89552/1‑149 (MQ=255)                                                                                                                                             ttGTTCACAAGCGCGGCACAACCTATGGCGGCTATGCGAAATCCGGGCGTATCCGGAAACTCGAAATGACTGCCACAAGCATCAGGGCTGAATCCTTCGAATTTGAAGGGGAGGATGATTTGCAGCGACCTTTCAGGGTGGTGGGTGGc      >  2:176589/1‑149 (MQ=255)                                                                                                                                               gTTCACAAGCGCGGCACAACCTATGGCGGCTATGCGAAATCCGGGCGTATCCGGAAACTCGAAATGACTGCCACAAGCATCAGGGCTGAATCCTTCGAATTTGAAGGGGAGGATGATTTGCAGCGACCTTTCAGGGTGGTGGGTGGcc     >  4:144171/1‑148 (MQ=255)                                                                                                                                               gTTCACAAGCGCGGCACAACCTATGGCGGCTATGCGAAATCCGGGCGTATCCGGAAACTCGAAATGACTGCCACAAGCATCAGGGCTGAATCCTTCGAATTTGAAGGGGAGGATGATTTGCAGCGACCTTTCAGGGTGGTGGGTGGcc     >  3:29035/1‑148 (MQ=255)                                                                                                                                                 tCACAAGCGCGGCACAACCTATGGCGGCTATGCGAAATCCGGGCGTATCCGGAAACTCGAAATGACTGCCACAAGCATCAGGGCTGAATCCTTCGAATTTGAAGGGGAGGATGATTTGCAGCGACCTTTCAGGGTGGTGGGTGGCCCAt  >  7:341130/1‑149 (MQ=255)                                                                                                                                                  |                                                                                                                                                    CTACGGGAGCGCCGAAAAATCAGCGTGGTTGTACGCGTACAAGGGGCAGGCCGAAGATGCGGACGTCTACATTGACTTCGATCTGGAACTACAGGAGGGCGTGCGGGATGTCGTCATCGGGGGTGAGGCGAACCGGGCGCTTGTTCACAAGCGCGGCACAACCTATGGCGGCTATGCGAAATCCGGGCGTATCCGGAAACTCGAAATGACTGCCACAAGCATCAGGGCTGAATCCTTCGAATTTGAAGGGGAGGATGATTTGCAGCGACCTTTCAGGGTGGTGGGTGGCCCAT  >  NC_002947/1457773‑1458065

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