Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,798,964	G→T	5.9%	P300Q (CCA→CAA)	PP_5077 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,798,964	0	G	T	5.9%	275.9 / 5.5	85	P300Q (CCA→CAA)	PP_5077	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (43/37);  new base T (2/3);  total (45/40)
Fisher's exact test for biased strand distribution p-value = 6.63e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.21e-01

TCCCTCTTCACCCTCGCCCATGCCTGCGGCTTGGGCCAGCGGCTCGCGCATCACCGGCTGCGACTGGCCAACCAGTGGCAGCGGCATTGGCTGGGACGAGCCAGAGAACTCGATGGCAGCACCCTGCTTGCCGTCGCCCAGTGGCAACTGGGCCTGGGCAGCCGGC‑GCTTCGGCA‑GGGGCCTTGTCGCCTTTCTTGGGCATCAATACCGCGGCCGCAACGGCAACCACGACCACAGCGGACAGCGCGAGCACATGTTTTTTAGGCATCTTGAACCCCATGGTTGGTCGCT  >  NC_002947/5798822‑5799111                                                                                                                                               |                                                                                                                                                      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                                                                                                                  <  6:265481/135‑1 (MQ=255)                                            tCGCGCATCACCGGCTGCGACTGGCCAACCAGTGGCAGCGGCATTGGCTGGGACGAGCCAGAGAACTCGATGGCAGCACCCTGCTTGCCGTCGCCCAGTGGCAACTGGGCCTGGGCAGCCGGC‑GCTTCGGCA‑ggg                                                                                                                  >  5:265481/1‑135 (MQ=255)                                              gcgcATCACCGGCTGCGACTGGCCAACCAGTGGCAGCGGCATTGGCTGGGACGAGCCAGAGAACTCGATGGCAGCACCCTGCTTGCCGTCGCCCAGTGGCAACTGGGCCTGGGCAGCCGGC‑GCTTCGGCA‑GGGGCCTTGTCGCCTTTCt                                                                                                  >  2:320652/1‑149 (MQ=255)                                              gcgcATCACCGGCTGCGACTGGCAAACCAGTGGCAGCGGCATTGGCTGGGACGAGCCAGAGAACTCAATGGCAGCACCCTGCTTGCCGTCGCCCAGTGGCAACTGGGCCTGGGCAGCCGGC‑GCTTCGGCA‑GGGGCCTTGTCGCCTTTCt                                                                                                  >  4:236568/1‑149 (MQ=255)                                               cgcATCACCGGCTGCGACTGGCCAACCAGTGGCAGCGGCATTGGCTGGGACGAGCCAGAGAACTCGATGGCAGCACCCTGCTTGCCGTCGCCCAGTGGCAACTGGGCCTGGGCAGCCGGC‑GCTTCGGCA‑GGGGCCTTGTCGCCTTTCtt                                                                                                 <  7:363918/149‑1 (MQ=255)                                               cgcATCACCGGCTGCGACTGGCCAACCAGTGGCAGCGGCATTGGCTGGGACGAGCCAGAGAACTCGATGGCAGCACCCTGCTTGCCGTCGCCCAGTGGCAACTGGGCCTGGGCAGCCGGC‑GCTTCGGCA‑GGGGCCTTGTCGCCTTTCtt                                                                                                 <  2:110268/149‑1 (MQ=255)                                                gcATCACCGGCTGCGACTGGCCAACCAGTGGCAGCGGCATTGGCTGGGACGAGCCAGAGAACTCGATGGCAGCACCCTGCTTGCCGTCGCCCAGTGGCAACTGGGCCTGGGCAGCCGGC‑GCTTCGGCA‑GGGGCCTTGTCGCCTTTCTTg                                                                                                >  4:80762/1‑149 (MQ=255)                                                gcATCACCGGCTGCGACTGGCCAACCAGTGGCAGCGGCATTGGCTGGGACGAGCCAGAGAACTCGATGGCAGCACCCTGCTTGCCGTCGCCCAGTGGCAACTGGGCCTGGGCAGCCGGC‑GCTTCGGCA‑GGGGCCTTGTCGCCTTTCTTg                                                                                                <  5:264589/149‑1 (MQ=255)                                                 cATCACCGGCTGCGACTGGCCAACCAGTGGCAGCGGCATTGGCTGGGACGAGCCAGAGAACTCGATGGCAGCACCCTGCTTGCCGTCGCCCAGTGGCAACTGGGCCTGGGCAGCCGGC‑GCTTCGGCA‑GGGGCCTTGTCGCCTTTCTTgg                                                                                               >  8:336512/1‑149 (MQ=255)                                                  aTCACCGGCTGCGACTGGCCAACCAGTGGCAGCGGCATTGGCTGGGACGAGCCAGAGAACTCGATGGCAGCACCCTGCTTGCCGTCGCCCAGTGGCAACTGGGCCTGGGCAGCCGGC‑GCTTCGGCA‑GGGGCCTTGTCGCCTTTCTTggg                                                                                              <  6:139218/149‑1 (MQ=255)                                                  aTCACCGGCTGCGACTGGCCAACCAGTGGCAGCGGCATTGGCTGGGACGAGCCAGAGAACTCGATGGCAGCACCCTGCTTGCCGTCGCCCAGTGGCAACTGGGCCTGGGCAGCCGGC‑GCTTCGGCA‑GGGGCCTTGTCGCCTTTCTTggg                                                                                              >  5:139218/1‑149 (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                                >  4:171787/1‑149 (MQ=255)                                                                                                                gATGGCAGCACCCTGCTTGCCGTCGCCCAGTGGCAACTGGGCCTGGGCAGCCGGC‑GCTTCGGCA‑GGGGCCTTGTCGCCTTTCTTGGGCATCAATACCGCGGCCGCAACGGCAACCACGACCACAGCGGACAGCGCGAGCACATGttttt                                >  2:171938/1‑149 (MQ=255)                                                                                                                gATGGCAGCACCCTGCTTGCCGTCGCCCAGTGGCAACTGGGCCTGGGCAGCCGGC‑GCTTCGGCA‑GGGGCCTTGTCGCCTTTCTTGGGCATCAATACCGCGGCCGCAACGGCAACCACGACCACAGCGGACAGCGCGAGCACATGtttt                                 <  2:216824/148‑1 (MQ=255)                                                                                                                gATGGCAGCACCCTGCTTGCCGTCGCCCAGTGGCAACTGGGCCTGGGCAGCCGGC‑GCTTCGGCA‑GGGGCCTTGTCGCCTTTCTTGGGCATCAATACCGCGGCCGCAACGGCAACCACGACCACAGCGGACAGCGCGAGCACATGtttt                                 <  6:37509/148‑1 (MQ=255)                                                                                                                 aTGGCAGCACCCTGCTTGCCGTCGCCCAGTGGCAACTGGGCCTGGGCAGCCGGC‑GCtt                                                                                                                           >  6:56236/1‑58 (MQ=255)                                                                                                                 aTGGCAGCACCCTGCTTGCCGTCGCCCAGTGGCAACTGGGCCTGGGCAGCCGGC‑GCtt                                                                                                                           <  5:56236/58‑1 (MQ=255)                                                                                                                 aTGGCAGCACCCTGCTTGCCGTCGCCCAGTGGCAACTGGGCCTGGGCAGCCGGC‑GCTTCGGCA‑GGGGCCTTGTCGCCTTTCTTGGGCATCAATACCGCGGCCGCAACGGCAACCACGACCACAGCGGACAGCGCGAGCACATGtttttt                               <  2:368686/149‑1 (MQ=255)                                                                                                                    gcagcaCCCTGCTTGCCGTCGCCCAGTGGCAACTGGGCCTGGGCAGCCGGC‑GCTTCGGCA‑GGGGCCTTGTCGCCTTTCTTGGGCATCAATACCGCGGCCGCAACGGCAACCACGACCACAGCGGACAGCGCGAGCACATGTTTTTTAgg                            >  2:329655/1‑149 (MQ=255)                                                                                                                    gcagcaCCCTGCTTGCCGTCGCCCAGTGGCAACTGGGCCTGGGCAGCCGGC‑GCTTCGGCA‑GGGGCCTTGTCGCCTTTCTTGGGCATCAATACCGCGGCCGCAACGGCAACCACGACCACAGCGGACAGCGCGAGCACATGTTTTTTAg                             <  3:80762/148‑1 (MQ=255)                                                                                                                    gcagcaCCCTGCTTGCCGTCGCCCAGTGGCAACTGGGCCTGGGCAGCCGGC‑GCTTCGGCA‑GGGGCCTTGTCGCCTTTCTTGGGCATCAATACCGCGGCCGCAACGGCAACCACGACCACAGCGGACAGCGCGAGCACATGTTTTTTAg                             <  7:190068/148‑1 (MQ=255)                                                                                                                    gcagcaCCCTGCTTGCCGTCGCCCAGTGGCAACTGGGCCTGGGCAGCCGGC‑GCTTCGGCA‑GGGGCCTTGTCGCCTTTCTTGGGCATCAATACCGCGGCCGCAACGGCAACCACGACCACAGCGGACAGCGCGAGCACATGTTTTTTAg                             <  3:236568/148‑1 (MQ=255)                                                                                                                          cccTGCTTGCCGTCGCCCAGTGGCAACTGGGCCTGGGCAGCCGGC‑GCTTCGGCA‑GGGGCCTTGTCGCCTTTCTTGGGCATCAATACCGCGGCCGCAACGGCAACCACGACCACAGCGGACAGCGCGAGCACATGTTTTTTAGGCATCt                       >  4:84457/1‑148 (MQ=255)                                                                                                                          cccTGCTTGCCGTCGCCCAGTGGCAACTGGGCCTGGGCAGCCGGC‑GCTTCGGCA‑GGGGCCTTGTCGCCTTTCTTGGGCATCAATACCGCGGCCGCAACGGCAACCACGACCACAGCGGACAGCGCGAGCACATGTTTTTTAGGCATCt                       <  3:84457/148‑1 (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            >  4:186641/5‑149 (MQ=255)                                                                                                                                    ccaggCCCAGTTGCCACTGGGCCTGGGCAGCCGGC‑GCTTCGGCA‑GGGGCCTTGTCGCCTTTCTTGGGCATCAATACCGCGGCCGCAACGGCAACCACGACCACAGCGGACAGCGCGAGCACATGTTTTTTAGGCATCTTGAACCCCAtg            >  7:249730/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)                                                                                                                                               |                                                                                                                                                      TCCCTCTTCACCCTCGCCCATGCCTGCGGCTTGGGCCAGCGGCTCGCGCATCACCGGCTGCGACTGGCCAACCAGTGGCAGCGGCATTGGCTGGGACGAGCCAGAGAACTCGATGGCAGCACCCTGCTTGCCGTCGCCCAGTGGCAACTGGGCCTGGGCAGCCGGC‑GCTTCGGCA‑GGGGCCTTGTCGCCTTTCTTGGGCATCAATACCGCGGCCGCAACGGCAACCACGACCACAGCGGACAGCGCGAGCACATGTTTTTTAGGCATCTTGAACCCCATGGTTGGTCGCT  >  NC_002947/5798822‑5799111

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