Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	807,888	A→T	6.1%	L49Q (CTG→CAG)	PP_0694 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	807,888	0	A	T	6.1%	196.6 / 2.7	65	L49Q (CTG→CAG)	PP_0694	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (33/28);  new base T (2/2);  total (35/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 = 7.83e-01

TCGGGTCGCGGCCCTGGTCAGCAAGGTGCTCGACCAGGATGAAGGTGAACAGCGACAGGCGCTGCGCGGTCTTGTTGACCTGCGCGGCGGTCATGTCCGGCACGAAGTAGTAGAAACCACGGGTATCGCACACCAGTTCGAAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCTGGAAGTTCGACAGCTGGGCGTACAGCTCAGGGTCGCGCCGGCTGACGTGGAAGCCTTTGAACAGCTCGCGGAAGATCGGCGCGAGCTGGGACAGTTCGGAAAGATCAAGATGCATGG  >  NC_002947/807742‑808035                                                                                                                                                   |                                                                                                                                                     tCGGGTCGCGGCCCTGGTCAGCAAGGTGCTCGACCAGGATGAAGGTGAACAGCGACAGGCGCTGCGCGGTCTTGTTGACCTGCGCGGCGGTCATGTCCGGCACGAAGTAGTAGAAACCACGGGTATCGCACACCAGTTCGAAGCCCAgg                                                                                                                                                    <  5:312559/149‑1 (MQ=255)   gggTCGCGGCCCTGGTCAGCAAGGTGCTCGACCAGGATGAAGGTGAACAGCGACAGGCGCTGCGCGGTCTTGTTGACCTGCGCGGCGGTCATGTCCGGCACGAAGTAGTAGAAACCACGGGTATCGCACACCAGTTCGAAGCCCAGGGc                                                                                                                                                  >  1:234612/1‑149 (MQ=255)    ggTCGCGGCCCTGGTCAGCAAGGTGCTCGACCAGGATGAAGGTGAACAGCGACAGGCGCTGCGCGGTCTTGTTGACCTGCGCGGCGGTCATGTCCGGCACGAAGTAGTAGAAACCACGGGTATCGCACACCAGTTCGAAGCCCAGGGcc                                                                                                                                                 <  3:72642/149‑1 (MQ=255)               tGGTCAGCAAGGTGCTCGACCAGGATGAAGGTGAACAGCGACAGGCGCTGCGCGGTCTTGTTGACCTGCGCGGCGGTCATGTCCGGCACGAAGTAGTAGAAACCACGGGTATCGCACACCAGTTCGAAGCCCAGGGCCTTGAACAGG‑Gt                                                                                                                                     >  1:228482/1‑149 (MQ=255)               tGGTCAGCAAGGTGCTCGACCAGGATGAAGGTGAACAGCGACAGGCGCTGCGCGGTCTTGTTGACCTGCGCGGCGGTCATGTCCGGCACGAAGTAGTAGAAACCACGGGTATCGCACACCAGTTCGAAGCCCAGGGCCTTGAACAGG‑Gt                                                                                                                                     >  2:324463/1‑149 (MQ=255)                    aGCAAGGTGCTCGACCAGGATGAAGGTGAACAGCGACAGGCGCTGCGCGGTCTTGTTGACCTGCGCGGCGGTCATGTCCGGCACGAAGTAGTAGAAACCACGGGTATCGCACACCAGTTCGAAGCCCAGGGCCTTGAACAGG‑GTGCgg                                                                                                                                 <  1:173696/148‑1 (MQ=255)                    aGCAAGGTGCTCGACCAGGATGAAGGTGAACAGCGACAGGCGCTGCGCGGTCTTGTTGACCTGCGCGGCGGTCATGTCCGGCACGAAGTAGTAGAAACCACGGGTATCGCACACCAGTTCGAAGCCCAGGGCCTTGAACAGG‑GTGCgg                                                                                                                                 <  2:254562/148‑1 (MQ=255)                      cAAGGTGCTCGACCAGGATGAAGGTGAACAGCGACAGGCGCTGCGCGGTCTTGTTGACCTGCGCGGCGGTCATGTCCGGCACGAAGTAGTAGAAACCACGGGTATCGCACACCAGTTCGAAGCCCAGGGCCTTg                                                                                                                                              >  8:217271/1‑134 (MQ=255)                      cAAGGTGCTCGACCAGGATGAAGGTGAACAGCGACAGGCGCTGCGCGGTCTTGTTGACCTGCGCGGCGGTCATGTCCGGCACGAAGTAGTAGAAACCACGGGTATCGCACACCAGTTCGAAGCCCAGGGCCTTGAACAGG‑GTGCGGTAc                                                                                                                              >  3:65299/1‑149 (MQ=255)                      cAAGGTGCTCGACCAGGATGAAGGTGAACAGCGACAGGCGCTGCGCGGTCGTGTTGACCTGCGCGGCGGTCATGTCCGGCACGAAGTAGTAGAAACCACGGGTATCGCACACCAGTTCGAAGCCCAGGGCCTTg                                                                                                                                              <  7:217271/134‑1 (MQ=255)                       aaGGTGCTCGACCAGGATGAAGGTGAACAGCGACAGGCGCTGCGCGGTCTTGTTGACCTGCGCGGCGGTCATGTCCGGCACGAAGTAGTAGAAACCACGGGTATCGCACACCAGTTCGAAGCCCAGGGCCTTGAACAGG‑GTGCGGTACt                                                                                                                             <  6:264093/149‑1 (MQ=255)                        aGGTGCTCGACCAGGATGAAGGTGAACAGCGACAGGCGCTGCGCGGTCTTGTTGACCTGCGCGGCGGTCATGTCCGGCACGAAGTAGTAGAAACCACGGGTATCGCACACCAGTTCGAAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTg                                                                                                                            <  2:228482/149‑1 (MQ=255)                          gTGCTCGACCAGGATGAAGGTGAACAGCGACAGGCGCTGCGCGGTCTTGTTGACCTGCGCGGCGGTCATGTCCGGCACGAAGTAGTAGAAACCACGGGTATCGCACACCAGTTCGAAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGt                                                                                                                          >  1:40205/1‑149 (MQ=255)                          gTGCTCGACCAGGATGAAGGTGAACAGCGACAGGCGCTGCGCGGTCTTGTTGACCTGCGCGGCGGTCATGTCCGGCACGAAGTAGTAGAAACCACGGGTATCGCACACCAGTTCGAAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGt                                                                                                                          >  8:207/1‑149 (MQ=255)                             cTCGACCAGGATGAAGGTGAACAGCGACAGGCGCTGCGCGGTCTTGTTGACCTGCGCGGCGGTCATGTCCGGCACGAAGTAGTAGAAACCACGGGTATCGCACACCAGTTCGAAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCt                                                                                                                       >  4:245153/1‑149 (MQ=255)                              tCGACCAGGATGAAGGTGAACAGCGACAGGCGCTGCGCGGTCTTGTTGACCTGCGCGGCGGTCATGTCCGGCACGAAGTAGTAGAAACCACGGGTATCGCACACCAGTTCGAAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCTg                                                                                                                      >  3:324278/1‑149 (MQ=255)                                 aCCAGGATGAAGGTGAACAGCGACAGGCGCTGCGCGGTCTTGTTGACCTGCGCGGCGGTCATGTCCGGCACGAAGTAGTAGAAACCACGGGTATCGCACACCAGTTCGAAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCTGGaa                                                                                                                   <  1:188009/149‑1 (MQ=255)                                    aGGATGAAGGTGAACAGCGACAGGCGCTGCGCGGTCTTGTTGACCTGCGCGGCGGTCATGTCCGGCACGAAGTAGTAGAAACCACGGGTATCGCACACCAGTTCGAAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCTGGAAGtt                                                                                                                <  4:324278/149‑1 (MQ=255)                                          aaGGTGAACAGCGACAGGCGCTGCGCGGTCTTGTTGACCTGCGCGGCGGTCATGTCCGGCACGAAGTAGTAGAAACCACGGGTATCGCACACCAGTTCGAAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCTGGAAGTTCGACAg                                                                                                          >  6:242516/1‑149 (MQ=255)                                              tGAACAGCGACAGGCGCTGCGCGGTCTTGTTGACCTGCGCGGCGGTCATGTCCGGCACGAAGTAGTAGAAACCACGGGTATCGCACACCAGTTCGAAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCTGGAAGTTCGACAGCTgg                                                                                                      <  3:245153/149‑1 (MQ=255)                                               gAACAGCGACAGGCGCTGCGCGGTCTTGTTGACCTGCGCGGCGGTCATGTCCGGCACGAAGTAGTAGAAACCACGGGTATCGCACACCAGTTCGAAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCTGGa                                                                                                                    <  7:168479/134‑1 (MQ=255)                                               gAACAGCGACAGGCGCTGCGCGGTCTTGTTGACCTGCGCGGCGGTCATGTCCGGCACGAAGTAGTAGAAACCACGGGTATCGCACACCAGTTCGAAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCTGGa                                                                                                                    >  8:168479/1‑134 (MQ=255)                                               gAACAGCGACAGGCGCTGCGCGGTCTTGTTGACCTGCGCGGCGGTCATGTCCGGCACGAAGTAGTAGAAACCACGGGTATCGCACACCAGTTCGAAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCTGGAAGTTCGACAGCTggg                                                                                                     >  4:173303/1‑149 (MQ=255)                                                           gcgcTGCGCGGTCTTGTTGACCTGCGCGGCGGTCATGTCCGGCACGAAGTAGTAGAAACCACGGGTATCGCACACCAGTTCGAAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCTGGAAGTTCGACAGCTGGGCGTACAGCTCAg                                                                                         >  5:123857/1‑149 (MQ=255)                                                            cgcTGCGCGGTCTTGTTGACCTGCGCGGCGGTCATGTCCGGCACGAAGTAGTAGAAACCACGGGTATCGCACACCAGTTCGAAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCTGGAAGTTCGACAGCTGGGCGTACAGCTCAgg                                                                                        <  2:40205/149‑1 (MQ=255)                                                                      tCTTGTTGACCTGCGCGGCGGTCATGTCCGGCACGAAGTAGTAGAAACCACGGGTATCGCACACCAGTTCGAAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCTGGAAGTTCGACAGCTGGGCGTACAGCTCAGGGTCGCGCCgg                                                                              <  2:234612/149‑1 (MQ=255)                                                                        ttgttgACCTGCGCGGCGGTCATGTCCGGCACGAAGTAGTAGAAACCACGGGTATCGCACACCAGTTCGAAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCTGGAAGTTCGACAGCTGGGCGTACAGCTCAGGGTCGCGCCGGCt                                                                            >  6:181838/1‑149 (MQ=255)                                                                              aCCTGCGCGGCGGTCATGTCCGGCACGAAGTAGTAGAAACCACGGGTATCGCACACCAGTTCGAAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCTGGAAGTTCGACAGCTGGGCGTACAGCTCAGGGTCGCGCCGGCTGACGt                                                                       <  1:242729/148‑1 (MQ=255)                                                                              aCCTGCGCGGCGGTCATGTCCGGCACGAAGTAGTAGAAACCACGGGTATCGCACACCAGTTCGAAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCTGGAAGTTCGACAGCTGGGCGTACAGCTCAGGGTCGCGCCGGCTGACGt                                                                       <  4:147269/148‑1 (MQ=255)                                                                                cTGCGCGGCGGTCATGTCCGGCACGAAGTAGTAGAAACCACGGGTATCGCACACCAGTTCGAAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCTGGAAGTTCGACAGCTGGGCGTACAGCTCAGGGTCGCGCCGGCTGACGTgg                                                                     >  4:125892/1‑148 (MQ=255)                                                                                cTGCGCGGCGGTCATGTCCGGCACGAAGTAGTAGAAACCACGGGTATCGCACACCAGTTCGAAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCTGGAAGTTCGACAGCTGGGCGTACAGCTCAGGGTCGCGCCGGCTGACGTgg                                                                     >  6:228628/1‑148 (MQ=255)                                                                                 tGCGCGGCGGTCATGTCCGGCGCGATGTAGTAGCAACGACGGGTATCGCCAACCTGTTCGGTGGGCAGGTGCGTGATGAGG‑GTGGGGTACTGGTCCTGGGTGATAGACAGCTGGGGGTGGAGGTGGGGGTCGGGGCGGCTGACGTgg                                                                     >  2:223027/1‑147 (MQ=255)                                                                                      ggcggTCATGTCCGGCACGATGTAGTAGAAACCACGGGTATCGCACACCAGTTCGAAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCTGGAAGTTCGACAGCTGGGCGTACAGCTCAGGGTCGCGCCGGCTGACGTGGAAGCCtt                                                              >  2:122645/1‑149 (MQ=255)                                                                                            cATGTCCGGCACGAAGTAGTAGAAACCACGGGTATCGCACACCAGTTCGAAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCTGGAAGTTCGACAGCTGGGCGTACAGCTCAGGGTCGCGCCGGCTGACGTGGAAGCCTTTGaa                                                          >  7:8241/1‑147 (MQ=255)                                                                                              tGTCCGGCACGAAGTAGTAGAAACCACGGGTATCGCACACCAGTTCGAAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCTGGAAGTTCGACAGCTGGGCGTACAGCTCAGGGTCGCGCCGGCTGACGTGGAAGCCTTTGAACAGc                                                      <  1:324463/149‑1 (MQ=255)                                                                                                tCCGGCACGAAGTAGTAGAAACCACGGGTATCGCACACCAGTTCGAAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCTGGAAGTTCGACAGCTGGGCGTACAGCTCAGGGTCGCGCCGGCTGACGTGGAAGCCTTTGAACAGCTc                                                    >  5:92366/1‑149 (MQ=255)                                                                                                          agtagtAGAAACCACGGGTATCGCACACCAGTTCGAAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCTGGAAGTTCGACAGCTGGGCGTACAGCTCAGGGTCGCGCCGGCTGACGTGGAAGCCTTTGAACAGCTCGCgg                                                >  4:61533/1‑143 (MQ=255)                                                                                                             agtagAAACCACGGGTATCGCACACCAGTTCGAAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCTGGAAGTTCGACAGCTGGGCGTACAGCTCAGGGTCGCGCCGGCTGACGTGGAAGCCTTTGAACAGCTCGCGGAAGATCGgc                                       >  3:217843/1‑149 (MQ=255)                                                                                                              gtagAAACCACGGGTATCGCACACCAGTTCGAAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCTGGAAGTTCGACAGCTGGGCGTACAGCTCAGGGTCGCGCCGGCTGACGTGGAAGCCTTTGAAc                                                         <  3:127179/130‑1 (MQ=255)                                                                                                              gtagAAACCACGGGTATCGCACACCAGTTCGAAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCTGGAAGTTCGACAGCTGGGCGTACAGCTCAGGGTCGCGCCGGCTGACGTGGAAGCCTTTGAAc                                                         >  4:127179/1‑130 (MQ=255)                                                                                                              gtagAAACCACGGGTATCGCACACCAGTTCGAAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCTGGAAGTTCGACAGCTGGGCGTACAGCTCAGGGTCGCGCCGGCTGACGTGGAAGCCTTTGAACAGCTCGCGGAAGATCGgcg                                      <  6:92366/149‑1 (MQ=255)                                                                                                                agAAACCACGGGTATCGCACACCAGTTCGAAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCTGGAAGTTCGACAGCTGGGCGTACAGCTCAGGGTCGCGCCGGCTGACGTGGAAGCCTTTGAACAGCTCGCGGAAGATCGgcgc                                     >  5:255597/1‑148 (MQ=255)                                                                                                                 gAAACCACGGGTATCGCACACCAGTTCGAAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCTGGAAGTTCGACAGCTGGGCGTACAGCTCAGGGTCGCGCCGGCTGACGTGGAAGCCTTTGAACAGCTCGCGGAAGATCGgcgc                                     >  3:145406/1‑147 (MQ=255)                                                                                                                    aCCACGGGTATCGCACACCAGTTCGAAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCTGGAAGTTCGACTGCTGGGCGTACAGCTCAGGGTCGCGCCGGCTGACGTGGAAGCCTTTGAACAGCTCGCGGAAGATCGGCGCGAGCt                                <  4:65299/149‑1 (MQ=255)                                                                                                                    aCCACGGGTATCGCACACCAGTTCGAAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCTGGAAGTTCGACAGCTGGGCGTACAGCTCAGGGTCGCGCCGGCTGACGTGGAAGCCTTTGAACAGCTCGCGGAAGATCGGCGCGAGCt                                >  5:23810/1‑149 (MQ=255)                                                                                                                      cACGGGTATCGCACACCAGTTCGAAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCTGGAAGTTCGACAGCTGGGCGTACAGCTCAGGGTCGCGCCGGCTGACGTGGAAGCCTTTGAACAGCTCGCGGAAGATCGGCGCGAGCTgg                              <  3:173303/149‑1 (MQ=255)                                                                                                                      cACGGGTATCGCACACCAGTTCGAAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGTACTGGAAGTTCGACAGCTGGGCGTACAGCTCAGGGTCGCGCCGGCTGACGTGGAAGCCTTTGAACAGCTCGCGGAATATCGGCGCGAGCTgg                              >  6:261871/1‑149 (MQ=255)                                                                                                                         gggTATCGCACACCAGTTCGAAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCTGGAAGTTCGACAGCTGGGCGTACAGCTCAGGGTCGCGCCGGCTGACGTGGAAGCCTTTGAACAGCTCGCGGAAGATCGGCGCGAGCTGGGa                            >  7:224855/1‑148 (MQ=255)                                                                                                                         gggTATCGCACACCAGTTCGAAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCTGGAAGTTCGACAGCTGGGCGTACAGCTCAGGGTCGCGCCGGCTGACGTGGAAGCCTTTGAACAGCTCGCGGAAGATCGGCGCGAGCTGGGAc                           <  6:74687/149‑1 (MQ=255)                                                                                                                         cagTACCGCAC‑CCTGTTCAAGGCCCTGGGCTTCGAACTGGTGTGCGGTACTGGTCCTGGAAGTTCGACAGCTGGGCGTACAGCTCAGGGTCGCGCCGGCTGACGTGGAAGCCTTTGAACAGCTCGCGGAAGATCGGcgcg                                    >  8:10509/3‑140 (MQ=25)                                                                                                                         cagTACCGCAC‑CCTGTTCAAGGCCCTGGGCTTCGAACTGGTGTGCGGTACTGGTCCTGGAAGTTCGACAGCTGGGCGTACAGCTCAGGGTCGCGCCGGCTGACGTGGAAGCCTTTGAACAGCTCGCGGAAGATCGGcgcg                                    <  7:10509/138‑1 (MQ=25)                                                                                                                         cagTACCGCAC‑CCTGTTCAAGGCCCTGGGCTTCGAACTGGTGTGCGGTACTGGTCCTGGAAGTTCGACAGCTGGGCGTACAGCTCAGGGTCGCGCCGGCTGACGTGGAAGCCTTTGAACAGCTCGCGGAAGATCGGCGCGAGCTGGGa                            >  3:140471/3‑148 (MQ=25)                                                                                                                                gCACACCAGTTCGAAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCTGGAAGTTCGACAGCTGGGCGTACAGCTCAGGGTCGCGCCGGCTGACGTGGAAGCCTTTGAACAGCTCGCGGAAGATCGGCGCGAGCTGGGACAGTTCgg                    <  4:217843/149‑1 (MQ=255)                                                                                                                                 cacaCCAGTTCGAAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCTGGAAGTTCGACAGCTGGGCGTACAGCTCAGGGTCGCGCCGGCTGACGTGGAAGCCTTTGAACAGCTCGCGGAAGATCGGCGCGAGCTGGGACAGTTCgg                    >  2:110660/1‑148 (MQ=255)                                                                                                                                 cacaCCAGTTCGAAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCTGGAAGTTCGACAGCTGGGCGTACAGCTCAGGGTCGCGCCGGCTGACGTGGAAGCCTTTGAACAGCTCGCGGAAGATCGGCGCGAGCTGGGACAGTTCGGa                   <  6:23810/149‑1 (MQ=255)                                                                                                                                  acacCAGTTCGAAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCTGGAAGTTCGACAGCTGGGCGTACAGCTCAGGGTCGCGCCGGCTGACGTGGAAGCCTTTGAACAGCTCGCGGAAGATCGGCGCGAGCTGGGACAGTTCGGaa                  <  5:181838/149‑1 (MQ=255)                                                                                                                                  acacCAGTTCGAAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCTGGAAGTTCGACAGCTGGGCGTACAGCTCAGGGTCGCGCCGGCTGACGTGGAAGCCTTTGAACAGCTCGCGGAAGATCGGCGCCAGCTGGGACACTACGGa                   >  6:236010/1‑148 (MQ=255)                                                                                                                                   cacCAGTTCGAAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCTGGAAGTTCGACAGCTGGGCGTACAGCTCAGGGTCGCGCCGGCTGACGTGGAAGCCTTTGAACAGCTCGCGGAAGATCGGCGCGAGCTGGGACAGTTCGGaaa                 <  8:8241/149‑1 (MQ=255)                                                                                                                                     ccAGTTCGAAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCTGGAAGTTCGACAGCTGGGCGTACAGCTCAGGGTCGCGCCGGCTGACGTGGAAGCCTTTGAACAGCTCGCGGAAGATCGGCGCGAGCTgg                              <  7:25003/134‑1 (MQ=255)                                                                                                                                     ccAGTTCGAAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCTGGAAGTTCGACAGCTGGGCGTACAGCTCAGGGTCGCGCCGGCTGACGTGGAAGCCTTTGAACAGCTCGCGGAAGATCGGCGCGAGCTgg                              >  8:25003/1‑134 (MQ=255)                                                                                                                                            gAAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCTGGAAGTTCGACAGCTGGGCGTACAGCTCAGGGTCGCGCCGGCTGACGTGGAAGCCTTTGAACAGCTCGCGGAAGATCGGCGCGAGCTGGGACAGTTCGGAAAGATCAAGAt        >  7:173536/1‑149 (MQ=255)                                                                                                                                             aggCCCTGGGCTTCGAACTGGTGTGCGGTACTGGTCCTGGAAGTTCGACAGCTGGGCGTACAGCTCAGGGTCGCGCCGGCTGACGTGGAAGCCTTTGAACAGCTCGCGGAAGATCGGCGCGAGCTGGGACAGTTCGGAAAGATCAAGAt        <  4:140471/147‑1 (MQ=255)                                                                                                                                                cccAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCTGGAAGTTCGTCAGCTGTGCGTACAGCTCAGGGTCGCGCCGGCTGACGTGGAAGCCTTTGAACAGCTCGCGGAAGATCGGCGCTAGCTGGGACAGTTCGGAAAGATCAAGATGCAt    <  7:207/149‑1 (MQ=255)                                                                                                                                                cccAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCTGGAAGTTCGACAGCTGGGCGTACAGCTCAGGGTCGCGCCGGCTGACGTGGAAGCCTTTGAACAGCTCGCGGAAGATCGGCGCGAGCTGGGACAGTTCGGAAAGATCAAGATGCAt    <  3:61533/149‑1 (MQ=255)                                                                                                                                                  cAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCTGGAAGTTCGACAGCTGGGCGTACAGCTCAGGGTCGCGCCGGCTGACGTGGAAGCCTTTGAACAGCTCGCGGAAGATCGGCGCGAGCTGGGACAGTTCGGAAAGATCAAGATGCAtgg  >  4:14780/1‑149 (MQ=255)                                                                                                                                                   |                                                                                                                                                     TCGGGTCGCGGCCCTGGTCAGCAAGGTGCTCGACCAGGATGAAGGTGAACAGCGACAGGCGCTGCGCGGTCTTGTTGACCTGCGCGGCGGTCATGTCCGGCACGAAGTAGTAGAAACCACGGGTATCGCACACCAGTTCGAAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCTGGAAGTTCGACAGCTGGGCGTACAGCTCAGGGTCGCGCCGGCTGACGTGGAAGCCTTTGAACAGCTCGCGGAAGATCGGCGCGAGCTGGGACAGTTCGGAAAGATCAAGATGCATGG  >  NC_002947/807742‑808035

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