Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,182,077	G→T	6.8%	intergenic (‑42/‑113)	PP_1935 ← / → PP_5487	Cro/CI family transcriptional regulator/hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,182,077	0	G	T	6.8%	167.9 / 2.4	59	intergenic (‑42/‑113)	PP_1935/PP_5487	Cro/CI family transcriptional regulator/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (29/26);  new base T (2/2);  total (31/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.17e-01

GATCAAAAGCTTCCTGCGTAAGGCCTCGCTGAACACGCAGTTGGCGGAATACGCGGCCAAAGGCAATTTTGTAATTCATCACGCGTACGCGTCCAAAAGGAAGCATTCAACCTATTGCGTCGAGGCACGAACAGGCAGTATGCTGCCGGTTACGTCATTTTC‑TGCGTTTTCGCCTTTCCGAAATCTGAGTGCTCCCGCCTTTAGAGAGTCCTCGGGAACGACTTCCTGCACCCAAAAATTAGGATCGCAAGAGTGTAAATATGATAAACGCGAGTGAGTACTTTATTAAAAATTAT  >  NC_002947/2181930‑2182225                                                                                                                                                    |                                                                                                                                                      gATCAAAAGCTTCCTGCGTAAGGCCTCGCTGAACACGCAGTTGGCGGAATACGCGGCCAAAGGCAATTTTGTAATTCATCACGCGTACGCGTCCAAAAGGAAGCATTCAACCTATTGCGTCGAGGCACGAACAGGCAGTATGCTGCCgg                                                                                                                                                      >  3:29580/1‑149 (MQ=255)  aTCAAAAGCTTCCTGCGTAAGGCCTCGCTGAACACGCAGTTGGCGGAATACGCGGCCAAAGGCAATTTTGTAATTCATCACGCGTACGCGTCCAAAAGGAAGCATTCAACCTATTGCGTCGAGGCACGAACAGGCAGCATACTGCCTGt                                                                                                                                                     <  1:42553/149‑1 (MQ=255)      aaaGCTTCCTGCGTAAGGCCTCGCTGAACACGCAGTTGGCGGAATACGCGGCCAAAGGCAATTTTGTAATTCATCACGCGTACGCGTCCAAAAGGAAGCATTCAACCTATTGCGTCGAGGCACGAACAGGCAGTATGCTGCCGGTTACg                                                                                                                                                 >  8:55528/1‑149 (MQ=255)        aGCTTCCTGCGTAAGGCCTCGCTGAACACGCAGTTGGCGGAATACTCGGCCTAAGGCATTTTTGAAATTCATCACGCGTACGCGTCCAAAAGGAAGCAGTCAACCTATTGCGTCGAGGCAGGAACAGGCAGTATGCTGCCGGTTACGTc                                                                                                                                               <  1:183470/149‑1 (MQ=255)               tGCGTAAGGCCTCGCTGAACACGCAGTTGGCGGAATACGCGGCCAAAGGCAATTTTGTAATTCATCACGCGTACGCGTCCAAAAGGAAGCATTCAACCTATTGCGTCGAGGCACGAACAGGCAGTATGCTGCCGGTTACGTCATTTTC‑t                                                                                                                                       >  2:230663/1‑149 (MQ=255)                 cGTAAGGCCTCGCTGAACACGCAGTTGGCGGAATACGCGGCCAAAGGCAATTTTGTAATTCATCACGCGTACGCGTCCAAAAGGAAGCATTCAACCTATTGCGTCGAGACACGAACAGGCAGTATGCTGCCGGTTACGTCATTTTC‑t                                                                                                                                       <  1:61497/147‑1 (MQ=255)                    aaGGCCTCGCTGAACACGCAGTTGGCGGAATACGCGGCCAAAGGCAATTTTGTAATTCATCACGCGTACGCGTCCAAAAGGAAGCATTCAACCTATTGCGTCGAGGCACGAACAGGCAGTATGCTGCCGGTTACGTCAtttt                                                                                                                                          >  5:20077/1‑142 (MQ=255)                    aaGGCCTCGCTGAACACGCAGTTGGCGGAATACGCGGCCAAAGGCAATTTTGTAATTCATCACGCGTACGCGTCCAAAAGGAAGCATTCAACCTATTGCGTCGAGGCACGAACAGGCAGTATGCTGCCGGTTACGTCAtttt                                                                                                                                          <  6:20077/142‑1 (MQ=255)                     aGGCCTCGCTGAACACGCAGTTGGCGGAATACGCGGCCAAAGGCAATTTTGTAATTCATCACGCGTACGCGTCCAAAAGGAAGCATTCAACCTATTGCGTCGAGGCACGAACAGGCAGTATGCTGCCGGTTACGTCATTTTC‑TGCGttt                                                                                                                                 <  7:55528/149‑1 (MQ=255)                       gCCTCGCTGAACACGCAGTTGGCGGAATACGCGGCCAAAGGCAATTTTGTAATTCATCACGCGTACGCGTCCAAAAGGAAGCATTCAACCTATTGCGTCGAGGCACGAACAGGCAGTATGCTGCCGGTTACGTCATTTTC‑TGCGTTTTc                                                                                                                               >  7:197314/1‑149 (MQ=255)                           cGCTGAACACGCAGTTGGCGGAATACGCGGCCAAAGGCAATTTTGTAATTCATCACGCGTACGCGTCCAAAAGGAAGCATTCAACCTATTGCGTCGAGGCACGAACAGGCAGTATGCTGCCGGTTACGTCATTTTC‑TGCGTTTTCGCCt                                                                                                                           >  3:304167/1‑149 (MQ=255)                            gCTGAACACGCAGTTGGCGGAATACGCGGCCAAAGGCAATTTTGTAATTCATCACGCGTACGCGTCCAAAAGGAAGCATTCAACCTATTGCGTCGAGGCACGAACAGGCAGTATGCTGCCGGTTACGTCATTTTC‑TGCGTTTTCGCCtt                                                                                                                          >  4:212643/1‑149 (MQ=255)                                aaCACGCAGTTGGCGGAATACGCGGCCAAAGGCAATTTTGTAATTCATCACGCGTACGCGTCCAAAAGGAAGCATTCAACCTATTGCGTCGAGGCACGAACAGGCAGTATGCTGCCGGTTACGTCATTTTC‑TGCGTTTTCGCCTTTcc                                                                                                                       >  6:278526/1‑148 (MQ=255)                                aaCACGCAGTTGGCGGAATACGCGGCCAAAGGCAATTTTGTAATTCATCACGCGTACGCGTCCAAAAGGAAGCATTCAACCTATTGCGTCGAGGCACGAACAGGCAGTATGCTGCCGGTTACGTCATTTTC‑TGCGTTTTCGCCTTTc                                                                                                                        <  3:212643/147‑1 (MQ=255)                                 acacGCAGTTGGCGGAATACGCGGCCAAAGGCAATTTTGTAATTCATCACGCGTACGCGTCCAAAAGGAAGCATTCAACCTATTGCGTCGAGGCACGAACAGGCAGTATGCTGCCGGTTACGTCATTTTC‑TGCGTTTTCGCCTTTCCg                                                                                                                      <  4:29580/148‑1 (MQ=255)                                  cacGCAGTTGGCGGAATACGCGGCCAAAGGCAATTTTGTAATTCATCACGCGTACGCGTCCAAAAGGAAGCATTCAACCTATTGCGTCGAGGCACGAACAGGCAGTATGCTGCCGGTTACGTCATTTTC‑TGCGTTTTCGCCTTTCCGa                                                                                                                     >  8:292840/1‑148 (MQ=255)                                   acGCAGTTGGTGGAATACGCGGCCAAAGGCAATTTTGTAATTCATCACGCGTACGCGTCCAAAAGGAAGCATTCAACCTATTGCGTCGAGGCACGAACAGGCAGTATGCTGCCGGTTACGTCATTTTC‑TGCGTTTTCGCCTTTCCGaa                                                                                                                    >  8:7380/1‑148 (MQ=255)                                     gCAGTTGGCGGAATACGCGGCCAAAGGCAATTTTGTAATTCATCACGCGTACGCGTCCAAAAGGAAGCATTCAACCTATTGCGTCGAGGCAAGAACAGGCAGTATGCTGCCGGTTACGTCATTTTC‑TGCGTTTTCGCCTTTCCGAAATc                                                                                                                 >  7:300739/1‑149 (MQ=255)                                       aGTTGGCGGAATACGCGGCCAAAGGCAATTTTGTAATTCATCACGCGTACGCGTCCAAAAGGAAGCATTCAACCTATTGCGTCGAGGCACGAACAGGCAGTATGCTGCCTGTTACGTCATTTTC‑TGCGTTTTCGCCTTTCCGAAATCt                                                                                                                >  4:194938/1‑148 (MQ=255)                                       aGTTGGCGGAATACGCGGCCAAAGGCAATTTTGTAATTCATCAAGCGTACGCGTCCAAAAGGAAGCATTCAACCTATTGCGTCGAGGCACGAACAGGCAGTATGCTGCCGGTTACGTCATTTTC‑TGCGTTTTCGCCTTTCCGAAATCt                                                                                                                >  1:63536/1‑148 (MQ=255)                                        gTTGGCGGAATACGCGGCCAAAGGCAATTTTGTAATTCATCACGCGTACGCGTCCAAAAGGAAGCATTCAACCTATTGCGTCGAGGCACGAACAGGCAGTATGCTGCCGGTTACGTCATTTTC‑TGCGTTTTCGCCTTTCCGAAATCTGa                                                                                                              <  8:15047/149‑1 (MQ=255)                                              ggAATACGCGGCCAAAGGCAATTTTGTAATTCATCACGCGTACGCGTCCAAAAGGAAGCATTCAACCTATTGCGTCGAGGCACGAACAGGCAGTATGCTGCCGGTTACGTCATTTTC‑TGCGTTTTCGCCTTTCCGAAATCTGAGTGCTc                                                                                                        >  8:296022/1‑149 (MQ=255)                                                       ggCCAAAGGCAATTTTGTAATTCATCACGCGTACGCGTCCAAAAGGAAGCATTCAACCTATTGCGTCGAGGCACGAACAGGCAGTATGCTGCCGGTTACGTCATTTTC‑TGCGTTTTCGCCTTTCCGAAATCTGAGTGCTCCCGCCTTTa                                                                                               <  1:110672/149‑1 (MQ=255)                                                         ccAAAGGCAATTTTGTAATTCATCACGCGTACGCGTCCAAAAGGAAGCATTCAACCTATTGCGTCGAGGCACGAACAGGCAGTATGCTGCCGGTTACGTCATTTTC‑TGCGTTTTCGCCTTTCCGAAATCTGAGTGCTCCCGCCtt                                                                                                 >  3:299119/1‑145 (MQ=255)                                                         ccAAAGGCAATTTTGTAATTCATCACGCGTACGCGTCCAAAAGGAAGCATTCAACCTATTGCGTCGAGGCACGAACAGGCAGTATGCTGCCGGTTACGTCATTTTC‑TGCGTTTTCGCCTTTCCGAAATCTGAGTGCTCCCGCCTTTaga                                                                                             <  7:7380/149‑1 (MQ=255)                                                         ccAAAGGCAATTTTGTAATTCATCACGCGTACGCGTCCAAAAGGAAGCATTCAACCTATTGCGTCGAGGCACGAACAGGCAGTATGCTGCCGGTTACGTCATTTTC‑TGCGTTTTCGCCTTTCCGAAATCTGAGTGCTCCCGCCTTTaga                                                                                             <  5:235591/149‑1 (MQ=255)                                                            aatgCAATTTTGTAATTCATCAACGGTAAGGGTCCAAAATGTATCATTGAACCTATTGCGTCGCGGCACGAACATGCACTATGCTTCCGGTTAAGTCATTTTCACGCGTTTTCGCCTTTCCGAAATCTGAGTGCTCCCGCCTTTagaga                                                                                           <  1:230663/146‑1 (MQ=255)                                                            aaGGCAATTTTGTAATTCATCACGCGTACGCGTCCAAAAGGAAGCATTCAACCTATTGCGTCGAGGCACGAACAGGCAGTATGCTGCCGGTTACGTCATTTTC‑TGCGTTTTCGCCTTTCCGAAATCTGAGTGCTCCCGCCTTTagagag                                                                                          >  8:203144/1‑149 (MQ=255)                                                            aaGGCAATTTTGTAATTCATCACGCGTACGCGTCCAAAAGGAAGCATTCAACCTATTGCGTCGAGGCACGAACAGGCAGTATGCTGCCGGTTACGTCATTTTC‑TGCGTTTTCGCCTTTCCGAAATCTGAGTGCTACCGCCTTTagagag                                                                                          >  5:240353/1‑149 (MQ=255)                                                                              atcaCGCGTACGCGTCCAAAAGGAAGCATTCAACCTATTGCGTCGAGGCACGAACAGGCAGTATGCTGCCGGTTACGTCATTTTC‑TGCGTTTTCGCCTTTCCGAAATCTGAGTGCTCCCGCCTTTAGAGAGTCCTCGGGAACGACTTc                                                                         <  7:296022/148‑1 (MQ=255)                                                                                 acgcgtacgcgtCCAAAAGGAAGCATTCAACCTATTGCGTCGAGGCACGAACAGGCAGTATGCTGCCGGTTACGTCATTTTC‑TGCGTTTTCGCCTTTCCGAAATCTGAGTGCTCCCGCCTTTAGAGAGTCCTCGGGAACGACTTCCTg                                                                      <  4:287637/148‑1 (MQ=255)                                                                                  cgagtacgcgtCCAAAAGGAAGCATTCAACCTATTGCGTCGAGGCACGAACAGGCAGTATGCTGCCGGTTACGTCATTTTC‑TGCGTTTTCGCCTTTCCGAAATCTGAGTGCTCCCACCTTTAGAGAGTCCTCGGGAACGACTTCCTGCa                                                                    <  8:300739/146‑1 (MQ=255)                                                                                         gcgtCCAAAAGGAAGCATTCAACCTATTGCGTCGAGGCACGAACAGGCAGTATGCTGCCGGTTACGTCATTTTC‑TGCGTTTTCGCCTTTCCGAAATCTGAGTGCTCCCGCCTTTagagag                                                                                          <  3:250365/120‑1 (MQ=255)                                                                                         gcgtCCAAAAGGAAGCATTCAACCTATTGCGTCGAGGCACGAACAGGCAGTATGCTGCCGGTTACGTCATTTTC‑TGCGTTTTCGCCTTTCCGAAATCTGAGTGCTCCCGCCTTTagagag                                                                                          >  4:250365/1‑120 (MQ=255)                                                                                          cgtCCAAAAGGAAGCATTCAACCTATTGCGTCGAGGCACGAACAGGCAGTATGCTGCCTGTTACGTCATTTTC‑TGCGTTTTCGCCTTTCCGAAATCTGAGTGCTCCCGCCTTTAGAGAGTCCTCGGGAACGACTTCCTGCACCCaaaaa                                                            <  3:194938/149‑1 (MQ=255)                                                                                            tCCAAAAGGAAGCATTCAACCTATTGCGTCGAGGCACGAACAGGCAGTATGCTGCCGGTTACGTCATTTTC‑TGCGTTTTCGCCTTTCCGAAATCTGAGTGCTCCCGCCTTTAGAGAGTCCTCGGGAACGACt                                                                           <  2:25740/132‑1 (MQ=255)                                                                                            tCCAAAAGGAAGCATTCAACCTATTGCGTCGAGGCACGAACAGGCAGTATGCTGCCGGTTACGTCATTTTC‑TGCGTTTTCGCCTTTCCGAAATCTGAGTGCTCCCGCCTTTAGAGAGTCCTCGGGAACGACt                                                                           >  1:25740/1‑132 (MQ=255)                                                                                              cAAAAGGAAGCATTCATCCTATTGCGTCGAGGCACGAACATGCAGTATTCTGCCGGTTACGTCATTTTC‑TTCGTTTTCGCCTTTCCGAAATCTGATTGCTCCCGCCTTTAGATAGTCCTCGGGTACTACTTCCTGCACCCAAAAATTAg                                                        <  7:203144/149‑1 (MQ=255)                                                                                               aaaaGGAAGCATTCAACCTATTGCGTCGAGGCACGAACAGGCAGTATGCTGCCGGTTACGTCATTTTC‑TGCGTTTTCGCCTTTCCGAAATCTGAGTGCTCCCGCCTTTAGAGAGTCCTCGGGAACGACTTCCTGCACCCAAAAATTAg                                                        >  5:18515/1‑148 (MQ=255)                                                                                                     aaGCATTCAACCTATTGCGTCGAGGCACGAACAGGCAGTATGCTGCCGGTTACGTCATTTTC‑TGCGTTTTCGCCTTTCCGAAATCTGAGTGCTCCCGCCTTTagaga                                                                                           >  8:24288/1‑107 (MQ=255)                                                                                                     aaGCATTCAACCTATTGCGTCGAGGCACGAACAGGCAGTATGCTGCCGGTTACGTCATTTTC‑TGCGTTTTCGCCTTTCCGAAATCTGAGTGCTCCCGCCTTTagaga                                                                                           <  7:24288/107‑1 (MQ=255)                                                                                                     aaGCATTCAACCTATTGCGTCGAGGCACGAACAGGCAGTATGCTGCCGGTTACGTCATTTTC‑TGCGTTTTCGCCTTTCCGAAATCTGAGTGCTCCCGCCTTTAGAGAGTCCTCGGGAACGACTTCCTGCACCCAAAAATTAGGATCGCa                                                 <  4:304167/149‑1 (MQ=255)                                                                                                             aaCCTATTGCGTCGAGGCACGAACAGGCAGTATGCTGCCGGTTACGTCATTTTC‑TGCGTTTTCGCCTTTCCGAAATCTGAGTGCTCCCGCCTTTAGAGAGTCCTCGGGAACGACTTCCTGCACCCAAAAATTAGGATCGCAAGAgtgt                                          <  2:75927/148‑1 (MQ=255)                                                                                                                  aTTGCGTCGAGGCACGAACAGGCAGTATGCTGCCGGTTACGTCATTTTC‑TGCGTTTTCGCCTTTCCGAAATCTGAGTGCTCCCGCCTTTAGAGAGTCCTCGGGAACGACTTCCTGCACCCAAAAATTAGGATCGCAAGAGTGTAAatat                                    >  2:165720/1‑149 (MQ=255)                                                                                                                     gCGTCGAGGCACGAACAGGCAGTATGCTGCCGGTTACGTCATTTTC‑TGCGTTTTCGCCTTTCCGAAATCTGAGTGCTCCCGCCTTTAGAGAGTCCTCGGGAACGACTTCCTGCACCCAAAAATTAGGATCGCAAGAGTGTAAATATGa                                  <  1:165720/148‑1 (MQ=255)                                                                                                                       gTCGAGGCACGAACAGGCAGTATGCTGCCGGTTACGTCATTTTC‑TGCGTTTTCGCCTTTCCGAAATCTGAGTGCTCCCGCCTTTAGAGAGTCCTCGGGAACGACTTCCTGCACCCAAAAATTAGGATCGCAAGAGTGTAAATATGATa                                >  8:40938/1‑148 (MQ=255)                                                                                                                            tgacgtaaCCGGCAGCATACTGCCTGTTACGTCATTTTC‑TGCGTTTTCGCCTTTCCGAAATCTGAGTGCTCCCGCCTTTAGAGAGTCCTCGGGAACGACTTCCTGCACCCAAAAATTAGGATCGCAAGAGTGTAAATATGATAAAcgc                           >  8:45398/7‑148 (MQ=255)                                                                                                                             gCACGAACAGGCAGTATGCTGCCGGTTACGTCATTTTC‑TGCGTTTTCGCCTTTCCGAAATCTGAGTGCTCCCGCCTTTAGAGAGTCCTCGGGAACGACTTCCTGCACCCAAAAATTAGGATCGCAAGAGTGTAAATATGATAAAcgc                           >  3:75999/1‑147 (MQ=255)                                                                                                                                 gAACAGGCAGTATGCTGCCGGTTACGTCATTTTC‑TGCGTTTTCGCCTTTCCGAAATCTGAGTGCTCCCGCCTTTAGAGAGTCCTCGGGAACGACTTCCTGCACCCAAAAATTAGGATCGCAAGAGTGTAAATATGATAAACGCgagtga                     >  7:227630/1‑149 (MQ=255)                                                                                                                                 gAACAGGCAGTATGCTGCCGGTTACGTCATTTTC‑TGCGTTTTCGCCTTTCCGAAATCTGAGTGCTCCCGCCTTTAGAGAGTCCTCGGGAACGACTTCCTGCACCCAAAAATTAGGATCGCAAGAGTGTAAATATGATAAACGCgagt                       >  6:151621/1‑147 (MQ=255)                                                                                                                                  aaCCGGCAGCATACTGCCGGTTACGTCATTTTC‑TGCGTTTTCGCCTTTCCGAAATCTGAGTGCTCCCGCCTTTAGAGAGTCCTCGGGAACGACTTCCTGCACCCAAAAATTAGGATCGCAAGAGTGTAAatat                                    >  6:204665/1‑133 (MQ=255)                                                                                                                                  aaCCGGCAGCATACTGCCGGTTACGTCATTTTC‑TGCGTTTTCGCCTTTCCGAAATCTGAGTGCTCCCGCCTTTAGAGAGTCCTCGGGAACGACTTCCTGCACCCAAAAATTAGGATCGCAAGAGTGTAAatat                                    <  5:204665/133‑1 (MQ=255)                                                                                                                                   aCAGGCAGTATGCTGCCGGTTACGTCATTTTC‑TGCGTTTTCGCCTTTCCGAAATCTGAGTGCTCCCGCCTTTAGAGAGTCCTCGGGAACGACTTCCTGCACCCAAAAATTAg                                                        <  2:231834/112‑1 (MQ=255)                                                                                                                                   aCAGGCAGTATGCTGCCGGTTACGTCATTTTC‑TGCGTTTTCGCCTTTCCGAAATCTGAGTGCTCCCGCCTTTAGAGAGTCCTCGGGAACGACTTCCTGCACCCAAAAATTAg                                                        >  1:231834/1‑112 (MQ=255)                                                                                                                                       gCAGTATGCTGCCGGTTACGTCATTTTC‑TGCGTTTTCGCCTTTCCGAAATCTGAGTGCTCCCGCCTTTAGAGAGTCCTCGGGAACGACTTCCTGCACCCAAAAATTAGGATCGCAAGAGTGTAAATATGATAAACGCGAGTGAGTACtt               >  1:292317/1‑149 (MQ=255)                                                                                                                                           tATGCTGCCGGTTACGTCATTTTC‑TGCGTTTTCGCCTTTCCGAAATCTGAGTGCTCCCGCCTTTAGAGAGTCCTCGGGAACGACTTCCTGCACCCAAAAATTAGGATCGCAAGAGTGTAAATATGATAAACGCGAGTGAGTACTttatt           <  7:166639/149‑1 (MQ=255)                                                                                                                                               ctgcCGGTTACGTCATTTTC‑TGCGTTTTCGCCTTTCCGAAATCTGAGTGCTCCCGCCTTTAGAGAGTCCTCGGGAACGACTTCCTGCACCCAAAAATTAGGATCGCAAGAGTGTAAATATGATAAACGCGAGTGAGTACTTTATTaaa        >  2:107718/1‑148 (MQ=255)                                                                                                                                               ctgcCGGTTACGTCATTTTC‑TGCGTTTTCGCCTTTCCGAAATATGAGTGCTCCCGCCTTTAGATAGTCCTCGGGAACGACTTCCTGCACCCAAAAATTAGGATCGCAAGAGTGTAAATTTGATTAATGCGAGTGAGTACTTATTTaaa        <  8:47697/148‑1 (MQ=255)                                                                                                                                                 gcCGGTTACGTCATTTTC‑TGCGTTTTCGCCTTTCCGAAATCTGAGTGCTCCCGCCTTTAGAGAGTCCTCGGGAACGACTTCCTGCACCCAAAAATTAGGATCGCAAGAGTGTAAATATGATAAACGCGAGTGAGTACTTTATTAAAAAt     >  4:232572/1‑149 (MQ=255)                                                                                                                                                  ccGGTTACGTCATTTTC‑TGCGTTTTCGCCTTTCCGAAATCTGAGTGCTCCCGCCTTTAGAGAGTCCTCGGGAACGACTTCCTGCACCCAAAAATTAGGATCGCAAGAGTGTAAATATGATAAACGCGAGTGAGTACTTTATTaaaaa      <  6:18515/147‑1 (MQ=255)                                                                                                                                                    ggTTACGTCATTTTC‑TGCGTTTTCGCCTTTCCGAAATCTGAGTGCTCCCGCCTTTAGAGAGTCCTCGGGAACGACTTCCTGCACCCAAAAATTAGGATCGCAAGAGTGTAAATATGATAAACGCGAGTGAGTACTTTATTAAAAATTAt  <  5:151621/149‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                                      GATCAAAAGCTTCCTGCGTAAGGCCTCGCTGAACACGCAGTTGGCGGAATACGCGGCCAAAGGCAATTTTGTAATTCATCACGCGTACGCGTCCAAAAGGAAGCATTCAACCTATTGCGTCGAGGCACGAACAGGCAGTATGCTGCCGGTTACGTCATTTTC‑TGCGTTTTCGCCTTTCCGAAATCTGAGTGCTCCCGCCTTTAGAGAGTCCTCGGGAACGACTTCCTGCACCCAAAAATTAGGATCGCAAGAGTGTAAATATGATAAACGCGAGTGAGTACTTTATTAAAAATTAT  >  NC_002947/2181930‑2182225

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