Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,483,015	T→A	5.8%	Q970L (CAG→CTG)	PP_3091 ←	membrane protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,483,015	0	T	A	5.8%	212.9 / 2.2	69	Q970L (CAG→CTG)	PP_3091	membrane protein
Reads supporting (aligned to +/- strand):  ref base T (32/32);  new base A (2/2);  total (34/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.97e-01

AGCCCCTGCAGCGAGCGCGACAGGCCATCCTGGCTGGTGTCGACGCTGCTTTCCACG‑TAGTTGCGCACGGTGGTGATCTCGCTCGGCTGGCCTTCCAGGGTCTCGCTGATCAACTGTTTGCTGCTCTTGCCCACATCCTGCGAGCGCTGGATGTTGTTCATGCGCACCTTGAGTTTGCGCAAGGCGGCCAGGTAGCGATCCATGATGGTGCTGTCGGCGCCTTCGGCGTTGTTTTCGGCGAACA‑CCCGCGCCACCGGCTCGAAGCGCTTGGCCAGGCTGCCGTCGTCCACCGCC  >  NC_002947/3482868‑3483161                                                                                                                                                     |                                                                                                                                                    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                                                                                                                               >  5:352809/1‑146 (MQ=255)                         ccATCCTGGCTGGTGTCGACGCTGCTTTCCACG‑TAGTTGCGCACGGTGGTGATCTCGCTCGGCTGGCCTTCCAGGGTCTCGCTGATCAACTGTTTGCTGCTCTTGCCCACATCCTGCGAGCGCTGGATGTTGTTCATGCGCACCtt                                                                                                                               <  6:352809/146‑1 (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                                                                                          >  3:74748/1‑145 (MQ=255)                                                                gcgcACGGTGGTGATCTCGCTCGGCTGGCCTTCCAGGGTCTCGCTGATCAACTGTTTGCTGCTCTTGCCCACATCCTGCGAGCGCTGGATGTTGTTCATGCGCACCTTGAGTTTGCGCAAGGCGGCCAGGTAGCGATCCatgatg                                                                                          <  4:74748/145‑1 (MQ=255)                                                                 cgcACGGTGGTGATCTCGCTCGGCTGGCCTTCCAGGGTCTCGCTGATCAACTGTTTGCTGCTCTTGCCCACATCCTGCGAGCGCTGGATGTTGTTCATGCGCACCTTGAGTTTGCGCAAGGCGGCCAGGTAGCGATCCATGATGGTGCt                                                                                     <  4:197554/149‑1 (MQ=255)                                                                            gATCTCGCTCGGCTGGCCTTCCAGGGTCTCGCTGATCAACTGTTTGCTGCTCTTGCCCACATCCTGCGAGCGCTGGATGTTGTTCATGCGCACCTTGAGTTTGCGCAAGGCGGCCAgg                                                                                                         <  6:255968/118‑1 (MQ=255)                                                                            gATCTCGCTCGGCTGGCCTTCCAGGGTCTCGCTGATCAACTGTTTGCTGCTCTTGCCCACATCCTGCGAGCGCTGGATGTTGTTCATGCGCACCTTGAGTTTGCGCAAGGCGGCCAgg                                                                                                         >  5:255968/1‑118 (MQ=255)                                                                               ctcgctcgGCTGGCCTTCCAGGGTCTCGCTGATCAACTGTTTGCTGCTCTTGCCCACATCCTGCGAGCGCTGGATGTTGTTCATGCGCACCTTGAGTTTGCGCAAGGCGGCCAGGTAGCGATCCATGATGGTGCTGTCGGCGCCTTCgg                                                                       >  5:356608/1‑149 (MQ=255)                                                                                  gctcgGCTGGCCTTCCAGGGTCTCGCTGATCAACTGTTTGCTGCTCTTGCCCACATCCTGCGAGCGCTGGATGTTGTTCATGCGCACCTTGAGTTTGCGCAAGGCGGCCAGGTAGCGATCCATGATGGTg                                                                                       >  6:257568/1‑130 (MQ=255)                                                                                  gctcgGCTGGCCTTCCAGGGTCTCGCTGATCAACTGTTTGCTGCTCTTGCCCACATCCTGCGAGCGCTGGATGTTGTTCATGCGCACCTTGAGTTTGCGCAAGGCGGCCAGGTAGCGATCCATGATGGTg                                                                                       <  5:257568/130‑1 (MQ=255)                                                                                    tcgGCTGGCCTTCCAGGGTCTCGCTGATCAACTGTTTGCTGCTCTTGCCCACATCCTGCGAGCGCTGGATGTTGTTCATGCGCACCTTGAGTTTGCGCAAGGCGGCCAGGTAGCGATCCATGATGGTGCt                                                                                     >  4:80414/1‑130 (MQ=255)                                                                                    tcgGCTGGCCTTCCAGGGTCTCGCTGATCAACTGTTTGCTGCTCTTGCCCACATCCTGCGAGCGCTGGATGTTGTTCATGCGCACCTTGAGTTTGCGCAAGGCGGCCAGGTAGCGATCCATGATGGTGCt                                                                                     <  3:80414/130‑1 (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                                                            <  3:90141/128‑1 (MQ=255)                                                                                                               tCAACTGTTTGCTGCTCTTGCCCACATCCTGCGAGCGCTGGATGTTGTTCATGCGCACCTTGAGTTTGCGCAAGGCGGCCAGGTAGCGATCCATGATGGTGCTGTCGGCGCCTTCGGCGTTGTTTTCg                                                            >  4:90141/1‑128 (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                                                        <  1:265851/109‑1 (MQ=255)                                                                                                                                   acaacaTCCAGCGCTCGCAGGATGTTGTTCATGCGCACCTTGAGTTTGCGCAAGGCGGCCAGGTAGCGATCCATGATGGTGCTGTCGGCGCCTTCGGCGTTGTTTTCGGCGa                                                        >  2:265851/4‑112 (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)                                                                                                                                                   gcTGGATGTTGTTCATGCGCACCTTGAGTTTGCGCAAGGCGGCCAGGTAGCGATCCATGATGGTGCTGTCGGCGCCTTCGGCGTTGTTTTCGGCGAACA‑CCCGCGCCACCGGCTCGAAGCGCTTGGCCAGGCTGCCGTCGTCCAccgcc  <  8:114473/149‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                    AGCCCCTGCAGCGAGCGCGACAGGCCATCCTGGCTGGTGTCGACGCTGCTTTCCACG‑TAGTTGCGCACGGTGGTGATCTCGCTCGGCTGGCCTTCCAGGGTCTCGCTGATCAACTGTTTGCTGCTCTTGCCCACATCCTGCGAGCGCTGGATGTTGTTCATGCGCACCTTGAGTTTGCGCAAGGCGGCCAGGTAGCGATCCATGATGGTGCTGTCGGCGCCTTCGGCGTTGTTTTCGGCGAACA‑CCCGCGCCACCGGCTCGAAGCGCTTGGCCAGGCTGCCGTCGTCCACCGCC  >  NC_002947/3482868‑3483161

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