Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,466,896	C→G	7.9%	intergenic (‑28/+96)	PP_3081 ← / ← PP_5560	hypothetical protein/hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,466,896	0	C	G	7.9%	216.5 / 5.7	72	intergenic (‑28/+96)	PP_3081/PP_5560	hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (30/34);  new base G (3/3);  total (35/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTTGTCCCAGCTCAGTTGCGGGGTCCGCTCTTCACCGGCACCGCATGCCGGGTCCTGGCCGACCTGTTGCAGGGCAATGACGCCGCCATCAGTGAAAGACGAGGAACTGACCGACATCGCGCCTTGGCCGCCTGC‑GGCGGCCAACGCTGTGCAGGGCAGC‑AAAGCAGCGAACAGCCAGTGTCTGAGGTGCATGGTCAATCCTTGGCAGGGCGTCACGGCTCAGTCTAGCCCAGGTCTGCTTCAGCGGCTGTCCTGGTCCTCGGCAGTGCCTCCGCTGGCTCCGGCACC  >  NC_002947/3466752‑3467039                                                                                                                                                  |                                                                                                                                                 gTTGTCCCAGCTCAGTTGCGGGGTCCGCTCTTCACCGGCACCGCATGCCGGGTCCTGGCCGACCTGTTGCAGGGCAATGACGCCGCCATCAGTGAAAGACGAGGAACTGACCGACATCGCGCCTTGGCCGCCTGC‑GGCGGCCAACGCtgt                                                                                                                                             <  1:526904/150‑1 (MQ=255) gTTGTCCCAGCTCAGTTGCGGGGTCCGCTCTTCACCGGCACCGCATGCCGGGTCCTGGCCGACCTGTTGCAGGGCAATGACGCCGCCATCAGTGAAAGACGAGGAACTGACCGACATCGCGCCTTGGCCGCCTGC‑GGCGGCCAACGCtgt                                                                                                                                             >  1:823393/1‑150 (MQ=255) gTTGTCCCAGCTCAGTTGCGGGGTCCGCTCTTCACCGGCACCGCATGCCGGGTCCTGGCCGACCTGTTGCAGGGCAATGACGCCGCCATCAGTGAAAGACGAGGAACTGACCGACATCGCGCCTTGGCCGCCTGC‑GGCGGCCAACGCtgt                                                                                                                                             >  2:482800/1‑150 (MQ=255)     tCCCAGCTCAGTTGCGGGGTCCGCTCTTCACCGGCACCGCATGCCGGGTCCTGGCCGACCTGTTGCAGGGCAATGACGCCGCCATCAGTGAAAGACGAGGAACTGACCGACATCGCGCCTTGGCCGCCTGC‑GGCGGCCAACGCtgttcgc                                                                                                                                         >  2:400045/1‑146 (MQ=255)     tCCCAGCTCAGTTGCGGGGTCCGCTCTTCACCGGCACCGCATGCCGGGTCCTGGCCGACCTGTTGCAGGGCAATGACGCCGCCATCAGTGAAAGACGAGGAACTGACCGACATCGCGCCTTGGCCGCCTGC‑GGCGGCCAACGCTGTGCAg                                                                                                                                         >  1:1021849/1‑150 (MQ=255)       ccAGCTCAGTTGCGGGGTACGCTCTTCACCGGCACCGCATGCCGGGTCCTGGCCGACCTGTTGCAGGGCAATGACGCCGCCATCAGTGAAAGACGAGGAACTGACCGACATCGCGCCTTGGCCGCCTGC‑CGGGGACAAAGCTGTGCCggg                                                                                                                                       >  2:166771/1‑150 (MQ=255)         ctcTCAGTTGGGGGCTCCGCTCTTCCCCGGCACCGCATGCCGGGTCCTGGCCGACGTGTTGCAGGGCAATGACGCCGCCATCAGTGAAAGACGAGGAACTGAACGAACTCGCGCCTTGGCCGCCTGC‑GGCGGCGAACGCTGTGCAGGgca                                                                                                                                     <  2:170551/148‑1 (MQ=255)         aGCTCAGTTGCGGGGTCCGCTCTTCACCGGCACCGCATGCCGGGTCCTGGCCGACCTGTTGCAGGGCAATGACGCCGCCATCAGTGAAAGACGAGGAACTGACCGACATCGCGCCTTGGCCGCCTGC‑GGCGGCCAACGCTGTGCAGGgca                                                                                                                                     <  1:639631/150‑1 (MQ=255)         aGCTCAGTTGCGGGGTCCGCTCTTCACCGGCACCGCATGCCGGGTCCTGGCCGACCTGTTGCAGGGCAATGACGCCGCCATCAGTGAAAGACGAGGAACTGACCGACATCGCGCCTTGGCCGCCTGC‑GGCGGCCAACGCTGTGCAGGgca                                                                                                                                     <  1:827047/150‑1 (MQ=255)         aGCTCAGTTGCGGGGTCCGCTCTTCACCGGCACCGCATGCCGGGTCCTGGCCGACCTGTTGCAGGGCAATGACGCCGCCATCAGTGAAAGACGAGGAACTGACCGACATCGCGCCTTGGCCGCCTGC‑GGCGGCCAACGCTGTGCAGGgca                                                                                                                                     <  1:430151/150‑1 (MQ=255)            tCAGTTGCGGGGTCCGCTCTTCACCGGCACCGCATGCCGGGTCCTGGCCGACCTGTTGCAGGGCAATGACGCCGCCATCAGTGAAAGACGAGGAACTGACCGACATCGCGCCTTGGCCGCCTGC‑GGGGCACAACGCTGTCCAGGgcagag                                                                                                                                  >  1:170551/1‑148 (MQ=255)               gTTGCGGGGTCCGCTCTTCACCGGCACCGCATGCCGGGTCCTGGCCGACCTGTTGCAGGGCAATGACGCCGCCATCAGTGAAAGACGAGGAACTGACCGACATCGCGCCTTGGCCGCCTGC‑GGCGGCCAACGCTGTGCAGGGCAGC‑AAag                                                                                                                              >  2:211312/1‑150 (MQ=255)                   cGGGGTCCGCTCTTCACCGGCACCGCATGCCGGGTCCTGGCCGACCTGTTGCAGGGCAATGACGCCGCCATCAGTGAAAGACGAGGAACTGACCGACATCGCGCCTTGGCCGCCTGC‑GGCGGCCAACGCTGTGCAGGGCAGC‑AAagcagc                                                                                                                          <  1:1179814/150‑1 (MQ=255)                     gggTCCGCTCTTCACCGGCACCGCATGCCGGGTCCTGGCCGACCTGTTGCAGGGCAATGACGCCGCCATCAGTGAAAGACGAGGAACTGACCGACATCGCGCCTTGGCCGCCTGC‑GGCGGCCAACGCTGTGCAGGGCAGC‑AAAGCAGCGa                                                                                                                        >  2:775861/1‑150 (MQ=255)                           gCTCTTCACCGGCACCGCATGCCGGGTCCTGGCCGACCTGTTGCAGGGCAATGACGCCGCCATCAGTGAAAGACGAGGAACTGACCGACATCGCGCCTTGGCCGCCGCA‑GGCGGCCAAGGc                                                                                                                                                >  1:1126404/1‑121 (MQ=255)                           gCTCTTCACCGGCACCGCATGCCGGGTCCTGGCCGACCTGTTGCAGGGCAATGACGCCGCCATCAGTGAAAGACGAGGAACTGACCGACATCGCGCCTTGGCCGCCGCA‑GGCGGCCAAGGc                                                                                                                                                <  2:1126404/121‑1 (MQ=255)                            ctctTCACCGGCACCGCATGCCGGGTCCTGGCCGACCTGTTGCAGGGCAATGACGCCGCCATCAGTGAAAGACGAGGAACTGACCGACATCGCGCCTTGGCCGCCTGC‑GGCGGCCAACGCtg                                                                                                                                              >  1:1057232/1‑122 (MQ=255)                            ctctTCACCGGCACCGCATGCCGGGTCCTGGCCGACCTGTTGCAGGGCAATGACGCCGCCATCAGTGAAAGACGAGGAACTGACCGACATCGCGCCTTGGCCGCCTGC‑GGCGGCCAACGCtg                                                                                                                                              <  2:1057232/122‑1 (MQ=255)                               ttCACCGGCACCGCATGCCGGGTCCTGGCCGACCTGTTGCAGGGCAATGACGCCGCCATCAGTGAAAGACGAGGAACTGACCGACATCGCGCCTTGGCCGCCTGC‑GGCGGCCAACGCTGTGCAGGGCAGC‑AAAGCAGCGAACAGCCAgtg                                                                                                              >  1:259027/1‑150 (MQ=255)                               ttCACCGGCACCGCATGCCGGGTCCTGGCCGACCTGTTGCAGGGCAATGACGCCGCCATCAGTGAAAGACGAGGAACTGACCGACATCGCGCCTTGGCCGCCGCA‑GGCGGCCAAGGc                                                                                                                                                <  2:42619/117‑1 (MQ=255)                                  aCCGGCACCGCATGCCGGGTCCTGGCCGACCTGTTGCAGGGCAATGACGCCGCCATCAGTGAAAGACGAGGAACTGACCGACATCGCGCCTTGGCCGCCTGC‑GGCGGCCAACGCTGTGCAGGGCAGC‑AAAGCAGCGAACAGCCAGTGTCt                                                                                                           >  1:864368/1‑150 (MQ=255)                                      gCACCGCATGCCGGGTCCTGGCCGACCTGTTGCAGGGCAATGACGCCGCCATCAGTGAAAGACGAGGAACTGACCGACATCGCGCCTTGGCCGCCTGC‑GGCGGCCAACGCTGTGCAGGGCAGC‑AAAGCAGCGAACAGCCAGTGTCTGAgg                                                                                                       >  2:702225/1‑150 (MQ=255)                                      gCACCGCATGCCGGGTCCTGGCCGACCTGTTGCAGGGCAATGACGCCGCCATCAGTGAAAGACGAGGAACTGACCGACATCGCGCCTTGGCCGCCTGC‑GGCGGCCAACGCTGTGCAGGGCAGC‑AAAGCAGCGAACAGCCAGTGTCTGAgg                                                                                                       >  1:1041542/1‑150 (MQ=255)                                         ccGCATGCCGGGTCCTGGCCGACCTGTTGCAGGGCAATGACGCCGCCATCAGTGAAAGACGAGGAACTGACCGACATCGCGCCTTGGCCGCCTGC‑GGCGGCCAACGCTGTGCAGGGCAGC‑AAAGCAGCGAACAGCCAGTGTCTGAGGTGc                                                                                                    <  2:964811/150‑1 (MQ=255)                                             aTGCCGGGTCCTGGCCGACCTGTTGCAGGGCAATGACGCCGCCATCAGTGAAAGACGAGGAACTGACCGACATCGCGCCTTGGCCGCCTGC‑GGCGGCCAACGCTGTGCAGGGCAGC‑AAAGCAGCGAACAGCCAGTGTCTGAGGTGCATgg                                                                                                <  1:702225/150‑1 (MQ=255)                                                 cGGGTCCTGGCCGACCTGTTGCAGGGCAATGACGCCGCCATCAGTGAAAGACGAGGAACTGACCGACATCGCGCCTTGGCCGCCGCA‑GGCGGCCAAGGcg                                                                                                                                               >  1:412361/1‑99 (MQ=255)                                                 cGGGTCCTGGCCGACCTGTTGCAGGGCAATGACGCCGCCATCAGTGAAAGACGAGGAACTGACCGACATCGCGCCTTGGCCGCCGCA‑GGCGGCAAAGGcg                                                                                                                                               <  2:412361/100‑2 (MQ=255)                                                   ggTCCTGGCCGACCTGTTGCAGGGCAATGACGCCGCCATCAGTGAAAGACGAGGAACTGACCGACATCGCGCCTTGGCCGCCTGC‑TGCGGCCAACGCTGTGCAGGGCAGC‑AAAGCAGCGAACAGCCAGTGTCTGAGGTGCATGGTCAATc                                                                                          <  2:1041542/150‑1 (MQ=255)                                                   ggTCCTGGCCGACCTGTTGCAGGGCAATGACGCCGCCATCAGTGAAAGACGAGGAACTGACCGACATCGCGCCTTGGCCGCCTGC‑GGCGGCCAACGCTGTGCAGGGCAGC‑AAAGCAGCGAACAGCCAGTGTCTGAGGTGCATGGTCAATc                                                                                          <  1:166771/150‑1 (MQ=255)                                                          gCCGACCTGTTGCAGGGCAATGACGCCGCCATCAGTGAAAGACGAGGAACTGACCGACATCGCGCCTTGGCCGCCTGC‑GGCGGCCAACGCTGTGCAg                                                                                                                                         >  1:1096802/1‑97 (MQ=255)                                                          gCCGACCTGTTGCAGGGCAATGACGCCGCCATCAGTGAAAGACGAGGAACTGACCGACATCGCGCCTTGGCCGCCTGC‑GGCGGCCAACGCTGTGCAg                                                                                                                                         <  2:1096802/97‑1 (MQ=255)                                                          gCCGACCTGTTGCAGGGCAATGACGCCGCCATCAGTGAAAGACGAGGAACTGACCGACATCGCGCCTTGGCCGCCTGC‑GGCGGCCAACGCTGTGCAGGGCAGC‑AAAGCAGCGAACAGCCAGTGTCTGAGGTGCATGGTCAATCCTTGGCa                                                                                   <  1:211312/150‑1 (MQ=255)                                                           ccGACCTGTTGCAGGGCAATGACGCCGCCATCAGTGAAAGACGAGGAACTGACCGACATCGCGCCTTGGCCGCCTGC‑GGCGGCCAACGCTGTGCAGGGCAGC‑AAAGCAGCGAACAGCCAGTGTCTGAGGTGCATGGTCAATCCTTGGCAg                                                                                  <  1:775861/150‑1 (MQ=255)                                                             gACCTGTTGCAGGGCAATGACGCCGCCATCAGTGAAAGACGAGGAACTGACCGACATCGCGCCTTGGCCGCCTGC‑GGCGGCCAACGCTGTGCAGGGCAGC‑AAAGCAGCGAACAGCCAGTGTCTGAGGTGCATGGTCAATCCTTGGCAggg                                                                                <  1:56397/150‑1 (MQ=255)                                                                 tGTTGCAGGGCAATGACGCCGCCATCAGTGAAAGACGAGGAACTGACCGACATCGCGCCTTGGCCGCCTGC‑GGCGGCCAACGCTGTGCAGGGCAGC‑AAAGCAGCGAACAGCCAGTGTCTGa                                                                                                         <  2:323271/121‑1 (MQ=255)                                                                 tGTTGCAGGGCAATGACGCCGCCATCAGTGAAAGACGAGGAACTGACCGACATCGCGCCTTGGCCGCCTGC‑GGCGGCCAACGCTGTGCAGGGCAGC‑AAAGCAGCGAACAGCCAGTGTCTGa                                                                                                         >  1:323271/1‑121 (MQ=255)                                                                         ggCAATGACGCCGCCATCAGTGAAAGACGAGGAACTGACCGACATCGCGCCTTGGCCGCCTGC‑GGCGGCCAACGCTGTGCAGGGCAGC‑AAAGCAGCGAACAGCCAGTGTCTGAGGTGCATGGTCAATCCTTGGCAGGGCGTCACGGCTCa                                                                    >  1:257303/1‑150 (MQ=255)                                                                              tGACGCCGCCATCAGTGAAAGACGAGGAACTGACCGACATCGCGCCTTGGCCGCCTGC‑GGCGGCCAACGCTGTGCAGGGCAGC‑AAagc                                                                                                                             <  1:430886/88‑1 (MQ=255)                                                                              tGACGCCGCCATCAGTGAAAGACGAGGAACTGACCGACATCGCGCCTTGGCCGCCTGC‑GGCGGCCAACGCTGTGCAGGGCAGC‑AAagc                                                                                                                             >  2:430886/1‑88 (MQ=255)                                                                                    cgccATCAGTGAAAGACGAGGAACTGACCGACATCGCGCCTCGGCCGCCTGC‑GGCGGCCAACGCTGTGCAGGGCAGC‑AAAGCAGCGAACAGCCAGTGTCTGAGGTGCATGGTCAATCCTTGGCAGGGCGTCACGGCTCAGTCTAGCCCAg                                                         <  2:823393/150‑1 (MQ=255)                                                                                     gccATCAGTGAAAGACGAGGAACTGACCGACATCGCGCCTTGGCCGCCCGC‑GGCAGCGAGAACT‑CGCGTGGCGCCGTAGGCGGGGAAGAGCCAGTGTCTGAGGTGCATGGTCAATCCTTGGCAGGGCGTCACGGCTCAGTCTAGCCCAgg                                                        >  1:1052491/1‑150 (MQ=255)                                                                                        aTCAGTGAAAGACGAGGAACTGACCGACATCGCGCCTTGGCCGCCTGC‑GGCGGCCAACGCTGTGCAGGGCAGC‑AAAGCAGCGAACAGCCAGTGTCTGAGGTGCATGGTCAATCCTTGGCAGGGCGTCACGGCTCAGTCTAGCCCAGGTCt                                                     >  1:827476/1‑150 (MQ=255)                                                                                             tGAAAGACGAGGAACTGACCGACATCGCGCCTTGGCCGCCTGC‑GGCGGCCAACGCTGTGCAGGGCAGC‑AAAGCAGCGAACAGCCAGTGTCTGAGGTGCATGGTCAATCCTTGGCAGGGCGTCACGGCTCAGTCTAGCCCAGGTCTGCTTc                                                >  2:206762/1‑150 (MQ=255)                                                                                                 aGACGAGGAACTGACCGACATCGCGCCTTGGCCGCCTGC‑GGCGGCCAACGCTGTGCAGGGCAGC‑AAAGCAGCGaa                                                                                                                       <  2:370810/75‑1 (MQ=255)                                                                                                 aGACGAGGAACTGACCGACATCGCGCCTTGGCCGCCTGC‑GGCGGCCAACGCTGTGCAGGGCAGC‑AAAGCAGCGaa                                                                                                                       >  1:370810/1‑75 (MQ=255)                                                                                                 aGACGAGGAACTGACCGACATCGCGCCTTGGCCGCCTGC‑GGCGGCCAACGCTGTGCAGGGCAGC‑AAAGCAGCGAACAGCCAGTGTCTGAGGTGCATGGTCAATCCTTGGCAGGGCGTCACGGCTCAGTCTAGCCCAGGTCTGCTTCAGCg                                            >  1:844554/1‑150 (MQ=255)                                                                                                  gACGAGGAACTGACCGACATCGCGCCTTGGCCGCCTGC‑GGCGGCCAACGCTGTGCAGGGCAGC‑AAAGCAGCGAACAGCCAGTGTCTGAGGTGCATGGTCAATCCTTGGCAGGGCGTCACGGCTCAGTCTAGCCCAGGTCTGCTTCAGCgg                                           >  2:809515/1‑150 (MQ=255)                                                                                                     gAGGAACTGACCGACATCGCGCCTTGGCCGCCTGC‑GGCGGCCAACGCTGTGCAGGGCAGC‑AAAGCAGCGAACAGCCAGTGTCTGAGGTGCATGGTCAATCCTTgg                                                                                     >  1:665638/1‑105 (MQ=255)                                                                                                     gAGGAACTGACCGACATCGCGCCTTGGCCGCCTGC‑GGCGGCCAACGCTGTGCAGGGCAGC‑AAAGCAGCGAACAGCCAGTGTCTGAGGTGCATGGTCAATCCTTgg                                                                                     <  2:665638/105‑1 (MQ=255)                                                                                                      aGGAACTGACCGACATCGCGCCTTGGCCGCCTGC‑GGCGGCCAACGCTGTGCAGGGCAGC‑AAAGCAGCGAACAGCCAGTGTCTGAGGTGCATGGTCAATCCTTGGCAGGGCGTCACGGCTCAGTCTAGCCCAGGTCTGCTTCAGCGGCTgg                                       <  2:827476/150‑2 (MQ=255)                                                                                                        gAACTGACCGACATCGCGCCTTGGCCGCCTGC‑GGCGGCCAACGCTGTGCAGGGCAGC‑AAAGCAGCGAACAGCCAGTGTCTGAGGTGCATGGTCAATCCTTGGCAGGGCGTCACGGCTCAGTCTAGCCCAGGTCTGCTTCAGCGGCTGTcc                                     >  2:913919/1‑150 (MQ=255)                                                                                                                 gACATCGCGCCTTGGCCGCCTGC‑GGCGGCCAACGCTGTGCAGGGCAGC‑AAAGCAGCGAACAGCCAGTGTCTGAGGTGCATGGTCAATCCTTGGCAGGGCGTCACGGCTCAGTCTAGCCCAGGTCTGCTTCAGCGGCTGTCCTGGTCCTCg                            <  1:267409/150‑1 (MQ=255)                                                                                                                  aCATCGCGCCTTGGCCGCCTGC‑GGCGGCCAACGCTGTGCAGGGCAGC‑AAAGCAGCGAACAGCCAGTGTCTGAGGTGCATGGTCAATCCTTGGCAGGGCGTCACGGCTCAGTCTAGCCCAGGTCTGCTTCAGCGGCTGTCCTGGTCCTCgg                           <  1:482800/150‑1 (MQ=255)                                                                                                                      cgcgCCTTGGCCGCCTGC‑GGCGGCCAACGCTGTGCAGGGCAGC‑AAAGCAGCGAACAGCCAGTGTCTGAGGTGCATGGTCAATCCTTGGCAGGGCGTCACGGCTCAGTCTAGCCCAGGTCTGCTTCAGCGGCTGTCCTGGTCCTCGGCAGt                       <  1:809515/150‑1 (MQ=255)                                                                                                                       gcgcCTTGGCCGCCTGC‑GGCGGCCAACGCTGTGCAGGGCAGC‑AAAGCAGCGAACAGc                                                                                                                   <  1:263277/57‑1 (MQ=255)                                                                                                                       gcgcCTTGGCCGCCTGC‑GGCGGCCAACGCTGTGCAGGGCAGC‑AAAGCAGCGAACAGc                                                                                                                   >  2:263277/1‑57 (MQ=255)                                                                                                                        cgcCTTGGCCGCCTGC‑GGCGGCCAACGCTGTGCAGGGCAGC‑AAAGCAGCGAACAGCCAGTGTCTGAGGTGCATGGTCAATCCTTGGCAGGGCGTCACGGCTCAGTCTAGCCCAGGTCTGCTTCAGCGGCTGTCCTGGTCCTCGGCAGTGc                     <  2:844554/150‑1 (MQ=255)                                                                                                                         gcCTTGGCCGCCTGC‑GGCGGCCAACGCTGTGCAGGGCAGC‑AAAGCAGCGAACAGCCAGTGTCTGAGGTGCATGGTCAATCCTTGGCAGGGCGTCACGGCTCAGTCTAGCCCAGGTCTGCTTCAGCGGCTGTCCTGGTCCTCGGCAGTGcc                    <  2:442895/150‑1 (MQ=255)                                                                                                                           gttGGCCGCC‑GCAGGCGGCCAACGCTGTGCAGGGCAGC‑AAAGCAGCGAACAGCCAGTGTCTGAGGTGCATGGTCAATCCTTGGCAGGGCGTCACGGCTCAGTCTAGCCCAGGTCTGCTTCAGCGGCTGTCCTGGTCCTCGGCAGTGCCTc                  >  2:204376/2‑150 (MQ=255)                                                                                                                           cTTGGCCGCCTGC‑GGCGGCCAACGCTGTGCAGGGCAGC‑AAAGCAGCGAACAGCCAGTGTCTGAGGTGCATGGTCAATCCTTGGCAGGGCGTCACGGCTCAGTCTAGCCCAGGTCTg                                                    <  1:460321/116‑1 (MQ=255)                                                                                                                           cTTGGCCGCCTGC‑GGCGGCCAACGCTGTGCAGGGCAGC‑AAAGCAGCGAACAGCCAGTGTCTGAGGTGCATGGTCAATCCTTGGCAGGGCGTCACGGCTCAGTCTAGCCCAGGTCTg                                                    >  2:460321/1‑116 (MQ=255)                                                                                                                           cTTGGCCGCCTGC‑GGCGGCCAACGCTGTGCAGGGCAGC‑AAAGCAGCGAACAGCCAGTGTCTGAGGTGCATGGTCAATCCTTGGCAGGGCGTCACGGCTCAGTCTAGCCCAGGTCTGCTTCAGCGGCTGTCCTGGTCCTCGGCAGTGCCTc                  >  1:541100/1‑150 (MQ=255)                                                                                                                            ttGGCCGCCTGC‑GGCGGCCAAGGCTGTGCAGGGCAGC‑AAAGCAGCGAACAGCCAGTGTCTGAGGTGCATGGTCAATCCTTGGCAGGGCGTCACGGCTCAGTCTAGCCCAGGTCTGCTTCAGCGGCTGTCCTGGTCATCGGCAGTGCCTcc                 >  1:678581/1‑150 (MQ=255)                                                                                                                            ttGGCCGCCTGC‑GGCGGCCAACGCTGTGCAGGGCAGC‑AAAGCAGCGAACAGCCAGTGTCTGAGGTGCATGGTCAATCCTTGGCAGGGCGTCACGGCTCAGTCTAGCCCAGGTCTGCTTCAGCGGCTGTCCTGGTCCTCGGCAGTGCCTcc                 >  1:426307/1‑150 (MQ=255)                                                                                                                            ttGGCCGCCTGC‑GGCGGCCAACGCTGTGCAGGGCAGC‑AAAGCAGCGAACAGCCAGTGTCTGAGGTGCATGGTCAATCCTTGGCAGGGCGTCACGGCTCAGTCTAGCCCAGGTCTGCTTCAGCGGCTGTCCTGGTCCTCGGCAGTGCCTcc                 >  2:1085819/1‑150 (MQ=255)                                                                                                                              ggCCGCCTGC‑GGCGGCCAACGCTGTGCAGGGCAGC‑AAAGCAGCGAACAGCCAGTGTCTGAGGTGCATGGTCAATCCTTGGCAGGGCGTCACGGCTCAGTCTAGCCCAGGTCTGCTTCAGCGGCTGTCCTGGTCCTCGGCAGTGCCTCCGc               <  1:206762/150‑1 (MQ=255)                                                                                                                               gccgccTGC‑GGCGGCCAACGCTGTGCAGGGCAGC‑AAAGCAGCGAACAGCCAGTGTCTGAGGTGCATGGTCAATCCTTGGCAGGGCGTCACGGCTCAGTCTAGCCCAGGTCTGCTTCAGCGGCTGTCCTGGTCCTCGGCAGTGCCTCCGCt              <  1:913919/150‑1 (MQ=255)                                                                                                                                    cTGC‑GGCGGCCAACGCTGTGCAGGGCAGC‑AACGCAGCGAACAGCCAGTGTCTGAGGTGCATGGTCAATCCTTGGCAGGGCGTCACGGCTCAGTCTAGCCCAGGTCTGCTTCAGCGGCTGTCCTGGTCCTCGGCAGTGCCTCCGCTGGCTc         <  2:96135/150‑1 (MQ=255)                                                                                                                                     tGC‑GGCGGCCAACGCTGTGCAGGGCAGC‑AAAGCAGCGAACAGCCAGTGTCTGAGGTGCATGGTCAATCCTTGGCAGGGCGTCACGGCTCAGTCTAGCCCAGGTCTGCTTCAGCGGCTGTCCTGGTCCTCGGCAGTGCCTCCGCTGGCTcc        <  1:64408/150‑1 (MQ=255)                                                                                                                                            ggcCAACGCTGTGCAGGGCAGC‑AAAGCAGCGAACAGCCAGTGTCTGAGGTGCATGGTCAATCCTTGGCAGGGCGTCACGGCTCAGTCTAGCCCAGGTCTGCTTCAGCGGCTGTACTGGTCCTCGGCAGTGCCTCCGCTGGCTCCGGctct  <  2:1052491/150‑4 (MQ=255)                                                                                                                                             gcCAACGCTGTGCAGGGCAGC‑AAAGCAGCGAACAGCCAGTGTCTGAGGTGCATGGTCAATCCTTGGCAGGGCGTCACGGCTCAGTCTAGCCCAGGTCTGCTTCAGCGGCTGTCCTGGTCCTCGGCa                         >  2:674715/1‑126 (MQ=255)                                                                                                                                             gcCAACGCTGTGCAGGGCAGC‑AAAGCAGCGAACAGCCAGTGTCTGAGGTGCATGGTCAATCCTTGGCAGGGCGTCACGGCTCAGTCTAGCCCAGGTCTGCTTCAGCGGCTGTCCTGGTCCTCGGCa                         <  1:674715/126‑1 (MQ=255)                                                                                                                                                  |                                                                                                                                                 GTTGTCCCAGCTCAGTTGCGGGGTCCGCTCTTCACCGGCACCGCATGCCGGGTCCTGGCCGACCTGTTGCAGGGCAATGACGCCGCCATCAGTGAAAGACGAGGAACTGACCGACATCGCGCCTTGGCCGCCTGC‑GGCGGCCAACGCTGTGCAGGGCAGC‑AAAGCAGCGAACAGCCAGTGTCTGAGGTGCATGGTCAATCCTTGGCAGGGCGTCACGGCTCAGTCTAGCCCAGGTCTGCTTCAGCGGCTGTCCTGGTCCTCGGCAGTGCCTCCGCTGGCTCCGGCACC  >  NC_002947/3466752‑3467039

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