Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,550,238	T→C	5.3%	S29S (AGT→AGC)	PP_3134 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,550,238	0	T	C	5.3%	225.8 / 2.1	75	S29S (AGT→AGC)	PP_3134	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (34/37);  new base C (2/2);  total (36/39)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.05e-01

CACCTCGATAACAAGATGAAAGCGGCCAGCGGCGCCGCGGCGTGATGAAGGTTGTACACGGGTGAGTGCGATGCGGTCAGCACCTAGAAGAACGAGACCGGAACGCCTCAGGAGACTCCAGATTGTTCAGATCGCGAGGCCAGGCAGTATGAAAGGATGGATACGCAAGGCAGTGGCGCGCTACGCCCATGCAACCGGGCGAGGTGGCAAATTCTACAGGAGGTTCTGCAACCCCTCTGGTCTTGAATGGGGGGCCTATCTGGCCCGATGGGGCAACTTTCATTCGGTCGGCAG  >  NC_002947/3550091‑3550384                                                                                                                                                    |                                                                                                                                                   caccTCGATAACAAGATGAAAGCGGCCAGCGGCGCCGCGGCGTGATGAAGGTTGTACACGGGTGAGTGCGATGCGGTCAGCACCTAGAAGAACGAGACCGGAACGCCTCAGGAGACTCCAGATTGTTCAGATCGCGAGGCCAGGCAGTa                                                                                                                                                   >  8:35802/1‑149 (MQ=255) caccTCGATAACAAGATGAAAGCGGCCAGCGGCGCCGCGGCGTGATGAAGGTTGTACACGGGTGAGTGCGATGCGGTCAGCACCTAGAAGAACGAGACCGGAACGCCTCAGGAGACTCCAGATTGTTCAGATCGCGAGGCCAGGCAGTa                                                                                                                                                   <  7:317901/149‑1 (MQ=255)       gATAACAAGATGAAAGCGGCCAGCGGCGCCGCGGCGTGATGAAGGTTGTACACGGGTGAGTGCGATGCGGTCAGCACCTAGAAGAACGAGACCGGAACGCCTCAGGAGACTCCAGATTGTTCAGATCGCGAGGCCAGGCAGTATGaaa                                                                                                                                              <  1:231110/148‑1 (MQ=255)       gATAACAAGATGAAAGCGGCCAGCGGCGCCGCGGCGTGATGAAGGTTGTACACGGGTGAGTGCGATGCGGTCAGCACCTAGAAGAACGAGACCGGAACGCCTCAGGAGACTCCAGATTGTTCAGATCGCGAGGCCAGGCAGTATGaa                                                                                                                                               <  6:225134/147‑1 (MQ=255)       gATAACAAGATGAAAGCGGCCAGCGGCGCCGCGGCGTGATGAAGGTTGTACACGGGTGAGTGCGATGCGGTCAGCACCTAGAAGAACGAGACCGGAACGCCTCAGGAGACTCCAGATTGTTCAGATCGCGAGGCCAGGCAGTATGAAAg                                                                                                                                             >  1:164306/1‑149 (MQ=255)       gATAACAAGATGAAAGCGGCCAGCGGCGCCGCGGCGTGATGAAGGTTGTACACGGGTGAGTGCGATGCGGTCAGCACCTAGAAGAACGAGACCGGAACGCCTCAGGAGACTCCAGATTGTTCAGATCGCGAGGCCAGGCAGTATGAAAg                                                                                                                                             >  6:206252/1‑149 (MQ=255)         tAACAAGATGAAAGCGGCCAGCGGCGCCGCGGCGTGATGAAGGTTGTACACGGGTGAGTGCGATGCGGTCAGCACCTAGAAGAACGAGACCGGAACGCCTCAGGAGACTCCAGATTGTTCAGATCGCGAGGCCAGGCAGTATGAAAgga                                                                                                                                           <  2:271572/149‑1 (MQ=255)             aaGATGAAAGCGGCCAGCGGCGCCGCGGCGTGATGAAGGTTGTACACGGGTGAGTGCGATGCGGTCAGCACCTAGAAGAACGAGACAGGAACGCCTCAGGAGACTCCAGAATGTTCAGATCGCGAGGCCAGGCAGTATGAAAggatgga                                                                                                                                       >  5:108711/1‑149 (MQ=255)               gATGAAAGCGGCCAGCGGCGCCGCGGCGTGATGAAGGTTGTACACGGGTGAGTGCGATGCGGTCAGCACCTAGAAGAACGAGACCGGAACGCCTCAGGAGACTCCAGATTGTTCAGATCGCGAGGCCAGGCAGTATGAAAGGATGGATa                                                                                                                                     >  8:233385/1‑149 (MQ=255)               gATGAAAGCGGCCAGCGGCGCCGCGGCGTGATGAAGGTTGTACACGGGTGAGTGCGATGCGGTCAGCACCTAGAAGAACGAGACCGGAACGCCTCAGGAGACTCCAGATTGTTCAGATCGCGAGGCCAGGCAGTATGAAAGGATGGATa                                                                                                                                     >  7:164562/1‑149 (MQ=255)                    aaGCGGCCAGCGGCGCCGCGGCGTGATGAAGGTTGTACACGGGTGAGTGCGATGCGGTCAGCACCTAGAAGAACGAGACCGGAACGCCTCAGGAGACTCCAGATTGTTCAGATCGCGAGGCCAGGCAGTATGAAAggatggat                                                                                                                                      >  8:17159/1‑143 (MQ=255)                    aaGCGGCCAGCGGCGCCGCGGCGTGATGAAGGTTGTACACGGGTGAGTGCGATGCGGTCAGCACCTAGAAGAACGAGACCGGAACGCCTCAGGAGACTCCAGATTGTTCAGATCGCGAGGCCAGGCAGTATGAAAggatggat                                                                                                                                      <  7:17159/143‑1 (MQ=255)                     aGCGTCCAGCGGCGACGCGGCGTTATGAAGGTTGTACACGGGTGAGTGCGATGCGGTCAGCACCTAGAAGAACGAGACCGGAACGCCTCAGGAGACTCCAGATTTTTCAGATCGCGAGGCCAGGCAGTATGAAAGGATGGATACGCAAg                                                                                                                               <  3:216034/149‑1 (MQ=255)                      gCGGCCAGCGGCGCCGCGGCGTGATGAAGGTTGTACACGGGTGAGTGCGATGCGGTCAGCACCTAGAAGAACGAGACCGGAACGCCTCAGGAGACTCCAGATTGTTCAGATCGCGAGGCCAGGCAGTATGAAAGGATGGATACGCAAgg                                                                                                                              <  3:190111/149‑1 (MQ=255)                        ggCCAGCGGCGCCGCGGCGTGATGAAGGTTGTACACGGGTGAGTGCGATGCGGTCAGCACCTAGAAGAACGAGACCGGAACGCCTCAGGAGACTCCAGATTGTTCAGATCGCGAGGCCAGGCAGTATGAAAGGATGGATACGCAAgg                                                                                                                              >  5:156026/1‑147 (MQ=255)                          ccAGCGGCGCCGCGGCGTGATGAAGGTTGTACACGGGTGAGTGCGATGCGGTCAGCACCTAGAAGAACGAGACCGGAACGCCTCAGGAGACTCCAGATTGTTCAGATCGCGAGGCCAGGCAGTATGAAAGGATGGATACGCAAGGCAGt                                                                                                                          >  1:386669/1‑149 (MQ=255)                           cAGCGGCGCCGCGGCGTGATGAAGGTTGTACACGGGTGAGTGCGATGCGGTCAGCACCTAGAAGAACGAGACCGGAACGCCTCAGGAGACTCCAGATTGTTCAGATCGCGAGGCCAGGCAGTATGAAAGGATGGATACGCAAGGCAGTg                                                                                                                         <  3:93642/149‑1 (MQ=255)                             gcggcgCCGCGGCGTGATGAAGGTTGTACACGGGTGAGTGCGATGCGGTCAGCACCTAGAAGAACGAGACCGGAACGCCTCAGGAGACTCCAGATTGTTCAGATCGCGAGGCCAGGCAGTATGAAAGGATGGATACGCAAGGCAGTGgc                                                                                                                       >  2:260719/1‑149 (MQ=255)                               ggcgCCGCGGCGTGATGAAGGTTGTACACGGGTGAGTGCGATGCGGTCAGCACCTAGAAGAACGAGACCGGAACGCCTCAGGAGACTCCAGATTGTTCAGATCGCGAGGCCAGGCAGTATGAAAGGATGGATACGCAAGGCAGTGgcgc                                                                                                                     <  5:26172/149‑1 (MQ=255)                                    cgcgGCGTGATGAAGGTTGTACACGGGTGAGTGCGATGCGGTCAGCACCTAGAAGAACGAGACCGGAACGCCTCAGGAGACTCCAGATTGTTCAGATCGCGAGGCCAGGCAGTATGAAAGGATGGATACGCAAGGCAGTGGCGCGCTAc                                                                                                                >  2:298477/1‑149 (MQ=255)                                     gcggcgTGATGAAGGTTGTACACGGGTGAGTGCGATGCGGTCAGCACCTAGAAGAACGAGACCGGAACGCCTCAGGAGACTCCAGATTGTTCAGATCGCGAGGCCAGGCAGTATGAAAGGATGGATACGCAAGGCAGTGGCGCGCTACg                                                                                                               >  8:312728/1‑149 (MQ=255)                                         cgTGATGAAGGTTGTACACGGGTGAGTGCGATGCGGTCAGCACCTAGAAGAACGAGACCGGAACGCCTCAGGAGACTCCAGATTGTTCAGATCGCGAGGCCAGGCAGTATGAAAGGATGGATACGCAAGGCAGTGGCGCGCTACGCCCa                                                                                                           <  7:236238/149‑1 (MQ=255)                                         cgTGATGAAGGTTGTACACGGGTGAGTGCGATGCGGTCAGCACCTAGAAGAACGAGACCGGAACGCCTCAGGAGACTCCAGATTGTTCAGATCGCGAGGCCAGGCAGTATGAAAGGATGGATACGCAAGGCAGTGGCGCGCTACGCCCa                                                                                                           <  4:31186/149‑1 (MQ=255)                                           tgatgaAGGTTGTACACGGGTGAGTGCGATGCGGTCAGCACCTAGAAGAACGAGACCGGAACGCCTCAGGAGACTCCAGATTGTTCAGATCGCGAGGCCAGGCAGTATGAAAGGATGGATACGCAAGGCAGTGGCGCGCTACGCCCATg                                                                                                         <  7:312728/149‑1 (MQ=255)                                            gatgaAGGTTGTACACGGGTGAGTGCGATGCGGTCAGCACCTAGAAGAACGAGACCGGAACGCCTCAGGAGACTCCAGATTGTTCAGATCGCGAGGCCAGGCAGTATGAAAGGATGGATACGCAAGGCAGTGGCGCGCTACGCCCATg                                                                                                         >  4:116600/1‑148 (MQ=255)                                             atgaAGGTTGTACACGGGTGAGTGCGATGCGGTCAGCACCTAGAAGAACGAGACCGGAACGCCTCAGGAGACTCCAGATTGTTCAGATCGCGAGGCCAGGCAGTATGAAAGGATGGATACGCAAGGCAGTGGCGCGCTACGCCCATGc                                                                                                        >  3:374382/1‑148 (MQ=255)                                                aaGGTTGTACACGGGTGAGTGCGATGCGGTCAGCACCTAGAAGAACGAGACCGGAACGCCTCAGGAGACTCCAGATTGTTCAGATCGCGAGGCCAGGCAGTATGAAAGGATGGATACGCAAGGCAGTGGCGCGCTACGCCCATGCAAcc                                                                                                    <  2:94647/149‑1 (MQ=255)                                                aaGGTTGTACACGGGTGAGTGCGATGCGGTCAGCACCTAGAAGAACGAGACCGGAACGCCTCAGGAGACTCCAGATTGTTCAGATCGCGAGGCCAGGCAGTATGAAAGGATGGATACGCAAGGCAGTGGCGCGCTACGCCCATGCAAc                                                                                                     <  1:158297/148‑1 (MQ=255)                                                    ttGTACACGGGTGAGTGCGATGCGGTCAGCACCTAGAAGAACGAGACCGGAACGCCTCAGGAGACTCCAGATTGTTCAGATCGCGAGGCCAGGCAGTATGAAAGGATGGATACGCAAGGCAGTGGCGCGCTACGCCCATg                                                                                                         >  3:225119/1‑140 (MQ=255)                                                    ttGTACACGGGTGAGTGCGATGCGGTCAGCACCTAGAAGAACGAGACCGGAACGCCTCAGGAGACTCCAGATTGTTCAGATCGCGAGGCCAGGCAGTATGAAAGGATGGATACGCAAGGCAGTGGCGCGCTACGCCCATg                                                                                                         <  4:225119/140‑1 (MQ=255)                                                         cacGGGTGAGTGCGATGCGGTCAGCACCTAGAAGAACGAGACCGGAACGCCTCAGGAGACTCCAGATTGTTCAGATCGCGAGGCCAGGCAGCATGAAAGGATGGATACGCAAGGCAGTGGCGCGCTACGCCCATGCAACCGGGCGAGGt                                                                                           >  4:357506/1‑149 (MQ=255)                                                         cacGGGTGAGTGCGATGCGGTCAGCACCTAGAAGAACGAGACCGGAACGCCTCAGGAGACTCCAGATTGTTCAGATCGCGAGGCCAGGCAGCATGAAAGGATGGATACGCAAGGCAGTGGCGCGCTACGCCCATGCAACCGGGCGAGGt                                                                                           >  4:363330/1‑149 (MQ=255)                                                          acGGGTGAGTGCGATGCGGTCAGCACCTAGAAGAACGAGACCGGAACGCCTCAGGAGACTCCAGATTGTTCAGATCGCGAGGCCAGGCAGTATGAAAGGATGGATACGCAAGGCAGTGGCGCGCTACGCCCATGCAACCGGGCGAGGTg                                                                                          <  1:375860/149‑1 (MQ=255)                                                           cGGGTGAGTGCGATGCGGTCAGCACCTAGAAGAACGAGACCGGAACGCCTCAGGAGACTCCAGATTGTTCAGATCGCGAGGCCAGGCAGTATGAAAGGATGGATACGCAAGGCAGTGGCGCGCTACGCCCATGCAACCGGGCGAGGTgg                                                                                         <  5:206252/149‑1 (MQ=255)                                                              gTGAGTGCGATGCGGTCAGCACCTAGAAGAACGAGACCGGAACGCCTCAGGAGACTCCAGATTGTTCAGATCGCGAGGCCAGGCAGTATGAAAGGATGGATACGCAAGGCAGTGGCGCGCTACGCCCATGCAACCGGGCGAGGTGGCa                                                                                       >  4:148899/1‑148 (MQ=255)                                                                 aGTGCGATGCGGTCAGCACCTAGAAGAACGAGACCGGAACGCCTCAGGAGAATCCAGATTGTTCAGATCGCGAGGACAGGCAGTATGAAAGCATGGATACGCAATGCAGTGGCGCCCTACTCCCATGCAACCGGGCGAGGTGGCAAAt                                                                                    >  3:259963/1‑148 (MQ=255)                                                                     cGATGCGGTCAGCACCTAGAAGAACGAGACCGGAACGCCTCAGGAGACTCCAGATTGTTCAGATCGCGAGGCCAGGCAGTATGAAAGGATGGATACGCAAGGCAGTGGCGCGCTACGCCCATGCAACCGGGCGAGGTGGCAAATTCTAc                                                                               <  1:260719/149‑1 (MQ=255)                                                                       aTGCGGTCAGCACCTAGAAGAACGAGACCGGAACGCCTCAGGAGACTCCAGATTGTTCAGATCGCGAGGCCAGGCAGTATGAAAGGATGGATACGCAAGGCAGTGGCGCGCTACGCCCATGCAACCGGGCGAGGTGGCAAATTCTACag                                                                             <  3:116600/149‑1 (MQ=255)                                                                         gCGGTCAGCACCTAGAAGAACGAGACCGGAACGCCTCAGGAGACTCCAGATTGTTCAGATCGCGAGGCCAGGCAGTATGAAAGGATGGATACGCAAGGCAGTGGCGCGCTACGCCCATGCAACCGGGCGAGGTGGCAAATTCTACag                                                                             >  8:214532/1‑147 (MQ=255)                                                                          cGGTCAGCACCTAGAAGAACGAGACCGGAACGCCTCAGGAGACTCCAGATTGTTCAGATCGCGAGGCCAGGCAGTATGAAAGGATGGATACGCAAGGCAGTGGCGCGCTACGCCCATGCAACCGGGCGAGGTGGCAAATTCAACagga                                                                           <  8:164562/148‑1 (MQ=255)                                                                            gTCAGCACCTAGAAGAACGAGACCGGAACGCCTCAGGAGACTCCAGATTGTTCAGATCGCGAGGCCAGGCAGTATGAAAGGATGGATACGCAAGGCAGTGGCGCGCTACGCCCATGCAACCGGGCGAGGTGGCAAATTCTACAGGAGGt                                                                        >  8:307334/1‑149 (MQ=255)                                                                               aGCACCTAGAAGAACGAGACCGGAACGCCTCAGGAGACTCCAGATTGTTCAGATCGCGAGGCCAGGCAGTATGAAAGGATGGATACGCAAGGCAGTGGCGCGCTACGCCCATGCAACCGGGCGAGGTGGCAAATTCTACAGGAGGTTc                                                                      <  7:214532/148‑1 (MQ=255)                                                                                    cTAGAAGAACGAGACCGGAACGCCTCAGGAGACTCCAGATTGTTCAGATCGCGAGGCCAGGCAGTATGAAAGGATGGATACGCAAGGCAGTGGCGCGCTACGCCCATGCAACCGGGCGAGGTGGCAAATTCTACAGGAGGTTCTGCAAc                                                                >  5:288760/1‑149 (MQ=255)                                                                                         agaaCGAGACCGGAACGCCTCAGGAGACTCCAGATTGTTCAGATCGCGAGGCCAGGCAGTATGAAAGGATGGATACGCAAGGCAGTGGCGCGCTACGCCCATGCAACCGGGCGAGGTGGCAAATTCTACAGGAGGTTCTGCAACCCctc                                                           >  8:255851/1‑149 (MQ=255)                                                                                          gaaCGAGACCGGAACGCCTCAGGAGACTCCAGATTGTTCAGATCGCGAGGCCAGGCAGTATGAAAGGATGGATACGCAAGGCAGTGGCGCGCTACGCCCATGCAACCGGGCGAGGTGGCAAATTCTACAGGAGGTTCTGCAACCCctct                                                          <  4:374382/149‑1 (MQ=255)                                                                                           aaCGAGACCGGAACGCCTCAGGAGACTCCAGATTGTTCAGATCGCGAGGCCAGGCAGTATGAAAGGATGGATACGCAAGGCAGTGGCGCGCTACGCCCATGCAACCGGGCGAGGTGGCAAATTCTACAGGAGGTTCTGCAACCCCTCTg                                                         <  1:298477/149‑1 (MQ=255)                                                                                              gagaCCGGAACGCCTCAGGAGACTCCAGATTGTTCAGATCGCGAGGCCAGGCAGTATGAAAGGATGGATACGCAAGGCAGTGGCGCGCTACGCCCATGCAACCGGGCGAGGTGGCAAATTCTACAGGAGGTTCTGGAACCCCTCTGGTc                                                      >  6:270412/1‑149 (MQ=255)                                                                                                 cccGGAAAGCCTCAGGAGACTCCAGATTGTTCAGATCGCGAGGCCAGGCAGTATGAAAGGATGGATACGCAAGGCAGTGGCGCGCTACGCCAATGCAACCGGGCGAGGTGGCAAATTCTACAGGAGGTTCTGCAACCCCTCTGGTCtt                                                    >  6:308640/2‑148 (MQ=255)                                                                                                   cGGAACGCCTCAGGAGACTCCAGATTGTTCAGATCGCGAGGCCAGGCAGTATGAAAGGATGGATACGCAAGGCAGTGGCGCGCTACGCCCATGCAACCGGGCGAGGTGGCAAATTCTACAGGAGGTTCTGCAACCCCTCTGGTCTTGaa                                                 >  3:376324/1‑149 (MQ=255)                                                                                                   cGGAACGCCTCAGGAGACTCCAGATTGTTCAGATCGCGAGGCCAGGCAGTATGAAAGGATGGATACGCAAGGCAGTGGCGCGCTACGCCCATGCAACCGGGCGAGGTGGCAAATTCTACAGGAGGTTCTGCAACCCCTCTGGTCTTGaa                                                 >  1:227124/1‑149 (MQ=255)                                                                                                     gAACGCCTCAGGAGACTCCAGATTGTTCAGATCGCGAGGCCAGGCAGTATGAAAGGATGGATACGCAAGGCAGTGGCGCGCTACGCCCATGCAACCGGGCGAGGTGGCAAATTCTACAGGAGGTTCTGCAACCCCTCTGGTCTTGAATg                                               <  6:288760/149‑1 (MQ=255)                                                                                                         gCCTCAGGAGACTCCAGATTGTTCAGATCGCGAGGCCAGGCAGTATGAAAGGATGGATACGCAAGGCAGTGGCGCGCTACGCCCATGCAACCGGGCGAGGTGGCAAATTCTACAGGAGGTTCTGCAACCCCTCTGGTCTTGAATggggg                                           <  7:35802/149‑1 (MQ=255)                                                                                                          ccTCAGGAGACTCCAGATTGTTCAGATCGCGAGGCCAGGCAGTAt                                                                                                                                                  >  6:79198/1‑45 (MQ=255)                                                                                                          ccTCAGGAGACTCCAGATTGTTCAGATCGCGAGGCCAGGCAGTATGAAAGGATGGATACGCAAGGCAGTGGCGCGCTACGCCCATGCAACCGGGCGAGGTGGCAAATTCTACAGGAGGTTCTGCAACCCCTCTGGTCTTGAATgggggg                                          >  5:91269/1‑149 (MQ=255)                                                                                                                   aCTCCAGATTGTTCAGATCGCGAGGCCAGGCAGTATGAAAGGATGGATACGCAAGGCAGTGGCGCGCTACGCCCATGCAACCGGGCGAGGTGGCAAATTCTACAGGAGGTTCTGCAACCCCTCTGGTCTTGAATGGGGGGCCTATCTgg                                 >  7:10105/1‑149 (MQ=255)                                                                                                                     tCCAGTTTGTTCAGATCGCGAGGGCATGCATTTTGTAAGGATTGATACGCAACGGAGTTGCGCGGTACGCCCAAGCAACCGGGCGTGGTGGCAAATTCAACAGGAGGTTCTGCAACCCCTCTGGTCTTGAATGGGGGGCCTATCTGGcc                               <  2:164306/149‑1 (MQ=255)                                                                                                                     tCCAGATTGTTCAGATCGCGAGGCCAGGCAGTATGAAAGGATGGATACGCAAGGCAGTGGCGCGCTACGCCCATGCAACCGGGCGAGGTGGCAAATTCTACAGGAGGTTCTGCAACCCCTCTGGTCTTGAATGGGGGGCCTATCTggtc                               >  2:271946/1‑147 (MQ=255)                                                                                                                     tCCAGATTGTTCAGATCGCGAGGCCAGGCAGTATGAAAGGATGGATACGCAAGGCAGTGGCGCGCTACGCCCATGCAACCGGGCGAGGTGGCAAATTCTACAGGAGGTTCTGCAACCCCTCTGGTCTTGAATGGGGGGCCTATCTGGcc                               >  1:25325/1‑149 (MQ=255)                                                                                                                     tCCAGATTGTTCAGATCGCGAGGCCAGGCAGTATGAAAGGATGGATACGCAAGGCAGTGGCGCGCTACGCCCATGCAACCGGGCGAGGTGGCAAATTCTACAGGAGGTTCTGCAACCCCTCTGGTCTTGAATGGGGGGCCTATCTGGcc                               >  6:357507/1‑149 (MQ=255)                                                                                                                         gATTGTTCAGATCGCGAGGCCAGGCAGTATGAAAGGATGGATACGCAAGGCAGTGGCGCGCTACGCCCATGCAACCGGGCGAGGTGGCAAATTCTACAGGAGGTTCTGCAACCCCTCTGGTCTTGAATGGGGGGCCTATCTGGCCCGAt                           <  7:233385/149‑1 (MQ=255)                                                                                                                         gATTGTTCAGATCGCGAGGCCAGGCAGTATGAAAGGATGGATACGCAAGGCAGTGGCGCGCTACGCCCATGCAACCGGGCGAGGTGGCAAATTCTACAGGAGGTTCTGCAACCCCTCTGGTCTTGAATGGGGGGCCTATCTGGCCCGAt                           >  2:117718/1‑149 (MQ=255)                                                                                                                          aTTGTTCAGATCGCGAGGCCAGGCAGTATGAAAGGATGGATACGCAAGGCAGTGGCGCGCTACGCCCATGCAACCGGGCGAGGTGGCAAATTCTACAGGAGGTTCTGCAACCCCTCTGGTCTTGAATGGGGGGCCTATCTGGCCCGATg                          >  3:155351/1‑149 (MQ=255)                                                                                                                           ttGTTCAGATCGCGAGGCCAGGCAGTATGAAAGGATGGATACGCAAGGCAGTGGCGCGCTACGCCCATGCAACCGGGCGAGGTGGCAAATTCTACAGGAGGTTCTGCAACCCCTCTGGTCTTGAATGGGGGGCCTATCTGGCCCGATgg                         <  2:386669/149‑1 (MQ=255)                                                                                                                           ttGTTCAGATCGCGAGGCCAGGCAGTATGAAAGGATGGATACGCAAGGCAGTGGCGCGCTACGCCCATGCAACCGGGCGAGGTGGCAAATTCTACAGGAGGTTCTGCAACCCCTCTGGTCTTGAATGGGGGGCCTATCTGGCCCGATg                          <  3:148899/148‑1 (MQ=255)                                                                                                                             gTTCAGATCGCGAGGCCAGGCAGTATGAAAGGATGGATACGCAAGGCAGTGGCGCGCTACGCCCATGCAACCGGGCGAGGTGGCAAATTCTACAGGAGGTTCTGCAACCCCTCTGGTCTTGAATGGGGGGCCTATCTGGCCCGATggg                        <  7:307334/148‑1 (MQ=255)                                                                                                                             gTTCAGATCGCGAGGCCAGGCAGTATGAAAGGATGGATACGCAAGGCAGTGGCGCGCTACGCCCATGCAACCGGGCGAGGTGGCAAATTCTACAGGAGGTTCTGCAACCCCTCTGGTCTTGAATGGGGGGCCTATCTGGCCCGATggg                        <  6:156026/148‑1 (MQ=255)                                                                                                                                         aGGCCAGGCAGTATGAAAGGATGGATACGCAAGGCAGTGGCGCGCTACGCCCATGCAACCGGGCGAGGTGGCAAATTCTACAGGAGGTTCTGCAACCCCTCTGGTCTTGAATGGGGGGCCTATCTGGCCCGATGGGGCAACTTTCATtc           >  4:80631/1‑149 (MQ=255)                                                                                                                                            ccAGGCAGTATGAAAGGATGGATACGCAAGGCAGTGGCGCGCTACGCCCATGCAACCGGGCGAGGTGGCAAATTCTACAGGAGGTTCTGCAACCCCTCTGGTCTTGAATGGGGGGCCTATCTGGCCCGATGGGGCAACTTTCATtcggt        <  5:308640/149‑1 (MQ=255)                                                                                                                                             cAGGCAGCATGAAAGGATGGATACGCAAGGCAGTGGCGCGCTACGCCCATGCAACCGGGCGAGGTGGCAAATTCTACAGGAGGTTCTGCAACCCCTCTGGTCTTGAATGGGGGGCCTATCTGGCCCGATGGGGCAACTTTCATtcggtc       <  3:363330/149‑1 (MQ=255)                                                                                                                                             cAGGCAGCATGAAAGGATGGATACGCAAGGCAGTGGCGCGCTACGCCCATGCAACCGGGCGAGGTGGCAAATTCTACAGGAGGTTCTGCAACCCCTCTGGTCTTGAATGGGGGGCCTATCTGGCCCGATGGGGCAACTTTCATtcggtc       <  3:357506/149‑1 (MQ=255)                                                                                                                                                gCAGTATGAAAGGATGGATACGCAAGGCAGTGGCGCGCTACGCCCATGCAACCGGGCGAGGTGGCAAATTCTACAGGAGGTTCTGCAACCCCTCTGGTCTTGAATGGGGGGCCTATCTGGCCCGATGGGGCAACTTTCATTCGGTCGgc    <  5:139714/149‑1 (MQ=255)                                                                                                                                                 cAGTATGAAAGGATGGATACGCAAGGCAGTGGCGCGCTACGCCCATGCAACCGGGCGAGGTGGCAAATTCTACAGGAGGTTCTGCAACCCCTCTGGTCTTGAATGGGGGGCCTATCTGGCCCGATGGGGCAACTTTCATTCGGTCGgca   <  2:227124/149‑1 (MQ=255)                                                                                                                                                 cAGTATGAAAGGATGGATACGCAAGGCAGTGGCGCGCTACGCCCATGCAACCGGGCGAGGTGGCAAATTCTACAGGAGGTTCTGCAACCCCTCTGGTCTTGAATGGGGGGCCTATCTGGCCCGATGGGGCAACTTTCATTCGGTCGgca   <  4:376324/149‑1 (MQ=255)                                                                                                                                                 cAGTATGAAAGGATGGATACGCAAGGCAGGGGCGCGCTACGCCCATGCAACCGGGCGAGGTGGCAAATTCTACAGGAGGTTCTGCAACCCCTCTGGTCTTGAATGGGGGGCCTATCTGGCCCGATGGGGCAACTTTCATTCGGTCGgca   <  8:286058/149‑1 (MQ=255)                                                                                                                                                  aGTATGAAAGGATGGATACGCAAGGCAGTGGCGCGCTACGCCCATGCAACCGGGCGAGGTGGCAAATTCTACAGGAGGTTCTGCAACCCCTCTGGTCTTGAATGGGGGGCCTATCTGGCCCGATGGGGCAACTTTCATTCGGTCGgcag  <  8:10105/149‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                                   CACCTCGATAACAAGATGAAAGCGGCCAGCGGCGCCGCGGCGTGATGAAGGTTGTACACGGGTGAGTGCGATGCGGTCAGCACCTAGAAGAACGAGACCGGAACGCCTCAGGAGACTCCAGATTGTTCAGATCGCGAGGCCAGGCAGTATGAAAGGATGGATACGCAAGGCAGTGGCGCGCTACGCCCATGCAACCGGGCGAGGTGGCAAATTCTACAGGAGGTTCTGCAACCCCTCTGGTCTTGAATGGGGGGCCTATCTGGCCCGATGGGGCAACTTTCATTCGGTCGGCAG  >  NC_002947/3550091‑3550384

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑