Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,798,967	A→C	5.7%	L299W (TTG→TGG)	PP_5077 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,798,967	0	A	C	5.7%	281.1 / 6.9	88	L299W (TTG→TGG)	PP_5077	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (44/39);  new base C (2/3);  total (46/42)
Fisher's exact test for biased strand distribution p-value = 6.66e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.42e-01

TCCCTCTTCACCCTCGCCCATGCCTGCGGCTTGGGCCAGCGGCTCGCGCATCACCGGCTGCGACTGGCCAACCAGTGGCAGCGGCATTGGCTGGGACGAGCCAGAGAACTCGATGGCAGCACCCTGCTTGCCGTCGCCCAGTGGCAACTGGGCCTGGGCAGCCGGC‑GCTTCGGCA‑GGGGCCTTGTCGCCTTTCTTGGGCATCAATACCGCGGCCGCAACGGCAACCACGACCACAGCGGACAGCGCGAGCACATGTTTTTTAGGCATCTTGAACCCCATGGTTGGTCGCTTGG  >  NC_002947/5798822‑5799114                                                                                                                                                  |                                                                                                                                                      tccctcTTCACCCTCGCCCATGCCTGCGGCTTGGGCCAGCGGCTCGCGCATCACCGGCTGCGACTGGCCAACCAGTGGCAGCGGCATTGGCTGGGACGAGCCAGAGAACTCGATGGCAGCACCCTGCTTGCCGTCGCCCAGTGGCAAct                                                                                                                                                    >  4:209874/1‑149 (MQ=255)       ttCACCCTCGCCCATGCCTGCGGCTTGGGCCAGCGGCTCGCGCATCACCGGCTGCGACTGGCCAACCAGTGGCAGCGGCATTGGCTGGGACGAGCCAGAGAACTCGATGGCAGCACCCTGCTTGCCGTCGCCCAGTGGCAActgggcct                                                                                                                                              >  7:225791/1‑149 (MQ=255)         cACCCTCGCCCATGCCTGCGGCTTGGGCCAGCGGCTCGCGCATCACCGGCTGCGACTGGCCAACCAGTGGCAGCGGCATTGGCTGGGACGAGCCAGAGAACTCGATGGCAGCACCCTGCTTGCCGTCGCCCAGTTGCCActgggcga                                                                                                                                              <  1:222656/147‑3 (MQ=255)         cACCCTCGCCCATGCCTGCGGCTTGGGCCAGCGGCTCGCGCATCACCGGCTGCGACTGGCCAACCAGTGGCAGCGGCATTGGCTGGGACGAGCCAGAGAACTCGATGGCAGCACCCTGCTTGCCGTCGCCCAGTGGCAActgggcctgg                                                                                                                                            >  2:225774/1‑149 (MQ=255)          aCCCTCGCCCATGCCTGCGGCTTGGGCCAGCGGCTCGCGCATCACCGGCTGCGACTGGCCAACCAGTGGCAGCGGCATTGGCTGGGACGAGCCAGAGAACTCGATGGCAGCACCCTGCTTGCCGTCGCCCAGTGGCAActgggcctggg                                                                                                                                           >  2:49013/1‑149 (MQ=255)            ccTCGCCCATGCCTGCGGCTTGGGCCAGCGGCTCGCGCATCACCGGCTGCGACTGGCCAACCAGTGGCAGCGGCATTGGCTGGGACGAGCCAGAGAACTCGATGGCAGCACCCTGCTTGCCGTCGCCCAGTGGCAACTGGGCCTGGGCa                                                                                                                                         <  8:225791/149‑1 (MQ=255)              tCGCCCATGCCTGCGGCTTGGGCCAGCGGCTCGCGCATCACCGGCTGCGACTGGCCAACCAGTGGCAGCGGCATTGGCTGGGACGAGCCAGAGAACTCGATGGCAGCACCCTGCTTGCCGTCGCCCAGTGGCAACTGGGCCTGGGCAGc                                                                                                                                       >  1:110268/1‑149 (MQ=255)                  ccATGCCTGCGGCTTGGGCCAGCGGCTCGCGCATCACCGGCTGCGACTGGCCAACCAGTGGCAGCGGCATTGGCTGGGACGAGCCAGAGAACTCGATGGCAGCACCCTGCTTGCCGTCGCCCAGTGGCAACTGGGCCTGGGCAGCCGgc                                                                                                                                   >  5:37509/1‑149 (MQ=255)                  ccATGCCTGCGGCTTGGGCCAGCGGCTCGCGCATCACCGGCTGCGACTGGCCAACCAGTGGCAGCGGCATTGGCTGGGACGAGCCAGAGAACTCGATGGCAGCACCCTGCTTGCCGTCGCCCAGTGGCAACTGGGCCTGGGCAGCCGgc                                                                                                                                   >  5:97380/1‑149 (MQ=255)                       ccTGCGGCTTGGGCCAGCGGCTCGCGCATCACCGGCTGCGACTGGCCAACCAGTGGCAGCGGCATTGGCTGGGACGAGCCAGAGAACTCGATGGCAGCACCCTGCTTGCCGTCGCCCAGTGGCAACTGGGCCTGGGCAGCCGGC‑GCTTc                                                                                                                             >  3:355270/1‑149 (MQ=255)                        cTGCGGCTTGGGCCAGCGGCTCGCGCATCACCGGCTGCGACTGGCCAACCAGTGGCAGCGGCATTGGCTGGGACGAGCCAGAGAACTCGATGGCAGCACCCTGCTTGCCGTCGCCCAGTGGCAACTGGGCCTGGGCAGCCGGC‑GCTTCg                                                                                                                            >  3:318119/1‑149 (MQ=255)                              cTTGGGCCAGCGGCTCGCGCATCACCGGCTGCGACTGGCCAACCAGTGGCAGCGGCATTGGCTGGGACGAGCCAGAGAACTCGATGGCAGCACCCTGCTTGCCGTCGCCCAGTGGCAACTGGGCCTGGGCAGCCGGC‑GCTTCGGCA‑ggg                                                                                                                     <  2:374267/149‑1 (MQ=255)                              cTTGGGCCAGCGGCTCGCGCATCACCGGCTGCGACTGGCCAACCAGTGGCAGCGGCATTGGCTGGGACGAGCCAGAGAACTCGATGGCAGCACCCTGCTTGCCGTCGCCCAGTGGCAACTGGGCCTGGGCAGCCGGC‑GCTTCGGCA‑ggg                                                                                                                     <  8:212457/149‑1 (MQ=255)                                tGGGCCAGCGGCTCGCGCATCACCGGCTGCGACTGGCCAACCAGTGGCAGCGGCATTGGCTGGGACGAGCCAGAGAACTCGATGGCAGCACCCTGCTTGCCGTCGCCCAGTGGCAACTGGGCCTGGGCAGCCGGC‑GCTTCGGCA‑GGGGc                                                                                                                   <  3:87072/149‑1 (MQ=255)                                 tggCCAGCGGCTCGCGCATCACCGGCTGCGACTGGCCAACCAGTGGCAGCGGCATTGGCTGGGACGAGCCAGAGAACTCGATGGCAGCACCCTGCTTGCCGTCGCCCAGTGGCAACTGGGCCTGGGCAGCCGGC‑GCTTCGGCA‑GGGGcc                                                                                                                  >  4:10599/2‑149 (MQ=255)                                   gCCAGCGGCTCGCGCATCACCGGCTGCGACTGGCCAACCAGTGGCAGCGGCATTGGCTGGGACGAGCCAGAGAACTCGATGGCAGCACCCTGCTTGCCGTCGCCCAGTGGCAACTGGGCCTGGGCAGCCGGC‑GCTTCGGCA‑GGGGCCtt                                                                                                                <  6:97380/149‑1 (MQ=255)                                       gCGGCTCGCGCATCACCGGCTGCGACTGGCCAACCAGTGGCAGCGGCATTGGCTGGGACGAGCCAGAGAACTCGATGGCAGCACCCTGCTTGCCGTCGCCCAGTGGCAACTGGGCCTGGGCAGCCGGC‑GCTTCGGCA‑GGGGCCTTGTCg                                                                                                            <  3:224564/149‑1 (MQ=255)                                            tCGCGCATCACCGGCTGCGACTGGCCAACCAGTGGCAGCGGCATTGGCTGGGACGAGCCAGAGACCTCGATGGCAGCACCCTGCTTGCCGTCGCCCAGTGGCAACTGGGCCTGGGCAGCCGGC‑GCTTCGGCA‑GGGGCCTTGTCGCCtgt                                                                                                       >  7:345594/1‑147 (MQ=255)                                            tCGCGCATCACCGGCTGCGACTGGCCAACCAGTGGCAGCGGCATTGGCTGGGACGAGCCAGAGAACTCGATGGCAGCACCCTGCTTGCCGTCGCCCAGTGGCAACTGGGCCTGGGCAGCCGGC‑GCTTCGGCA‑ggg                                                                                                                     >  5:265481/1‑135 (MQ=255)                                            tCGCGCATCACCGGCTGCGACTGGCCAACCAGTGGCAGCGGCATTGGCTGGGACGAGCCAGAGAACTCGATGGCAGCACCCTGCTTGCCGTCGCCCAGTGGCAACTGGGCCTGGGCAGCCGGC‑GCTTCGGCA‑ggg                                                                                                                     <  6:265481/135‑1 (MQ=255)                                              gcgcATCACCGGCTGCGACTGGCCAACCAGTGGCAGCGGCATTGGCTGGGACGAGCCAGAGAACTCGATGGCAGCACCCTGCTTGCCGTCGCCCAGTGGCAACTGGGCCTGGGCAGCCGGC‑GCTTCGGCA‑GGGGCCTTGTCGCCTTTCt                                                                                                     >  2:320652/1‑149 (MQ=255)                                              gcgcATCACCGGCTGCGACTGGCAAACCAGTGGCAGCGGCATTGGCTGGGACGAGCCAGAGAACTCAATGGCAGCACCCTGCTTGCCGTCGCCCAGTGGCAACTGGGCCTGGGCAGCCGGC‑GCTTCGGCA‑GGGGCCTTGTCGCCTTTCt                                                                                                     >  4:236568/1‑149 (MQ=255)                                               cgcATCACCGGCTGCGACTGGCCAACCAGTGGCAGCGGCATTGGCTGGGACGAGCCAGAGAACTCGATGGCAGCACCCTGCTTGCCGTCGCCCAGTGGCAACTGGGCCTGGGCAGCCGGC‑GCTTCGGCA‑GGGGCCTTGTCGCCTTTCtt                                                                                                    <  2:110268/149‑1 (MQ=255)                                               cgcATCACCGGCTGCGACTGGCCAACCAGTGGCAGCGGCATTGGCTGGGACGAGCCAGAGAACTCGATGGCAGCACCCTGCTTGCCGTCGCCCAGTGGCAACTGGGCCTGGGCAGCCGGC‑GCTTCGGCA‑GGGGCCTTGTCGCCTTTCtt                                                                                                    <  7:363918/149‑1 (MQ=255)                                                gcATCACCGGCTGCGACTGGCCAACCAGTGGCAGCGGCATTGGCTGGGACGAGCCAGAGAACTCGATGGCAGCACCCTGCTTGCCGTCGCCCAGTGGCAACTGGGCCTGGGCAGCCGGC‑GCTTCGGCA‑GGGGCCTTGTCGCCTTTCTTg                                                                                                   <  5:264589/149‑1 (MQ=255)                                                gcATCACCGGCTGCGACTGGCCAACCAGTGGCAGCGGCATTGGCTGGGACGAGCCAGAGAACTCGATGGCAGCACCCTGCTTGCCGTCGCCCAGTGGCAACTGGGCCTGGGCAGCCGGC‑GCTTCGGCA‑GGGGCCTTGTCGCCTTTCTTg                                                                                                   >  4:80762/1‑149 (MQ=255)                                                 cATCACCGGCTGCGACTGGCCAACCAGTGGCAGCGGCATTGGCTGGGACGAGCCAGAGAACTCGATGGCAGCACCCTGCTTGCCGTCGCCCAGTGGCAACTGGGCCTGGGCAGCCGGC‑GCTTCGGCA‑GGGGCCTTGTCGCCTTTCTTgg                                                                                                  >  8:336512/1‑149 (MQ=255)                                                  aTCACCGGCTGCGACTGGCCAACCAGTGGCAGCGGCATTGGCTGGGACGAGCCAGAGAACTCGATGGCAGCACCCTGCTTGCCGTCGCCCAGTGGCAACTGGGCCTGGGCAGCCGGC‑GCTTCGGCA‑GGGGCCTTGTCGCCTTTCTTggg                                                                                                 >  5:139218/1‑149 (MQ=255)                                                  aTCACCGGCTGCGACTGGCCAACCAGTGGCAGCGGCATTGGCTGGGACGAGCCAGAGAACTCGATGGCAGCACCCTGCTTGCCGTCGCCCAGTGGCAACTGGGCCTGGGCAGCCGGC‑GCTTCGGCA‑GGGGCCTTGTCGCCTTTCTTggg                                                                                                 <  6:139218/149‑1 (MQ=255)                                                   tCACCGGCTGCGACTGGCCAACCAGTGGCAGCGGCATTGGCTGGGACGAGCCAGAGAACTCGATGGCAGCACCCTGCTTGCCGTCGCCCAGTGGCAACTGGGCCTGGGCAGCCGGC‑GCTTCGGCA‑GGGGCCTTGTCGCCTTTCTTggg                                                                                                 >  2:94995/1‑148 (MQ=255)                                                   tCACCGGCTGCGACTGGCCAACCAGTGGCAGCGGCATTGGCTGGGACGAGCCAGAGAACTCGATGGCAGCACCCTGCTTGCCGTCGCCCAGTGGCAACTGGGCCTGGGCAGCCGGC‑GCTTCGGCA‑GGGGCCTTGTCGCCTTTCTTGGGc                                                                                                <  4:50808/149‑1 (MQ=255)                                                      acGGCTGCGACTGGCCAACCAGGGGCAGCGGCATTGGCTGGGACGAGCCAGAGAACTCGATGGCAGCACCCTGCTTGCCGTCGCCCAGTGGCAACTGGGCCTGGGCAGCCGGC‑GCTTCGGCA‑GGGGCCTTGTCGCCTTTCTTGGGCATc                                                                                             <  4:27422/148‑1 (MQ=255)                                                       cGGCTGCGACTGGCCAACCAGTGGCAGCGGCATTGGCTGGGACGAGCCAGAGAACTCGATGGCAGCACCCTGCTTGCCGTCGCCCAGTGGCAACTGGGCCTGGGCAGCCGGC‑GCTTCGGCACGGGGCCTTGTCGCCTTTCTTGGGCATc                                                                                             <  1:225774/149‑1 (MQ=255)                                                       cGGCTGCGACTGGCCAACCAGTGGCAGCGGCATTGGCTGGGACGAGCCAGAGAACTCGATGGCAGCACCCTGCTTGCCGTCGCCCAGTGGCAACTGGGCCTGGGCAGCCGGC‑GCTTCGGCA‑GGGGCCTTGTCACCTTTCTTGGGCATCa                                                                                            <  3:209874/149‑1 (MQ=255)                                                         gCTGCGACTGGCCAACCAGTGGCAGCGGCATTGGCTGGGACGAGCCAGAGAACTCGATGGCAGCACCCTGCTTGCCGTCGCCCAGTGGCAACTGGGCCTGGGCAGCCGGC‑GCTTCGGCA‑GGGGCCTTGTCGCCTTTCTTGGGCATCaa                                                                                           >  3:55243/1‑148 (MQ=255)                                                          cTGCGACTGGCCAACCAGTGGCAGCGGCATTGGCTGGGACGAGCCAGAGAACTCGATGGCAGCACCCTGCTTGCCGTCGCCCAGTGGCAACTGGGCCTGGGCAGCCGGC‑GCTTCGGCA‑GGGGCCTTGTCGCCTTTCTTGGGCATCAATa                                                                                         <  5:173044/149‑1 (MQ=255)                                                              gACTTGCCAACCAGTGGCAGCGGCATTGGCTGGGACGAGCCAGAGAACTCGATGGCAGCACCCTGCTTGCCGTCGCCCAGTTGCCActgggcgac                                                                                                                                             <  7:342085/95‑4 (MQ=255)                                                              gACTGGCCAACCAGTGGCAGCGGCATTGGCTGGGACGAGCCAGAGAACTCGATGGCAGCACCCTGCTTGCCGTCGCCCAGTTGCCActgggcgac                                                                                                                                             <  7:336339/95‑4 (MQ=255)                                                              gACTGGCCAACCAGTGGCAGCGGCATTGGCTGGGACGAGCCAGAGAACTCGATGGCAGCACCCTGCTTGCCGTCGCCCAGTGGCAACTGGGCCTGGGCAGCCGGC‑GCTTCGGCA‑GGGGCCTTGTCGCCTTTCTTGGGCATCAATACcgc                                                                                     >  7:24675/1‑149 (MQ=255)                                                              gACTGGCCAACCAGTGGCAGCGGCATTGGCTGGGACGAGCCAGAGAACTCGATGGCAGCACCCTGCTTGCCGTCGCCCAGTGGCAACTGGGCCTGGGCAGCCGGC‑GCTTCGGCA‑GGGGCCTTGTCGCCTTTCTTGGGCATCAATACcgc                                                                                     >  5:231197/1‑149 (MQ=255)                                                                 tGGCCAACCAGTGGCAGCGGCATTGGCTGGGACGAGCCAGAGAACTCGATGGCAGCACCCTGCTTGCCGTCGCCCAGTGGCAACTGGGCCTGGGCAGCCGGC‑GCTTCGGCA‑GGGGCCTTGTCGCCGTTCTTGGGCATCAATTCCGCGGc                                                                                  >  6:289615/1‑149 (MQ=255)                                                                   gCCAACCAGTGGCAGCGGCATTGGCTGGGACGAGCCAGAGAACTCGATGGCAGCACCCTGCTTGCCGTCGCCCAGTGGCAACTGGGCCTGGGCAGCCGGC‑GCTTCGGCA‑GGGGCCTTGTCGCCTTTCTTGGGCATCAATACCGCGGc                                                                                  <  8:24675/147‑1 (MQ=255)                                                                      aaCCAGTGGCAGCGGCATTGGCTGGGACGAGCCAGAGAACTCGATGGCAGCACCCTGCTTGCCGTCGCCCAGTGGCAACTGGGCCTGGGCAGCTGGC‑GCTTCGGCA‑GGGGCCTTGTCGCCTTTCTTGGGCATCAATACCGCGGCCGCaa                                                                             >  3:216/1‑149 (MQ=255)                                                                         cAGTGGCAGCGGCATTGGCTGGGACGAGCCAGAGAACTCGATGGCAGCACCCTGCTTGCCGTCGCCCAGTGGCAACTGGGCCTGGGCAGCCGGC‑GCTTCGGCA‑GGGGCCTTGTCGCCTTTCTTGGGCATc                                                                                             <  6:403044/130‑1 (MQ=255)                                                                              gCAGCGGCATTGGCTGGGACGAGCCAGAGAACTCGATGGCAGCACCCTGCTTGCCGTCGCCCAGTGGCAACTGGGCCTGGGCAGCCGGC‑GCTTCGGCA‑GGGGCCtt                                                                                                                <  6:39140/106‑1 (MQ=255)                                                                              gCAGCGGCATTGGCTGGGACGAGCCAGAGAACTCGATGGCAGCACCCTGCTTGCCGTCGCCCAGTGGCAACTGGGCCTGGGCAGCCGGC‑GCTTCGGCA‑GGGGCCtt                                                                                                                >  5:39140/1‑106 (MQ=255)                                                                                  cGGCATTGGCTGGGACGAGCCAGAGAACTCGATGGCAGCACCCTGCTTGCCGTCGCCCAGTGGCAACTGGGCCTGGGCAGCCGGC‑GCTTCGGCA‑GGGGCCTTGTCGCCTTTCTTGTGCATTAATACCGCGGCCGCAACTGCAACCACGa                                                                 <  6:259781/149‑1 (MQ=255)                                                                                    gCATTGGCTGGGACGAGCCAGAGAACTCGATGGCAGCACCCTGCTTGCCGTCGCCCAGTGGCAACTGGGCCTGGGCAGCCGGC‑GCTTCGGCA‑GGGGCCTTGTCGCCTTTCTTGGGCATCAATACCGCGGCCGCAACGGCAACCACGAc                                                                >  3:246564/1‑148 (MQ=255)                                                                                          gCTGGGACGAGCCAGAGAACTCGATGGCAGCACCCTGCTTGCCGTCGCCCAGTGGCAACTGGGCCTGGGCAGCCGGC‑TCTTCGGCA‑GGGGCGTTGTCGCCTTTCTTGGGCATCAATACCGCGGCCGCAACGGCAACCACGACCACAGCg                                                         <  5:289615/149‑1 (MQ=255)                                                                                                    gCCAGAGAACTCGATGGCAGCACCCTGCTTGCCGTCGCCCAGTGGCAACTGGGCCTGGGCAGCCGGC‑GCTTCGGCA‑GGGGCCTTGTCGCCTTTCTTGGGCATCAATACCGCGGCCGCAACGGCAACCACGACcaca                                                            >  7:303907/1‑136 (MQ=255)                                                                                                    gCCAGAGAACTCGATGGCAGCACCCTGCTTGCCGTCGCCCAGTGGCAACTGGGCCTGGGCAGCCGGC‑GCTTCGGCA‑GGGGCCTTGTCGCCTTTCTTGGGCATCAATACCGCGGCCGCAACGGCAACCACGACcaca                                                            <  8:303907/136‑1 (MQ=255)                                                                                                    gCCAGAGAACTCGATGGCAGCACCCTGCTTGCCGTCGCCCAGTGGCAACTGGGCCTGGGCAGCCGGC‑GCTTCGGCA‑GGGGCCTTGTCGCCTTTCTTGGGCATCAATACCGCGGCCGCAACGGCAACCACGACCACAGCGGACAGCGCGa                                               <  4:246564/149‑1 (MQ=255)                                                                                                     ccAGAGAACTCGATGGCAGCACCCTGCTTGCCGTCGCCCAGTGGCAACTGGGCCTGGGCAGCCGGC‑GCTTCGGCA‑GGGGCCTTGTCGCCTTTCTTGGGCATCAATACCGCGGCCGCAACGGCAACCACGACCACAGCGGACAGCGCGAg                                              >  1:216824/1‑149 (MQ=255)                                                                                                       agagAACTCGATGGCAGCACCCTGCTTGCCGTCGCCCAGTGGCAACTGGGCCTGGGCAGCCGGC‑GCTTCGGCA‑GGGGCCTTGTCGCCTTTCTTGGGCATCAATACCGCGGCCGCAACGGCAACCACGACCACAGCGGACAGCGCGAGca                                            >  1:368686/1‑149 (MQ=255)                                                                                                             cTCGATGGCAGCACCCTGCTTGCCGTCGCCCAGTGGCAACTGGGCCTGGGCAGCCGGC‑GCTTCGGCA‑GGGGCCTTGTCGCCTTTCTTGGGCATCAATACCGCGGCCGCAACGGCAACCACGACCACAGCGGACAGCGCGAGCACATGtt                                      <  1:49013/149‑1 (MQ=255)                                                                                                             cTCGATGGCAGCACCCTGCTTGCCGTCGCCCAGTGGCAACTGGGCCTGGGCAGCCGGC‑GCTTCGGCA‑GGGGCCTTGTCGCCTTTCTTGGGCATCAATACCGCGGCCGCAACGGCAACCACGACCACAGCGGACAGCGCGAGCACATGtt                                      >  6:178550/1‑149 (MQ=255)                                                                                                              tCGATGGCAGCACCCTGCTTGCCGTCGCCCAGTGGCAACTGGGCCTGGGCAGCCGGC‑GCTTCGGCA‑GGGGCCTTGTCGCCTTTCTTGGGCATCAATACCGCGGCCGCAACGGCAACCACGACCACAGCGGACAGCGCGAGCACATGttt                                     >  4:23433/1‑149 (MQ=255)                                                                                                                gATGGCAGCACCCTGCTTGCCGTCGCCCAGTGGCAACTGGGCCTGGGCAGCCGGC‑GCTTCGGCA‑GGGGCCTTGTCGCCTTTCTTGGGCATCAATACCGCGGCCGCAACGGCAACCACGACCACAGCGGACAGCGCGAGCACATGttttt                                   >  2:171938/1‑149 (MQ=255)                                                                                                                gATGGCAGCACCCTGCTTGCCGTCGCCCAGTGGCAACTGGGCCTGGGCAGCCGGC‑GCTTCGGCA‑GGGGCCTTGTCGCCTTTCTTGGGCATCAATACCGCGGCCGCAACGGCAACCACGACCACAGCGGACAGCGCGAGCACATGttttt                                   >  4:171787/1‑149 (MQ=255)                                                                                                                gATGGCAGCACCCTGCTTGCCGTCGCCCAGTGGCAACTGGGCCTGGGCAGCCGGC‑GCTTCGGCA‑GGGGCCTTGTCGCCTTTCTTGGGCATCAATACCGCGGCCGCAACGGCAACCACGACCACAGCGGACAGCGCGAGCACATGtttt                                    <  2:216824/148‑1 (MQ=255)                                                                                                                gATGGCAGCACCCTGCTTGCCGTCGCCCAGTGGCAACTGGGCCTGGGCAGCCGGC‑GCTTCGGCA‑GGGGCCTTGTCGCCTTTCTTGGGCATCAATACCGCGGCCGCAACGGCAACCACGACCACAGCGGACAGCGCGAGCACATGtttt                                    <  6:37509/148‑1 (MQ=255)                                                                                                                 aTGGCAGCACCCTGCTTGCCGTCGCCCAGTGGCAACTGGGCCTGGGCAGCCGGC‑GCtt                                                                                                                              <  5:56236/58‑1 (MQ=255)                                                                                                                 aTGGCAGCACCCTGCTTGCCGTCGCCCAGTGGCAACTGGGCCTGGGCAGCCGGC‑GCtt                                                                                                                              >  6:56236/1‑58 (MQ=255)                                                                                                                 aTGGCAGCACCCTGCTTGCCGTCGCCCAGTGGCAACTGGGCCTGGGCAGCCGGC‑GCTTCGGCA‑GGGGCCTTGTCGCCTTTCTTGGGCATCAATACCGCGGCCGCAACGGCAACCACGACCACAGCGGACAGCGCGAGCACATGtttttt                                  <  2:368686/149‑1 (MQ=255)                                                                                                                    gcagcaCCCTGCTTGCCGTCGCCCAGTGGCAACTGGGCCTGGGCAGCCGGC‑GCTTCGGCA‑GGGGCCTTGTCGCCTTTCTTGGGCATCAATACCGCGGCCGCAACGGCAACCACGACCACAGCGGACAGCGCGAGCACATGTTTTTTAgg                               >  2:329655/1‑149 (MQ=255)                                                                                                                    gcagcaCCCTGCTTGCCGTCGCCCAGTGGCAACTGGGCCTGGGCAGCCGGC‑GCTTCGGCA‑GGGGCCTTGTCGCCTTTCTTGGGCATCAATACCGCGGCCGCAACGGCAACCACGACCACAGCGGACAGCGCGAGCACATGTTTTTTAg                                <  7:190068/148‑1 (MQ=255)                                                                                                                    gcagcaCCCTGCTTGCCGTCGCCCAGTGGCAACTGGGCCTGGGCAGCCGGC‑GCTTCGGCA‑GGGGCCTTGTCGCCTTTCTTGGGCATCAATACCGCGGCCGCAACGGCAACCACGACCACAGCGGACAGCGCGAGCACATGTTTTTTAg                                <  3:236568/148‑1 (MQ=255)                                                                                                                    gcagcaCCCTGCTTGCCGTCGCCCAGTGGCAACTGGGCCTGGGCAGCCGGC‑GCTTCGGCA‑GGGGCCTTGTCGCCTTTCTTGGGCATCAATACCGCGGCCGCAACGGCAACCACGACCACAGCGGACAGCGCGAGCACATGTTTTTTAg                                <  3:80762/148‑1 (MQ=255)                                                                                                                          cccTGCTTGCCGTCGCCCAGTGGCAACTGGGCCTGGGCAGCCGGC‑GCTTCGGCA‑GGGGCCTTGTCGCCTTTCTTGGGCATCAATACCGCGGCCGCAACGGCAACCACGACCACAGCGGACAGCGCGAGCACATGTTTTTTAGGCATCt                          <  3:84457/148‑1 (MQ=255)                                                                                                                          cccTGCTTGCCGTCGCCCAGTGGCAACTGGGCCTGGGCAGCCGGC‑GCTTCGGCA‑GGGGCCTTGTCGCCTTTCTTGGGCATCAATACCGCGGCCGCAACGGCAACCACGACCACAGCGGACAGCGCGAGCACATGTTTTTTAGGCATCt                          >  4:84457/1‑148 (MQ=255)                                                                                                                             tGCTTGCCGTCGCCCAGTGGCAACTGGGCCTGGGCAGCCGGC‑GCTTCGGCA‑GGGGCCTTGTCGCCTTTCTTGGGCATCAATACCGCGGCCGCAACGGCAACCACGACCACAGCGGACAGCGCGAGCACATGTTTTTTAGGCATCTTGaa                      >  5:96083/1‑149 (MQ=255)                                                                                                                              gCTTGCCGTCGCCCAGTGGCAACTGGGCCTGGGCAGCCGGC‑GCTTCGGCA‑GGGGCCTTGTCGCCTTTCTTGGGCATCAATACCGCGGCCGCAACGGCAACCACGACCACAGCGGACAGCGCGAGCACATGTTTTTTAGGCATCTTGAAc                     >  2:228964/1‑149 (MQ=255)                                                                                                                                ttGCCGTCGCCCAGTGGCAACTGGGCCTGGGCAGCCGGC‑GCTTCGGCA‑GGGGCCTTGTCGCCTTTCTTGGGCATCAATACCGCGGCCGCAACGGCAACCACGACCACAGCGGACAGCGCGAGCACATGtttt                                    <  8:269436/132‑1 (MQ=255)                                                                                                                                ttGCCGTCGCCCAGTGGCAACTGGGCCTGGGCAGCCGGC‑GCTTCGGCA‑GGGGCCTTGTCGCCTTTCTTGGGCATCAATACCGCGGCCGCAACGGCAACCACGACCACAGCGGACAGCGCGAGCACATGtttt                                    >  7:269436/1‑132 (MQ=255)                                                                                                                                   ccGTCGCCCAGTGGCAACTGGGCCTGGGCAGCCGGC‑GCTTCGGCA‑GGGGCCTTGTCGCCTTTCTTGGGCATCAATACCGCGGCCGCAACGGCAACCACGACCACAGCGGACAGCGCGAGCACATGTTTTTTAGGCATCTTGAACCCCAt                <  3:23433/149‑1 (MQ=255)                                                                                                                                    ccaggCCCAGTTGCCACTGGGCCTGGGCAGCCGGC‑GCTTCGGCA‑GGGGCCTTGTCGCCTTTCTTGGGCATCAATACCGCGGCCGCAACGGCAACCACGACCACAGCGGACAGCGCGAGCACATGTTTTTTAGGCATCTTGAACCCCAtg               >  7:249730/5‑149 (MQ=255)                                                                                                                                    ccaggCCCAGTTGCCACTGGGCCTGGGCAGCCGGC‑GCTTCGGCA‑GGGGCCTTGTCGCCTTTCTTGGGCATCAATACCGCGGCCGCAACGGCAACCACGACCACAGCGGACAGCGCGAGCACATGTTTTTTAGGCATCTTGAACCCCAtg               >  4:186641/5‑149 (MQ=255)                                                                                                                                    cGTCGCCCAGTGGCAACTGGGCCTGGGCAGCCGGC‑GCTTCGGCA‑GGGGCCTTGTCGCCTTTCTTGGGCATCAATACCGCGGCCGCAACGGCAACCACGACCACAGCGGACAGCGCGAGCACATGTTTTTTAGGCATCTTGAACCCCAtg               >  7:218889/1‑149 (MQ=255)                                                                                                                                     gTCGCCCAGTGGCAACTGGGCCTGGGCAGCCGGC‑GCTTCGGCA‑GGGGCCTTGTCGCCTTTCTTGGGCATCAATACCGCGGCCGCAACGGCAACCACGACCACAGCGGACAGCGCGAGCACATGTTTTTTAGGCATCTTGAACCCCAtgg              >  3:295361/1‑149 (MQ=255)                                                                                                                                      tCGCCCAGTGGCAACTGGGCCTGGGCAGCCGGC‑GCTTCGGCA‑GGGGCCTTGTCGCCTTTCTTGGGCATCAATACCGCGGCCGCAACGGCAACCACGACCACAGCGGACAGCGCGAGCACATGTTTTTTAGGCATCTTGAACCCCAtggt             <  1:320652/149‑1 (MQ=255)                                                                                                                                        gCCCAGTGGCAACTGGGCCTGGGCAGCCGGC‑GCTTCGGCA‑GGGGCCTTGTCGCCTTTCTTGGGCATCAATACCGCGGCCGCAACGGCAACCACGACCACAGCGGACAGCGCGAGCACATGTTTTTTAGGCATCTTGAACCCCAtggttg           <  1:171938/149‑1 (MQ=255)                                                                                                                                         cccAGTGGCAACTGGGCCTGGGCAGCCGGC‑GCTTCGGCA‑GGGGCCTTGTCGCCTTTCTTGGGCATCAATACCGCGGCCGCAACGGCAACCACGACCACAGCGGACAGCGCGAGCACATGTTTTTTAGGCATCTTGAACCCCAtggttgg          <  6:96083/149‑1 (MQ=255)                                                                                                                                            aGTGGCAACTGGGCCTGGGCAGCCGGCGGCTTCGGCA‑GGGGCCTTGTCGCCTTTCTTGGGCATCAATACCGCGGCCGCAACGGCAACCACGACCACAGCGGACAGCGCGAGCACATGTTTTTTAGGCATCTTGAACCCCATGGTTGGTc        >  1:225426/1‑149 (MQ=255)                                                                                                                                             gTGGCAACTGGGCCTGGGCAGCCGGC‑GCTTCGGCA‑GGGGCCTTGTCGCCTTTCTTGGGCATCAATACCGCGGCCGCAACGGCAACCACGACCACAGCGGACAGCGCGAGCACATGTTTTTTAGGCATCTTGAACCCCATGGTTGGTCGc      >  6:213065/1‑149 (MQ=255)                                                                                                                                              tGGCAACTGGGCCTGGGCAGCCGGC‑GCTTCGGCA‑GGGGCCTTGTCGCCTTTCTTGGGCATCAATACCGCGGCCGCAACGGCAACCACGACCACAGCGGACAGCGCGAGCACATGTTTTTTAGGCATCTTGAACCCCATGGTTGGTCGCt     >  6:140883/1‑149 (MQ=255)                                                                                                                                                gCAACTGGGCCTGGGCAGCCGGC‑GCTTCGGCA‑GGGGCCTTGTCGCCTTTCTTGGGCATCAATACCGCGGCCGCAACGGCAACCACGACCACAGCGGACAGCGCGAGCACATGTTTTTTAGGCATCTTGAACCCCATGGTTGGTCGCtt    <  1:94995/148‑1 (MQ=255)                                                                                                                                                 cAACTGGGCCTGGGCAGCCGGC‑GCTTCGGCA‑GGGGCCTTGTCGCCTTTCTTGGGCATCAATACCGCGGCCGCAACGGCAACCACGACCACAGCGGACAGCGCGAGCACATGTTTTTTAGGCATCTTGAACCCCATGGTTGGTCGCTTgg  >  2:182827/1‑149 (MQ=255)                                                                                                                                                 cAACTGGGCCTGGGCAGCCGGC‑GCTTCGGCA‑GGGGCCTTGTCGCCTTTCTTGGGCATCAATACCGCGGCCGCAACGGCAACCACGACCACAGCGGACAGCGCGAGCACATGTTTTTTAGGCATCTTGAACCCCATGGTTGGTCGCTTgg  <  3:10599/149‑1 (MQ=255)                                                                                                                                                  |                                                                                                                                                      TCCCTCTTCACCCTCGCCCATGCCTGCGGCTTGGGCCAGCGGCTCGCGCATCACCGGCTGCGACTGGCCAACCAGTGGCAGCGGCATTGGCTGGGACGAGCCAGAGAACTCGATGGCAGCACCCTGCTTGCCGTCGCCCAGTGGCAACTGGGCCTGGGCAGCCGGC‑GCTTCGGCA‑GGGGCCTTGTCGCCTTTCTTGGGCATCAATACCGCGGCCGCAACGGCAACCACGACCACAGCGGACAGCGCGAGCACATGTTTTTTAGGCATCTTGAACCCCATGGTTGGTCGCTTGG  >  NC_002947/5798822‑5799114

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