Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	807,893	C→T	6.4%	K47K (AAG→AAA)	PP_0694 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	807,893	0	C	T	6.4%	184.4 / 2.5	62	K47K (AAG→AAA)	PP_0694	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (32/26);  new base T (2/2);  total (34/28)
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.64e-01

GGTCGCGGCCCTGGTCAGCAAGGTGCTCGACCAGGATGAAGGTGAACAGCGACAGGCGCTGCGCGGTCTTGTTGACCTGCGCGGCGGTCATGTCCGGCACGAAGTAGTAGAAACCACGGGTATCGCACACCAGTTCGAAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCTGGAAGTTCGACAGCTGGGCGTACAGCTCAGGGTCGCGCCGGCTGACGTGGAAGCCTTTGAACAGCTCGCGGAAGATCGGCGCGAGCTGGGACAGTTCGGAAAGATCAAGATGCATGG  >  NC_002947/807745‑808035                                                                                                                                                     |                                                                                                                                                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                                                                                                                          >  8:207/1‑149 (MQ=255)                       gTGCTCGACCAGGATGAAGGTGAACAGCGACAGGCGCTGCGCGGTCTTGTTGACCTGCGCGGCGGTCATGTCCGGCACGAAGTAGTAGAAACCACGGGTATCGCACACCAGTTCGAAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGt                                                                                                                          >  1:40205/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                                                                                                                    >  8:168479/1‑134 (MQ=255)                                            gAACAGCGACAGGCGCTGCGCGGTCTTGTTGACCTGCGCGGCGGTCATGTCCGGCACGAAGTAGTAGAAACCACGGGTATCGCACACCAGTTCGAAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCTGGa                                                                                                                    <  7:168479/134‑1 (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                                                                       <  4:147269/148‑1 (MQ=255)                                                                           aCCTGCGCGGCGGTCATGTCCGGCACGAAGTAGTAGAAACCACGGGTATCGCACACCAGTTCGAAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCTGGAAGTTCGACAGCTGGGCGTACAGCTCAGGGTCGCGCCGGCTGACGt                                                                       <  1:242729/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                                                         >  4:127179/1‑130 (MQ=255)                                                                                                           gtagAAACCACGGGTATCGCACACCAGTTCGAAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCTGGAAGTTCGACAGCTGGGCGTACAGCTCAGGGTCGCGCCGGCTGACGTGGAAGCCTTTGAAc                                                         <  3:127179/130‑1 (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                                    <  7:10509/138‑1 (MQ=25)                                                                                                                      cagTACCGCAC‑CCTGTTCAAGGCCCTGGGCTTCGAACTGGTGTGCGGTACTGGTCCTGGAAGTTCGACAGCTGGGCGTACAGCTCAGGGTCGCGCCGGCTGACGTGGAAGCCTTTGAACAGCTCGCGGAAGATCGGcgcg                                    >  8:10509/3‑140 (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)                                                                                                                                                     |                                                                                                                                                GGTCGCGGCCCTGGTCAGCAAGGTGCTCGACCAGGATGAAGGTGAACAGCGACAGGCGCTGCGCGGTCTTGTTGACCTGCGCGGCGGTCATGTCCGGCACGAAGTAGTAGAAACCACGGGTATCGCACACCAGTTCGAAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCTGGAAGTTCGACAGCTGGGCGTACAGCTCAGGGTCGCGCCGGCTGACGTGGAAGCCTTTGAACAGCTCGCGGAAGATCGGCGCGAGCTGGGACAGTTCGGAAAGATCAAGATGCATGG  >  NC_002947/807745‑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: ‑