Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,845,658	Δ1 bp	6.1%	intergenic (+27/‑75)	coaD → / → yfhL	phosphopantetheine adenylyltransferase/4Fe‑4S cluster‑containing protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,845,658	0	T	.	6.1%	187.7 / 7.1	66	intergenic (+27/‑75)	coaD/yfhL	phosphopantetheine adenylyltransferase/4Fe‑4S cluster‑containing protein
Reads supporting (aligned to +/- strand):  ref base T (33/26);  new base . (2/2);  total (36/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.60e-01

GTTATTCGTTCATTTCCTCGACCCTGGTCCGGGAGATTGCAGCGCTGGGCGGTGATATCAGCAAGTTCGTCCACCCGGTGGTGGCCGACGCACTGACTGAACGCTTCAAGAAGCAAGCCTGACGCCTGATTGTC‑GCGCCCGC‑GTGCACTGCGGGCGCGAATGCGGCACAATTGTGCCCATTGCGTTGAACATGCCCGGGCTGTGAGCCCCGGCCGGAGTCCCCATGTCCCTGATCATCACCGACGATTGCATCAATTGCGACGTCTGCGAACCCGAGTGCCCGAACGAGGCT  >  NC_002947/5845510‑5845801                                                                                                                                                       |                                                                                                                                                gTTATTCGTTCATTTCCTCGACCCTGGTCCGGGAGATTGCAGCGCTGGGCGGTGATATCAGCAAGTTCGTCCACCCGGTGGTGGCCGACGCACTGACTGAACGCTTCAAGAAGCAAGCCTGACGCCTGATTGTC‑GCGCCCGC‑GTGCACTg                                                                                                                                                >  2:1401563/1‑150 (MQ=255) gTTATTCGTTCATTTCCTCGACCCTGGTCCGGGAGATTGCAGCGCTGGGCGGTGATATCAGCAAGTTCGTCCACCCGGTGGTGGCCGACGCACTGACTGAACGCTTCAAGAAGCAAGCCTGACGCCTGATTGTC‑GCGCCCGC‑GTGCACTg                                                                                                                                                <  2:194094/150‑1 (MQ=255)   tATTCGTTCATTTCCTCGACCCTGGTCCGGGAGATTGCAGCGCTGGGCGGTGATATCAGCAAGTTCGTCCACCCGGTGGTGGCCGACGCACTGACTGAACGCTTCAAGAAGCAAGCCTGACGCCTGATTGTC‑GCGCCCGC‑GTGCACTGCg                                                                                                                                              <  2:1334353/150‑1 (MQ=255)    aTTCGTTCATTTCCTCGACCCTGGTCCGGGAGATTGCAGCGCTGGGCGGTCATATCAGCAAGTTCTTCCACCCGGTGGTGGCCGACGCACTGACTGAACGCTTCAAGAAGCAAGCCTGACGCCTGATTGTC‑GCGCCCGC‑GTGCACTGCgg                                                                                                                                             <  1:1325541/150‑1 (MQ=255)        tttCATTTCCTCGACCCTGGTCCGGGAGATTGCAGCGCTGGGCGGTGATATCAGCAAGTTCGTCCACTCGGGGGTGGCCGACGCACTGACTGAACGCTTCAAGAAGCAAGCCTGACGCCTGATTGTC‑GCGCCCGCAGTGCAC‑GCGGgcgc                                                                                                                                         <  1:1096468/149‑1 (MQ=255)        gTTCATTTCCTCGACCCTGGTCCGGGAGATTGCAGCGCTGGGCGGTGATATCAGCAAGTTCGTCCACCCGGTGGTGGCCGACGCACTGACTGAACGCTTCAAGAAGCAAGCCTGACGCCTGATTGTC‑GCGCCCGC‑GTGCACTGCGGgcgc                                                                                                                                         <  1:862780/150‑1 (MQ=255)         ttCATTTCCTCGACCCTGGTCCGGGAGATTGCAGCGCTGGGCGGTGATATCAGCAAGTTCGTCCGCCCGGTGGTGGCCGACGCACTGACTGAACGCTTCAAGAAGCAAGCCTGACGCCTGATTGTC‑GCGCCCGC‑GTGCACTGCGGGcgcg                                                                                                                                        >  1:2028571/1‑150 (MQ=255)              ttCCTCGACCCTGGTCCGGGAGATTGCAGCGCTGGGCGGTGATATCAGCAAGTTCGTCCACCCGGTGGTGGCCGACGCACTGACTGAACGCTTCAAGAAGCAAGCCTGACGCCTGATTGTC‑GCGCCCGC‑GTGCACTGCGGGCGCGAATGc                                                                                                                                   <  2:2111618/150‑1 (MQ=255)                  tCGACCCTGGTCCGGGAGATTGCAGCGCTGGGCGGTGATATCAGCAAGTTCGTCCACCCGGTGGTGGCCGACGCACTGACTGAACGCTTCAAGAAGCAAGCCTGACGCCTGATTGTC‑GCGCCCGC‑GTGCACTGCGGGCGCGAATTCGGca                                                                                                                               >  1:2168488/1‑150 (MQ=255)                       ccTGGTCCGGGAGATTGCAGCGCTGGGCGGTGATATCAGCAAGTTCGTCCACCCGGTGGTGGCCGACGCACTGACTGAACGCTTCAAGAAGCAAGCCTGACGCCTGATTGTC‑GCGCCCGC‑GTGCACTGCGGGCGCGAATGCGGCACAAtt                                                                                                                          >  1:424491/1‑150 (MQ=255)                       ccTGGTCCGGGAGATTGCAGCGCTGGGCGGTGATATCAGCAAGTTCGTCCACCCGGTGGTGGCCGACGCACTGACTGAACGCTTCAAGAAGCAAGCCTGACGCCTGATTGTC‑GCACCCGC‑GTGCACTGCGGGCGCGAATGCGGCACAAtt                                                                                                                          <  2:547786/150‑1 (MQ=255)                       ccTGGTCCGGGAGATTGCAGCGCTGGGCGGTGATATCAGCAAGTTCGTCCAACCGGTGGTGGCCGACGCACTGACTGAACGCTTCAAGAAGCAAGCCTGACGCCTGATTGTC‑GCGCCCGC‑GTGCACTGCGGGCGCGAATGCGGCACAAtt                                                                                                                          >  2:571331/1‑150 (MQ=255)                         tGGTCCGGGAGATTGCAGCGCTGGGCGGTGATATCAGCAAGTTCGTCCACCCGGTGGTGGCCGACGCACTGACTGAACGCTTCAAGAAGCAAGCCTGACGCCTGATTGTC‑GCGCCCGC‑GTGCACTGCGGGCGCGAATGCGGCACAATtgt                                                                                                                        <  1:2022048/150‑1 (MQ=255)                              cGGGAGATTGCAGCGCTGGGCGGTGATATCAGCAAGTTCGTCCAAACGGAGGTGGCCGACGCACTGACTGAACGCTTCAAGAAGCAAGCCTGACGCCTGATTGTC‑GCGCCCGC‑GTGCACTGCGGGCGCGAATGCGGCACAATTGTGCCCa                                                                                                                   >  1:426682/1‑150 (MQ=255)                                  agaTTGCAGCGCTGGGCGGTGATATCAGCAAGTTCGTCCACCCGGTGGTGGCCGACGCACTGACTGAACGCTTCAAGAAGCAAGCCTGACGCCTGATTGTC‑GCGCCCGC‑GTGCACTGCGGGCGCGAATGCGGCACAATTGTGCCCATTGc                                                                                                               <  1:571331/150‑1 (MQ=255)                                      tGCAGCGCTGGGCGGTGATATCAGCAAGTTCGTCCACCCGGTGGTGGCCGACGCACTGACTGAACGCTTCAAGAAGCAAGCCTGACGCCTGATTGAC‑GCGAACGC‑CTGCACCGCCCACGCCAGTGCGGCACGATTGTGCCCATTGCGTTg                                                                                                           <  2:1194987/150‑1 (MQ=255)                                         aGCGCTGGGCGGTGATATCAGCAAGTTCGTCCACCCGGTGGTGGCCGACGCACTGACTGAACGCTTCAAGAAGCAAGCCTGACGCCTGATTGTC‑GCGCCCGC‑GTGCACTGCGGGCGCGAATGCGGCACAATTGTGCCCATTGCGTTGAAc                                                                                                        >  2:896952/1‑150 (MQ=255)                                          gcgcTGGGCGGTGATATCAGCAAGTTCGTCCACCCGGTGGTGGCCGACGCACTGACTGAACGCTTCAAGAAGCAAGCCTGACGCCTGATTGTC‑GCGCCCGC‑GTGCACTGCGGGCGCGAATACGGCACAATTGTGCCCATTGCGTTGAACa                                                                                                       >  1:937524/1‑150 (MQ=255)                                                   ggTGATATCAGCAAGTTCGTCCACCCGGTGGTGGCCCACGCAATCACTGAACGTTTCAACAAGCAAGCCTGACGCCTGATTGTC‑GCGCCCGC‑GTGCACTGCGGGCGCGAATGCGGCACAATTGTGCCCATTGCGTTGAAGATGCGCGGGc                                                                                              >  1:942432/1‑150 (MQ=255)                                                         atCAGCAAGTTCGTCCACCCGGTGGTGGCCGACGCACTGACTGCCCGCATCACAATGCACGCCTGGATGATG‑CGGTCGGTGTTGTC‑GTTCTTGTGGGGCTGGAAGGCGGCACAATTGTGCCCATTGCGTTGAACATGCCCGGGCTGTGAg                                                                                        <  2:2015928/150‑1 (MQ=255)                                                          ccagtccgTTCGTCCACCCGTTGGTGGCCGACGCACTGACTGAACGCTTCAAGAAGCAAGCCTGACGCCTGATTGTC‑GCGCCCGC‑GTGCACTGCGGGCGCGAATGCGGCACAATTGTGCCCATTGCGTTGAACATGCCCGGGCTGTGAGc                                                                                       <  2:1885861/143‑1 (MQ=255)                                                            aGCAAGTTCGTCCACCCGGTGGTGGCCGACGCACTGACTGAACGCTTCAAGAAGCAAGCCTGACGCCTGATTGTC‑GCGCCCGC‑GTGCACTGCGGGCGCGAATGCGGCACAATTGTGCCCATTGCGTTGAACATGCCCGGGCTGTGAGccc                                                                                     >  2:554034/1‑150 (MQ=255)                                                            aGCAAGTTCGTCCACCCGGTGGTGGCCGACGCACTGACTGAACGCTTCAAGAAGCAAGCCTGACGCCTGATTGTC‑GCGCCCGC‑GTGCACTGCGGGCGCGAATGCGGCACAATTGTGCCCATTGCGTTGAACATGCCCGGGCTGTGAGccc                                                                                     >  1:2072002/1‑150 (MQ=255)                                                            aGCAAGTTCGTCCACCCGGTGGTGGCCGACGCACTGACTGAACGCTTCAAGAAGCAAGCCTGACGCCTGATTGTC‑GCGCCCGC‑GTGCACTGCGGGCGCGAATGCGGCACAATTGTGCCCATTGCGTTGAACATGCCCGGGCTGTGAGccc                                                                                     >  1:77404/1‑150 (MQ=255)                                                             gCAAGTTCGTCCACCCGGTGGTGGCCGACGCACTGACTGAACGCTTCAAGAAGCAAGCCTGACGCCTGATTGTC‑GCGCCCGCAGTGCAC‑GCGGGCGCGa                                                                                                                                       >  1:1368738/1‑99 (MQ=39)                                                             gCAAGTTCGTCCACCCGGTGGTGGCCGACGCACTGACTGAACGCTTCAAGAAGCAAGCCTGACGCCTGATTGTC‑GCGCCCGCAGTGCAC‑GCGGGCGCGa                                                                                                                                       <  2:1368738/99‑1 (MQ=39)                                                             gCAAGTTCGTCCACCCGGTGGTGGCCGACGCACTGACTGAACGCTTCAAGAAGCAAGCCTGACGCCTGATTGTC‑GCGCCCGC‑GTGCACTGCGGGCGCGAATGCGGCACAATTGTGCCCATTGCGTTGAACATGCCCGGGCTGTGAGcccc                                                                                    <  1:554034/150‑1 (MQ=255)                                                             gCAAGTTCGTCCACCCGGTGGTGGCCGACGCACTGACTGAACGCTTCAAGAAGCAAGCCTGACGCCTGATTGTC‑GCGCCCGC‑GTGCACTGCGGGCGCGAATGCGGCACAATTGTGCCCATTGCGTTGAACATGCCCGGGCTGTGAGcccc                                                                                    >  1:1226266/1‑150 (MQ=255)                                                               aaGTTCGTCCACCCGGTGGTGGCCGACGCACTGACTGAACGCTTCAAGAAGCAAGCCTGACGCCTGATTGTC‑GCGCCCGC‑GTGCACTGCGGGCGCGAATGCGGCACAATTGTGcc                                                                                                                     <  2:494212/115‑1 (MQ=255)                                                               aaGTTCGTCCACCCGGTGGTGGCCGACGCACTGACTGAACGCTTCAAGAAGCAAGCCTGACGCCTGATTGTC‑GCGCCCGC‑GTGCACTGCGGGCGCGAATGCGGCACAATTGTGcc                                                                                                                     >  1:494212/1‑115 (MQ=255)                                                                    cGTCCACCCGGTGGTGGCCGACGCACTGACTGAACGCTTCAAGAAGCAAGCCTGACGCCTGATTGTC‑GCGCCCGC‑GTGCACTGCGGGCGCGAATGCGGCACAATTGTGCCCATTGCGTTGAACATGCCCGGGCTGTGAGCCCCGGCCGGa                                                                             <  2:424491/150‑1 (MQ=255)                                                                     gTCCACCCGGTGGTGGCCGACGCACTGACTGAACGCTTCAAGAAGCAAGCCTGACGCCTGATTGTC‑GCGCCCGC‑GTGCACTGCGGGCGCGAATTCGGCACAATTGTGCCCATTGCGTTGAACATGCCCGGGCTGTGAGCCCCGGCCGGAg                                                                            <  2:2168488/150‑1 (MQ=255)                                                                      tCCACCCGGTGGTGGCCGACGCACTGACTGAACGCTTCAAGAAGCAAGCCTGACGCCTGATTGTC‑GCGCCCGC‑GTGCACTGCGGGCGCGAATGCGGCACAATTGTGCCCATTGCGTTGAACATGCCCGGGCTGTGAGCCCCGGCCGGAGt                                                                           >  2:1316954/1‑150 (MQ=255)                                                                         aCCCGGTGGTGGCCGACGCACTGACTGAACGCTTCAAGAAGCAAGCCTGACGCCTGATTGTC‑GCGCCCGC‑GTGCACTGCGGGCGCGAATGCGGCACAATTGTGCCCATTGCGTTGAACATGCCCGGGCTGTGATCCCCGGCCGGAggctc                                                                        >  2:1052826/1‑146 (MQ=255)                                                                         aCCCGGTGGTAGCCGACGCACTGACTGAACGCTTCAAGAAGCAAGCCTGACGCCTGATTGTC‑GCGCCCGC‑GTGCACTGCGGGCGCGAATGCGGCACAATTGTGCCCATTGCGTTGAACATGCCCGGGCTGTGAGCCCCGGCCGGAGTccc                                                                        >  2:11355/1‑150 (MQ=255)                                                                            cGGTGGTGGCCGACGCACTGACTGAACGCTTCAAGAAGCAAGCCTGACGCCTGATTGTC‑GCGCCCGC‑GTGCACTGCGGGCGCGAATGCGGCACAATTGTGCCCATTGCGTTGAACATGCCCGGGCTGTGAGCCCCGGCCGGAGTCCCCAt                                                                     <  1:896952/150‑1 (MQ=255)                                                                                          gCACTGACTGAACGCTTCAAGAAGCAAGCCTGACGCCTGATTGTC‑GCGCCCGC‑GTGCACTGCGGGCGCGAATGCGGCACAATTGTGCCCATTGCGTTGAACATGCCCGGGCTGTGAGCCCCGGCCGGAGTCCCCATGTCCCTGAtcatca                                                       <  2:2028571/150‑1 (MQ=255)                                                                                            actgactgAACGCTTCAAGAAGCAAGCCTGACGCCTGATTGTC‑GCGCCCGC‑GTGCACTGCGGGCGCGAATGCGGCACAATTGTGCCCATTGCGTTGAACATGCCCGGGCTGTGAGCCCCGGCCGGa                                                                             <  2:1350996/126‑1 (MQ=255)                                                                                            actgactgAACGCTTCAAGAAGCAAGCCTGACGCCTGATTGTC‑GCGCCCGC‑GTGCACTGCGGGCGCGAATGCGGCACAATTGTGCCCATTGCGTTGAACATGCCCGGGCTGTGAGCCCCGGCCGGa                                                                             >  1:1350996/1‑126 (MQ=255)                                                                                            actgactgAACGCTTCAAGAAGCAAGCCTGACGCCTGATTGTC‑GCGCCCGC‑GTGCACTGCGGGCGCGAATGCGGCACAATTGTGCCCATTGCGTTGAACATGCCCGGGCTGTGAGCCCCGGCCGGAGTCCCCATGTCCCTGATCATCAcc                                                     >  2:920319/1‑150 (MQ=255)                                                                                                actgaACGCTTCAAGAAGCAAGCCTGACGCCTGATTGTC‑GCGCCCGC‑GTGCGCGGCGGGCGCGAATGCGGCACAATTGTGCCCATTGCGTTGAACATGCCCGGGCTGTGAGCCCCGGCCGGAGTCCCCATGTCCCTGATCATCACcgacg                                                 >  1:640288/1‑150 (MQ=255)                                                                                                   gaACGCTTCAAGAAGCAAGCCTGACGCCTGATTGTC‑GCGCCCGC‑GTGCACTGCGGGCGCGAATGCGGCACAATTGTGCCCATTGCGTTGAACATGCCCGGGCTGTGAGCCCCGGCCGGAGTCCCCACCCGCCCGATGATCACCGACTAtt                                              <  2:2322774/150‑1 (MQ=255)                                                                                                         ttCAAGAAGCAAGCCTGACGCCTGATTGTC‑GCGCCCGC‑GTGCACTGCGGGCGCGAATGCGGCACAATTGTGCCCATTGCGTTGAACATGCCCGGGCTGTGAGCCCCGGCCGGAGTCCCCATGTCCCTGATCATCACCGACGATTGCATCa                                        >  1:1538329/1‑150 (MQ=255)                                                                                                         ttCAAGAAGCAAGCCTGACGCCTGATTGTC‑GCGCCCGC‑GTGCACTGCGGGCGCGAATGCGGCACAATTGTGCCCATTGCGTTGAACATGCCCGGGCTGTGAGCCCCGGCCGGAGTCCCCATGTCCCTGATCATCACCGACGATTGCATCa                                        >  2:171305/1‑150 (MQ=255)                                                                                                             agaagCAAGCCTGACGCCTGATTGTC‑GCGCCCGC‑GTGCACTGCGGGCGCGAATGCGGCACAATTGTGCCCATTGCGTTGAACATGCCCGGGCTGTGAGCCCCGGCCGGAGTCCCCATGTCCCTGATCATCACCGACGATTGCATCAATTg                                    >  1:1992209/1‑150 (MQ=255)                                                                                                                     gcCTGACGCCTGAGTGGC‑GCGTCCGC‑GTGCACTGCGGGCGCGAATGCGGCACAATTGTGCCCATTGCGTTGAACATGCCCGGGCTGTGAGCCCCGGCCGGAGTCCCCATGTCCCTGATCATCACCGACGATTGCATCAATTGCGACGTCt                            <  2:426682/150‑1 (MQ=255)                                                                                                                       cTGACGCCTGATTGTC‑GCGCCCGC‑GTGCACTGCGGGCGCGAATGCGGCACAATTGTGCCCATTGCGTTGAACATGCCCGGGCTGTGAGCCCCGGCCGGAGTCCCCATGTCCCTGATCATCACCGACGATTGCATCAATTGCGACGTCTGc                          <  1:1401563/150‑1 (MQ=255)                                                                                                                          aCGCCTGATTGTC‑GCGCCCGC‑GTGCACTGCGGGCGCGAATGCGGCACAATTGTGCCCATTGCGTTGAACATGCCCGGGCTGTGAGCCCCGGCCGGAGTCCCCATGTCCCTGATCATCACCGACGATTGCATCAATTGCGACGTCTGCGaa                       <  1:287771/150‑1 (MQ=255)                                                                                                                            gCCTGATTGTC‑GCGCCCGC‑GTGCACTGCGGGCGCGAATGCGGCACAATTGTGCCCATTACGTTGAACATGCCCGGGCTGTGAGCCCCGGCCGGAGTCCCCATGTCCCTGATCATCACCGACGATCTCATCAATTGCGACGTCTGCGAAcc                     >  2:745990/1‑150 (MQ=255)                                                                                                                                 aTTGTC‑GCGCCGGC‑GTGCACTGCGGGCGCGAATGCGGCACAATTGTGCCCATTGCGTTGAACATGCCCGGGCTGTGAGCCCCGGCCGGAGTCCCCATGTCCCTGATCATCACCGACGATTGCATCAATTGCGACGTCTGCGAACCCGAGt                <  2:999207/150‑1 (MQ=255)                                                                                                                                 aTTGTC‑GCGCCCGC‑GTGCACTGCGGGCGCGAATGCGGCAGAATTGTGCCCATTGCGTTGAACATGCCCGGGCTGTGAGCCCCGGCCGGAGGCCCCATGTCCCTGATCATCACCGACGATTGCATCAATTGCGACGTCTGCGAACCCGAGt                >  2:2083013/1‑150 (MQ=255)                                                                                                                                 aTTGAC‑CCGCCCGC‑GTCCACTGCGGGCGCGAATGCGGCACAATTGTGCCCATTGCGTTGAACATGCCCGGGCTGTGAGCCCCGGCCGGAGTCCCCATGTCCCTGATCATCACCGACGATTGCATCAATTGCGACGTCTGCGAACCCGAGt                <  1:860501/150‑1 (MQ=255)                                                                                                                                  ttGTC‑GCGCCCGC‑GTGCACTGCGGGCGCGAATGCGGCACAATTGTGCCCATTGCGTTGAACATGCCCGGGCTGTGAGCCCCGGCCGGAGTCCCCATGTCCCTGATCATCACCGACGATTGCATCAATTGCGACGTCTGCGAACCCGAGTg               >  2:1728843/1‑150 (MQ=255)                                                                                                                                        gcgcCCGCAGTGCAC‑GCGGGCGCGAATGCGGCACAATTGTGCCCATTGCGTTGAACATGCCCGGGCTGTGAGCCCCGGCCGGAGTCCCCATGTCCCTGATCATCACCGACGATTGCATCAATTGCGACGTCTGCGAACCCGAGTGCCCGa          >  1:718186/1‑150 (MQ=255)                                                                                                                                          gcCCGC‑GTGCACTGCGGGCGCGAATGCGGCACAATTGTGCCCATTGCGTTGAACATGCCCGGGCTGTGAGCCCCGGCCGGAGTCCCCATGTCCCTGATCATCACCGACGATTGCATCAATTGCGACGTCTGCGAACCCGAGTGCCCGAAc        >  1:2304290/1‑150 (MQ=255)                                                                                                                                            ccGC‑GTGCACTGCGGGCGCGAATGCGGCACAATTGTGCCCATTGCGTTGAAGATGCCCGGGCTGTGAGCCCCGGCCGGAGTCCCCATGTCCATGATCAGAACCGACGATTGCATCAATTGCGACGTCTGCGAACCCGAGTGCCCGAACGa      >  1:1813865/1‑150 (MQ=255)                                                                                                                                            ccGC‑GTGCACTGCGGGCGCGAATGCGGCACAATTGTGCCCATTGCGTTGAACATGCCCGGGCTGTGAGCCCCGGCCGGAGTCCCCATGTCCCTGATCATCACCGACGATTGCATCAATTGCGACGTCTGCGAACCCGAGTGCCCGAACGa      <  2:1698103/150‑1 (MQ=255)                                                                                                                                              gc‑gTGCACTGCGGGCGCGAATGCGGCACAATTGTGCCCATTGCGTTGAACATGCCCGGGCTGTGAGCCCCGGCCGGAGTCCCCATGTCCCTGATCATCACCGACGATTGCATCAATTGCGACGTCTGCGAACCCGAGTGCCCGAACGAg     <  2:2084568/149‑1 (MQ=255)                                                                                                                                              gc‑gTGCACTGCGGGCGCGAATGCGGCACAATTGTGCCCATTGCGTTGAACATGCCCGGGCTGTGAGCCCCGGCCGGAGTCCCCATGTCCCTGATCATCACCGACGATTGCATCAATTGCGACGTCTGCGAACCCGAGTGCCCGAACGAg     >  1:2084568/1‑149 (MQ=255)                                                                                                                                              cc‑gTGCACTGCGGGCGCGAATGCGGCACAATTGTGCCCATTGCGTTGAACATGCCCGGGCTGTGAGCCCCGGCCGGAGTCCCCATGTCCCTGATCATCACCGACGATTGCATCAATTGCGACGTCTGCGAACCCGAGTGCCCGAACGAgg    <  1:1052826/149‑1 (MQ=255)                                                                                                                                               c‑gTGCACTGCGGGCGCGAATGCGGCACAATTGTGCCCATTGCGTTGAACATGCCCGGGCTGTGAGCCCCGGCCGGAGTCCCCATGTCCCTGATCATCACCGACGATTGCATCAATTGCGACGTCTGCGAACCCGAGTGCCCGAACGAGGc   >  1:1703003/1‑150 (MQ=255)                                                                                                                                                 gTGCACTGCGGGCGCGAATGCGGCACAATTGTGCCCATTGCGTTGAACATGCCCGGGCTGTGAGCCCCGGCCGGAGTCCCCATGTCCCTGATCATCACCGACGATTGCATCAATTGCGACGTCTGCGAACCCGAGTGCCCGAACGAGGCt  >  2:2048218/1‑150 (MQ=255)                                                                                                                                                 gTGCACTGCGGGCGCGAATGCGGCACAATTGTGCCCATTGCGTTGAACATGCCCGGGCTGTGAGCCCCGGCCGGAGTCCCAATGTCCCTGATCATCACCGACGATTGCATCAATTGCGACGTCTGCGAACCCGAGTGCCCGAACGAGGCt  >  1:2121320/1‑150 (MQ=255)                                                                                                                                                 gTGCACTGCGGGCGCGAATGCGGCACAATTGTGCCCATTGCGTTGAACATGCCCGGGCTGTGAGCCCCGGACGGAGTCCCCATGTCCATGATCATCACCGAACATTTCATCTGTTTCGACTTCTTCTAACCCGAGTGCCCGACCGAGGCt  >  2:1550621/1‑150 (MQ=255)                                                                                                                                                    cACTGCGGGCGCGAATGCGGCACAATTGTGCCCATTGCGTTGAACATGCCCGGGCTGTGAGCCCCGGCCGGAGTCCCCATGTCCCTGATCATCACCGACGATTGCATCAATTGCGACGTCTGCGAACCCGAGTGCCCGaa         >  1:1427019/1‑140 (MQ=255)                                                                                                                                                    cACTGCGGGCGCGAATGCGGCACAATTGTGCCCATTGCGTTGAACATGCCCGGGCTGTGAGCCCCGGCCGGAGTCCCCATGTCCCTGATCATCACCGACGATTGCATCAATTGCGACGTCTGCGAACCCGAGTGCCCGaa         <  2:1427019/140‑1 (MQ=255)                                                                                                                                                       |                                                                                                                                                GTTATTCGTTCATTTCCTCGACCCTGGTCCGGGAGATTGCAGCGCTGGGCGGTGATATCAGCAAGTTCGTCCACCCGGTGGTGGCCGACGCACTGACTGAACGCTTCAAGAAGCAAGCCTGACGCCTGATTGTC‑GCGCCCGC‑GTGCACTGCGGGCGCGAATGCGGCACAATTGTGCCCATTGCGTTGAACATGCCCGGGCTGTGAGCCCCGGCCGGAGTCCCCATGTCCCTGATCATCACCGACGATTGCATCAATTGCGACGTCTGCGAACCCGAGTGCCCGAACGAGGCT  >  NC_002947/5845510‑5845801

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑