Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,625,580	T→A	5.6%	intergenic (‑50/+27)	PP_3195 ← / ← PP_3196	hypothetical protein/nodulation protein d1

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,625,580	0	T	A	5.6%	218.4 / 4.0	72	intergenic (‑50/+27)	PP_3195/PP_3196	hypothetical protein/nodulation protein d1
Reads supporting (aligned to +/- strand):  ref base T (37/30);  new base A (2/2);  total (39/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.85e-01

TTGTCCTTCTGCGGTCGTGGAGTCCGGCGCCCACTTGAGCGACAGGCAGGCAGGGCGCGGTTTGAAGGCAAGGTACGCAGGCAGGGCAGGGTGGGGCCAATCGTGATCAGGGATGCGAGGTATCGCGAAACGCGATGCCCCAG‑CACCGCTGCGGTGTCTGGGGCGATGGCCGCCGCTCAGGCGGTGGCGCGCAGAGCATCTGGCTGCGTGCTGGCGCACTGAATGACATCGCTGATGCGCTCGCGAAGCCATCTGTGCATGGGGTCGCCGTCCATGCTGCGGTGCCACAGCATGAA  >  NC_002947/3625431‑3625726                                                                                                                                                       |                                                                                                                                                   ttgtCCTTCTGCGGTCGTGGAGTCCGGCGCCCACTTGAGCGACAGGCAGGCAGGGCGCGGTTTGAAGGCAAGGTACGCAGGCAGGGCAGGGTGGGGCCAATCGTGATCAGGGATGCGAGGTATCGCGAAACGCGATGCCCCAG‑CACCGCt                                                                                                                                                    >  1:1795908/1‑150 (MQ=255)      cTTCTGCGGTCGTGGAGTCCGGCGCCCACTTGAGCGACAGGCAGGCAGGGCGCGGTTTGAAGGCAAGGTACGCAGGCAGGGCAGGGTGGGGCCAATCGTGATCAGGGATGCGAGGTATCGCGAAACGCGATGCCCCAG‑CACCGCTGCGgt                                                                                                                                               >  1:846660/1‑150 (MQ=255)         cTGCGGTCGTGGAGGCCTGCGCCGACATGCGGGACAGGCAGGCAGTGGGCGGTTTGAAGGCAAGGTACGCAGACAGGGCCGGGTGGGGCCAATCGTGATCAGGTATGCGAGGTCTCGCGAAACGGGATGACCGAG‑CATCGCTACGGCGTc                                                                                                                                            >  2:114780/1‑150 (MQ=255)          tGCGGTCGTGGAGTCCGGCGCCCACTTGAGCGACAGGCAGGCAGGGCGCGGTTTGAAGGCAAGGTACGCAGGCAGGGCAGGGTGGGGCCAATCGTGATCAGGGATGCGAGGTATCGCGAAACGCGATGCCCCAG‑CACCGCTGCGGTGTCt                                                                                                                                           >  2:1901086/1‑150 (MQ=255)           gCGGGCGTGGAGTCCGGCGCCCACTTGAGCGACAGGCAGGCAGGGCGCGGTTTCAAGGCAAGGTACGCAGGCAGGGCAGGGTGGGGCCAATCGTGATCAGGGATGCGAGCTATCGCGAAACGCGATGCCCCAG‑CACCGCTGCGGTGTCTg                                                                                                                                          <  1:1854358/150‑1 (MQ=255)              gTCGTGGAGTCCGGCGCCCACTTGAGCGACAGGCAGGCAGGGCGCGGTTTGAAGGCAAGGTACGCAGGCAGGGCAGGGTGGGGCCAACCGGTTTTTGGGATGCGAGGTATCGCGAAACGCGATGCCCCAG‑CACCGCTGCGGTGTCTgggg                                                                                                                                       <  2:264673/150‑1 (MQ=255)               tCGTGGAGTCCGGCGCCCACTTGAGCGACAGGCAGGCAGGGCGCGGTTTGAAGGCAAGGTACGCAGGCAGGGCAGGGTGGGGCCAATCGTGATCAGGGATGCGAGGTATCGCGAAACGCGATGCCCCAG‑CACCGCTGCGGTGTCTGGGGc                                                                                                                                      >  2:1513120/1‑150 (MQ=255)                      gTCCGGCGCCCACTTGAGCGACAGGCAGGCAGGGCGCGGTTTGAAGGCAAGGTACGCAGGCAGGGCAGGGTGGGGCCAATCGTGATCAGGGATGCGAGGTATCGCGAAACGCGATGCCCCAG‑CACCGCTGCGGTGTCTGGGGCGATGgcc                                                                                                                               <  2:1018812/150‑1 (MQ=255)                      gTCCGGCGCCCACTTGAGCGACAGGCAGGCAGGGCGCGGTTTGAAGGCAAGGTACGCAGGCAGGGCAGGGTGGGGCCAATCGTGATCAGGGATGCGAGGTATCGCGAAACGCGATGCCCCAG‑CACCGCTGCGGTGTCCGGGGCGATGgcc                                                                                                                               >  2:1739106/1‑150 (MQ=255)                       tCCGGCGCCCACTTGAGCGACAGGCAGGCAGGGCGCGGTTTGAAGGCAAGGTACTCAGGCAGGGCAGGGTGGGGCCAATCGTGATCAGGGATGCGAGGTATCGCGAAACGCGATGCCCCAG‑CACCGCTGCGGTGTCTGGGGCGATGgccg                                                                                                                              >  1:732092/1‑150 (MQ=255)                        ccGGCGCCCACTTGAGCGAAAGGCAGGCAGGGCGCGGTTTGAAGGCAAGGTACGCAGGCAGGGCAGGGTGGGGCCAATCGTGATCAGGGATGCGAGGTATCGCGAAACGCGATGCCCCAG‑CACCGCTGCGGTGTCTGGGGCGATGgccgc                                                                                                                             <  1:836490/150‑1 (MQ=255)                                       gCGACAGGCAGGCAGGGCGCGGTTTGAAGGCAAGGTACGCAGGCAGGGCAGGGTGGGGCCAATCGTGATCAGGGATGCGAGGTATCGCGAAACGCGATGCCCCAG‑CACCGCTGCGGTGTCTGGGGCGATGGCCGCCGCTCAGGCGGTGgc                                                                                                              >  1:924048/1‑150 (MQ=255)                                            aggcaggcaggGCGCGGTTTGAAGGCAAGGTACGCAGGCAGGGCAGGGTGGGGCCAATCGTGATCAGGGATGCGAGGTATCGCGAAACGCGATGCCCCAG‑CACCGGTGCGGTGTGTGGGGCGATGGCCGCCGCTCAGGCGGTGGCGCGCa                                                                                                         >  1:1317081/1‑150 (MQ=255)                                            aggcaggcaggGCGCGGTTTGAAGGCAAGGTACGCAGGCAGGGCAGGGTGGGGCCAATCGTGATCAGGGATGCGAGGTATCGCGAAACGCGATGCCCCAG‑CACCGCTGCGGTGTCTGGGGCGATGGCCGCCGCTCAGGCGGTGGCGCGCa                                                                                                         <  1:1892330/150‑1 (MQ=255)                                             ggcaggcaggGCGCGGTTTGAAGGCAAGGTACGCAGGCAGGGCAGGGTGGGGCCAATCGTGATTAGGGATGCGAGGTATCGCGAAACGCGATGCCCCAG‑CACCGCTGc                                                                                                                                                  <  2:1557197/108‑1 (MQ=255)                                             ggcaggcaggGCGCGGTTTGAAGGCAAGGTACGCAGGCAGGGCAGGGTGGGGCCAATCGTGATTAGGGATGCGAGGTATCGCGAAACGCGATGCCCCAG‑CACCGCTGc                                                                                                                                                  >  1:1557197/1‑108 (MQ=255)                                             ggcaggcaggGCGCGGTTTGAAGGCAAGGTACGCAGGCAGGGCAGGGTGGGGCCAATCGTGATCCGGGATGCGAGGTATCGCGAAACGCGATGCCCCAG‑CACCGCTGCGGTGTCTGGGGCGATGGCCGCCGCTCAGGCGGTGGCGCGCag                                                                                                        >  2:1030374/1‑150 (MQ=255)                                                aggcaggGCGCGGTTTGAAGGCAAGGTACGCAGGCAGGGCAGGGTGGGGCCAATCGTGATCAGGGATGCGAGGTATCGCGAAACGCGACGCCCCAGACACCGCAGCg                                                                                                                                                 <  1:998443/107‑1 (MQ=255)                                                aggcaggGCGCGGTTTGAAGGCAAGGTACGCAGGCAGGGCAGGGTGGGGCCAATCGTGATCAGGGATGCGAGGTATCGCGAAACGCGACGCCCCAGACACCGCAGCg                                                                                                                                                 >  2:998443/1‑107 (MQ=255)                                                 ggcaggGCGCGGTTTGAAGGCAAGGTACGCAGGCAGGGCAGGGTGGGGCCAATCGTGATCAGGGATGCGAGGTATCGCGAAACGCGATGCCCCAG‑CACCGCTGCGGTGTCTGGGGCGATGGCCGCCGCTCAGGCGGTGGCGCGCAGAGCa                                                                                                    <  1:922829/150‑1 (MQ=255)                                                  gcaggGCGCGGTTTGAAGGCAAGGTACGCAGGCAGGGCAGGGTGGGGCCAATCGTGATCAGGGATGCGAGGTATCGCGAAACGCGATGCCCCAG‑CACCGCTGCGGTGTCTGGGGCGATGGCCGCCGCTCAGGCGGTGGCGCGCAGAGCAt                                                                                                   >  1:855020/1‑150 (MQ=255)                                                  gcaggGCGCGGTTTGAAGGCAAGGTACGCAGGCAGGGCAGGGTGGGGCCAATCGTGATCAGGGATGCGAGGTATCGCGAAACGCGATGCCCCAG‑CACCGCTGCGGTGTCTGGGGCGATGGCCGCCGCTCAGGCGGTGGCGCGCAGAGCAt                                                                                                   >  1:700036/1‑150 (MQ=255)                                                     gggCGCGGTCTGAAGGCAAGGTACGCAGGCAGGGCAGGGTGGGGCCAATCGTGATCAGGGATGCGAGGTATCGCGAAACGCGATGCCCCAG‑CACCGCTGCGGTGTCTGGGGCGATGGCCGCCGCTCAGGCGGTGGCGCGCAGAGCATCTg                                                                                                <  1:1444566/150‑1 (MQ=255)                                                            gTTTGAAGGCAAGGTACGCAGGCAGGGCAGGGTGGGGCCAATCGTGATCAGGGATGCGAGGTATCGCGAAACGCGATGCCCCAG‑CACCGCTGCGGTGTCTGGGGCGATGGCCGCCGCTCAGGCGGt                                                                                                                 >  1:1779273/1‑126 (MQ=255)                                                            gTTTGAAGGCAAGGTACGCAGGCAGGGCAGGGTGGGGCCAATCGTGATCAGGGATGCGAGGTATCGCGAAACGCGATGCCCCAG‑CACCGCTGCGGTGTCTGGGGCGATGGCCGCCGCTCAGGCGGt                                                                                                                 <  2:1779273/126‑1 (MQ=255)                                                             tttGAAGGCAAGGTACGCAGGCAGGGCAGGTTGGGGCCAATCGTGATCAGGGATGCGAGGTATCGCGAAACGCGATGCCCCAG‑CACCGCTGCGGTGTCTGGGGCGATGGCCGCCGCTCAGGCGGTGGCGCGCAGAGCATCTGGCTGCGTg                                                                                        >  2:395655/1‑150 (MQ=255)                                                             tttGAAGGCAAGGTACGCAGGCAGGGCAGGGTGGGGCCAATCGTGATCAGGGATGCGAGGTATCGCGAAACGCGATGCCCCAG‑CACCGCTGCGGTGTCTGGGGCGATGGCCGCCGCTCAGGCGGTGGCGCGCAGAGCATCTGGCTGCGTg                                                                                        >  2:363371/1‑150 (MQ=255)                                                                  gggCAAGGTGCGCGGGCAGGGGGGGGTGGGGGCGATCGGGAGCAGGGATGCGAGGTATCGGGAAACGCGGTGCGCCAG‑CACGGCGGCGGTGTGTGGGGGGATGGCCGCCGCTCAGGGGGTGGCGCGCAGAGCATCTGGCTGCGTGCTGgc                                                                                   <  1:1030374/149‑1 (MQ=255)                                                                   ggCAAGGTACGCAGGCAGGGCAGGGTGGGGCCAATCGTGATCAGGGATGCGAGGTATCGCGAAACGCGATGCCCCAG‑CACCGCTGCGGTGTCTGGGGCGATGGCCGCCGCTCAGGCGGTGGCGCGCAGAGCATCTGGCTGCGTGCTGgcg                                                                                  >  1:1075294/1‑150 (MQ=255)                                                                       aGGTACGCAGGCAGGGCAGGGTGGGGCCAATCGTGATCAGGGATGCGAGGTATCGCGAAACGCGATGCCCCAG‑CACCGCTGCGGTGTCTGGGGCGATGGCCGCCGCTCAGGCGGTGGCGCGCAGAGCATCTGGCTGCGTGCTGGCGCACt                                                                              <  1:1901086/150‑1 (MQ=255)                                                                        ggTACGCAGGCAGGGCAGGGTGGGGCCAATCGTGATCAGGGATGCGAGGTATCGCGAAACGCGATGCCCCAG‑CACCGCTGCGGTGTCTGGGGCGa                                                                                                                                    <  1:1763810/95‑1 (MQ=255)                                                                        ggTACGCAGGCAGGGCAGGGTGGGGCCAATCGTGATCAGGGATGCGAGGTATCGCGAAACGCGATGCCCCAG‑CACCGCTGCGGTGTCTGGGGCGa                                                                                                                                    >  2:1763810/1‑95 (MQ=255)                                                                             gcaggcagGGCAGGGTGGGGCCAATCGTGATCAGGGATGCGAGGTATCGCGAAACGCGATGCCCCAG‑CACCGCTGCGGTGACTGGGGCGATGGCCGCCGCTCAGGCGGTGGCGCGCAGAGCATCTGGCTGCGTGCTGGCGCACTGAATGa                                                                        >  2:1516709/1‑150 (MQ=255)                                                                                   agggcagggTGGGGCCAATCGTGATCAGGGATGCGAGGTATCGCGAAACGCGATGCCCCAG‑CACCGCTGCGGTGTCTGGGGCGATGGCCGCCGCTCAGGCGGTGTCGCGCAGAGCATCTGGCTGCGTGCTGGCGCACTGAATGACATCGc                                                                  >  1:1667951/1‑150 (MQ=255)                                                                                     ggcagggTGGGGCCAATCGTGATCAGGGATGCGAGGTATCGCGAAACGCGATGCCCCAG‑CACCGCTGCGGTGTCTGGGGCGATGGCCGCCGCTCAGGCGGTGGCGCGCAGAGCATCTGGCTGCGTGCTGGCGCACTGAATGACATCGCTg                                                                >  1:1424885/1‑150 (MQ=255)                                                                                       cagggTGGGGCCAATCGTGATCAGGGATGCACTGGAGTTCGAAACTCTATGCCCCAG‑CACCGCTGCGGTGTCTGGGGCGATGGCCGCCGCTCAGGCGGTGGCGCGCAGAGCATCTGGCTGCGTGCTGGCGCACTGAATGACATCGCTGAt                                                              <  1:1397903/150‑1 (MQ=255)                                                                                        agggTGGGGCCAATCGTGATCAGGGATGCGAGGTATCGCGAAACGCGATGCCCCAG‑CACCGCTGCGGTGTCTGGGGCGATGGccgccg                                                                                                                           >  2:1542847/1‑88 (MQ=255)                                                                                        agggTGGGGCCAATCGTGATCAGGGATGCGAGGTATCGCGAAACGCGATGCCCCAG‑CACCGCTGCGGTGTCTGGGGCGATGGccgccg                                                                                                                           <  1:1542847/88‑1 (MQ=255)                                                                                         gggTGGGGCCAATCGTGATCAGGGATGCGAGGTATCGCGAAACGCGATGCCCCAG‑CACCGCTGCGGTGTCTGGGGCGATGGCCGCCGCTCAGGCGGTGGCGCGCAGAGCATCTGGCTGCGTGCTGGCGCACTGAATGACATCGCTGATgc                                                            <  2:1795908/150‑1 (MQ=255)                                                                                                  cAATCGTGATCAGGGATGCGAGGTATCGCGAAACGCGATGCCCCAG‑CACCGCTGCGGTGTCTGGGGCGATGGCCGCCGCTCAGGCGGTGGCGCGCAGAGCATCTGGCTGCGTGCTGGCGCACTGAATGACATCACTGATGCGCTCGCGaa                                                   >  1:1830781/1‑150 (MQ=255)                                                                                                     tCGTGATCAGGGATGCGAGGTATCGCGAAACGCGATGCCCCAG‑CACCGCTGCGGTGTCTGGGGCGATGGCCGCCGCTCAGGCGGTGGCGCGCAGAGCATCTGGCTGCGTGCTGGCGCACTGAATGACATCGCTGATGCGCTCGCGAAGcc                                                <  2:855020/150‑1 (MQ=255)                                                                                                       gTGATCAGGGATGCGAGGTATCACGAAACGCAATGCCCCAG‑CACCTCTTCGGTGTCTGGGGCGATGGCCGCCGCTCAGGCGGTGGCGCGCAGAGCATCTGGCTGCGTGCTGGCGCACTGAATGACATCGCTGATGCGCTCGCGAAGCCAt                                              <  2:1667951/150‑1 (MQ=255)                                                                                                         gATCAGGGATGCGAGGTATCGCGAAACGCGATGCCCCAG‑CACCGCTGCGGTGTCTGGGGCGATGGCCGCCGCTCAGGCGGTGGCGCGCAGAGCATCTGGCTGCGTGCTGGCGCACTGAATGACATCGCTGATGCGCTCGCGAAGCCATCt                                            <  2:1572001/150‑1 (MQ=255)                                                                                                           tCAGGGATGCGAGGTATCGCGAAACGCGATGCCCCAG‑CACCGCTGCGGTGTCTGGGGCGATGGCCGCCGCTCAGGCGGTGGCGCGCAGAGCATCTGGCTGCGTGCTGGCGCACTGAATGACATCGCTGATGCGCTCGCGAAGCCATCtgt                                          <  1:1513120/150‑1 (MQ=255)                                                                                                            cAGGGATGCGAGGTATCGCGAAACGCGATGCCCCAG‑CACCGCTGCGGTGTCTGGGGCGATGGCCGCCGCTCAGGCGGTGGCGCGCAGAGCATCTGGCTGCGTGCTGGCGCACTGAATGACATCGCTGATGCGCTCGCGAAGCCATCtgtg                                         <  1:1739106/150‑1 (MQ=255)                                                                                                             aGGGATGCGAGGTATCGCGAAACGCGATGCCCCAG‑CACCGCTGCGGTGTCTGGGGCGATGGCCGCCGCTCAGGCGGTGGCGCGCAGAGCATCTGGCTGCGTGCTGGCGCACTGAATGACATCGCTGATGCGCTCGCGAAGCCATCTGTGc                                        >  1:974380/1‑150 (MQ=255)                                                                                                               ggATGCGAGGTATCGCGAAACGCGATGCCCCAG‑CACCGCTGCGGTGTCTGGGGCGATGGCCGCCGCTCAGGCGGTGGCGCGCAGAGCATCTGGCTGCGTGCTGGCGCACTGAATGACATCGCTGATGCGCTCGCGAAGCCATCTGTGCAt                                      <  2:732092/150‑1 (MQ=255)                                                                                                                  tGCGAGGTATCGCGAAACGCGATGCCCCAG‑CACCGCTGCGGTGTCTGGGGCGATGGCCGCCGCTCAGGCGGTGGCGCGCAGAGCATCTGGc                                                                                              >  1:1196242/1‑91 (MQ=255)                                                                                                                  tGCGAGGTATCGCGAAACGCGATGCCCCAG‑CACCGCTGCGGTGTCTGGGGCGATGGCCGCCGCTCAGGCGGTGGCGCGCAGAGCATCTGGc                                                                                              <  2:1196242/91‑1 (MQ=255)                                                                                                                   tcGAGGTATCGCGAAACGCGATGCCCCAG‑CACCGCTGCGGTGTCTGGGGCGATGGCCGCCGCTCAGGCGGTGGCGCGCAGAGCATCTGGCTGCGTGCTGGCGCACTGAATGACATCGCTGATGCGCTCGCGAAGCCATCTGTGCATgggg                                  >  2:518656/2‑150 (MQ=255)                                                                                                                   gCGAGGTATCGCGAAACGCGATGCCCCAG‑CACCGCTGCGGTGTCTGGGGCGATGGCCGCCGCTCAGGCGGTGGCGCGCAGAGCATCTGGCTGCGTGCTGGCGCACTGAATGACATCGCTGATGCGCTCGCGAAGCCATCTGTGCATgggg                                  <  2:974380/150‑1 (MQ=255)                                                                                                                           tCGCGAAACGCGATGCCCCAG‑CACCGCTGCGGTGTCTGGGGCGATGGCCGCCGCTCAGGCGGTGGCGCGCAGAGCATCTGGCTGCGTGCTGGCGCACTGAATGACATCGCTGATGCGCTCGCGAAGCCATCTGTGCa                                       <  1:1072033/137‑1 (MQ=255)                                                                                                                           tCGCGAAACGCGATGCCCCAG‑CACCGCTGCGGTGTCTGGGGCGATGGCCGCCGCTCAGGCGGTGGCGCGCAGAGCATCTGGCTGCGTGCTGGCGCACTGAATGACATCGCTGATGCGCTCGCGAAGCCATCTGTGCa                                       >  2:1072033/1‑137 (MQ=255)                                                                                                                           tCGCGAAACGCGATGCCCCAG‑CACCGCTGCGGTGTCTGGGGCGATGGCCGCCGCTCAGGCGGTGGCGCGCAGAGCATCTGGCTGCGTGCTGGCGCACTGAATGACATCGCTGATGCGCTCGCGAAGCCATCTGTGCATGGGGTCGCCGTc                          >  1:471234/1‑150 (MQ=255)                                                                                                                               gAAACGCGATGCCCCAG‑CACCGCTGCGGTGTCTGGGGCGATGGCCGCCGCTCAGGCGGTGGCGCGCAGAGCATCTGGCTGCGTGCTGGCGCACTGAATGACATCGCTGATGCGCTCGCGAAGCCATCTGTGCATGGGGTCGCCGTCCAtt                      >  2:1093892/1‑149 (MQ=255)                                                                                                                                aaaCGCGATGCCCCAG‑CACCGCTGCGGTGTCTGGGGCGATGGCCGCCGCTCAGGCGGTGGCGCGCAGAGCATCTGGCTGCGTGCTGGCGCACTGAATGACATCGCTGATGCGCTCGCGAAGCCATCTGTGCATGGGGTCGCCGTCCAtgc                     >  2:570354/1‑150 (MQ=255)                                                                                                                                       aTGCCCCAG‑CACCGCTGCGGTGTCTGGGGCGATGGCCGCCGCTCAGGCGGTGGCGCGCAGAGCATCTGGCTGCGTGCTGGCGCACTGAATGACATCGCTGATGCGCTCGCGAAGCCATCTGTGCATGGGGTCGCCGTCCATGCTGCGGTg              >  2:1246323/1‑150 (MQ=255)                                                                                                                                        tGCCCCAG‑CACCGCTGCGGTGTCTGGGGAGATGGCCGCCGCTCAGGCGGTGGCGCGCAGAGCATCTGGCTGCGTGCTGGCGCACTGAATGACATCGCTGATGCGCTCGCGAAGCCATCTGTGCATGGGGTCGCCGTCCATGCTGCGGTGc             <  1:1426390/150‑1 (MQ=255)                                                                                                                                         gCCCCAGACACCGCAGCGGTG‑CTGGGGCGATGGCCGCCGCTCAGGCGGTGGCGCGCAGAGCATCTGGCTGCGTGCTGGCGCACTGAATGACATCGCTGATGCGCTCGCGAAGCCATCTGTGCATGGGGTCGCCGTCCATGCTGCGGTGcc            >  1:191310/1‑150 (MQ=255)                                                                                                                                           cccAG‑CACCGCTGCGGTGTCTGGGGCGATGGCCGCCGCTCAGGCGGTGGCGCGCAGAGCAt                                                                                                   <  2:1797566/61‑1 (MQ=255)                                                                                                                                           cccAG‑CACCGCTGCGGTGTCTGGGGCGATGGCCGCCGCTCAGGCGGTGGCGCGCAGAGCAt                                                                                                   >  1:1797566/1‑61 (MQ=255)                                                                                                                                            ccAG‑CACCGCTGCGGTGTCTGGGGCGATGGCCGCCGCTCAGGCGGTGGCGCGCAGAGCATCTGGCTGCGTGCTGGCGCACTGAATGACATCGCTGATGCGCTCGCGAAGCCATCTGTGCATGGGGTCGCCGTCCATGCTGCGGt               >  1:71967/1‑144 (MQ=255)                                                                                                                                            ccAG‑CACCGCTGCGGTGTCTGGGGCGATGGCCGCCGCTCAGGCGGTGGCGCGCAGAGCATCTGGCTGCGTGCTGGCGCACTGAATGACATCGCTGATGCGCTCGCGAAGCCATCTGTGCATGGGGTCGCCGTCCATGCTGCGGt               <  2:71967/144‑1 (MQ=255)                                                                                                                                            ccAG‑CACCGCTGCGGTGTCTGGGGCGATGGCCGCCGCTCAGGCGGTGGCGCGCAGAGCATCTGGCTGCGTGCTGGCGCACTGAATGACATCGCTGATGCGCTCGCGAAGCCATCTGTGCATGGGGTCGCCGTCCATGCTGCGGTGCcaca         >  1:1119013/1‑150 (MQ=255)                                                                                                                                            ccAG‑CACCGCTGCGGTGTCTGGGGCGATGGCCGCCGCTCAGGCGGTGGCGCGCAGAGCATCTGGCTGCGTGCTGGCGCACTGAATGACATCGCTGATGCGCTCGCGAAGCCATCTGTGCATGGGGTCGCCGTCCATGCTGCGGTGCcaca         >  1:482791/1‑150 (MQ=255)                                                                                                                                               g‑CACCGCTGCGGTGTCTGGGGCGATGGCCGCCGCTCAGGCGGTGGCGCGCAGAGCATCTGGCTGCGTGCTGGCGCACTGAATGACATCGCTGATGCGCTCGCGAAGCCATCTGTGCATGGGGTCGCCGTCCATGCTGCGGTGCCACAGCa      >  1:557399/1‑150 (MQ=255)                                                                                                                                                 cACCGCTGCGGTGTCTGGGGCGATGGCCGCCGCTCAGGCGGTGGCGCGCAGAGCATCTGGCTGCGTGCTGGCGCACTGAATGACATCACTGATGCGCTCGCGAAGCCATCTGTGCATGGGGTCGCCGTCCATGCTGCGGTGCCACAGCAt     <  2:1830781/150‑1 (MQ=255)                                                                                                                                                  aCCGCTGCGGTGTCTGGGGCGATGGCCGCCGCTCAGGCGGTGGCGCGCAGAGCATCTGGCTGCGTGCTGGCGCACTGAATGACATCGc                                                                  >  2:945287/1‑88 (MQ=255)                                                                                                                                                  aCCGCTGCGGTGTCTGGGGCGATGGCCGCCGCTCAGGCGGTGGCGCGCAGAGCATCTGGCTGCGTGCTGGCGCACTGAATGACATCGc                                                                  <  1:945287/88‑1 (MQ=255)                                                                                                                                                  aCCGCTGCGGTGTCTGGGGCGATGGCCGCCGCTCAGGCGGTGGCGCGCAGAGCATCTGGCTGCGTGCTGGCGCACTGAATGACATCGCTGATGCGCTCGCGAAGc                                                 <  2:665249/105‑1 (MQ=255)                                                                                                                                                  aCCGCTGCGGTGTCTGGGGCGATGGCCGCCGCTCAGGCGGTGGCGCGCAGAGCATCTGGCTGCGTGCTGGCGCACTGAATGACATCGCTGATGCGCTCGCGAAGc                                                 >  1:665249/1‑105 (MQ=255)                                                                                                                                                   ccGCTGCGGTGTCTGGGGCGATGGCCGCCGCTCAGGCGGTGGCGCGCAGAGCATCTGGCTGCGTGCTGGCGCACTGAATGACATCGCTGATGCGCTCGCGAAGCCATCTGTGCATGGGGTCGCCGTCCATGCTGCGGTGCCACAGCATGa   <  2:924048/150‑1 (MQ=255)                                                                                                                                                   ccGCAGCGGTG‑CTGGGGCGATGGCCGCCGCTCAGGCGGTGGCGCGCAGAGCATCTGGCTGCGTGCTGGCGCACTGAATGACATCGCTGATGCGCTCGCGAAGCCATCTGTGCATGGGGTCGCCGTCCATGCTGCGGTGCCACAGCATGaa  <  2:191310/150‑1 (MQ=255)                                                                                                                                                       |                                                                                                                                                   TTGTCCTTCTGCGGTCGTGGAGTCCGGCGCCCACTTGAGCGACAGGCAGGCAGGGCGCGGTTTGAAGGCAAGGTACGCAGGCAGGGCAGGGTGGGGCCAATCGTGATCAGGGATGCGAGGTATCGCGAAACGCGATGCCCCAG‑CACCGCTGCGGTGTCTGGGGCGATGGCCGCCGCTCAGGCGGTGGCGCGCAGAGCATCTGGCTGCGTGCTGGCGCACTGAATGACATCGCTGATGCGCTCGCGAAGCCATCTGTGCATGGGGTCGCCGTCCATGCTGCGGTGCCACAGCATGAA  >  NC_002947/3625431‑3625726

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