Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,649,273	G→A	7.0%	G547G (GGC→GGT)	PP_4113 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,649,273	0	G	A	7.0%	207.1 / 3.0	71	G547G (GGC→GGT)	PP_4113	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (34/31);  new base A (2/3);  total (37/34)
Fisher's exact test for biased strand distribution p-value = 6.69e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTGTGGTGCAGTCCTGACGGACGTTAGAAAAACGAGAGGCGCGATGTTACGGAATTTGTCTGACAGCGGCCACCGCTCAACGCACCCGTTCCAGGCTCTGCACGGCCAGGTCCAGCGCACGCTGCATGTCCAGGCGCGCC‑GAACGGCCCAGGTCCGGGTCGTTCAGGCGGCCACCACGCTCTGAGGGCAGGATCGGCGCCGTCAGGTCTTCGGCATCCATCGCTTCGAAATCCTGTTCGGCGCCCAGGGTCAAGGCACTGCGCCGCAACAACATGCGCACCTGTTCCGGGGT  >  NC_002947/4649128‑4649419                                                                                                                                                   |                                                                                                                                                   cTGTGGTGCAGTCCTGACGGACGTTAGAAAAACGAGAGGCGCGATGTTACGGAATTTGTCTGACAGCGGCCACCGCTCAACGCACCCGTTCCAGGCTCTGCACGGCCAGGTCCAGCGCACGCTGCATGTCCAGGCGCGCC‑GAACGGcc                                                                                                                                                  >  5:199349/1‑148 (MQ=255)  tgtgGTGCAGTCCTGACGGACGTTAGAAAAACGAGAGGCGCGATGTTACGGAATTTGTCTGACAGCGGCCACCGCTCAACGCACCCGTTCCAGGCTCTGCACGGCCAGGTCCAGCGCACGCTGCATGTCCAGGCGCGCC‑GAACGGcc                                                                                                                                                  >  1:140470/1‑147 (MQ=255)    tggtgCAGTCCTGACGGACGTTAGAAAAACGAGAGGCGCGATGTTACGGAATTTGTCTGACAGCGGCCACCGCTCAACGCACCCGTTCCAGGCTCTGCACGGCCAGGTCCAGCGCACGCTGCATGTCCAGGCGCGCC‑GAACGGCCCAgg                                                                                                                                              >  7:118116/1‑149 (MQ=255)    tggtgCAGTCCTGACGGACGTTAGAAAAACGAGAGGCGCGATGTTACGGAATTTGTCTGACAGCGGCCACCGCTCAACGCACCCGTTCCAGGCTCTGCACGGCCAGGTCCAGCGCACGCTGCATGTCCAGGCGCGCC‑GAACGGCCCAg                                                                                                                                               <  3:201289/148‑1 (MQ=255)      gtgCAGTCCTGACGGACGTTAGAAAAACGAGAGGCGCGATGTTACGGAATTTGTCTGACAGCGGCCACCGCTCAACGCACCCGTTCCAGGCTCTGCACGGCCAGGTCCAGCGCACGCTGCATGTCCAGGCGCGCC‑GAACGAcccggacc                                                                                                                                            >  6:358854/1‑144 (MQ=255)       tgCAGTCCTGACGGACGTTAGAAAAACGAGAGGCGCGATGTTACGGAATTTGTCTGACAGCGGCCACCGCTCAACGCACCCGTTCCAGGCTCTGCACGGCCAGGTCCAGCGCACGCTGCATGTCCAGGCGCGCC‑GAACGGCCCAGGTcc                                                                                                                                           <  1:40046/149‑1 (MQ=255)       tgCAGTCCTGACGGACGTTAGAAAAACGAGAGGCGCGATGTTACGGAATTTGTCTGACAGCGGCCACCGCTCAACGCACCCGTTCCAGGCTCTGCACGGCCAGGTCCAGCGCACGCTGCATGTCCAGGCGCGCC‑GAACGGCCCAGGTcc                                                                                                                                           <  5:258769/149‑1 (MQ=255)         cAGTCCTGACGGACGTTAGAAAAACGAGAGGCGCGATGTTACGGAATTTGTCTGACAGCGGCCACCGCTCAACGCACCCGTTCCAGGCTCTGCACGGCCAGGTCCAGCGCACGCTGCATGTCCAGGCGCGCC‑GAACGGCCCAGGTCCgg                                                                                                                                         <  6:287768/149‑1 (MQ=255)              cTGACGGACGTTAGAAAAACGAGAGGCGCGATGTTACGGAATTTGTCTGACAGCGGCCACCGCTCAACGCACCCGTTCCAGGCTCTGCACGGCCAGGTCCAGCGCACGCTGCATGTCCAGGCGCGCC‑GAACGGCCCAGGTCCGGGTCGt                                                                                                                                    >  5:81539/1‑149 (MQ=255)                gacgcacgTTAGAAAAACGAGAGTAGCCATGTAAATGAATTTGTCTGACATAGGCCACCGCTCAACCAACCCCATCCAGCCTATGCACGTCCAGGTCCAGCTCACGCTGCATGTCCAGACGCGCC‑GAACGGCCCAGGTCCGGGTCgat                                                                                                                                   >  6:295872/1‑146 (MQ=255)                gacgcacgTTAGAAAAACGAGAGGCGCGATGTTACGGAATTTGTCTGACAGCGCCCACCGCTCAACGCACCCGATCCAGCCTCTGCACGGCCAGGTCCAGCGCACGCTGCAAGTCCAGGCGCGCC‑GAACGGCCCAGGACCGGGTCCTTc                                                                                                                                  >  6:290333/1‑149 (MQ=255)                   ggacgTTAGAAAAACGAGAGGCGCGATGTTACGGAATTTGTCTGACAGCGGCCACCGCTCAACGCACCCGTTCCAGGCTCTGCACGGCCAGGTCCAGCGCACGCTGCATGTCCAGGCGCGCC‑GAACGACCCGGACCTGGGCCGTTcggc                                                                                                                               <  5:358854/149‑4 (MQ=255)                      cgTTAGAAAAACGAGAGGCGCGATGTTACGGAATTTGTCTGACAGCGGCCACCGCTCAACGCACCCGTTCCAGGCTCTGCACGGCCAGGTCCAGCGCACGCTGCATGTCCAGGCGCGCCTGAACGACCCGGACCTGGGCCGTTcggcgc                                                                                                                             <  6:16524/149‑6 (MQ=255)                       gatAGAAAAACGAGAGGCGCGATGTTACGTAATATGTCTGACAGCGGCCACCGCTCAACGCACCCGTTCCAGGCTCTGCACGGCCAGGTCCAGCGCACGCAGCATGTCCAGGCGCGCC‑GAACGGCACAGGTCCGGGACGTACAcgctgc                                                                                                                           >  2:40660/3‑149 (MQ=255)                        ttAGAAAAACGAGAGGCGCGATGTTACGGAATTTGTCTGACAGCGGCCACCGCTCAACGCACCCGTTCCAGGCTCTGCACGGCCAGGTCCAGCGCACGCTGCATGTCCAGGCGCGCC‑GAACGGCCCAGGTCCGGGTCGTTCAggcggc                                                                                                                           >  5:266454/1‑148 (MQ=255)                         tAGAAAAACGAGAGGCGAGATGTTACGGAATTTGTCTGACAGCGGCCACCGCTAAACGCACCAGTTCCAGGCTCTGCACGGCAAGGTCCAGCGCAAGCTGCAAGTCCAGGCGCGCC‑GAACCCCCCAGGACCGCGTCGTACAGGCGGcca                                                                                                                         >  2:343493/1‑149 (MQ=255)                          aGAAAAACGAGAGGCGCGATGTTACGGAATTTGTCTGACAGCGGCCACCGCTCAACGCACCCGTTCCAGGCTCTGCACGGCCAGGTCCAGCGCACGCTGCATGTCCAGGCGCGCCTGAACGACCCGGACCTGGGCCGTTcggc                                                                                                                               <  8:355446/143‑4 (MQ=255)                          aGAAAAACGAGAGGCGCGATGTTACGGAATTTGTCTGACAGCGGCCACCGCTCAACGCACCCGTTCCAGGCTCTGCACGGCCAGGTCCAGCGCACGCTGCATGTCCAGGCGCGCCTGAACGACCCGGACCTGGGCCGTTcggc                                                                                                                               >  7:355446/1‑140 (MQ=255)                            aaaaaCGAGAGGCGCGATGTTACGGAATTTGTCTGACAGCGGCCACCGCTCAACGCACCCGTTCCAGGCTCTGCACGGCCAGGTCCAGCGCACGCTGCATGTCCAGGCGCGCC‑GAACGGCCCAGGTCCGGGTCGTTCAggcgg                                                                                                                            >  8:16226/1‑143 (MQ=255)                            aaaaaCGAGAGGCGCGATGTTACGGAATTTGTCTGACAGCGGCCACCGCTCAACGCACCCGTTCCAGGCTCTGCACGGCCAGGTCCAGCGCACGCTGCATGTCCAGGCGCGCC‑GAACGGCCCAGGTCCGGGTCGTTCAggcgg                                                                                                                            <  7:16226/143‑1 (MQ=255)                                    gagGCGCGATGTTACGGAATTTGTCTGACAGCGGCCACCGCTCAACGCACCCGTTCCAGGCTCTGCACGGCCAGGTCCAGCGCACGCTGCATGTCCAGGCGCGCC‑GAACGGCCCAGGTCCGGGTCGTTCAGGCGGCCACCACGctct                                                                                                                >  5:86407/1‑147 (MQ=255)                                    gagGCGCGATGTTACGGAATTTGTCTGACAGCGGCCACCGCTCAACGCACCCGTTCCAGGCTCTGCACGGCCAGGTCCAGCGCACGCTGCATGTCCAGGCGCGCC‑GAACGGCCCAGGTCCGGGTCGTTCAGGCGGCCACCACGCTCTGa                                                                                                              <  6:81539/149‑1 (MQ=255)                                     agGCGCGATGTTACGGAATTTGTCTGACAGCGGCCACCGCTCAACGCACCCGTTCCAGGCTCTGCACGGCCAGGTCCAGCGCACGCTGCATGTCCAGGCGCGCC‑GAACGGCCCAGGTCCGGGTCGTTCAGGCGGCCACCACGCTCTGAg                                                                                                             <  6:199349/149‑1 (MQ=255)                                          cgATGTTACGGAATTTGTCTGACAGCGGCCACCGCTCAACGCACCCGTTCCAGGCTCTGCACGGCCAGGTCCAGCGCACGCTGCATGTCCAGGCGCGCC‑GAACGGCCCAGGTCCGGGTCGTTCAgg                                                                                                                               <  2:256650/126‑1 (MQ=255)                                          cgATGTTACGGAATTTGTCTGACAGCGGCCACCGCTCAACGCACCCGTTCCAGGCTCTGCACGGCCAGGTCCAGCGCACGCTGCATGTCCAGGCGCGCC‑GAACGGCCCAGGTCCGGGTCGTTCAgg                                                                                                                               >  1:256650/1‑126 (MQ=255)                                               ttACGGAATTTGTCTGACAGCGGCCACCGCTCAACGCACCCGTTCCAGGCTCTGCACGGCCAGGTCCAGCGCACGCTGCATGTCCAGGCGCGCC‑GAACGGCCCAGGTCCGGGTCGTTCAGGCGGCCACCACGCTCTGAGGGCAGGATc                                                                                                    <  3:151300/148‑1 (MQ=255)                                               ttACGGAATTTGTCTGACAGCGGCCACCGCTCAACGCACCCGTTCCAGGCTCTGCACGGCCAGGTCCAGCGCACGCTGCATGTCCAGGCGCGCC‑GAACGGCCCAGGTCCGGGTCGTTCAGGCGGCCACCACGCTCTGAGGGCAGGATCg                                                                                                   >  7:139762/1‑149 (MQ=255)                                                  cGGAATTTGTCTGACAGCGGCCACCGCTCAACGCACCCGTTCCAGGCTCTGCACGGCCAGGTCCAGCGCACGCTGCATGTCCAGGCGCGCC‑GAACGGCCCAGGTCCGGGTCGTTCAGGCGGCCACCACGCTCTGAGGGCAGGATCGgcg                                                                                                >  8:34432/1‑149 (MQ=255)                                                  cGGAATTTGTCTGACAGCGGCCACCGCTCAACGCACCCGTTCCAGGCTCTGCACGGCCAGGTCCAGCGCACGCTGCATGTCCAGGCGCGCC‑GAACGGCCCAGGTCCGGGTCGTTCAGGCGGCCACCACGCTCTGAGGGCAGGATCGgcg                                                                                                >  3:371836/1‑149 (MQ=255)                                                  cGGAATTTGTCTGACAGCGGCCACCGCTCAACGCACCCGTTCCAGGCTCTGCACGGCCAGGTCCAGCGCACGCTGCATGTCCAGGCGCGCC‑GAACGGCCCAGGTCCGGGTCGTTCAGGCGGCCACCACGCTCTGAGGGCAGGATCGgc                                                                                                 <  7:5671/148‑1 (MQ=255)                                                  cGGAATTTGTCTGACAGCGGCCACCGCTCAACGCACCCGTTCCAGGCTCTGCACGGCCAGGTCCAGCGCACGCTGCATGTCCAGGCGCGCC‑GAACGGCCCAGGTCCGGGTCGTTCAGGCGGCCACCACGCTCTGAGGGCAGGATCGgc                                                                                                 <  4:394622/148‑1 (MQ=255)                                                  cGGAATTTGTCTGACAGCGGCCACCGCTCAACGCACCCGTTCCAGGCTCTGCACGGCCAGGTCCAGCGCACGCTGCATGTCCAGGCGCGCC‑GAACGGCCCAGGTCCGGGTCGTTCAGGCGGCCACCACGCTCTGAGGGCAGGATCGgc                                                                                                 <  8:139762/148‑1 (MQ=255)                                                  cGGAATTTGTCTGACAGCGGCCACCGCTCAACGCACCCGTTCCAGGCTCTGCACGGCCAGGTCCAGCGCACGCTGCATGTCCAGGCGCGCC‑GAACGGCCCAGGTCCGGGTCGTTCAGGCGGCCACCACGCTCTGAGGGCAGGATCGgc                                                                                                 >  3:394622/1‑148 (MQ=255)                                                     aaTTTGTCTGACAGCGGCCACCGCTCAACGCACCCGTTCCAGGCTCTGCACGGCCAGGTCCAGCGCACGCTGCATGTCCAGGCGCGCC‑GAACGGCCCAGGTCCGGGTCGTTCAGGCGGCCACCACGCTCTGAGGGCAGGATCGGCGCCg                                                                                             >  1:212701/1‑149 (MQ=255)                                                       tttGTCTGACAGCGGCCACCGCTCAACGCACCCGTTCCAGGCTCTGCACGGCCAGGTCCAGCGCACGCTGCATGTCCAGGCGCGCC‑GAACGGCCCAGGTCCGGGTCGTTCAGGCGGCCACCACGCTCTGAGGGCAGGATCGGCGCCGTc                                                                                           <  8:118116/149‑1 (MQ=255)                                                         tGTCTGACAGCGGCCACCGCTCAACGCACCCGTTCCAGGCTCTGCACGGCCAGGTCCAGCGCACGCTGCATGTCCAGGCGCGCC‑GAACGGCCCAGGTCCGGGTCGTTCAGGCGGCCACCACGCTCTGAGGGCAGGATCGGCGCCGTCa                                                                                          <  2:140470/148‑1 (MQ=255)                                                               aCAGCGGCCACCGCTCAACGCACCCGTTCCAGGCTCTGCACGGCCAGGTCCAGCGCACGCTGCATGTCCAGGCGCGCC‑GAACGGCCCAGGTCCGGGTCGTTCAGGCGGCCACCACGCTCTGAGGGCAGGATCGgcg                                                                                                <  8:213954/136‑1 (MQ=255)                                                               aCAGCGGCCACCGCTCAACGCACCCGTTCCAGGCTCTGCACGGCCAGGTCCAGCGCACGCTGCATGTCCAGGCGCGCC‑GAACGGCCCAGGTCCGGGTCGTTCAGGCGGCCACCACGCTCTGAGGGCAGGATCGgcg                                                                                                >  7:213954/1‑136 (MQ=255)                                                                cAGCGGCCACCGCTCAACGCACCCGTTCGAGGCTCTGCACGGCCAGGTCCAGCGCACGCTGCATGTCCAGGCGCGCC‑GAACGGCCCAGGTCCGGGTCGTTCAGGCGGCCACCACGCTGTGAGGGCAGGATCGGCGCCGTCAGGTCTTCg                                                                                  <  1:40660/149‑1 (MQ=255)                                                                          cGCTCAACGCACCCGTTCCAGGCTCTGCACGGCCAGGTCCAGCGCACGCTGCATGTCCAGGCGCGCC‑GAACGGCCCAGGTCCGGGTCGTTCAGGCGGCCACCACGCTCTGAGGGCAGGATCGGCGCCGTCAGGTCTTCGGCATCCATCg                                                                        >  8:107525/1‑149 (MQ=255)                                                                               aaCGCACCCGTTCCAGGCTCTGCACGGCCAGGTCCAGCGCACGCTGCATGTCCAGGCGCGCC‑GAACGGCCCAGGTCCgg                                                                                                                                         >  7:20774/1‑79 (MQ=255)                                                                               aaCGCACCCGTTCCAGGCTCTGCACGGCCAGGTCCAGCGCACGCTGCATGTCCAGGCGCGCC‑GAACGGCCCAGGTCCgg                                                                                                                                         <  8:20774/79‑1 (MQ=255)                                                                               aaCGCACCCGTTCCAGGCTCTGCACGGCCAGGTCCAGCGCACGCTGCATGTCCAGGCGCGCC‑GAACGGCCCAGGTCCGGGTCGTTCAGGCGGCCACCACGCTCTGAGGGCAGGATCGGCGCCGTCAGGTCTTCGGCATCCATCGCtt                                                                     <  5:295872/147‑1 (MQ=255)                                                                               aaCGCACCCGTTCCAGGCTCTGCACGGCCAGGTCCAGCGCACGCTGCATGTCCAGGCGCGCC‑GAACGGCCCAGGTCCGGGTCGTTCAGGCGGCCACCACGCTCTGAGGGCAGGATCGGCGCCGTCAGGTCTTCGGCATCCATCGCTTc                                                                    <  6:266454/148‑1 (MQ=255)                                                                               aaCGCACCCGTTCCAGGCTCTGCACGGCCAGGTCCAGCGCACGCTGCATGTCCAGGCGCGCC‑GAACGGCCCAGGGCCGGGTCGTTCAGGCGGCCACCACGCTCTGAGGGCAGGATCGGCGCCGTCAGGTCTTCGGCATCCATCGCtt                                                                     <  5:290333/147‑1 (MQ=255)                                                                                    aCCCGTTCCAGGCTCTGCACGGCCAGGTCCAGCGCACGCTGCATGTCCAGGCGCGCC‑GAACGGCCCAGGTCCGGGTCGTTCAGGCGGCCACCACGCTCTGAGGGCAGGATCGGCGCCGTCAGGTCTTCGGCATCCATCGCTTCGAAAt                                                               >  1:9942/1‑148 (MQ=255)                                                                                    aCCCGTTCCAGGCTCTGCACGGCCAGGTCCAGCGCACGCTGCATGTCCAGGCGCGCC‑GAACGGCCCAGGTCCGGGTCGTTCAGGCGGCCACCACGCTCTGAGGGCAGGATCGGCGCCGTCAGGTCTTCGGCATCCATCGCTTCGAAAt                                                               <  2:9942/148‑1 (MQ=255)                                                                                    aCCCGTTCCAGGCTCTGCACGGCCAGGTCCAGCGCACGCTGCATGTCCAGGCGCGCC‑GAACGGCCCAGGTCCGGGTCGTTCAGGCGGCCACCACGCTCTGAGGGCAGGATCGGCGCCGTCAGGTCTTCGGCATCCATCGCTTCGAAATc                                                              >  4:109763/1‑149 (MQ=255)                                                                                       cGTTCCAGGCTCTGCACGGCCAGGTCCAGCGCACGCTGCATGTCCAGGCGCGCC‑GAACGGCCCAGGTCCGGGTCGTTCAGGCGGCCACCACGCTCTGAGGGCAGGATCGGCGCCGTCAGGTCTTCGGCATCCATCGCTTCGAAATCCTg                                                           <  6:86407/149‑1 (MQ=255)                                                                                         ttCCAGGCTCTGCACGGCCAGGTCCAGCGCACGCTGCATGTCCAGGCGCGCC‑GAACGGCCCAGGTCCGGGTCGTTCAGGCGGCCACCACGCTCTGAGGGCAGGATCGGCGCCGTCAGGTCTTCGGCATCCATCGCTTCGAAATCCTGtt                                                         <  4:371836/149‑1 (MQ=255)                                                                                            cAGGCTCTGCACGGCCAGGTCCAGCGCACGCTGCATGTCCAGGCGCGCC‑GAACGGCCCAGGTCCGGGTCGTTCAGGCGGCCACCACGCTCTGAGGGCAGGATCGGCGCCGTCAGGt                                                                                       >  2:286609/1‑116 (MQ=255)                                                                                            cAGGCTCTGCACGGCCAGGTCCAGCGCACGCTGCATGTCCAGGCGCGCC‑GAACGGCCCAGGTCCGGGTCGTTCAGGCGGCCACCACGCTCTGAGGGCAGGATCGGCGCCGTCAGGt                                                                                       <  1:286609/116‑1 (MQ=255)                                                                                                 tctGCACGGCCAGGTCCAGCGCACGCTGCATGTCCAGGCGCGCC‑GAACGGCCCAGGTCCGGGTCGTTCAGGCGGCCACCACGCTCTGAGGGCAGGATCGGCGCCGTCAGGTCTTCGGCATCCATCGCTTCGAAATCCTGTTCGGCGcc                                                  >  6:259501/1‑148 (MQ=255)                                                                                                    gCACGGCCAGGTCCAGCGCACGCTGCATGTCCAGGCGCGCC‑GAACGGCCCAGGTCCGGGTCGTTCAGGCGGCCACCACGCTCTGAGGGCAGGATCGGCGCCGTCAGGTCTTCGGCATCCATCGCTTCGAAATCCTGTTCGGCGCCCAgg                                              >  2:364495/1‑149 (MQ=255)                                                                                                     cACGGCCAGGTCCAGCGCACGCTGCATGTCCAGGCGCGCC‑GAACGGCCCAGGTCCGGGTCGTTCAGGCGGCCACCACGCTCTGAGGGCAGGATCGGCGCCGTCAGGTCTTCg                                                                                  <  4:133711/112‑1 (MQ=255)                                                                                                     cACGGCCAGGTCCAGCGCACGCTGCATGTCCAGGCGCGCC‑GAACGGCCCAGGTCCGGGTCGTTCAGGCGGCCACCACGCTCTGAGGGCAGGATCGGCGCCGTCAGGTCTTCg                                                                                  >  3:133711/1‑112 (MQ=255)                                                                                                       cGGCCAGGTCCAGCGCACGCTGCATGTCCAGGCGCGCC‑GAACGGCCCAGGTCCGGGTCGTTCAGGCGGCCACCACGCTCTGAGGGCAGGATCGGCGCCGTCAGGTCTTCGGCATCCATCGCTTCGAAATCCTGTTCGGCGCCCAGGGt                                            >  2:114781/1‑148 (MQ=255)                                                                                                         gCCAGGTCCAGCGCACGCTGCATGTCCAGGCGCGCC‑GAACGGCCCAGGTCCGGGTCGTTCAGGCGGCCACCACGCTCTGAGGGCAGGATCGGCGCCGTCAGGTCTTCGGCATCCATCGCTTCGAAATCCTGTTCGGCGCCCAGGGTCaa                                         <  2:212701/149‑1 (MQ=255)                                                                                                           cAGGTCCAGCGCACGCTGCATGTCCAGGCGCGCC‑GAACGGCCCAGGTCCGGGTCGTTCAGGCGGCCACCACGCTCTGAGGGCAGGATCGGCGCCGTCAGGTCTTCGGCATCCATCGCTTCGAAATc                                                              >  2:370687/1‑126 (MQ=255)                                                                                                           cAGGTCCAGCGCACGCTGCATGTCCAGGCGCGCC‑GAACGGCCCAGGTCCGGGTCGTTCAGGCGGCCACCACGCTCTGAGGGCAGGATCGGCGCCGTCAGGTCTTCGGCATCCATCGCTTCGAAATc                                                              <  1:370687/126‑1 (MQ=255)                                                                                                               tCCAGCGCACGCTGCATGTCCAGGCGCGCC‑GAACGGCCCAGGTCCGGGTCGTTCAGGCGGCCACCACGCTCTGAGGGCAGGATCGGCGCCGTCAGGTCTTCGGCATCCATCGCTTCGAAATCCTGTTCGGCGCCCAGGGTCAAGGCACt                                   >  3:268422/1‑149 (MQ=255)                                                                                                               tCCAGCACACGCTGCATGTCCAGGCGCGCC‑GAACGGCCCAGGTCCGGGTCGTTCAGGCGGCCACCACGCTCTGAGGGCAGGATCGGCGCCGTCAGGTCTTCGGCATCCATCGCTTCGAAATCCTGTTCGGCGCCCAGGGTCAAGGCACt                                   >  3:3571/1‑149 (MQ=255)                                                                                                                    cgcACGCTGCATGTCCAGGCGCGCC‑GAACGGCCCAGGTCCGGGTCGTTCAGGCGGCCACCACGCTCTGAGGGCAGGATCGGCGCCGTCAGGTCTTCGGCATCCATCGCTTCGAAATCCTGTTCGGCGCCCAGGGTCAAGGCACTGcgcc                              >  1:32463/1‑149 (MQ=255)                                                                                                                            gCATGTCCAGGCGCGCC‑GAACGGCCCAGGTCCGGGTCGTTCAGGCGGCCACCACGCTCTGAGGGCAGGATCGGCGCCGTCAGGTCTTCGGCATCCATCGCTTCGAAATCCTGTTCGGCGCCCAGGGTCAAGGCACTGCGCCGcaacaa                       >  6:230860/1‑148 (MQ=255)                                                                                                                               tGTCCAGGCGCGCC‑GAACGGCCCAGGTCCGGGTCGTTCAGGCGGCCACCACGCTCTGAGGGCAGGATCGGCGCCGTCAGGTCTTCGGCATCCATCGCTTCGAAATCCTGTTCGGCGCCCAGGGTCAAGGCACTGCGCCGCAACAACATg                   <  4:268422/149‑1 (MQ=255)                                                                                                                                      gcgcgcC‑GAACGGCCCAGGTCCGGGTCGTTCAGGCGGCCACCACGCTCTGAGGGCAGGATCGGCGCCGTCAGGTCTTCGGCATCCATCGCTTCGAAATCCTGTTCGGCGCCCAGGGTCAAGGCACTGCGCCGCAACAACATGCGCACCt            <  2:170154/149‑1 (MQ=255)                                                                                                                                          gcC‑GAACGGCCCAGGTCCGGGTCGTTCAGGCGGCCCCCACGCTCTGAGGGCAGGATCGGCGCCGTCAGGTCTTCGGCATCCATCGCTTCGAAATCCTGTTCGGCGCCCAGGGTCAAGGCACTGCGCCGCAACAACATGCGCACCTGTTc        >  2:202498/1‑149 (MQ=255)                                                                                                                                          gcC‑GAACGGCCCAGGTCCGGGTCGTTCAGGCGGCCACCACGCTCTGAGGGCAGGATCGGCGCCGTCAGGTCTTCGGCATCCATCGCTTCGAAATCCTGTTCGGCGCCCAgg                                              >  1:297187/1‑111 (MQ=255)                                                                                                                                          gcC‑GAACGGCCCAGGTCCGGGTCGTTCAGGCGGCCACCACGCTCTGAGGGCAGGATCGGCGCCGTCAGGTCTTCGGCATCCATCGCTTCGAAATCCTGTTCGGCGCCCAgg                                              <  2:297187/111‑1 (MQ=255)                                                                                                                                               aaCGGCCCAGGTCCGGGTCGTTCAGGCGGCCACCACGCTCTGAGGGCAGGATCGGCGCCGTCAGGTCTTCGGCATCCATCGCTTCGAAATCCTGTTCGGCGCCCAGGGTCAAGGCACTGCGCCGCAACAACATGCGCACCTGTTCCggg    <  1:114781/149‑1 (MQ=255)                                                                                                                                               aaCGGCCCAGGTCCGGGTCGTTCAGGCGGCCACCACGCTCTGAGGGCAGGATCGGCGCCGTCAGGTCTTCAGCATCCATCGCTTCGAAAACCTGATCGGCTCCCAGGGTCAAGGCACTGCGCCGGAACAACATGCGCACGTGTTCgggg    >  6:126101/1‑149 (MQ=255)                                                                                                                                                  ggCCCAGGTCCGGGTCGTTCAGGCGGCCACCACGCTCTGAGGGCAGGATCGGCGCCGTCAGGTCTTCGGCATCCATCGCTTCGAAATCCTGTTCGGCGCCCAGGGTCAAGGCACTGCGCCGCAACAACATGCGCACCTGTTCCGGGGt  >  6:360218/1‑148 (MQ=255)                                                                                                                                                   gCCCAGGTCCGGGTCGTTCAGGCGGCCACCACGCTCTGAGGGCAGGATCGGCGCCGTCAGGTCTTCGGCATCCATCGCTTCGAAATCCTGTTCGGCGCCCAGGGTCAAGGCACTGCGCCGCAACAACATGCGCACCTGTTCCGGGGt  >  2:356826/1‑147 (MQ=255)                                                                                                                                                   |                                                                                                                                                   CTGTGGTGCAGTCCTGACGGACGTTAGAAAAACGAGAGGCGCGATGTTACGGAATTTGTCTGACAGCGGCCACCGCTCAACGCACCCGTTCCAGGCTCTGCACGGCCAGGTCCAGCGCACGCTGCATGTCCAGGCGCGCC‑GAACGGCCCAGGTCCGGGTCGTTCAGGCGGCCACCACGCTCTGAGGGCAGGATCGGCGCCGTCAGGTCTTCGGCATCCATCGCTTCGAAATCCTGTTCGGCGCCCAGGGTCAAGGCACTGCGCCGCAACAACATGCGCACCTGTTCCGGGGT  >  NC_002947/4649128‑4649419

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