Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,466,885	G→C	8.1%	intergenic (‑17/+107)	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,885	0	G	C	8.1%	219.0 / 2.9	71	intergenic (‑17/+107)	PP_3081/PP_5560	hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (32/33);  new base C (3/3);  total (35/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.77e-01

TACCTGCGGGCAGGTTGTCCCAGCTCAGTTGCGGGGTCCGCTCTTCACCGGCACCGCATGCCGGGTCCTGGCCGACCTGTTGCAGGGCAATGACGCCGCCATCAGTGAAAGACGAGGAACTGACCGACATCGCGCCTTGGCCGCCTGC‑GGCGGCCAACGCTGTGCAGGGCAGC‑AAAGCAGCGAACAGCCAGTGTCTGAGGTGCATGGTCAATCCTTGGCAGGGCGTCACGGCTCAGTCTAGCCCAGGTCTGCTTCAGCGGCTGTCCTGGTCCTCGGCAGTGCCTCCGCTGGCTCC  >  NC_002947/3466739‑3467033                                                                                                                                                   |                                                                                                                                                       tACCTGCGGGCAGGTTGTCCCAGCTCAGTTGCGGGGTCCGCTCTTCACCGGCACCGCATGCCGGGTCCTGGCCGACCTGTTGCAGGGCAATGACGCCGCCATCAGTGAAAGACGAGGAACTGACCGACATCGCGCCTTGGCCGCCTgc‑gg                                                                                                                                                    >  2:526904/1‑150 (MQ=255)     tGCGGGCAGGTTGTCCCAGCTCAGTTGCGGGGTCCGCTCTTCACCGGCACCGCATGCCGGGTCCTGGCCGACCTGTTGCAGGGCAATGACGCCGCCATCAGTGAAAGACGAGGAACTGACCGACATCGCGCCTTGGCCGCCTGC‑ggcggc                                                                                                                                                <  2:147717/150‑1 (MQ=255)              gTTGTCCCAGCTCAGTTGCGGGGTCCGCTCTTCACCGGCACCGCATGCCGGGTCCTGGCCGACCTGTTGCAGGGCAATGACGCCGCCATCAGTGAAAGACGAGGAACTGACCGACATCGCGCCTTGGCCGCCTGC‑GGCGGCCAACGCtgt                                                                                                                                       >  1:823393/1‑150 (MQ=255)              gTTGTCCCAGCTCAGTTGCGGGGTCCGCTCTTCACCGGCACCGCATGCCGGGTCCTGGCCGACCTGTTGCAGGGCAATGACGCCGCCATCAGTGAAAGACGAGGAACTGACCGACATCGCGCCTTGGCCGCCTGC‑GGCGGCCAACGCtgt                                                                                                                                       <  1:526904/150‑1 (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                                                                                                                                          <  2:1126404/121‑1 (MQ=255)                                        gCTCTTCACCGGCACCGCATGCCGGGTCCTGGCCGACCTGTTGCAGGGCAATGACGCCGCCATCAGTGAAAGACGAGGAACTGACCGACATCGCGCCTTGGCCGCCGCA‑GGCGGCCAAGGc                                                                                                                                          >  1:1126404/1‑121 (MQ=255)                                         ctctTCACCGGCACCGCATGCCGGGTCCTGGCCGACCTGTTGCAGGGCAATGACGCCGCCATCAGTGAAAGACGAGGAACTGACCGACATCGCGCCTTGGCCGCCTGC‑GGCGGCCAACGCtg                                                                                                                                        <  2:1057232/122‑1 (MQ=255)                                         ctctTCACCGGCACCGCATGCCGGGTCCTGGCCGACCTGTTGCAGGGCAATGACGCCGCCATCAGTGAAAGACGAGGAACTGACCGACATCGCGCCTTGGCCGCCTGC‑GGCGGCCAACGCtg                                                                                                                                        >  1:1057232/1‑122 (MQ=255)                                            ttCACCGGCACCGCATGCCGGGTCCTGGCCGACCTGTTGCAGGGCAATGACGCCGCCATCAGTGAAAGACGAGGAACTGACCGACATCGCGCCTTGGCCGCCTGC‑GGCGGCCAACGCTGTGCAGGGCAGC‑AAAGCAGCGAACAGCCAgtg                                                                                                        >  1:259027/1‑150 (MQ=255)                                            ttCACCGGCACCGCATGCCGGGTCCTGGCCGACCTGTTGCAGGGCAATGACGCCGCCATCAGTGAAAGACGAGGAACTGACCGACATCGCGCCTTGGCCGCCGCA‑GGCGGCCaaggc                                                                                                                                          >  1:42619/1‑114 (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                                                                                                                                   <  2:1096802/97‑1 (MQ=255)                                                                       gCCGACCTGTTGCAGGGCAATGACGCCGCCATCAGTGAAAGACGAGGAACTGACCGACATCGCGCCTTGGCCGCCTGC‑GGCGGCCAACGCTGTGCAg                                                                                                                                   >  1:1096802/1‑97 (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                                                                                                   >  1:323271/1‑121 (MQ=255)                                                                              tGTTGCAGGGCAATGACGCCGCCATCAGTGAAAGACGAGGAACTGACCGACATCGCGCCTTGGCCGCCTGC‑GGCGGCCAACGCTGTGCAGGGCAGC‑AAAGCAGCGAACAGCCAGTGTCTGa                                                                                                   <  2:323271/121‑1 (MQ=255)                                                                                      ggCAATGACGCCGCCATCAGTGAAAGACGAGGAACTGACCGACATCGCGCCTTGGCCGCCTGC‑GGCGGCCAACGCTGTGCAGGGCAGC‑AAAGCAGCGAACAGCCAGTGTCTGAGGTGCATGGTCAATCCTTGGCAGGGCGTCACGGCTCa                                                              >  1:257303/1‑150 (MQ=255)                                                                                           tGACGCCGCCATCAGTGAAAGACGAGGAACTGACCGACATCGCGCCTTGGCCGCCTGC‑GGCGGCCAACGCTGTGCAGGGCAGC‑AAagc                                                                                                                       >  2:430886/1‑88 (MQ=255)                                                                                           tGACGCCGCCATCAGTGAAAGACGAGGAACTGACCGACATCGCGCCTTGGCCGCCTGC‑GGCGGCCAACGCTGTGCAGGGCAGC‑AAagc                                                                                                                       <  1:430886/88‑1 (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                                                                                                                 >  1:370810/1‑75 (MQ=255)                                                                                                              aGACGAGGAACTGACCGACATCGCGCCTTGGCCGCCTGC‑GGCGGCCAACGCTGTGCAGGGCAGC‑AAAGCAGCGaa                                                                                                                 <  2:370810/75‑1 (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)                                                                                                                                                   |                                                                                                                                                       TACCTGCGGGCAGGTTGTCCCAGCTCAGTTGCGGGGTCCGCTCTTCACCGGCACCGCATGCCGGGTCCTGGCCGACCTGTTGCAGGGCAATGACGCCGCCATCAGTGAAAGACGAGGAACTGACCGACATCGCGCCTTGGCCGCCTGC‑GGCGGCCAACGCTGTGCAGGGCAGC‑AAAGCAGCGAACAGCCAGTGTCTGAGGTGCATGGTCAATCCTTGGCAGGGCGTCACGGCTCAGTCTAGCCCAGGTCTGCTTCAGCGGCTGTCCTGGTCCTCGGCAGTGCCTCCGCTGGCTCC  >  NC_002947/3466739‑3467033

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