Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,385,686	G→A	7.0%	A125A (GCC→GCT)	PP_2987 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,385,686	0	G	A	7.0%	164.5 / 2.4	58	A125A (GCC→GCT)	PP_2987	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (26/27);  new base A (2/2);  total (28/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 = 1.00e+00

TGTTGACGGTTGATCCAGCATCACCCTGATCACTGAACAGCACCTCACCGGCCGCCGACCAAGGCTGCAGAGGGCGTGTCAACTTTCCTCAGGTGTGACAATGCGCA‑GGGGGCGCGTGCCGATCTTCGCCCCCGTAATCTCGTCCACGGCGACAGCTTTGATTTCTGTGCCTGTTTGCGCATCGAAGAAGCGCGTCAGCTTACCACCGCCTCGGTATTGCCGGCCCCATGCACCAATCATGAACAGCACAGGCAGAAAATCGCGACCCGCAGCGGTCAGCAGGTACTCCTCCC  >  NC_002947/3385539‑3385831                                                                                                                                                     |                                                                                                                                                  tGTTGACGGTTGATCCAGCATCACCCTGATCACTGAACAGCACCTCACCGGCCGCCGACCAAGGCTGCAGAGGGCGTGTCAACTTTCCTCAGGTGTGACAATGCGCA‑GGGGGCGCGTGCCGATCTTCGCCCCCGTAATCTCGTCCACgg                                                                                                                                                  >  8:7350/1‑149 (MQ=255)     gACGGTTGATCCAGCATCACCCTGATCACTGAACAGCACCTCACCGGCCGCCGACCAAGGCTGCAGAGGGCGTGTCAACTTTCCTCAGGTGTGACAATGCGCA‑GGGGGCGCGTGCCGATCTTCGCCCCCGTAATCTCGTCCACGGCGa                                                                                                                                               >  2:328826/1‑148 (MQ=255)         gTTGATCCAGCATCACCCTGATCACTGAACAGCACCTCACCGGCCGCCGACCAAGGCTGCAGAGGGCGTGTCAACTTTCCTCAGGTGTGACAATGCGCA‑TGGGGCGCGTGCCGATCTTCGCCCCCGTAATCTCGTCCACGGCGACAGCt                                                                                                                                          <  1:91815/149‑1 (MQ=255)         gTTGATCCAGCATCACCCTGATCACTGAACAGCACCTCACCGGCCGCCGACCAAGGCTGCAGAGGGCGTGTCAACTTTCCTCAGGTGTGACAATGCGCA‑GGGGGCGCGTGCCGATCTTCGCCCCCGTAATCTCGTCCACGGCGACAGCt                                                                                                                                          >  1:124114/1‑149 (MQ=255)              tCCAGCATCACCCTGATCACTGAACAGCACCTCACCGGCCGCCGACCAAGGCTGCAGAGGGCGTGTCAACTTTCCTCAGGTGTGACAATGCGCA‑GGGGGCGCGTGCCGATCTTCGCCCCCGTAATCTCGTCCACGGCGACAGCTTTGAt                                                                                                                                     <  1:1472/149‑1 (MQ=255)                cAGCATCACCCTGATCACTGAACAGCACCTCACCGGCCGCCGACCAAGGCTGCAGAGGGCGTGTCAACTTTCCTCAGGTGTGACAATGCGCA‑GGGGGCGCGTGCCGATCTTCGCCCCCGTAATCTCGTCCACGGCGACAGCTTTGAttt                                                                                                                                   >  1:283244/1‑149 (MQ=255)                  gCATCACCCTGATCACTGAACAGCACCTCACCGGCCGCCGACCAAGGCTGCAGAGGGCGTGTCAACTTTCCTCAGGTGTGACAATGCGCA‑GGGGGCGCGTGCCGATCTTCGCCCCCGTAATCTCGTCCACGGCGACAGCTTTGATTTCt                                                                                                                                 <  2:124114/149‑1 (MQ=255)                    aTCACCCTGATCACTGAACAGCACCTCACCGGCCGCCGACCAAGGCTGCAGAGGGCGTGTCAACTTTCCTCAGGTGTGACAATGCGCA‑GGGGGCGCGTGCCGATCTTCGCCCCCGTAATCTCGTCCACGGCGACAGCTTTGATTTCtgt                                                                                                                               <  8:218170/149‑1 (MQ=255)                        cccTGATCACTGAACAGCACCTCACCGGCCGCCGACCAAGGCTGCAGAGGGCGTGTCAACTTTCCTCAGGTGTGACAATGCGCA‑GGGGGCGCGTGCCGATCTTCGCCCCCGTAATCTCGTCCACGGCGACAGCTTTGATTTCTGTGCCt                                                                                                                           <  5:138525/149‑1 (MQ=255)                         ccTGATCACTGAACAGCACCTCACCGGCCGCCGACCAAGGCTGCAGAGGGCGTGTCAACTTTCCTCAGGTGTGACAATGCGCA‑GGGGGCGCGTGCCGATCTTCGCCCCCGTAATCTCGTCCACGGCGACAGCTTTGATTTCTGTGCCTg                                                                                                                          >  5:171059/1‑149 (MQ=255)                         ccTGATCACTGAACAGCACCTCACAGGCCGCCGACCAAGGCTGCAGAGGGCGTGTCAACTTTCCTCAGGTGTGACAATGCGCA‑GGGGGCGCGTGCCGATCTTCGCCCCCGTAATCTCGTCCAAGGCGACAGCTTTGATTTCTGTGCCTg                                                                                                                          >  7:54265/1‑149 (MQ=255)                          cTGATCACTGAACAGCACCTCACCGGCCGCCGACCAAGGCTGCAGAGGGCGTGTCAACTTTCCTCAGGTGTGACAATGCGCA‑GGGGGCGCGTGCCGATCTTCGCCCCCGTAATCTCGTCCACGGCGACAGCTTTGATTTCTGTGCCTGt                                                                                                                         >  4:292885/1‑149 (MQ=255)                          cTGATCACTGAACAGCACCTCACCGGCCGCCGACCAAGGCTGCAGAGGGCGTGTCAACTTTCCTCAGGTGTGACAATGCGCA‑GGGGGCGCGTGCCGATCTTCGCCCCCGTAATCTCGTCCACGGCGACAGCTTTGATTTCTGTGCCTGt                                                                                                                         >  6:271398/1‑149 (MQ=255)                            gATCACTGAACAGCACCTCACCGGCCGCCGACCAAGGCTGCAGAGGGCGTGTCAACTTTCCTCAGGTGTGACAATGCGCA‑GGGGGCGCGTGCCGATCTTCGCCCCCGTAATCTCGTCCACGGCGACAGCTTTGATTTCTGTGCCTGttt                                                                                                                       >  5:323359/1‑149 (MQ=255)                             aTCACTGAACAGCACCTCACCGGCCGCCGACCAAGGCTGCAGAGGGCGTGTCAACTTTCCTCAGGTGTGACAATGCGCA‑GGGGGCGCGTGCCGATCTTCGCCCCCGTAATCTCGTCCACGGCGACAGCTTTGATTTCTGTGCCTGttt                                                                                                                       <  5:371528/148‑1 (MQ=255)                                   gAACAGCACCTCACCGGCCGCCGACCAAGGCTGCAGAGGGCGTGTCAACTTTCCTCAGGTGTGACAATGCGCA‑GGGGGCGCGTGCCGATCTTCGCCCCCGTAATCTCGTCCACGGCGACAGCTTTGATTTCTGTGCCTGTTTGCGCATc                                                                                                                >  6:266718/1‑149 (MQ=255)                                      cAGCACCTCACCGGCCGCCGACCAAGGCTGCAGAGGGCGTGTCAACTTTCCTCAGGTGTGACAATGCGCA‑GGGGGCGCGTGCCGATCTTCGCCCCCGTAATCTCGTCCACGGCGACAGCTTTGATTTCTGTGCCTGTTTGCGCATCg                                                                                                               >  4:111056/1‑147 (MQ=255)                                               aCCGGCCGCCGACCAAGGCTGCAGAGGGCGTGTCAACTTTCCTCAGGTGTGACAATGCGCA‑GGGGGCGCGTGCCGATCTTCGCCCCCGTAATCTCGTCCACGGCGACAGCTTTGATTTCTGTGCCTGTTTGCGCATCGAAGAAGCGCGt                                                                                                    >  1:275264/1‑149 (MQ=255)                                                  ggCCGCCGACCAAGGCTGCAGAGGGCGTGTCAACTTTCCTCAGGTGTGACAATGCGCA‑GGGGGCGCGTGCCGATCTTCGCCCCCGTAATCTCGTCCACGGCGACAGCTTTGATTTCTGTGCCTGTTTGCGCATCGAAGAAGCGCGTCAg                                                                                                 >  3:128743/1‑149 (MQ=255)                                                   gcggcggACCAAGGCTGCAGAGGGCTTGTCAACTTTCCTCACGTGTGACAATGCGCA‑GGGTGCGCGTGCCGATCTTCGCCCCCGTAATCTCTTCCACGGCGACAGCTTTGATTTCTGTGCCTGTTTGCGCATCGAAGAAGCGCGTCAGc                                                                                                <  1:295749/143‑1 (MQ=255)                                                     cgccgACCAAGGCTGCAGAGGGCGTGTCAACTTTCCTCAGGTGTGACAATGCGCA‑GGGGGCGCGTGCCGATCTTCGCCCCCGTAATCTCGTCCACGGCGACAGCTTTGATTTCTGTGCCTGTTTGCGCATCGAAGAAGCGCGTCattgt                                                                                              >  6:285540/1‑145 (MQ=255)                                                      gccgACCAAGGCTGCAGAGGGCGTGTCAACTTTCCTCAGGTGTGACAATGCGCA‑GGGGGCGCGTGCCGATCTTCGCCCCCGTAATCTCGTCCACGGCGACAGCTTTGATTTCTGTGCCTGTTTGCGCATCGAAGAAGCGCGTCAGCTTa                                                                                             >  4:16532/1‑149 (MQ=255)                                                            cAAGGCTGCAGAGGGCGTGTCAACTTTCCTCAGGTGTGACAATGCGCA‑GGGGGCGCGTGCCGATCTTCGCCCCCGTAATCTCGTCCACGGCGACAGCTTTGATTTCTGTGCCTGTTTGCGCATCGAAGAAGCGCGTCAGCTTACCACCg                                                                                       <  1:328826/149‑1 (MQ=255)                                                             aaGGCTGCAGAGGGCGTGTCAACTTTCCTCAGGTGTGACAATGCGCA‑GGGGGCGCGTGCCGATCTTCGCCCCCGTAATCTCGTCCACGGCGACAGCTTTGATTTCTGTGCATGTTTGCGCATCGAAGAAGCGCGTCAGCTTACCACCGc                                                                                      >  2:344137/1‑149 (MQ=255)                                                               ggCTGCAGAGGGCGTGTCAACTTTCCTCAGGTGTGACAATGCGCA‑GGGGGCGCGTGCCGATCTTCGCCCCCGTAATCTCGTCCACGGCGACAGCTTTGATTTCTGTGCCTGTTTGCGCATCGAAGAAGCGCGTCAGCTTACCACCGCCt                                                                                    <  5:266718/149‑1 (MQ=255)                                                                     agagGGCGTGTCAACTTTCCTCAGGTGTGACAATGCGCA‑GGGGGCGCGTGCCGATCTTCGCCCCCGTAATCTCGTCCACGGCGACAGCTTTGATTTCTGTGCCTGTTTGCGCATCGAAGAAGCGCGTCAGCTTACCACCGCCTCGGTAt                                                                              <  7:7350/149‑1 (MQ=255)                                                                      gagGGCGTGTCAACTTTCCGCAGGTGTGACAATGCGCA‑GGGGGCGCGTGCCGATCTTCGCCCCCGTAATCTCGTCCACGGCGACAGCTTTGATTTCTGTGCCTGTTTGCGCATCGAAGAAGCGCGTCAGCTTACCACCGCCTCGGTAtt                                                                             >  7:323654/1‑149 (MQ=255)                                                                        gggCGTGTCAACTTTCCTCAGGTGTGACAATGCGCA‑GGGGGCGCGTGCCGATCTTCGCCCCCGTAATCTCGTCCACGGCGACAGCTTTGATTTCTGTGCCTGTTTGCGCATCGAAGAAGCGCGTCAGCTTACCACCGCCTCGGTATTGc                                                                           >  7:107394/1‑149 (MQ=255)                                                                         ggCGTGTCAACTTTCCTCAGGTGTGACAATGCGCA‑GGGGGCGCGTGCCGATCTTCGCCCCCGTAATCTCGTCCACGGCGACAGCTTTGATTTCTGTGCCTGTTTGCGCATCGAAGAAGCGCGTCAGCTTACCACCGCCTCGGTATTGcc                                                                          <  6:171059/149‑1 (MQ=255)                                                                         ggCGTGTCAACTTTCCTCAGGTGTGACAATGCGCA‑GGGGGCGCGTGCCGATCTTCGCCCCCGTAATCTCGTCCACGGCGACAGCTTTGATTTCTGTGCCTGTTTGCGCATCGAAGAAGCGCGTCAGCTTACCACCGCCTCGGTATTGcc                                                                          <  5:47513/149‑1 (MQ=255)                                                                         ggCGTGTCAACTTTCCTCAGGTGTGACAATGCGCA‑GGGGGCGCGTGCCGATCTTCGCCCCCGTAATCTCGTCCACGGCGACAGCTTTGATTTCTGTGCCTGTTTGCGCATCGAAGAAGCGCGTCAGCTTACCACCGCCTCGGTATTGcc                                                                          <  8:107394/149‑1 (MQ=255)                                                                                 aaCTTTCCTCAGGTGTGACAATGCGCC‑GGGGGCGCGTGCCGATCTTCGCCCCCGTAATCTCGTCCA‑AGCGACAGCTTTGATTTCTGTGCCTGTTTGCGCATCGAAGAAGCGCGTCAGCTTACCACCGCCTCGGTATTGCCGGCCCCATg                                                                 >  3:271216/1‑149 (MQ=255)                                                                                 aaCTTTCCTCAGGTGTGACAATGCGCA‑GGGGGCGCGTGCCGATCTTCGCCCCCGTAATCTCGTCCA‑AGCGACAGCTTTGAGTTCTGTGCCTGTTTGCGCATCGAAGAAGCGCGTCCGCTTTCCACCGCATCGGTATTGCCGGCGCCATg                                                                 >  3:266362/1‑149 (MQ=255)                                                                                      tCCTCAGGTGTGACAATGCGCA‑GGGGGCGCGTGCCGATCTTCGCCCCCGTAATCTCGTCCACGGCGACAGCTTTGATTTCTGTGCCTGTTTGCGCATCGAAGAAGCGCGTCAGCTTACCACCGCCTCGGTATTGCCGGCCCCATGCa                                                               <  5:271398/147‑1 (MQ=255)                                                                                      tCCTCAGGTGTGACAATGCGCA‑GGGGGCGCGTGCCGATCTTCGCCCCCGTAATCTCGTCCACGGCGACAGCTTTGATTTCTGTGCCTGTTTGCGCATCGAAGAAGCGCGTCAGCTTACCACCGCCTCGGTATTGCCGGCCCCATGCAc                                                              <  6:323359/148‑1 (MQ=255)                                                                                        cTCAGGTGTGACAATGCGCA‑GGGGGCGCGTGCCGATCTTCGCCCCCGTAATCTCGTCCACGGCGACAGCTTTGATTTCTGTGCCTGTTTGCGCATCGAAGAAGCGCGTCAGCTTACCACCGCCTCGGTATTGCCGGCCCCATGCACCaa                                                           >  6:74009/1‑149 (MQ=255)                                                                                           aGGTGTGACAATGCGCA‑GGGGGCGCGTGCCGATCTTCGCCCCCGTAATCTCGTCCACGGCGACAGCTTTGATTTCTGTGCCTGTTTGCGCATCGAAGAAGCGCGTCAGCTTACCACCGCCTCGGTATTGCCGGCCCCATGTACCAATc                                                         <  2:283244/148‑1 (MQ=255)                                                                                           aGGTGTGACAATGCGCA‑GGGGGCGCGTGCCGATCTTCGCCCCCGTAATCTCGTCCACGGCGACAGCTTTGATTTCTGTGCCTGTTTGCGCATCGAAGAAGCGCGTCAGCTTACCACCGCCTCGGTATTGCCGGCCCCATGCACCAATc                                                         <  5:799/148‑1 (MQ=255)                                                                                           aGGTGTGACAATGCGCA‑GGGGGCGCGTGCCGATCTTCGCCCCCGTAATCTCGTCCACGGCGACAGCTTTGATTTCTGTGCCTGTTTGCGCATCGAAGAAGCGCGTCAGCTTACCACCGCCTCGGTATTGCCGGCCCCATGCACCAATCa                                                        <  5:60302/149‑1 (MQ=255)                                                                                            ggTGTGACAATGCGCA‑GGGGGCGCGTGCCGATCTTCGCCCCCGTAATCTCGTCCACGGCGACAGCTTTGATTTCTGTGCCTGTTTGCGCATCGAAGAAGCGCGTCAGCTTACCACCGCCTCGGTATTGCCGGCCCCATGCACCAATCAt                                                       <  2:275264/149‑1 (MQ=255)                                                                                            ggTGTGACAATGCGCA‑GGGGGCGCGTGCCGATCTTCGCCCCCGTAATCTCGTCCACGGCGACAGCTTTGATTTCTGTGCCTGTTTGCGCATCGAAGAAGCGCGTCAGCTTACCACCGCCTCGGTATTGCCGGCCCCATGCACCAATCAt                                                       <  6:209545/149‑1 (MQ=255)                                                                                                tgACAATGCGCAGGGGGGCGCGTGCCGATCTTCGCCCCCGTAATCTCGTCCACGGCGACAGCTTTGATTTCTGTGCCTGTTTGCGCATCGAAGAAGCGCGTCAGCTTACCACCGCCTCGGTATTGCCGGCCCCATGCACCAATCATGaa                                                    >  6:13120/1‑149 (MQ=255)                                                                                                    aaTGCGCA‑GGGGG‑GCGTGCCGATCTTCGCCCCCGTAATCTCGTCCACGGCGACAGCTTTGATTTCTGTGCCTGTTTGCGCATCGAAGAAGCGCGTCAGCTTACCACCGCCTCGGTATTGCCGGCCCCATGCACCAATCATGAACAGca                                               >  1:338068/1‑148 (MQ=255)                                                                                                        cgcA‑GGGGGCGCGTGCCGATCTTCGCCCCCGTAATCTCGTCCACGGCGACAGCTTTGATTTCTGTGCCTGTTTGCGCATCGAAGAAGCGCGTCAGCTTACCACCGCCTCGGTATTGCCGGCCCCATGCACCAATCATGAACAGCACAgg                                           >  4:221457/1‑149 (MQ=255)                                                                                                        cgcA‑GGGGGCGCGTGCCGATCTTCGCCCCCGTAATCTCGTCCACGGCGACAGCTTTGATTTCTGTGCCTGTTTGCGCATCGAAGAAGCGCGTCAGCTTACCACCGCCTCGGTATTGCCGGCCCCATGCACCAATCATGAACAGCACAg                                            <  8:323654/148‑1 (MQ=255)                                                                                                         gcA‑GGGGGCGCGTGCCGATCTTCGCCCCCGTAATCTCGTCCACGGCGACAGCTTTGATTTCTGTGCCTGTTTGCGCATCGAAGAAGCGCGTCAGCTTACCACCGCCTCGGTATTGCCGGCCCCATGCAc                                                              <  4:32846/129‑1 (MQ=255)                                                                                                         gcA‑GGGGGCGCGTGCCGATCTTCGCCCCCGTAATCTCGTCCACGGCGACAGCTTTGATTTCTGTGCCTGTTTGCGCATCGAAGAAGCGCGTCAGCTTACCACCGCCTCGGTATTGCCGGCCCCATGCAc                                                              >  3:32846/1‑129 (MQ=255)                                                                                                                  cgcgTGCCGATCTTCGACGCCGTAATCTCGTCCACTGCTACATCTTTGATTTCTGTGCCTGTTTTAGCATAGAATAAGCGCGTCAGCTTACAACCGTCTCGTTATTGCCGGCCCCATGCATCAATCATGAACAGCACAGGCAGAATATc                                  <  8:54265/149‑1 (MQ=255)                                                                                                                    cgTGCCGATCTTCGCCCCCGTAATCTCGTCCACGGCGACAGCTTTGATTTCTGTGCCTGTTTGCGCATCGAAGAAGCGCGTCAGCTTACCACCGCCTCGGTATTGCCGGCCCCATGCACCAATCATGAACAGCACAGGCAGAAAATcgc                                <  3:292885/149‑1 (MQ=255)                                                                                                                                cGCCCCCGTAATCTCGTCCA‑AGCGCCATCTTTGATTTCTGAGCCTGTTTGCGCATCGAAGAAGCGCGTCAGCTTACCACCGCCTCGGTATTGCCGGCCCCATGCACCAATCATGAACAGCACAGGCAGAAAATCGCGACCCGCAGCg                     <  4:271216/147‑1 (MQ=255)                                                                                                                                cGCCCCCGTAATCTAGTCCA‑AGCGCCAGCTTTGATTTCTGTGCCGGTTTGCGCATCGAAGAAGCGCGTCAGCTTACCACCGCCTCGGTATTGCCGGCCCCATGCACCAATCATGAACAGCACAGGCAGAAAATCGCGACCCGCAGCg                     <  4:266362/147‑1 (MQ=255)                                                                                                                                     ccGTAATCTCGTCCACGGCGACAGCTTTGATTTCTGTGCCTGTTTGCGCATCGAAGAAGCGCGTCAGCTTACCACCGCCTCGGTATTGCCGGCCCCATGCACCAATCATGAACAGCACAGGCAGAAAATCGCGACCCGCAGCGGTcag                >  5:241516/1‑148 (MQ=255)                                                                                                                                        tAATCTCGTCCACGGCGACAGCTTTGATTTCTGTGCCTGTTTGCGCATCGAAGAAGCGCGTCAGCTTACCACCGCCTCGGTATTGCCGGCCCCATGCACCAATCATGAACAGCACAGGCAGAAAATCGCGACCCGCAGCGGTCAGCAgg            <  5:74009/149‑1 (MQ=255)                                                                                                                                        tAATCTCGTCCACGGCGACAGCTTTGATTTCTGTGCCTGTTTGCGCATCGAAGAAGCGCGTCAGCTTACCACCGCCTCGGTATTGCCGGCCCCATGCACCAATCATGAACAGCACAGGCAGAAAATCGCGACCCGCAGCGGTCAGCAgg            <  3:221457/149‑1 (MQ=255)                                                                                                                                        tAATCTCGTCCACGGCGACAGCTTTGATTTCTGTGCCTGTTTGCGCATCGAAGAAGCGCGTCAGCTTACCACCGCCTCGGTATTGCCGGCCCCATGCACCAATCATGAACAGCACAGGCAGAAAATCGCGACCCGCAGCGGTCAGCAgg            <  1:344137/149‑1 (MQ=255)                                                                                                                                             tcGTCCACGGCGACAGCTTTGATTTCTGTGCCTGTTTGCTCATCGAAGAAGCGCGTCAGCTTACCACCGCCTCGGTATTGCCTGCCCCATGCACCAATCATCAACAGCACAGGCAGAAAATCGCGACCCGCAGCGGTCAGCAGGTActc       >  8:269536/1‑149 (MQ=255)                                                                                                                                             tcGTCCACGGCGACAGCTTTGATTTCTGTGCCTGTTTGCGCATCGAAGAAGCGCGTCAGCTTACCACCGCCTCGGTATTGCCGGCCCCATGCACCAATCATGAACAGCACAGGCAGAAAATCGCGACCCGCAGCGGTCAGCAgg            >  2:294894/1‑144 (MQ=255)                                                                                                                                             tcGTCCACGGCGACAGCTTTGATTTCTGTGCCTGTTTGCGCATCGAAGAAGCGCGTCAGCTTACCACCGCCTCGGTATTGCCGGCCCCATGCACCAATCATGAACAGCACAGGCAGAAAATCGCGACCCGCAGCGGTCAGCAgg            <  1:294894/144‑1 (MQ=255)                                                                                                                                                   aCGGCGACAGCTTTGATTTCTGTGCCTGTTTGCGCATCGAAGAAGCGCGTCAGCTTACCACCGCCTCGGTATTGCCGGCCCCATGCACCAATCATGAACAGCACAGGCAGAAAATCGCGACCCGCAGCGGTCAGCAGGTACTCCTccc  >  4:107856/1‑148 (MQ=255)                                                                                                                                                     |                                                                                                                                                  TGTTGACGGTTGATCCAGCATCACCCTGATCACTGAACAGCACCTCACCGGCCGCCGACCAAGGCTGCAGAGGGCGTGTCAACTTTCCTCAGGTGTGACAATGCGCA‑GGGGGCGCGTGCCGATCTTCGCCCCCGTAATCTCGTCCACGGCGACAGCTTTGATTTCTGTGCCTGTTTGCGCATCGAAGAAGCGCGTCAGCTTACCACCGCCTCGGTATTGCCGGCCCCATGCACCAATCATGAACAGCACAGGCAGAAAATCGCGACCCGCAGCGGTCAGCAGGTACTCCTCCC  >  NC_002947/3385539‑3385831

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