Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,146,825	C→A	6.5%	S413I (AGC→ATC)	PP_1006 ←	heme receptor

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,146,825	0	C	A	6.5%	244.4 / 6.0	77	S413I (AGC→ATC)	PP_1006	heme receptor
Reads supporting (aligned to +/- strand):  ref base C (32/40);  new base A (2/3);  total (34/43)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.24e-01

GAATACCCACTGGCGATCCTTGTAGGTAGTCTGCCGGTCCCGAAACACTTGGCGTCCCGAGGCGACATACAGCTCATCGGTAC‑GTTGGTCGGTCTTCGCGATCTGGTAGTTGAAACTCCATTTGACGTGATCCGCGACCACGCTATCGAGGTCGAATTCATGATTGAGGCCGAAGCGTTCACGGGTGATCGTGTCGTTGCCTTGGCGCATGCGGTACGAGTTCATGCCACCGAAGCCCGGAATGAACGGCCCACCCACCGCGCTGAGGATGTTCTGGTCACGGTCATCC  >  NC_002947/1146683‑1146971                                                                                                                                                |                                                                                                                                                   gAATACCCACTGGCGATCCTTGTAGGTAGTCTGCCGGTCCCGAAACACTTGGCGTCCCGAGGCGACATACAGCTCATCGGTAC‑GTTGGTCGGTCTTCGCGATCTGGTAGTTGAAACTCCATTTGACGTGATCCGCGACCACGCTATc                                                                                                                                                >  5:178176/1‑147 (MQ=255)  aagaCCCACTGGCGATCCTTGTAGGTAGTCTGCCGGTCCCGAAACACTTGGCGTCCCGAGGCGACATACAGCTCATCGGTAC‑GTTGGTCGGTCTTCGCGATCTGGTAGTTGAAACTCCATTTGACGTGATCCGCGACCACGCTATCGAg                                                                                                                                             <  5:188851/146‑1 (MQ=255)   aTACCCACTGGCGATCCTTGTAGGTAGTCTGCCGGTCCCGAAACACTTGGCGTCCCGAGGCGACATACAGCTCATCGGTAC‑GTTGGTCGGTCTTCGCGATCTGGTAGTTGAAACTCCATTTGACGTGATCCGCGACCACGCTATCGAgg                                                                                                                                            >  2:43405/1‑149 (MQ=255)    tACCCACTGGCGATCCTTGTAGGTAGTCTGCCGGTCCCGAAACACTTGGCGTCCCGAGGCGACATACAGCTCATCGGTAC‑GTTGGTCGGTCTTCGCGATCTGGTAGTTGAAACTCCATTTGACGTGATCCGCGACCACGCTATCGAGGt                                                                                                                                           >  7:70098/1‑149 (MQ=255)        cACTGGCTATCCTTGTAGGTAGTCTGCCGGTCCCGAAACACTTGGCGTCCCGAGGCGACATACAGCTCATCGGTAC‑GTTGGTCGGTCTTCGCGATCTGGTAGTTGAAACTCCATTTGACGTGATCCGCGACCTCGATAGCGTGGTCgcg                                                                                                                                       <  1:38588/149‑3 (MQ=255)         aCTGGCGATCCTTGTAGGTAGTCTGCCGGTCCCGAAACACTTGGCGTCCCGAGGCGACATACAGCTCATCGGTAC‑GTTGGTCGGTCTTCGCGATCTGGTAGTTGAAACTCCATTTGACGTGATCCGCGACCACGCTATCGAGGTCGAAt                                                                                                                                      <  6:84866/149‑1 (MQ=255)           tGGCGATCCTTGTAGGTAGTCTGCCGGTCCCGAAACACTTGGCGTCCCGAGGCGACATACAGCTCATCGGTAC‑GTTGGTCGGTCTTCGCGATCTGGTAGTTGAAACTCCATTTGACGTGATCCGCGACCACGCTATCGAg                                                                                                                                             <  6:206812/140‑1 (MQ=255)           tGGCGATCCTTGTAGGTAGTCTGCCGGTCCCGAAACACTTGGCGTCCCGAGGCGACATACAGCTCATCGGTAC‑GTTGGTCGGTCTTCGCGATCTGGTAGTTGAAACTCCATTTGACGTGATCCGCGACCACGCTATCGAg                                                                                                                                             >  5:206812/1‑140 (MQ=255)               gATCCTTGTAGGTAGTCTGCCGTTCCCGAAACACTTGGCGTCCCGAGGCGACATACAGCTCATCGGTAC‑GTTGGTCGGTCTTCGCGATCTGGTAGTTGAAACTCCATTTGACGTGATCCGCGACCACGCTATCGAGGTCGAATTCATGa                                                                                                                                <  8:248260/149‑1 (MQ=255)               gATCCTTGTAGGTAGTCTGCCGGTCCCGAAACACTTGGCGTCCCGAGGCGACATACAGCTCATCGGTAC‑GTTGGTCGGTCTTCGCGATCTGGTAGTTGAAACTCCATTTGACGTGATCCGCGACCACGCTATCGAGGTCGAATTCATGa                                                                                                                                >  2:150339/1‑149 (MQ=255)               gATCCTTGTAGGTAGTCTGCCGGTCCCGAAACACTTGGCGTCCCGAGGCGACATACAGCTCATCGGTAC‑GTTGGTCGGTCTTCGCGATCTGGTAGTTGAAACTCCATTTGACGTGATCCGCGACCACGCTATCGAGGTCGAATTCATGa                                                                                                                                <  3:251017/149‑1 (MQ=255)                 tCCTTGTAGGTAGTCTGCCGGTCCCGAAACACTTTGCGTCCCGAGGCGACATACAGCTCATCGGTAC‑GTTGGTCGGTCTTCGCGATCTGGTAGTTGAAACTCCATTTGACGTGATCCGCGACCACGCTATCGAGGTCGAATTCATGAtt                                                                                                                              <  8:6838/149‑1 (MQ=255)                  ccTTGTAGGTAGTCTGCCGGTCCCGAAACACTTGGCGTCCCGAGGCGACATACAGCTCATCGGTAC‑GTTGGTCGGTCTTCGCGATCTGGTAGTTGAAACTCCATTTGACGTGATCCGCGACCACGCTATc                                                                                                                                                <  4:83098/130‑1 (MQ=255)                  ccTTGTAGGTAGTCTGCCGGTCCCGAAACACTTGGCGTCCCGAGGCGACATACAGCTCATCGGTAC‑GTTGGTCGGTCTTCGCGATCTGGTAGTTGAAACTCCATTTGACGTGATCCGCGACCACGCTATc                                                                                                                                                >  3:83098/1‑130 (MQ=255)                   cTTGTAGGTAGTCTGCCGGTCCCGAAACACTTGGCGTCCCGAGGCGACATACAGCTCATCGGTAC‑GTTGGTCGGTCTTCGCGATCTGGTAGTTGAAACTCCATTTGACGTGATCCGCGACCACGCTATCGAGGTCGAATTCATGATTGa                                                                                                                            <  3:143984/149‑1 (MQ=255)                   cTTGTAGGTAGTCTGCCGGTCCCGAAACACTTGGCGTCACGAGGCGACATACAGCTCATCGGTAC‑GTTGGTCGGTCTTCGCGATCTGGTAGTTGAAACTCCATTTGACGTGATCCGCGACCACGCTATCGAGGTCGAATTGATGAtt                                                                                                                              >  1:284114/1‑147 (MQ=255)                    ttGTAGGTAGTCTGCCGGTCCCGAAACACTTGGCGTCCCGAGGCGACATACAGCTCATCGGTAC‑GTTGGTCGGTCTTCGCGATCTGGTAGTTGAAACTCCATTTGACGTGATCCGTGACCACGCTATCGAGGTCGAATTCATGAtt                                                                                                                              >  6:36658/1‑146 (MQ=255)                    ttGTAGGTAGTCTGCCGGTCCCGAAACACTTGGCGTCCCGAGGCGACATACAGCTCATCGGTAC‑GTTGGTCGGTCTTCGCGATCTGGTAGTTGAAACTCCATTTGACGTGATCCGCGACCACGCTATCGAGGTCGAATTCATGATTGAg                                                                                                                           >  6:368441/1‑149 (MQ=255)                      gtaggtagTCTGCCGGTCCCGAAACACTTGGCGTCCCGAGGCGACATACAGCTCATCGGTAC‑GTTGGTCGGTCTTCGCGATCTGGTAGTTGAAACTCCATTTGACGTGATCCGTGACCACGCTATCGAGGTCGAATTCATGATTGa                                                                                                                            <  5:36658/146‑1 (MQ=255)                          gtagTCTGCCGGTCCCGAAACACTTGGCGTCCCGAGGCGACATACAGCTCATCGGTAC‑GTTGGTCGGTCTTCGCGATCTGGTAGTTGAAACTCCATTTGACGTGATCCGCGACCACGCTATCGAGGTCGAATTCATGATTGAGGCCGaa                                                                                                                     <  1:43405/149‑1 (MQ=255)                              tCTGCCGGTCCCGAAACACTTGGCGTCCCGAGGCGACATACAGCTCATCGGTAC‑GTTGGTCGGTCTTCGCGATCTGGTAGTTGAAACTCCATTTGACGTGATCCGCGACCACGCTATCGAGGTCGAATTCATGATTGAg                                                                                                                           >  5:274761/1‑139 (MQ=255)                              tCTGCCGGTCCCGAAACACTTGGCGTCCCGAGGCGACATACAGCTCATCGGTAC‑GTTGGTCGGTCTTCGCGATCTGGTAGTTGAAACTCCATTTGACGTGATCCGCGACCACGCTATCGAGGTCGAATTCATGATTGAg                                                                                                                           <  6:274761/139‑1 (MQ=255)                              tCTGCCGGTCCCGAAACACTTGGCGTCCCGAGGCGACATACAGCTCATCGGTAC‑GTTGGTCGGTCTTCGCGATCTGGTAGTTGAAACTCCATTTGACGTGATCCGCGACCACGCTATCGAGGTCGAATTCATGATTGAGGCCGAAGCGt                                                                                                                 >  8:322312/1‑149 (MQ=255)                                tGCCGGTCACGAAACACTTGGCGTCCCTAGGCGACATACAGCTCATCGGTAC‑GTTGGTCGGTCTTCGCGATCTGGTAGTTGAAACTCCATTTGACGTGATCCGCGACCACGCTATCGAGGTCGAATTCATGATTGAGGCCGAAGCGtt                                                                                                                >  7:106422/1‑148 (MQ=255)                                 gCCGGTCCCGAAACACTTGGCGTCCCGAGGCGACATACAGCTCATCGGTAC‑GTTGGTCGGTCTTCGCGATCTGGTAGTTGAAACTCCATTTGACGTGATCCGCGACCACGCTATCGAGGTCGAATTCATGATTGAGGCCGAAGCGt                                                                                                                 <  7:116831/146‑1 (MQ=255)                                 gCCGGTCCCGAAACACTTGGCGTCCCGAGGCGACATACAGCTAATCGGTAC‑GTTGGTCGGTCTTCGCGATCTGGTAGTTGAAACTCCATTTGACGTGATCCGCGACCACGCTATCGAGGTCGAATTCATGATTGAGGCCGAAGCGt                                                                                                                 >  8:116831/1‑146 (MQ=255)                                     gTCCCGAAACACTTGGCGTCCCGAGGCGACATACAGCTCATCGGTAC‑GTTGGTCGGTCTTCGCGATCTGGTAGTTGAAACTCCATTTGACGTGATCCGCGACCACGCTATCGAGGTCGAATTCATGATTGAGGCCGAAGCGTTCACgg                                                                                                           >  5:189480/1‑148 (MQ=255)                                     gTCCCGAAACACTTGGCGTCCCGAGGCGACATACAGCTCATCGGTAC‑GTTGGTCGGTCTTCGCGATCTGGTAGTTGAAACTCCATTTGACGTGATCCGCGACCACGCTATCGAGGTCGAATTCATGATTGAGGCCGAAGCGTTCACgg                                                                                                           <  6:189480/148‑1 (MQ=255)                                       cccGAACCACTTGGCGTCCCGAGGCGACATACAGCTCATCGGTAC‑GTTTGTCGGTCTTCGCGATCTGGTAGTTGAAACTCCATTTGACGTGATCCGCGACCACGCTATCGAGGTCGAATTCATGATTGAGGCCGAAGCGTTCACGGGTg                                                                                                        <  2:284114/149‑1 (MQ=255)                                       cccGAAACACTTGGCGTCCCGAGGCGACATACAGCTCATCGGTAC‑GTTGGTCGGTCTTCGCGATCTGGTAGTTGAAACTCCATTTGACGTGATCCGCGACCACGCTATCGAGGTCGAATTCATGATTGAGGCCGAAGCGTTCACGGGt                                                                                                         >  6:355795/1‑148 (MQ=255)                                          gAAACACTTGGCGTCCCGAGGCGACATACAGCTCATCGGTAC‑GTTGGTCGGTCTTCGCGATCTGGTAGTTGAAACTCCATTTGACGTGATCCGCGACCACGCTATCGAGGTCGAATTCATGATTGAGGCCGAAGCGTTCACGGGTGATc                                                                                                     <  6:178176/149‑1 (MQ=255)                                          gAAACACTTGGCGTCCCGAGGCGACATACAGCTCATCGGTAC‑GTTGGTCGGTCTTCGCGATCTGGTAGTTGAAACTCCATTTGACGTGATCCGCGACCACGCTATCGAGGTCGAATTCATGATTGAGGCCGAAGCGTTCACGGGTGATc                                                                                                     >  4:314137/1‑149 (MQ=255)                                          gAAACACTTGGCGTCCCGAGGCGACATACAGCTCATCGGTAC‑GCGGGGCGGTCTTCGCGATCTGGTAGTTGAAACTCCATTTGACGTGATCCGCGACCACGCTATCGAGGTCGAATTCATGATTGAGGCCGAAGCGTTCACGGGTGATc                                                                                                     <  1:302798/149‑1 (MQ=255)                                           aaaaACATGGCGTCACGAGGCTACAAACACCTCAAACGTAC‑CTTCGTCCGTCTACCCGATCAATTAGTTCAAACACCATTTGACGTGATCCGAGACCACGCTATCGAGGTCGAATTCAAGATTCACTCACAAGCTTTCACgggagctc                                                                                                     >  4:286669/1‑143 (MQ=255)                                            aaCACTTGGCGTCCCGAGGCGACATACAGCTCATCGGTAC‑GTTGGTCGGTCTTCGCGATCTGGTAGTTGAAACTCCATTTGACGTGATCCGCGACCACGCTATCGAGGTCGAATTCATGATTGAGGCCGAAGCGTTCACGGGTGATCgt                                                                                                   <  3:236961/149‑1 (MQ=255)                                             acacTTGGCGTCCCGAGGCGACATACAGCTCATCGGTAC‑GTTGGTCGGTCTTCGCGATCTGGTAGTTGAAACTCCATTTGACGTGATCCGCGACCACGCTATCGAGGTCGAATTCATGATTGAGGCCGAAGCGTTCACGGGTGATCgtg                                                                                                  <  5:368441/149‑1 (MQ=255)                                                 ttGGCGTCCCGAGGCGACATACAGCTCATCGGTAC‑GTTGGTCGGTCTTCGCGATCTGGTAGTTGAAACTCCATTTGACGTGATCCGCGACCACGCTATCGAGGTCGAATTCATGAtt                                                                                                                              <  4:47826/117‑1 (MQ=255)                                                 ttGGCGTCCCGAGGCGACATACAGCTCATCGGTAC‑GTTGGTCGGTCTTCGCGATCTGGTAGTTGAAACTCCATTTGACGTGATCCGCGACCACGCTATCGAGGTCGAATTCATGAtt                                                                                                                              >  3:47826/1‑117 (MQ=255)                                                    gCGTCCCGAGGCGACATACAGCTCATCGGTAC‑GTTGGTCGGTCTTCGCGATCTGGTAGTTGAAACTCCATTTGACGTGATCCGCGACCACGCTATCGAGGTCGAATTCATGATTGAGGCCGAAGCGTTCACGGGTGATCGTGTCGTTGc                                                                                           <  1:150339/149‑1 (MQ=255)                                                       tCCCGAGGCGACATACAGCTCATCGGTACGGTTGGTCGGTCTTCGCGATCTGGTAGTTGAAACTCCATTTGACGTGATCCGCGACCACGCTATCGAGGTCGAATTCATGATTGAGGCCGAAGCGTTCACGGGTGATCgtgt                                                                                                 >  1:215893/1‑141 (MQ=255)                                                       tCCCGAGGCGACATACAGCTCATCGGTAC‑GTTGGTCGGTCTTCGCGATCTGGTAGTTGAAACTCCATTTGACGTGATCCGCGACCACGCTATCGAGGTCGAATTCATGATTGAGGCCGAAGCGTTCACGGGTGATCgtgt                                                                                                 <  2:215893/140‑1 (MQ=255)                                                       tCCCGAGGCGACATACAGCTCATCGGTAC‑GTTGGTCGGTCTTCGCGATCTGGTAGTTGAAACTCCATTTGACGTGATCCGCGACCACGCTATCGAGGTCGAATTCATGATTGAGGCCGAAGCGTTCACGGGTGATCGTGTCGTTGCCtt                                                                                        <  2:231290/149‑1 (MQ=255)                                                         ccGAGGCGACATACAGCTCATCGGTAC‑GTTGGTCGGTCTTCGCGATCTGGTAGTTGAAACTCCATTTGACGTGATCCGCGACCACGCTATCGAGGTCGAATTCATGATTGAGGCCGAAGCGTTCACGGGTGATCGTGTCGTTGCCTTgg                                                                                      >  3:201245/1‑149 (MQ=255)                                                              gCGACATACAGCTCATCGGTAC‑GTTGGTCGGTCTTCGCGATCTGGTAGTTGAAACTCCATTTGACGTGATCCGCGACCTCGATAGCGTGGTCgcgg                                                                                                                                      >  7:234792/1‑93 (MQ=255)                                                              gCGACATACAGCTCATCGGTAC‑GTTGGTCGGTCTTCGCGATCTGGTAGTTGAAACTCCATTTGACGTGATCCGCGACCTCGATAGCGTGGTCgcgg                                                                                                                                      <  8:234792/96‑4 (MQ=255)                                                                    tACAGCTCATCGGTAA‑GTTGGTCGGTCTTCGCGATCTGGTAGTTGAAACTCCATTTGACGTGATCCGCGACCACGCTATCGAGGTCTAATTCATGATTGAGGCCGAAGCGTTCACGGGTGATCGTGTCGTTGCCTTGGCGCATGCGGTa                                                                           <  5:219921/149‑1 (MQ=255)                                                                      cAGCTCATCGGTAC‑GTTGGTCGGTCTTCGAGATCTGGTAGTTGAAACTCCATTTGACGTGATCCGCGACCACGCTATCGAGGTCGAATTCATGATTGAGGCCGAAGCGTTCACGGGTGATCGTGTCGTTGCCTTGGCGCATGCGGTAc                                                                          <  3:314137/148‑1 (MQ=255)                                                                         cTCATCGGTAC‑GTTGGTCGGTCTTCGCGATCTGGTAGTTGAAACTCCATTTGACGTGATCCGCGACCACGCTATCGAGGTCGAATTCATGATTGAGGCCGAAGCGTTCACGGGTGATCGTGTCGTTGCCTTg                                                                                       <  2:150035/132‑1 (MQ=255)                                                                         cTCATCGGTAC‑GTTGGTCGGTCTTCGCGATCTGGTAGTTGAAACTCCATTTGACGTGATCCGCGACCACGCTATCGAGGTCGAATTCATGATTGAGGCCGAAGCGTTCACGGGTGATCGTGTCGTTGCCTTg                                                                                       >  1:150035/1‑132 (MQ=255)                                                                            aTCGGTAC‑GTTGGTCGGTCTTCGCGATCTGGTAGTTGAAACTCCATTTGACGTGATCCGCGACCACGCTATCGAGGTCGAATTCATGATTGAGGCCGAAGCGTTCACGGGTGATCGTGTCGTTGCCTTGGCGCATGCGGTACGAGTTCa                                                                   <  7:322312/149‑1 (MQ=255)                                                                             tCGGTAC‑GTTGGTCGGTCTTCGCGATCTGGTAGTTGAAACTCCATTTGACGTGATCCGCGACCACGCTATCGAGGTCGAATTCATGATTGAGGCCGAAGCGTTCACGGGTGATCGTGTCGTTGCCTTGGCGCATGCGGTACGAGTTCAt                                                                  <  8:70098/149‑1 (MQ=255)                                                                              cGGTAC‑GTTGGTCGGTCTTCGCGATCTGGTAGTTGAAACTCCATTTGACGTGATCCGCGACCACGCTATCGAGGTCGAATTCATGATTGAGGCCGAAGCGTTCACGGGTGATCGTGTCGTTGCCTTGGCGCATGCGGTACGAGTTCAt                                                                  <  8:106422/148‑1 (MQ=255)                                                                               ggTAC‑GTTGGTCGGTCTTCGCGATCTGGTAGTTGAAACTCCATTTGACGTGATCCGCGACCACGCTATCGAGGTCGAATTCATGATTGAGGCCGAAGCGTTCACGGGTGATCGTGTCGTTGCCTTGGCGCATGCGGTACGAGTTCATGc                                                                <  5:355795/149‑1 (MQ=255)                                                                                gTAC‑GTTGGTCGGTCTTCGCGATCTGGTAGTTGAAACTCCATTTGACGTGATCCGCGACCTCGATAGCGTGGTCgc                                                                                                                                        <  6:11986/76‑2 (MQ=38)                                                                                gTAC‑GTTGGTCGGTCTTCGCGATCTGGTAGTTGAAACTCCATTTGACGTGATCCGCGACCTCGATAGCGTGGTCgc                                                                                                                                        >  5:11986/1‑75 (MQ=38)                                                                                         gtcggtcTTCGCGATCTGGTAGTTGAAACTCCATTTGACGTGATCCGCGACCACGCTATCGAGGTCGAATTCATGATTGAGGCCGAAGCGTTCACGGGTGATCGTGTCGTTGCCTTGGCGCATGCGGTACGAGTTCATGCCACCGAAGc                                                       >  1:83369/1‑149 (MQ=255)                                                                                          tcggtcTTCGCGATCTGGTAGTTGAAACTCCATTTGACGTGATCCGCGACCACGCTATCGAGGTCGAATTCATGATTGAGGCCGAAGCGTTCACGGGTGATCGTGTCGTTGCCTTGGCGCATGCGGTACGAGTTCATGCCACCGAAGcc                                                      >  4:64506/1‑149 (MQ=255)                                                                                            ggtcTTCGCGATCTGGTAGTTGAAACTCCATTTGACGTGATCCGCGACCACGCTATCGAGGTCGAATTCATGATTGAGGCCGAAGCGTTCACGGGTGATCGTGTCGTTGCCTTGGCGCATGCGGTACGAGTTCATGCCACCGAAGCCCg                                                    <  2:83369/149‑1 (MQ=255)                                                                                              tcTTCGCGATCTGGTAGTTGAAACTCCATTTGACGTGATCCGCGACCACGCTATCGAGGTCGAATTCATGATTGAGGCCGAAGCGTTCACGGGTGATCGTGTCGTTGCCTTGGCGCATGCGGTACGAGTTCATGCCACCGAAGCCCgg                                                   >  5:215729/1‑148 (MQ=255)                                                                                                ttCGCGATCTGGTAGTTGAAACTCCATTTGACGTGATCCGCGACCACGCTATCGAGGTCGAATTCATGATTGAGGCCGAAGCGTTCACGGGTGATCGTGTCGTTGCCTTGGCGCATg                                                                                >  8:329557/1‑117 (MQ=255)                                                                                                ttCGCGATCTGGTAGTTGAAACTCCATTTGACGTGATCCGCGACCACGCTATCGAGGTCGAATTCATGATTGAGGCCGAAGCGTTCACGGGTGATCGTGTCGTTGCCTTGGCGCATg                                                                                <  7:329557/117‑1 (MQ=255)                                                                                                ttCGCGATCTGGTAGTTGAAACTCCATTTGACGTGATCCGCGACCACGCTATCGAGGTCGAATTCATGATTGAGGCCGAAGCGTTCACGGGTGATCGTGTCGTTGCCTTGGCGCATGCGGTACGAGTTCATGCCACCGAAGCCCGGa                                                  >  4:151799/1‑147 (MQ=255)                                                                                                ttCGCGATCTGGTAGTTGAAACTCCATTTGACGTGATCCGCGACCACGCTATCGAGGTCGAATTCATGATTGAGGCCGAAGCGTTCACGGGTGATCGTGTCGTTGCCTTGGCGCATGCGGTACGAGTTCATGCCACCGAAGCCCGGa                                                  <  3:151799/147‑1 (MQ=255)                                                                                                   gcgATCTGGTAGTTGAAACTCCATTTGACGTGATCCGCGACCACGCTATCGAGGTCGAATTCATGATTGAGGCCGAAGCGTTCACGGGTGATCGTGTCGTTGCCTTGGCGCATGCGGTACGAGTTCATGCCACCGAAGCCCGGAATGaa                                             <  2:180237/149‑1 (MQ=255)                                                                                                          ggTAGTTGAAACTCCATTTGACGTGATCCGCGACCACGCTATCGAGGTCGAATTCATGATTGAGGCCGAAGCGTTCACGGGTGATCGTGTCGTTGGCTTGGCGCATGCGGTACGAGTTCATGCCACCGAAGCCCGGAATGAACGGccca                                      >  2:317522/1‑149 (MQ=255)                                                                                                           gTAGTTGAAACTCCATTTGACGTGATCCGCGACCACGCTATCGAGGTCGAATTCATGATTGAGGCCGAAGCGTTCACGGGTGATCGTGTCGTTGCCTTGGCGCATGCGGTACGAGTTCATGCCACCGAAGCCCGGAATGAACGGccca                                      <  3:64506/148‑1 (MQ=255)                                                                                                                     cTCCATTTGACGTGATCCGCGACCACGCTATCGAGGTCGAATTCATGATTGAGGCCGAAGCGTTCACGGGTGATCGTGTCGTTGCCTTGGCGCATGCGGTACGAGTTCATGCCACCGAAGCCCGGAATGAACGGCCCACCCACCGCGCt                           >  5:155151/1‑149 (MQ=255)                                                                                                                              aCGTGATCCGCGACCACGCTATCGAGGTCGAATTCATGATTGAGGCCGAAGCGTTCACGGGTGATCGTGTCGTTGCCTTGGCGCATGCGGTACGAGTTCATGCCACCGAAGCCCGGAATGAACGGCCCACCCACCGCGCTGAGGATGtt                  <  3:286669/149‑1 (MQ=255)                                                                                                                                   aTCCGCGACCACGCTATCGAGGTCGAATTCATGATTGAGGCCGAAGCGTTCACGGGTGATCGTGTCGTTGCCTTGGCGCATGCGGTACGAGTTCATGCCACCGAAGCCCGGAATGAACGGCCCACCCACCGCGCTGAGGATGTTCTGGt             >  4:49756/1‑149 (MQ=255)                                                                                                                                     ccGCGACCACGCTATCGAGGTCGAATTCATGATTGAGGCCGAAGCGTTCACGGGTGATCGTGTCGTTGCCTTGGCGCATGCGGTACGAGTTCATGCCACCGAAGCCCGGAATGAACGGCCCACCCACCGCGCTGAGGATGTTCTGGTCa           >  3:359135/1‑149 (MQ=255)                                                                                                                                      gcggatcACGCTATCGAGGTCGAATTCATGATTGAGGCCGAAGCGTTCACGGGTGATCGTGTCGTTGCCTTGGCGCATGCGGTACGAGTTCATGCCACCGAAGCCCGGAATGAACGGCCCACCCACCGCGCTGAGGATGTTc                 <  5:299591/136‑1 (MQ=255)                                                                                                                                      gcggatcACGCTATCGAGGTCGAATTCATGATTGAGGCCGAAGCGTTCACGGGTGATCGTGTCGTTGCCTTGGCGCATGCGGTACGAGTTCATGCCACCGAAGCCCGGAATGAACGGCCCACCCACCGCGCTGAGGATGTTc                 >  6:299591/7‑142 (MQ=255)                                                                                                                                       gcgACCACGCTATCGAGGTCGAATTCATGATTGAGGCCGAAGCGTTCACGGGTGATCGTGTCGTTGCCTTGGCGCATGCGGTACGAGTTCATGCCACCGAAGCCCGGAATGAACGGCCCACCCACCGCGCTGAGGATGTTCTGGTCACg         >  2:149985/1‑149 (MQ=255)                                                                                                                                       gcgACCACGCTATCGAGGTCGAATTCATGATTGAGGCCGAAGCGTTCACGGGTGATCGTGTCGTTGCCTTGGCGCATGCGGTACGAGTTCATGCCACCGAAGCCCGGAATGAACGGCCCACCCACCGCGCTGAGGATGTTCTGGTCACg         <  4:359135/149‑1 (MQ=255)                                                                                                                                         gACCACGCTATCGAGGTCGAATTCATGATTGAGGCCGAAGCGTTCACGGGTGATCGTGTCGTTGGCTTGGCGCATGCGGTACGAGTTCATGCCACCGAAGCCCGGAATGAACGGCCCACCCACCGCGCTGAGGATGTTCTGGTCACGGt       <  1:317522/149‑1 (MQ=255)                                                                                                                                             aCGCTATCGAGGTCGAATTCATGATTGAGGCCGAAGCGTTCACGGGTGATCGTGTCGTTGCCTTGGCGCATGCGGTACGAGTTCATGCCACCGAAGCCCGGAATGAACGGCCCACCCACCGCGCTGAGGATGTTCTGGTCACGGTCATc   <  1:149985/149‑1 (MQ=255)                                                                                                                                              cGCTATCGAGGGCGAATTCATGATTGAGGCCGAAGCGTTCACGGGTGATCGTGTCGTTGCCTTGGCGCATGCGGTACGAGTTCATGCCACCGAAGCCCGGAATGAACGGCCCACCCACCGCGCTGAGGATGTTCTGGTCACGGTCATcc  <  4:201245/149‑1 (MQ=255)                                                                                                                                                |                                                                                                                                                   GAATACCCACTGGCGATCCTTGTAGGTAGTCTGCCGGTCCCGAAACACTTGGCGTCCCGAGGCGACATACAGCTCATCGGTAC‑GTTGGTCGGTCTTCGCGATCTGGTAGTTGAAACTCCATTTGACGTGATCCGCGACCACGCTATCGAGGTCGAATTCATGATTGAGGCCGAAGCGTTCACGGGTGATCGTGTCGTTGCCTTGGCGCATGCGGTACGAGTTCATGCCACCGAAGCCCGGAATGAACGGCCCACCCACCGCGCTGAGGATGTTCTGGTCACGGTCATCC  >  NC_002947/1146683‑1146971

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: ‑