Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,125,301	A→T	9.2%	intergenic (+30/+51)	PP_3628 → / ← PP_3629	choline/carnitine/betaine transporter family protein/hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,125,301	0	A	T	9.2%	186.5 / 8.5	65	intergenic (+30/+51)	PP_3628/PP_3629	choline/carnitine/betaine transporter family protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (31/28);  new base T (2/4);  total (33/32)
Fisher's exact test for biased strand distribution p-value = 4.27e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.95e-01

ATCGAGAAGGCCGTACTGTATGGCGCGGTGAAGCATGGGGAGGAGTTCTACATTGTGGTGCAGGAGAAGAAGGCATCGGCCAAGACAGCGGCAGAGGCCGAAGTAACCCGCAACTGACGGCGTCAACCCGAAAACCACGGGCTGCTAAGCAGCCCGTTCGCG‑ATGCGAGACCGCGCCTGCTGAAACATGGATCTGCTTCAACCTTCCAACGCCGCCGCTATCTTGGCCACATCGATCTTGCCCATCTGCATCATGGCCTCGAACGCACGCTTGGCCGCCGCACGGTCGGGGTG  >  NC_002947/4125155‑4125447                                                                                                                                                   |                                                                                                                                                    aTCGAGAAGGCCGTACTGTATGGCGCGGTGAAGCATGGGGAGGAGTTCTACATTGTGGTGCAGGAGAAGAAGGCATCGGCCAAGACAGCGGCAGAGGCCGAAGTAACCCGCAACTGACGGCGTCAACCCGAAAACCACGGGCTGCTAag                                                                                                                                                   >  4:236349/1‑149 (MQ=255)    gtgaAGGCCGTACTGTATGGCGCGGTGAAGCATGGGGAGGAGCTCTACATTGTGGTGCAGGAGAAGAAGGCATCGGCCAAGACAGCGGCAGAGGCCGAAGTAACCCGCAACTGACGGCGTCAACCCGAAAACCACGGGCTGCTAaccag                                                                                                                                                <  2:61242/147‑1 (MQ=255)      gaAGGCCGTACTGTATGGCGCGGTGAAGCATGGGGAGGAGTTCTACATTGTGGTGCAGGAGAAGAAGGCATCGGCCAAGACAGCGGCAGAGGCCGAAGTAACCCGCAACTGACGGCGTCAACCCGAAAACCACGGGCTGCTAAGCAGcc                                                                                                                                              >  5:367620/1‑149 (MQ=255)       aaGGCCGTACTGTATGGCGCGTGGAAGCATGGGGAGGATTGCTACATTGTGGTGCAGGAGAAGAACGCATCGGCCAAGACAGCGGCAGAGGCCGAAGTAACCCGCAACTGACGGCGTCAACCCGAAAACCACGGGCTGCTAAGCAGccc                                                                                                                                             <  2:317756/149‑1 (MQ=255)       aaGGCCGTACTGTATGGCGCGGTGAAGCATGGGGAGGAGTTCTACATTGTGGTGCAGGAGAAGAAGGCATCGGCCAAGACAGCGGCAGAGGCCGAAGTAACCCGCAACTGACGGCGTCAACCCGAAAACCACGGGCTGCTTAGCAGccc                                                                                                                                             <  8:121337/149‑1 (MQ=255)       aaGGCCGTACTGTATGGCGCGGTGAAGCATGGGGAGGAGTTCTACATTGTGGTGCAGGAGAAGAAGGCATCGGCCAAGACAGCGGCAGAGGCCGAAGTAACCCGCAACTGACGGCGTCAACCCGAAAACCACGGGCTGCTTAGCAGccc                                                                                                                                             <  8:116271/149‑1 (MQ=255)        aGGCCGTACTGTATGGCGCGGTGAAGCATGGGGAGGAGTTCTACATTGTGGTGCAGGAGAAGAAGGCATCGGCCAAGACAGCGGCAGAGGCCGAAGTAACCCGCAACTGACGGCGTCAACCCGAAAACCACGGGCTGCTTAGCAGCCCg                                                                                                                                            <  2:108424/149‑1 (MQ=255)           ccGTACTGTATGGCGCGGTGAAGCATGGGGAGGAGTTCTACATTGTGGTGCAGGAGAAGAAGGCATCGGCCAAGACAGCGGCAGAGGCCGAAGTAACCCGCAACTGACGGCGTCAACCCGAAAACCACGGGCTGCTAAGCAGCCCGTTc                                                                                                                                         <  5:171405/149‑1 (MQ=255)                cTGTATGGCGCGGTGAAGCATGGGGAGGAGTTCTACATTGTGGTGCAGGAGAAGAAGGCATCGGCCAAGACAGCGGCAGAGGCCGAAGTAACCCGCAACTGACGGCGTCAACCCGAAAACCACGGGCTGCTAAGCAGCCCGTTCGCG‑At                                                                                                                                   <  6:145002/149‑1 (MQ=255)                       gcgcGGCGAAGCATGGGGAGGAGTTCTACATTGTGGTGCAGGAGAAGAAGGCATCGGCCAAGACAGCGGCAGAGGCCGAAGTAACCCGCAACTGACGGCGTCAACCCGAAAACCACGGGCTGCTAagcagc                                                                                                                                               >  4:83561/1‑131 (MQ=255)                       gcgcGGCGAAGCATGGGGAGGAGTTCTACATTGTGGTGCAGGAGAAGAAGGCATCGGCCAAGACAGCGGCAGAGGCCGAAGTAACCCGCAACTGACGGCGTCAACCCGAAAACCACGGGCTGCTAagcagc                                                                                                                                               <  3:83561/131‑1 (MQ=255)                            gTGAAGCATGGGGAGGAGTTCTACATTGTTGTGCAGGAGACGAAGGCATCGGCCATGACTGCGGCATTGGCCTAAGTAAGCGGCAACTGAGGGCGTCAAACGGAAAACCAGGGGCTGCTAAGCTGCCCGTTCGCT‑ATGCGAGACCgcgc                                                                                                                       <  2:358568/149‑1 (MQ=255)                            gTGAAGCATGGGGAGGAGTTCTACATTGTGGTGCAGGAGAAGAAGGCATCGGCCAAGACAGCGGCAGAGGCCGAAGTAACCCGCAACTGACGGCGTCAACCCGAAAACCACGGGCTGCTAAGCAGCCCGTTCGCG‑ATGCGAGACCgcgc                                                                                                                       >  1:225226/1‑149 (MQ=255)                               aaGCATGGGGAGGAGTTCTACATTGTGGTGCAGGAGAAGAAGGCATCGGCCAAGACAGCGGCAGAGGCCGAAGTAACCCGCAACTGACGGCGTCAACCCGAAAACCACGGGCTGCTAAGCAGCCCGTTCGCG‑ATGCGAGACCGCGCctg                                                                                                                    <  6:367620/149‑1 (MQ=255)                                 gCATGGGGAGGAGTTCTACATTGTGGTGCAGGAGAAGAAGGCATCGGCCAAGACAGCGGCAGAGGCCGAAGTAACCCGCAACTGACGGCGTCAACACGAAAACCACGGGCTGCTAAGCAGCCCGTTCGCG‑ATGCGAGACCGCGCctgct                                                                                                                  >  3:126481/1‑149 (MQ=255)                                              ttCTACATTGTGGTGCAGGAGAAGAAGGCATCGGCCAAGACAGCGGCAGAGGCCGAAGTAACCCGCAACTGACGGCGTCAACCCGAAAACCACGGGCTGCTAAGCAGCCCGTTCGCG‑ATGCGAGACCGCGCCTGCTGAAACATgg                                                                                                         >  4:17307/1‑145 (MQ=255)                                              ttCTACATTGTGGTGCAGGAGAAGAAGGCATCGGCCAAGACAGCGGCAGAGGCCGAAGTAACCCGCAACTGACGGCGTCAACCCGAAAACCACGGGCTGCTAAGCAGCCCGTTCGCG‑ATGCGAGACCGCGCCTGCTGAAACATgg                                                                                                         <  3:17307/145‑1 (MQ=255)                                              ttCTACATTGTGGTGCAGGAGAAGAAGGCATCGGCCAAGACAGCGGCAGAGGCCGAAGTAACCCGCAACTGACGGCGTCAACCCGAAAACCACGGGCTGCTAAGCAGCCCGTTCGCG‑ATGCGAGACCGCGCCTGCTGAAACATGGATCt                                                                                                     >  3:289089/1‑149 (MQ=255)                                              ttCTACATTGTGGTGCAGGAGAAGAAGGCATCGGCCAAGACAGCGGCAGAGGCCGAAGTAACCCGCAACTGACGGCGTCAACCCAAAAACCACGGGCTGCTTAGCAGCCCGtg                                                                                                                                          >  5:58886/1‑112 (MQ=255)                                              ttCTACATTGTGGTGCAGGAGAAGAAGGCATCGGCCAAGACAGCGGCAGAGGCCGAAGTAACCCGCAACTGACGGCGTCAACCCAAAAACCACGGGCTGCTTAGCAGCCCGtg                                                                                                                                          <  6:58886/113‑2 (MQ=255)                                                cTACATTGTGGTGCAGGAGAAGAAGGCATCGGCCAAGACAGCGGCAGAGGCCGAAGTAACCCGCAACTGACGGCGTCAACCCGAAAACCACGGGCTGCTAAGCAGCCCGTTCGCG‑ATGCga                                                                                                                               >  2:179598/1‑121 (MQ=255)                                                cTACATTGTGGTGCAGGAGAAGAAGGCATCGGCCAAGACAGCGGCAGAGGCCGAAGTAACCCGCAACTGACGGCGTCAACCCGAAAACCACGGGCTGCTAAGCAGCCCGTTCGCG‑ATGCga                                                                                                                               <  1:179598/121‑1 (MQ=255)                                                cTACATTGTGGTGCAGGAGAAGAAGGCATCGGCCAAGACAGCGGCAGAGGCCGAAGTAACCCGCAACTGACGGCGTCAACCCGAAAACCACGGGCTGCTAAGCAGCCCGTTCGCG‑ATGCGAGACCGCGCCTGCTGAAACATGGATCTGc                                                                                                   >  4:175331/1‑149 (MQ=255)                                                 tACATTGTGGTGCAGGAGAAGAAGGCATCGGCCAAGACAGCGGCAGAGGCCGAAGTAACCAGAAACTGACGGCGTCAACCCGAAAACCACGGGCAGCTAAGCAGCCCGATCGCG‑ATGCGAGACCGCGCCAGCTGAAACATCGATCTGCt                                                                                                  >  5:362074/1‑149 (MQ=255)                                                  aCATTGTGGTGCAGGAGAAGAAGGCATCGGCCAAGACAGCGGCAGAGGCCGAAGTAACCCGCAACTGACGGCGTCAACCCGAAAACCACGGGCTGCTAAGCAGCCCGTTCGCG‑ATGCGAGACCGCGCCTGCTGAAACATGGATCTGCtt                                                                                                 >  6:207952/1‑149 (MQ=255)                                                           tgCAGGAGAAGAAGGCATCGGCCAAGACAGCGGCAGAGGCCGAAGTAACCCGCAACTGACGGCGTCAACCCGAAAACCACGGGCTGCTAAGCAGCCCGTTCGCG‑ATGCGAGACCGCGCCTGCTGAAACATGGATCTGCTTCAACCTTcc                                                                                        <  4:289089/149‑1 (MQ=255)                                                              aGGAGAAGAAGGCATCGGCCAAGACAGCGGCAGAGGCCGAAGTAACCCGCAACTGACGGCGTCAACCCGAAAACCACGGGCTGCTAAGCAGCCCGTTCGCG‑ATGCGAGACCGCGCCTGCTGAAACATGGATCTGCTTCAACCTTCCAAc                                                                                     >  3:84484/1‑149 (MQ=255)                                                              aGGAGAAGAAGGCATCGGCCAAGACAGCGGCAGAGGCCGAAGTAACCCGCAACTGACGGCGTCAACCCGAAAACCACGGGCTGCTAAGCAGCCCGTTCGCG‑ATGCGAGACCGCGCCTGCTGAAACATGGATCTGCTTCAACCTTCCAAc                                                                                     <  2:225226/149‑1 (MQ=255)                                                              aGGAGAAGAAGGCATCGCCCAAGACAGCGGCAGAGGCCGAAGTAACCCGCAACTGACGGCGTCAACCCGAAAACCACGGGCTGCTAAGCAGCCCGTTCGCG‑ATGCGAGACCGCGCCTGCTGAAACATGGATCTGCTTCAACCTTCCAAc                                                                                     <  4:202101/149‑1 (MQ=255)                                                                gagaAGAAGGCATCGGCCAAGACAGCGGCAGAGGCCGAAGTAACCCGCAACTGACGGCGTCAACCCGAAAACCACGGGCTGCTAAGCAGCCCGTTCGCG‑ATGCGAGACCGCGCCTGCTGAAACATGGATCTGCTTCAACCTTCCAAcgc                                                                                   >  6:58422/1‑149 (MQ=255)                                                                  gaagaaGGCATCGGCCAAGACAGCGGCAGAGGCCGAAGTAACCCGCAACTGACGGCGTCAACCCGAAAACCACGGGCTGCTAAGCAGCCCGTTCGCG‑ATGCGAGACCGCGCCTGCTGAAACATGGATCTGCTTCAACCTTCCAAcgccg                                                                                 <  4:84484/149‑1 (MQ=255)                                                                    agaagGCATCGGCCAAGACAGCGGCAGAGGCCGAAGTAACCCGCAACTGACGGCGTCAACCCGAAAACCACGGGCTGCTAAGCAGCCCGTTCGCG‑ATGCGAGACCGCGCCTGCTGAAACATGGATCTGCTTCAACCTTCCAAcgccgcc                                                                               <  3:236349/149‑1 (MQ=255)                                                                     gaagGCATCGGCCAAGACAGCGGCAGAGGCCGAAGTAACCCGCAACTGACGGCGTCAACCCGAAAACCACGGGCTGCTAAGCAGCCCGTTcgc                                                                                                                                       >  7:29576/1‑93 (MQ=255)                                                                     gaagGCATCGGCCAAGACAGCGGCAGAGGCCGAAGTAACCCGCAACTGACGGCGTCAACCCGAAAACCACGGGCTGCTAAGCAGCCCGTTcgc                                                                                                                                       <  8:29576/93‑1 (MQ=255)                                                                              ggCCAAGACAGCGGCAGAGGCCGAAGTAACCCGCAACTGACGGCGTCAACCCGAAAACCACGGGCTGCTAAGCAGCCCGTTCGCG‑ATGCGAGACCGCGCCTGCTGAAACATGGATCTGCTTCAACCTTCCAACGCCGCCGCTATCTTgg                                                                     >  2:180754/1‑149 (MQ=255)                                                                                  aaGACAGCGGCAGAGGCCGAAGTAACCCGCAACTGACGGCGTCAACCCGAAAACCACGGGCTGCTa                                                                                                                                                     >  8:154493/1‑66 (MQ=255)                                                                                  aaGACAGCGGCAGAGGCCGAAGTAACCCGCAACTGACGGCGTCAACCCGAAAACCACGGGCTGCTa                                                                                                                                                     <  7:154493/66‑1 (MQ=255)                                                                                   aGACAGCGGCAGAGGCCGAAGTAACCCGCAACTGACGGCGTCAACCCGAAAACCACGGGCTGCTAAGCAGCCCGTTCGCG‑ATGCGAGACCGCGCCTGCTGAAACATGGATCTGCTTCAACCTTCCAACGCCGCCGCTATCTTGGCcaca                                                                >  3:111790/1‑149 (MQ=255)                                                                                    gACAGCGGCAGAGGCCGAAGTAACCCGCAACTGACGGCGTCAACCCGAAAACCACGGGCTGCTAAGCAGCCCGTTCGCG‑ATGCGAGACCGCGCCTGCTGAAACATGGATCTGCTTCAAcc                                                                                            >  8:8312/1‑120 (MQ=255)                                                                                    gACAGCGGCAGAGGCCGAAGTAACCCGCAACTGACGGCGTCAACCCGAAAACCACGGGCTGCTAAGCAGCCCGTTCGCG‑ATGCGAGACCGCGCCTGCTGAAACATGGATCTGCTTCAAcc                                                                                            <  7:8312/120‑1 (MQ=255)                                                                                       aGCGGCAGAGGCCGAAGTAACCCGCAACTGACGGCGTCAACCCGAAAACCACGGGCTGCTAAGCAGCCCGTTCGCG‑ATGCGAGACCGCGCCTGCTGAAACATGGATCTGCTTCAACCTTCCAACGCCGCCGCTATCTTGGCCACATCGa                                                            >  7:206571/1‑149 (MQ=255)                                                                                       aGCGGCAGAGGCCGAAGTAACCCGCAACTGACGGCGTCAACCCGAAAACCACGGGCTGCTAAGCAGCCCGTTCGCG‑ATGCGAGACCGCGCCTGCTGAAACATGGATCTGCTTCAACCTTCCAACGCCGCCGCTATCTTGGCCACATCGa                                                            >  8:76830/1‑149 (MQ=255)                                                                                            cAGAGGCCGAAGTAACCCGCAACTGACGGCGTCAACCCGAAAACCACGGGCTGCTAAGCAGCCCGTTCGCG‑ATGCGAGACCGCGCCTGCTGAAACATGGATCTGCTTCAACCTTCCAACGCCGCCGCTATCTTGGCCACATCGATCTTg                                                       <  7:76830/149‑1 (MQ=255)                                                                                             agagGCCGAAGTAACCCGCAACTGACGGCGTCAACCCGAAAACCACGGGCTGCTAAGCAGCCCGTTCGCG‑ATGCGAGACCGCGCCTGCTGAAACATGGATCTGCTTCAACCTTCCAACGCCGCCGCTATCTTGGCCACATCGAACTTGc                                                      <  4:111790/149‑1 (MQ=255)                                                                                                 gCCGAAGTAACCCGCAACTGACGGCGTCAACCCGAAAACCACGGGCTGCTAAGCAGCCCGTTCGCGAATTGGAGACCGCGCGTGCTGAAACATGGATCTGATTCAACCTGCCAACGCCGCCGCTATCTTGGCCACATCGAGCTTGCCCa                                                   >  1:290737/1‑149 (MQ=255)                                                                                                  ccGAAGTAACCCGCAACTGACGGCGTCAACCCGAAAACCACGGGCTGCTAAGCAGCCCGTTCGCG‑ATGCGAGACCGCGCCTGCTGAAACATGGATCTGCTTCAACCTTCCAACGCCGCCGCTATCTTGGCCACATCGATCTTGCCCATc                                                 >  7:161691/1‑149 (MQ=255)                                                                                                   cGAAGTAACCCGCAACTGACGGCGTCAACCCGAATACCACGGGCTGCTAAGCTGCCCGTTCGGG‑ATGCGAGACCGCGCCTGCTGAAACATGGATCTGCTTCAACCTTCCAACGCCGGCGCTATGTTGGCCACATCGATCTTGCCCATCt                                                <  6:362074/149‑1 (MQ=255)                                                                                                            ccGCAACTGACGGCGTCAACCCGAAAACCACGGGCTGCTAAGCAGCCCGTTCGCG‑ATGCGAGACCGCGCCTGCTGAAACATGGATCTGCTTCAACCTTCCAACGCCGCCGCTATCTTGGCCACATCGATCTTGCCCATCTGCATCATgg                                       >  5:134300/1‑149 (MQ=255)                                                                                                            ccGCAACTGACGGCGTCAACCCGAAAACCACGGGCTGCTAAGCAGCCCGTTCGCG‑ATGCGAGACCGCGCCTGCTGAAACATGGATCTGCTTCAACCTTCCAACGCCGCCGCTATCTTGGCCACATCGATCTTGCCCATCTGCATCATgg                                       >  2:165657/1‑149 (MQ=255)                                                                                                            ccGCAACTGACGGCGTCAACCCGAAAACCACGGGCTGCTAAGCAGCCCGTTCGCG‑ATGCGAGACCGCGCCTGCTGAAACATGGATCTGCTTCAACCTTCCAACGCCGCCGCTATCTTGGCCACATCGATCTTGCCCATCTGCATCATgg                                       >  8:66527/1‑149 (MQ=255)                                                                                                             cGCAACTGACGGCGTCAACCCGAAAACCACGGGCTGCTAAGCAGCCCGTTCGCG‑ATGCGAGACCGCGCCTGCTGAAACATGGATCTGCTTCAACCTTCCAACGCCGCCGCTATCTTGGCCACATCGATCTTGCCCATCTGCATCATGGc                                      >  7:197837/1‑149 (MQ=255)                                                                                                                aaCTGACGGCGTCAACCCGAAAACCACGGGCTGCTAAGCAGCCCGTTCGCG‑ATGCGAGACCGCGCCTGCTGAAACATGGATCTGCTTCAACCTTCCAACGCCGCCGCTATCTTGGCCACATCGATCTTGCCCATCTGCATCATGGCCTc                                   <  7:66527/149‑1 (MQ=255)                                                                                                                     aCGGCGTCAACCCGAAAACCACGGGCTGCTAAGCAGCCCGTTCGCG‑ATGCGAGACCGCGCCTGCTGAAACATGGATCTGCTTCAACCTTCCAACGCCGCCGCTATCTTGGCCACATCGATCTTGCCCATCTGCATCATGGCCTCGAac                               <  5:207952/148‑1 (MQ=255)                                                                                                                          gTCAACCCGAAAACCACGGGCTGCTAAGCAGCCCGTTCGCG‑ATGCGAGACCGCGCCTGCTGAAACATGGCTCTGCTGCAACCTTCCAACGCCGCCGCTAGCTTGGCCACATCGCTATGGCCCATCGGCATCATGGCCTCGAacgcccgc                         >  6:356525/1‑149 (MQ=255)                                                                                                                             aaCCCGAAATCCACGGGCTGCTAATCAGCCCGTTCGCG‑ATGCGAGACCGCGCCTGCTGAAACATGGATCTGCTTCAACCTTCCAACGCCGCCGCTATCTTGGCCACATCGATCTTGCCCATCTGCATCATGGCCTCGAACGCACGCTTg                      <  4:126481/149‑1 (MQ=255)                                                                                                                             aaCCCGAAAACCACGGGCTGCTAAGCAGCCCGTTCGCG‑ATGCGAGACCGCGCCTGCTGAAACATGGATCTGCTTCAACCTTCCAACGCCGCCGCTATCTTGGCCACATCGATCTTGCCCATCTGCATCATGGCCTCGAACGCACGCTTg                      <  5:58422/149‑1 (MQ=255)                                                                                                                              aCCCGAAAACCACGGGCTGCTAAGCAGCCCGTTCGCG‑ATGCGAGACCGCGCCTGCTGAAACATGGATCTGCTTCAACCTTCCAACGCCGCCGCTATCTTGGCCACATCGATCTTGCCCATCTGCATCATGGCCTCGAACGCACGCTTgg                     <  6:134300/149‑1 (MQ=255)                                                                                                                              aCCCGAAAACCACGGGCTGCTAAGCAGCCCGTTCGCG‑ATGCGAGACCGCGCCTGCTGAAACATGGATCTGCTTCAACCTTCCAACGCCGCCGCTATCTTGGCCACATCGATCTTGCCCATCTGCATCATGGCCTCGAACGCACGCTTgg                     <  8:206571/149‑1 (MQ=255)                                                                                                                              aCCCGAAAACCACGGGCTGCTAAGCAGCCCGTTCGCG‑ATGCGAGACCGCGCCTGCTGAAACATGGATCTGCTTCAACCTTCCAACGCCGCCGCTATCTTGGCCACATCGATCTTGCCCATCTGCATCATGGCCTCGAACGCACGCTTgg                     >  1:182104/1‑149 (MQ=255)                                                                                                                                     aaCCACGGGCTGCTAAGCAGCCCGTTCGCG‑ATGCGAGACCGCGCCTGCTGAAACATGGATCTGCTTCAACCTTCCAACGCCGCCGCTATCTTGGCCACATCGATCTTGCCCATCTGCATCATGGCCTCGAACGCACGCTTGGCcgccgc              >  8:307978/1‑149 (MQ=255)                                                                                                                                      cgaaCGGGCTGCTTAGCAGCCCGTTCGCG‑ATGCGAGACCGCGCCTGCTGAAACATGGATCTGCTTCAACCTTCCAACGCCGCCGCTATCTTGGCCACATCGATCTTGCCCATCTGCATCATGGCCTCGAACGCACGCTTGGCCGCCGCa             >  6:132505/4‑149 (MQ=255)                                                                                                                                        cACGGGCTGCTAAGCAGCCCGTTCGCG‑ATGCGAGACCGCGCCTGCTGAAACATGGATCTGCTTCAACCTTCCAACGCCGCCGCTATCTTGGCCACATCGATCTTg                                                       <  7:219547/105‑1 (MQ=255)                                                                                                                                        cACGGGCTGCTAAGCAGCCCGTTCGCG‑ATGCGAGACCGCGCCTGCTGAAACATGGATCTGCTTCAACCTTCCAACGCCGCCGCTATCTTGGCCACATCGATCTTg                                                       >  8:219547/1‑105 (MQ=255)                                                                                                                                           gggCTGCTAAGCAGCCCGTTCGCG‑ATGCGAGACCGCGCCTGCTGAAACATGGATCTGCTTCAACCTTCCAACGCCGCCGCTATCTTGGCCACATCGATCTTGCCCATCTGCATCATGGCCTCGAACGCACGCTTGGCCGCCGCACGGTc        >  4:45060/1‑149 (MQ=255)                                                                                                                                              ctgctAAGCAGCCCGTTCGCG‑ATGCGAGACCGCGCCTGCTGAAACATGGATCTGCTTCAACCTTCCAACGCCGCCGCTATCTTGGCCACATCGATCTTGCCCATCTGCATCATGGCCTCGAACGCACGCTTGGCCGCCGCACGGTCggg     <  1:180754/149‑1 (MQ=255)                                                                                                                                                 ctAAGCAGCCCGTTCGCG‑ATGCGAGACCGCGCCTGCTGAAACATGGATCTGCTTCAACCTTCCAACGCCGCCGCTATCTTGGCCACATCGATCTTGCCCATCTGCATCATGGCCTCGAACGCACGCTTGGCCGCCGCACGGTCGGGGt   <  3:175331/148‑1 (MQ=255)                                                                                                                                                 ctAAGCAGCCCGTTCGCG‑ATGCGAGACCGCGCCTGCTGAAACATGGATCTGCTTCAACCTTCCAACGCCGCCGCTATCTTGGCCACATCGATCTTGCCCATCTGCATCATGGCCTCGAACGCACGCTTGGCCGCCGCACGGTCGGGGTg  >  6:272706/1‑149 (MQ=255)                                                                                                                                                   |                                                                                                                                                    ATCGAGAAGGCCGTACTGTATGGCGCGGTGAAGCATGGGGAGGAGTTCTACATTGTGGTGCAGGAGAAGAAGGCATCGGCCAAGACAGCGGCAGAGGCCGAAGTAACCCGCAACTGACGGCGTCAACCCGAAAACCACGGGCTGCTAAGCAGCCCGTTCGCG‑ATGCGAGACCGCGCCTGCTGAAACATGGATCTGCTTCAACCTTCCAACGCCGCCGCTATCTTGGCCACATCGATCTTGCCCATCTGCATCATGGCCTCGAACGCACGCTTGGCCGCCGCACGGTCGGGGTG  >  NC_002947/4125155‑4125447

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