Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,483,011	G→T	6.0%	R971R (CGC→CGA)	PP_3091 ←	membrane protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,483,011	0	G	T	6.0%	219.7 / 2.6	69	R971R (CGC→CGA)	PP_3091	membrane protein
Reads supporting (aligned to +/- strand):  ref base G (33/31);  new base T (2/2);  total (35/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.36e-01

TGAACAGCCCCTGCAGCGAGCGCGACAGGCCATCCTGGCTGGTGTCGACGCTGCTTTCCACG‑TAGTTGCGCACGGTGGTGATCTCGCTCGGCTGGCCTTCCAGGGTCTCGCTGATCAACTGTTTGCTGCTCTTGCCCACATCCTGCGAGCGCTGGATGTTGTTCATGCGCACCTTGAGTTTGCGCAAGGCGGCCAGGTAGCGATCCATGATGGTGCTGTCGGCGCCTTCGGCGTTGTTTTCGGCGAACA‑CCCGCGCCACCGGCTCGAAGCGCTTGGCCAGGCTGCCGTCG  >  NC_002947/3482863‑3483152                                                                                                                                                      |                                                                                                                                               tGAACAGCCCCTGCAGCGAGCGCGACAGGTCATCCTGGCTGGTGTCGACGCTGCTTTCCACG‑TAGTTGCGCACGGTGGTGATCTCGCTCGGCTGGCCTTCCAGGGTCTCGCTGATCAACTGTTTGCTGCTCTTGCCCACATCCTGCGAg                                                                                                                                                >  5:166399/1‑149 (MQ=255) tGAACAGCCCCTGCAGCGAGCGCGACAGGCCATCCTGGCTGGTGTCGACGCTGCTTTCCACG‑TAGTTGCGCACGGTGGTGATCTCGCTCGGCTGGCCTTCCAGGGTCTCGCTGATCAACTGTTTGCTGCTCTTGCCCACATCCTGCGAg                                                                                                                                                >  1:196048/1‑149 (MQ=255) tGAACAGCCCCTGCAGCGAGCGCGACAGGCCATCCTGGCTGGTGTCGACGCTGCTTTCCACG‑TAGTTGCGCACGGTGGTGATCTCGCTCGGCTGGCCTTCCAGGGTCTCGCTGATCAACTGTTTGCTGCTCTTGCCCACATCCTGCGAg                                                                                                                                                >  3:223547/1‑149 (MQ=255)      aGCCCCTGCAGCGAGCGCGACAGCCCATCCTGGCAGCTGTCGACGCTCCTATCCACGAAAGTTCCGCACGGTGGTGATCTCGCACGGCTGGCATTCCAGGGTCACGCTGATCAACTGTTTGCTGCACTTGCCCACATCCTGCGAGCGCt                                                                                                                                            >  1:314451/1‑149 (MQ=255)       gCCCCTGCAGCGAGCGCGACAGGCCATCCTGGCTGGTGTCGACGCTGCTTTCCACG‑TAGTTGCGCTCGGTGGTGATCTCGCTCGGCTGGCCTTCCAGGGTCTCGCTGATCAACTGTTTGCTGCTCTTGCCCACATCCTGCGAGCGCTg                                                                                                                                           >  3:226152/1‑148 (MQ=255)        ccccTGCAGCGAGCGCGACAGGCCATCCTGGCTGGTGTCGACGCTGCTTTCCACG‑TAGTTGCGCACGGTGGTGATCTCGCTCGGCTGGCCTTCCAGGGTCTCGCTGATCAACTGTTTGCTGCTCTTGCCCACATCCTGCGAGCGCTGGa                                                                                                                                         >  2:201963/1‑149 (MQ=255)        ccccTGCAGCGAGCGCGACAGGCCATCCTGGCTGGTGTCGACGCTGCTTTCCACG‑TAGTTGCGCACGGTGGTGATCTCGCTCGGCTGGCCTTCCAGGGTCTCGCTGATCAACTGTTTGCTGCTCTTGCCCACATCCTGCGAGCGCTGGa                                                                                                                                         <  6:184022/149‑1 (MQ=255)             gCAGCGAGCGCGACAGGCCATCCTGGCTGGTGTCGACGCTGCTTTCCACG‑TAGTTGCGCACGGTGGTGATCTCGCTCGGCTGGCCTTCCAGGGTCTCGCTGATCAACTGTTTGCTGCTCTTGCCCACATCCTGCGAGCGCTGGAtgttg                                                                                                                                    >  7:329981/1‑149 (MQ=255)              cAGCGAGCGCGACAGGCCATCCTGGCTGGTGTCGACGCTGCTTTCCACG‑TAGTTGCGCACGGTGGTGATCTCGCTCGGCTGGCCTTCCAGGGTCTCGCTGATCAACTGTTTGCTGCTCTTGCCCACATCCTGCGAGCGCTGGAtgttg                                                                                                                                    <  2:391384/148‑1 (MQ=255)                    gcgcGACAGGCCATCCTGGCTGGTGTCGACGCTGCTTTCCACG‑TAGTTGCGCACGGTGGTGATCTCGCTCGGCTGGCCTTCCAGGGTCTCGCTGATCAACTGTTTGCTGCTCTTGCCCACATCCTGCGAGCGCTGGATGTTGTTCATgc                                                                                                                             >  8:292694/1‑149 (MQ=255)                         aCAGGCCATCCTGGCTGGTGTCGACGCTGCTTTCCACG‑TAGTTGCGCACGGTGGTGATCTCGCTCGGCTGGCCTTCCAGGGTCTCGCTGATCAACTGTTTGCTGCTCTTGCCCACATCCTGCGAGCGCTGGATGTTGTTCATGCGCAcc                                                                                                                        >  5:45101/1‑149 (MQ=255)                           aGGCCATCCTGGCTGGTGTCGACGCTGCTTTCCACG‑TAGTTGCGCACGGTGGTGATCTCGCTCGGCTGGCCTTCCAGGGTCTCGCTGATCAACTGTTTGCTGCTCTTGCCCACATCCTGCGAGCGCTGGATGTTGTTCATGCGCACCtt                                                                                                                      <  8:329981/149‑1 (MQ=255)                             gCCATCCTGGCTGGTGTCGACGCTGCTTTCCACG‑TAGTTGCGCACGGTGGTGATCTCGCTCGGCTGGCCTTCCAGGGTCTCGCTGATCAACTGTTTGCTGCTCTTGCCCACATCCTGCGAGCGCTGGATGTTGTTCATGCGCACCTTGa                                                                                                                    <  2:250422/149‑1 (MQ=255)                              ccATCCTGGCTGGTGTCGACGCTGCTTTCCACG‑TAGTTGCGCACGGTGGTGATCTCGCTCGGCTGGCCTTCCAGGGTCTCGCTGATCAACTGTTTGCTGCTCTTGCCCACATCCTGCGAGCGCTGGATGTTGTTCATGCGCACCtt                                                                                                                      <  6:352809/146‑1 (MQ=255)                              ccATCCTGGCTGGTGTCGACGCTGCTTTCCACG‑TAGTTGCGCACGGTGGTGATCTCGCTCGGCTGGCCTTCCAGGGTCTCGCTGATCAACTGTTTGCTGCTCTTGCCCACATCCTGCGAGCGCTGGATGTTGTTCATGCGCACCtt                                                                                                                      >  5:352809/1‑146 (MQ=255)                               cATCCTGGCTGGTGTCGACGCTGCTTTCCACG‑TAGTTGCGCACGGTGGTGATCTCGCTCGGCTGGCCTTCCAGGGTCTCGCTGATCAACTGTTTGCTGCTCTTGCCCACATCCTGCGAGCGCTGGATGTTGTTCATGCGCACCTGGAGt                                                                                                                  >  5:418475/1‑149 (MQ=255)                                aTCCTGGCGGGTGTCTAGGCTGCTTTCCACG‑TCGTTCCGCACGGTGGTCTCCTCCGCCGGCT‑GCCTTCCAGGGTCTCGCTGATCAACTGTTTGCTGCTCTTGCCCACATCCTGCGAGCGCTGGATGTTGTTCATGCGCACCTTGAGttt                                                                                                                <  3:101839/149‑1 (MQ=255)                                   ctggctggTGTCGACGCTGCTTTCCACG‑TAGTTGCGCACGGTGGTGATCTCGCTCGGCTGGCCTTCCAGGGTCTCGCTGATCAACTGTTTGCTGCTCTTGCCCACATCCTGCGAGCGCTGGATGTTGTTCATGCGCACCTTGAGTTTg                                                                                                               <  2:367506/148‑1 (MQ=255)                                        tggTGTCGACGCTGCTTTCCACG‑TAGTTGCGCACGGTGGTGATCTCGCTCGGCTGGCCTTCCAGGGTCTCGCTGATCAACTGTTTGCTGCTCTTGCCCACATCCTGCGAGCGCTGGATGTTGTTCATGCGCACCTTGAGTTTGCGCAAg                                                                                                         <  2:196048/149‑1 (MQ=255)                                               gACGCTGCTTTCCACG‑TAGTTGCGCACGGTGGTGATCTCGCTCGGCTGGCCTTCCAGGGTCTCGCTGATCAACTGTTTGCTGCTCTTGCCCACATCCTGCGAGCGCTGGATGTTGTTCATGCGCACCTTGAGTTTGCGCAAggcggc                                                                                                    >  8:276355/1‑147 (MQ=255)                                                aCGCTGCTTTCCACG‑TAGTTGCGCACGGTGGTGATCTCGCTCGGCTGGCCTTCCAGGGTCTCGCTGATCAACTGTTTGCTGCTCTTGCCCACATCCTGCGAGCGCTGGATGTTGTTCATGCGCa                                                                                                                          >  1:206767/1‑124 (MQ=255)                                                aCGCTGCTTTCCACG‑TAGTTGCGCACGGTGGTGATCTCGCTCGGCTGGCCTTCCAGGGTCTCGCTGATCAACTGTTTGCTGCTCTTGCCCACATCCTGCGAGCGCTGGATGTTGTTCATGCGCa                                                                                                                          <  2:206767/124‑1 (MQ=255)                                                aCGCTGCTTTCCACG‑TAGTTGCGCACGGTGGTGATCTCGCTCGGCTGGCCTTCCAGGGTCTCGCTGATCAACTGTTTGCTGCTCTTGCCCACATCCTGCGAGCGCTGGATGTTGTTCATGCGCACCTTGAGTTTGCGCAAGGCGGCCAg                                                                                                 <  6:166399/149‑1 (MQ=255)                                                   ctgctTTCCACG‑TAGTTGCGCACGGTGGTGATCTCGCTCGGCTGGCCTTCCAGGGTCTCGCTGATCAACTGTTTGCTGCTCTTGCCCACATCCTGCGAGCGCTGGATGTTGTTCATGCGCACCTTGAGTTTGCGCAAGGCGGCCAgg                                                                                                >  2:23712/1‑147 (MQ=255)                                                    tgctTTCCACG‑TAGTTGCGCACGGTGGTGATCTCGCTCGGCTGGCCTTCCAGGGTCTCGCTGATCAACTGTTTGCTGCTCTTGCCCACATCCTGCGAGCGCTGGATGTTGTTCATGCGCACCTTGAGTTTGCGCAAGGCGGCCa                                                                                                  <  6:78220/144‑1 (MQ=255)                                                    tgctTTCCACG‑TAGTTGCGCACGGTGGTGATCTCGCTCGGCTGGCCTTCCAGGGTCTCGCTGATCAACTGTTTGCTGCTCTTGCCCACATCCTGCGAGCGCTGGATGTTGTTCATGCGCACCTTGAGTTTGCGCAAGGCGGCCAgaga                                                                                              >  5:78220/1‑145 (MQ=255)                                                      ctTTCCACG‑TAGTTGCGCACGGTGGTGATCTCGCTCGGCTGGCCTTCCAGGGTCTCGCTGATCAACTGTTTGCTGCTCTTGCCCACATCCTGCGAGCGCTGGATGTTGTTCATGCGCACCTTGAGTTTGCGCAAGGCGGCCAGGTAg                                                                                             <  6:163461/147‑1 (MQ=255)                                                          ccACG‑TAGTTGCGCTCGGTGGTGATCTCGCTCGGCTGGCCTTCCAGGGTCTCGCTGATCAACTGTTTGCTGCTCTTGCCCACATCCTGCGAGCGCTGGATGTTGTTCATGCGCACCTTGAGTTTGCGCAAGGCGGCCAGGTAGCGATc                                                                                        <  4:226152/148‑1 (MQ=255)                                                           cACG‑TAGTTGCGCACGGTGGTGATCTCGCTCGGCTGGCCTTCCAGGGTCTCGCTGATCAACTGTTTGCTGCTCTTGCCCACATCCTGCGAGCGCTGGATGTTGTTCATGCGCACCTTGAGTTTGCGCAAGGCGGCCAGGTAGCGATCCa                                                                                      <  1:201963/149‑1 (MQ=255)                                                              g‑TAGTTGCGCACGGTGGTGATCTCGCTCGGCTGGCCTTCCAGGGTCTCGCTGATCAACTGTTTGCTGCTCTTGCCCACATCCTGCGAGCGCTGGATGTTGTTCATGCGCACCTTGAGTTTGCGCAAGGCGGCCAGGTAGCGATCCatga                                                                                   <  5:136594/149‑1 (MQ=255)                                                                  gTTGCGCACGGTGGTGATCTCGCTCGGCTGGCCTTCCAGGGTCTCGCTGATCAACTGTTTGCTGCTCTTGCCCACATCCTGCGAGCGCTGGATGTTGTTCATGCGCACCTTGAGTTTGCGCAAGGCGGCCAGGTAGCGATCCatgatg                                                                                 <  7:292694/148‑1 (MQ=255)                                                                   ttGCGCACGGTTGTGATCTCGCTCGGCTGGCCTTCCAGGGTCTGGCTGATCAACTGTTTGCTGCTCTTGCCCACATCCTGCGAGCGCTGGATGTTGTTCATGGGCACCTTGAGTTTGCGCAAGGCGGCCAGGTAGCGATCCATGATGGt                                                                               <  2:314451/149‑1 (MQ=255)                                                                     gcgcACGGTGGTGATCTCGCTCGGCTGGCCTTCCAGGGTCTCGCTGATCAACTGTTTGCTGCTCTTGCCCACATCCTGCGAGCGCTGGATGTTGTTCATGCGCACCTTGAGTTTGCGCAAGGCGGCCAGGTAGCGATCCatgatg                                                                                 <  4:74748/145‑1 (MQ=255)                                                                     gcgcACGGTGGTGATCTCGCTCGGCTGGCCTTCCAGGGTCTCGCTGATCAACTGTTTGCTGCTCTTGCCCACATCCTGCGAGCGCTGGATGTTGTTCATGCGCACCTTGAGTTTGCGCAAGGCGGCCAGGTAGCGATCCatgatg                                                                                 >  3:74748/1‑145 (MQ=255)                                                                      cgcACGGTGGTGATCTCGCTCGGCTGGCCTTCCAGGGTCTCGCTGATCAACTGTTTGCTGCTCTTGCCCACATCCTGCGAGCGCTGGATGTTGTTCATGCGCACCTTGAGTTTGCGCAAGGCGGCCAGGTAGCGATCCATGATGGTGCt                                                                            <  4:197554/149‑1 (MQ=255)                                                                                 gATCTCGCTCGGCTGGCCTTCCAGGGTCTCGCTGATCAACTGTTTGCTGCTCTTGCCCACATCCTGCGAGCGCTGGATGTTGTTCATGCGCACCTTGAGTTTGCGCAAGGCGGCCAgg                                                                                                >  5:255968/1‑118 (MQ=255)                                                                                 gATCTCGCTCGGCTGGCCTTCCAGGGTCTCGCTGATCAACTGTTTGCTGCTCTTGCCCACATCCTGCGAGCGCTGGATGTTGTTCATGCGCACCTTGAGTTTGCGCAAGGCGGCCAgg                                                                                                <  6:255968/118‑1 (MQ=255)                                                                                    ctcgctcgGCTGGCCTTCCAGGGTCTCGCTGATCAACTGTTTGCTGCTCTTGCCCACATCCTGCGAGCGCTGGATGTTGTTCATGCGCACCTTGAGTTTGCGCAAGGCGGCCAGGTAGCGATCCATGATGGTGCTGTCGGCGCCTTCgg                                                              >  5:356608/1‑149 (MQ=255)                                                                                       gctcgGCTGGCCTTCCAGGGTCTCGCTGATCAACTGTTTGCTGCTCTTGCCCACATCCTGCGAGCGCTGGATGTTGTTCATGCGCACCTTGAGTTTGCGCAAGGCGGCCAGGTAGCGATCCATGATGGTg                                                                              <  5:257568/130‑1 (MQ=255)                                                                                       gctcgGCTGGCCTTCCAGGGTCTCGCTGATCAACTGTTTGCTGCTCTTGCCCACATCCTGCGAGCGCTGGATGTTGTTCATGCGCACCTTGAGTTTGCGCAAGGCGGCCAGGTAGCGATCCATGATGGTg                                                                              >  6:257568/1‑130 (MQ=255)                                                                                         tcgGCTGGCCTTCCAGGGTCTCGCTGATCAACTGTTTGCTGCTCTTGCCCACATCCTGCGAGCGCTGGATGTTGTTCATGCGCACCTTGAGTTTGCGCAAGGCGGCCAGGTAGCGATCCATGATGGTGCt                                                                            <  3:80414/130‑1 (MQ=255)                                                                                         tcgGCTGGCCTTCCAGGGTCTCGCTGATCAACTGTTTGCTGCTCTTGCCCACATCCTGCGAGCGCTGGATGTTGTTCATGCGCACCTTGAGTTTGCGCAAGGCGGCCAGGTAGCGATCCATGATGGTGCt                                                                            >  4:80414/1‑130 (MQ=255)                                                                                                gCCTTCCAGGGTCTCGCTGATCAACTGTTTGCTGCTCTTGCCCACATCCTGCGAGCGCTGGATGTTGTTCATGCGCACCTTGAGTTTGCGCAAGGCGGCCAGGTAGCGATCCATGATGGTGCTGTCGGCGCCTTCGGCGTTGTTTTCg                                                   <  7:194369/148‑1 (MQ=255)                                                                                                   ttCCAGGGTCTCGCTGATCAACTGTTTGCTGCTCTTGCCCACATCCTGCGAGCGCTGGAAGTTGTTCATGCGCACCTTGAGTTTGCGCAAGGCGGGCAGGTAGCGATCCAGGATGGTGCTGTCGGCGCCTTCGGCGGTGgtttgtgcgg                                               >  4:299791/1‑143 (MQ=255)                                                                                                     ccAGGGTCTCGCTGATCAACTGTTTGCTGCTCTTGCCCACATCCTGCGAGCGCTGGATGTTGTTCATGCGCACCTTGAGTTTGCGCAAGGCGGCCAGGTAGCGATCCATGATGGTGCTGTCGGCGCCTTCGGCGTTGTTTTCGGCGaa                                              >  3:56393/1‑148 (MQ=255)                                                                                                       aGGGTCTCGCTGATCAACTGTTTGCTGCTCTTGCCCACATCCTGCGAGCGCTGGATGTTGTTCATGCGCACCTTGAGTTTGCGCAAGGCGGCCAGGTAGCGATCCATGATGGTGCTGTCGgcgc                                                                    >  6:405516/1‑124 (MQ=255)                                                                                                       aGGGTCTCGCTGATCAACTGTTTGCTGCTCTTGCCCACATCCTGCGAGCGCTGGATGTTGTTCATGCGCACCTTGAGTTTGCGCAAGGCGGCCAGGTAGCGATCCATGATGGTGCTGTCGgcgc                                                                    <  5:405516/124‑1 (MQ=255)                                                                                                          gTCTCGCTGATCAACTGTTTGCTGCTCTTGCCCACATCCTGCGAGCGCTGGATGTTGTTCATGCGCACCTTGAGTTTGCGCAAGGCGGCCAGGTAGCGATCCATGATGGTGCTGTCGGCGCCTTCGGCGTTGTTTTCGGCGAACA‑CCcg                                       >  5:148933/1‑149 (MQ=255)                                                                                                           tctcGCTGATCAACTGTTTGCTGCTCTTGCCCACATCCTGCGAGCGCTGGATGTTGTTCATGCGCACCTGGAGTTTGCGCAAGGCGGCCAGGTAGCGATCCATGATGGTGCTGTCGGCGCCTTCGGCGTTGTTTTCGGCGAACA‑CCcgc                                      <  6:418475/149‑1 (MQ=255)                                                                                                               gCTGATCAACTGTTTGCTGCTCTTGCCCACATCCTGCGAGCGCTGGATGTTGTTCATGCGCACCTTGAGTTTGCGCAAGGCGGCCAGGTAGCGATCCATGATGGTGCTGTCGGCGCCTTCGGCGTTGTTTTCGGCGAACA‑CCCGCGCCa                                  >  8:51244/1‑149 (MQ=255)                                                                                                               gCTGATCAACTGTTTGCTGCTCTTGCCCACATCCTGCGAGCGCTGGATGTTGTTCATGCGCACCTCGAGTTTGCGCAAGGCGGCCAGGTAGCGATCCATGATGGTGCTGTCGGCGCCGGGGGCGTTGTTTTGGGCGAACA‑CCCCCGccc                                  >  4:10144/1‑148 (MQ=255)                                                                                                                 gtaTCAACTGTTAGCTGCGCGTGCCCACATCCTGCGAGCGCGGGATGTTGTTCATGCGCACCTTGAGTTTGCGCAAGGCGGCCAGGTAGCGATCCATGATGGTGCTGTCGGCGCCTTCGGCGTTGTTTTCGGCGAACA‑CCCGCGCCAcc                                <  8:224888/147‑1 (MQ=255)                                                                                                                    tCAACTGTTTGCTGCTCTTGCCCACATCCTGCGAGCGCTGGATGTTGTTCATGCGCACCTTGAGTTTGCGCAAGGCGGCCAGGTAGCGATCCATGATGGTGCTGTCGGCGCCTTCGGCGTTGTTTTCg                                                   >  4:90141/1‑128 (MQ=255)                                                                                                                    tCAACTGTTTGCTGCTCTTGCCCACATCCTGCGAGCGCTGGATGTTGTTCATGCGCACCTTGAGTTTGCGCAAGGCGGCCAGGTAGCGATCCATGATGGTGCTGTCGGCGCCTTCGGCGTTGTTTTCg                                                   <  3:90141/128‑1 (MQ=255)                                                                                                                        cTGTTTGCTGCTCTTGCCCACATCCTGCGAGCGCTGGATGTTGTTCATGCGCACCTTGAGTTTGCGCAAGGCGGCCAGGTAGCGATCCATGATGGTGCTGTCGGCGCCTTCGGCGTTGTTTTCGGCGAACA‑CCCGCGCCACCGGCTCg                          >  5:216109/1‑148 (MQ=255)                                                                                                                           tttGCTGCTCTTGCCCACATCCTGCGAGCGCTGGATGTTGTTCATGCGCACCTTGAGTTTGCGCAAGGCGGCCAGGTAGCGATCCATGATGGTGCTGTCGGCGCCTTCGGCGTTGTTTTCGGCGAACA‑CCCGCGCCACCGGCTCGAAgc                      >  2:39667/1‑149 (MQ=255)                                                                                                                                 gctCTTGCCCACATCCTGCGAGCGCTGGATGTTGTTCATGCGCACCTTGAGTTTGCGCAAGGCGGCCAGGTAGCGATCCATGATGGTGCTGTCGGCGCCTTCGGCGTTGTTTTCGGCGAaca‑c                                          <  1:288475/123‑1 (MQ=255)                                                                                                                                 gctCTTGCCCACATCCTGCGAGCGCTGGATGTTGTTCATGCGCACCTTGAGTTTGCGCAAGGCGGCCAGGTAGCGATCCATGATGGTGCTGTCGGCGCCTTCGGCGTTGTTTTCGGCGAaca‑c                                          >  2:288475/1‑123 (MQ=255)                                                                                                                                  ctctTGCCCACATCCTGCGAGCGCTGGATGTTGTTCATGCGCACCTTGAGTTTGCGCAAGGCGGCCAGGTAGCGATCCATGATGGTGCTGTCGGCGCCTTCGGCGTTGTTTTCGGCGAACA‑CCCGCGCCACCGGCTCGAAGCGCTTg                 <  4:56393/147‑1 (MQ=255)                                                                                                                                  ctctTGCCCACATCCTGCGAGCGCTGGATGTTGTTCATGCGCACCTTGAGTTTGCGCAAGGCGGCCAGGTAGCGATCCATGATGGTGCTGTCGGCGCCTTCGGCGTTGTTTTCGGCGAACA‑CCCGCGCCACCGGCTCGAAGCGCTTGGc               >  6:87698/1‑149 (MQ=255)                                                                                                                                  ctctTGCCCACATCCTGCGAGCGCTGGATGTTGTTCATGCGCACCTTGAGTTTGCGCAAGGCGGCCAGGTAGCGATCCATGATGGTGCTGTCGGCGCCGTAGGCGTTGTTTTAGGCGAAAA‑CCCCCGCCACCGGCTCGAAGCGCTTGGc               >  2:228028/1‑149 (MQ=255)                                                                                                                                  ctctTGCCCAAATCATTCGTCCCCTTGATGTTTTTAATGCCCACCTTGCCTTTTCGCCAGGCGGCCATGTAGCGATCCATGATTGTGCTGTCGGCGCCTTCGGCGTTGTTTTCGGCGAACA‑CCCGCGCCACCGGCTCGAAGCGCTTg                 <  8:255579/147‑1 (MQ=255)                                                                                                                                     ttGCCCACATCCTGCGAGCGCTGGATGTTGTTCATGCGCACCTTGAGTTTGCGCAAGGCGGCCAGGTAGCGATCCATGATGGTGCTGTCGGCGCCTTAGGCGTTGTTTTGGGCGAACA‑CCCGCGCCACCGGCTCGAAGCGCTTGGCCAg            >  2:119614/1‑149 (MQ=255)                                                                                                                                       gCCCACATCCTGCGAGCGCTGGATGTTGTTCATGCGCACCTTGAGTTTGCGCAAGGCGGCCAGGTAGCGATCCATGATGGTGCTGTCGGCGCCTTCGGCGTTGTTTTCGGCGAACA‑CCCGCGCCACCGGCTCGAAGCGCTTGGCCAgg           <  1:228028/148‑1 (MQ=255)                                                                                                                                        acaacaTCCAGCGCTCGCAGGATGTTGTTCATGCGCACCTTGAGTTTGCGCAAGGCGGCCAGGTAGCGATCCATGATGGTGCTGTCGGCGCCTTCGGCGTTGTTTTCGGCGa                                               >  2:265851/4‑112 (MQ=255)                                                                                                                                        acaacaTCCAGCGCTCGCAGGATGTTGTTCATGCGCACCTTGAGTTTGCGCAAGGCGGCCAGGTAGCGATCCATGATGGTGCTGTCGGCGCCTTCGGCGTTGTTTTCGGCGa                                               <  1:265851/109‑1 (MQ=255)                                                                                                                                         ccACATCCTGCGAGCGCTGGATGTTGTTCATGCGCACCTTGAGTTTGCGCAAGGCGGCCAGGTAGCGATCCATGATGGTGCTGTCGGCGCCTTCGGCGTTGTTTTCGGCGAACA‑CCCGCGCCACCGGCTCGAAGCGCTTGGCCAGGCTg        >  1:105448/1‑149 (MQ=255)                                                                                                                                         ccACATCCTGCGAGCGCTGGATGTTGTTCATGCGCACCTTGAGTTTGCGCAAGGCGGCCAGGTAGCGATCCATGATGGTGCTGTCGGCGCCTTCGGCGTTGTTTTCGGCGAACA‑CCCGCGCCACCGGCTCGAAGCGCTTGGCCAGGCTg        >  5:197373/1‑149 (MQ=255)                                                                                                                                         caacaTCCAGCGCTCGCAGGATGTTGTTCATGCGCACCTTGAGTTTGCGCAAGGCGGCCAGGTAGCGATCCATGATGGTGCTGTCGGCGCCTTCGGCGTTGTTTTCGGCGAACA‑CCCGCGCCACCGGCTCGAAGCGCTTGGCCAGGCTg        >  2:66793/3‑149 (MQ=255)                                                                                                                                            caTCCAGCGCTCGCAGGATGTTGTTCATGCGCACCTTGAGTTTGCGCAAGGCGGCCAGGTAGCGATCCATGATGGTGCTGTCGGCGCCTTCGGCGTTGTTTTCGGCGAACACCCCGCGCCACCGGCTCGAAGCGCTTGGCCAGGCTGcc      <  1:66793/149‑1 (MQ=255)                                                                                                                                             aTCCTGCGAGCGCTGGATGTTGTTCATGCGCACCTTGAGTTTGCGCAAGGCGGCCAGGTAGCGATCCATGATGGTGCTGTCGGCGCCTTAGGCGTTGTTTTCGGCGAACC‑CCCCCGCCCCCGGCTCGAAGCGCGTGGCCAGGCTGCcgt    >  4:243218/1‑149 (MQ=255)                                                                                                                                                cTGCGAGCGCTGGATGTTGTTCATGCGCACCTTGAGTTTGCGCAAGGCGGCCAGGTAGCGATCCATGATGGTGCTGTCGGCGCCTTCGGCGTTGTTTTCGGCGAACA‑CCCGCGCCACCGGCTCGAAGCGCTTGGCCAGGCTGCcgtcg  <  7:276355/148‑1 (MQ=255)                                                                                                                                                      |                                                                                                                                               TGAACAGCCCCTGCAGCGAGCGCGACAGGCCATCCTGGCTGGTGTCGACGCTGCTTTCCACG‑TAGTTGCGCACGGTGGTGATCTCGCTCGGCTGGCCTTCCAGGGTCTCGCTGATCAACTGTTTGCTGCTCTTGCCCACATCCTGCGAGCGCTGGATGTTGTTCATGCGCACCTTGAGTTTGCGCAAGGCGGCCAGGTAGCGATCCATGATGGTGCTGTCGGCGCCTTCGGCGTTGTTTTCGGCGAACA‑CCCGCGCCACCGGCTCGAAGCGCTTGGCCAGGCTGCCGTCG  >  NC_002947/3482863‑3483152

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