Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,014,674	C→T	6.3%	I10I (ATC→ATT)	PP_3543 →	iron‑sulfur cluster‑binding protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,014,674	0	C	T	6.3%	189.5 / 2.1	63	I10I (ATC→ATT)	PP_3543	iron‑sulfur cluster‑binding protein
Reads supporting (aligned to +/- strand):  ref base C (30/29);  new base T (2/2);  total (32/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.71e-01

CTGCCGGGGCTATGCGGCAAATCTGATATGGTTTTTATCGTAAAATCTGAACCTGTGCTGATAACAGCTTCAGCGCCATAGTGCTCATGCCTGATCGAGGCCTGATGAGCACCGAACCATGAACGATAGAATCCCTTTTACAGAGATCGCCACCGCCCCAGCCAGCATCCAAGCACGACGCACCGACAGCGGCATCCATACCCGCAGCTTTACCGGGCTTTACCGTAACCTGCGCATCGGCTTTGCCGGCGCCCTGTTCGTGCTGTTCTTCGGCACGGCCTGGCTGAACTGG  >  NC_002947/4014527‑4014818                                                                                                                                                    |                                                                                                                                                 cTGCCGGGGCTATGCGGCAAATCTGATATGGTTTTTATCGTAAAATCTGAACCTGTGCTGATAACAGCTTCAGCGCCATAGTGCTCATGCCTGATCGAGGCCTGATGAGCACCGAACCATGAACGATAGAATCCCTTTTACAGAGATc                                                                                                                                                  <  2:349217/148‑1 (MQ=255) cTGCCGGGGCTATGCGGCAAATCTGATATGGTTTTTATCGTAAAATCTGAACCTGTGCTGATAACAGCTTCAGCGCCATAGTGCTCATGCCTGATCGAGGCCTGATGAGCACCGAACCATGAACGATAGAATCCCTTTTACAGAGATc                                                                                                                                                  >  1:349217/1‑148 (MQ=255)             tGCGGCAAATCTGATATGGTTTTTATCGTAAAATCTGAACCTGTGCTGATAACAGCTTCAGCGCCATAGTGCTCATGCCTGATCGAGGCCTGATGAGCACCGAACCATGAACGATAGAATCCCTTTTACAGAGATCGCCACCGCCcca                                                                                                                                      <  7:76459/148‑1 (MQ=255)                     ctCTGAGATGGTTTTTTTCGTAAACTGTGAACCTGTGCTGATAACAGCTTCAGTGAGTTAGTGCTCATGCGTGCTCGAGGCCTGATGAGCTCCGAACTATGAACGATAGAATCGCTTTTACAGAGATCTCCACCGCCCcagcgagctt                                                                                                                              <  8:172926/147‑3 (MQ=255)                        tGATATGGTTTTTATCGTAAAATCTGAACCTGTGCTGATAACAGCTTCAGCGCCATAGTGCTCATGACTGATCGAGGCGTGATGAGCACAGAACCATGAACGATAGAATCCCTTTTACAGAGATCGCCACCGCCCCAGCAAGCAtaaca                                                                                                                          >  7:172926/1‑145 (MQ=255)                          atatGGTTTTTTTCGTAAAATCTGAACCTGTGCTGATAACAGCTTCAGCGCCATAGTGCTCATGCCTGATCGAGGCCTGATGAGCACCGAACCATGAACGATAGAATCCCTTTTACAGAGATCGCCACCGCCCCAGCCAGCATCCaa                                                                                                                          <  5:356695/147‑1 (MQ=255)                          atatGGTTTTTATCGTAAAATCTGAACCTGTGCTGATAACAGCTTCAGCGCCATAGTGCTCATGCCTGATCGAGGCCTGATGAGCACCGAACCATGAACGATAGAATCCCTTTTACAGAGATCGCCACCGCCCCAGCCAGCATCCaa                                                                                                                          <  6:348291/147‑1 (MQ=255)                            atGGTTTTTATCGTAAAATCTGAACCTGTGCTGATAACAGCTTCAGCGCCATAGTGCTCATGCCTGATCGAGGCCTGATGAGCACCGAACCATGAACGATAGAATCCCTTTTACAGAGATCGCCACCGCCCCAGCCAGCATCCaa                                                                                                                          <  8:209497/145‑1 (MQ=255)                            atGGTTTTTATCGTAAAATCTGAACCTGTGCTGATAACAGCTTCAGCGCCATAGTGCTCATGCCTGATCGAGGCCTGATGAGCACCGAACCATGAACGATAGAATCCCTTTTACAGAGATCGCCACCGCCCCAGCCAGCATCCaa                                                                                                                          >  7:209497/1‑145 (MQ=255)                               gTTTTTATCGTAAAATCTGAACCTGTGCTGATAACAGCTTCAGCGCCATAGTGCTCATGCATGATCGAGGCCTGATGAGCACCGAACCATGAACGATAGAATCCCTTTTACAGAGATCGCCACCGCCCCAGCCAGCATCCAAGCacgac                                                                                                                   >  1:296366/1‑149 (MQ=255)                                tttttATCGTAAAATCTGAACCTGTGCTGATAACAGCTTCAGCGCCATAGTGCTCATGCCTGATCGAGGCCTGATGAGCACCGAACCATGAACGATAGAATCCCTTTTACAGAGATCGCCACCGCCCCAGCCAGCATCCAAGCacgacg                                                                                                                  >  1:5681/1‑149 (MQ=255)                                   ttATCGTAAAATCTGAACCTGTGCTGATAACAGCTTCAGCGCCATAGTGCTCATGCCTGATCGAGGCCTGATGAGCACCGAACCATGAACGATAGAATCCCTTTTACAGAGATCGCCACCGCCCCAGCCAGCATCCAAGCACGACGCa                                                                                                                <  2:341406/148‑1 (MQ=255)                                     aTCGTAAAATCTGAACCTGTGCTGATAACAGCTTCAGCGCCATAGTGCTCATGCCTGATCGAGGCCTGATGAGCACCGAACCATGAACGATAGAATCCCTTTTACAGAGATCGCCACCGCCCCAGCCAGCATCCAAGCACGACGCAcc                                                                                                              <  2:5681/148‑1 (MQ=255)                                         tAAAATCTGAACCTGTGCTGATAACAGCTTCAGCGCCATAGTGCTCATGCCTGATCGAGGCCTGATGAGCACCGAACCATGAACGATAGAATCCCTTTTACAGAGATTGCCACCGCCCCAGCCAGCATCCAAGc                                                                                                                        >  1:96724/1‑134 (MQ=255)                                         tAAAATCTGAACCTGTGCTGATAACAGCTTCAGCGCCATAGTGCTCATGCCTGATCGAGGCCTGATGAGCACCGAACCATGAACGATAGAATCCCTTTTACAGAGATTGCCACCGCCCCAGCCAGCATCCAAGc                                                                                                                        <  2:96724/134‑1 (MQ=255)                                         tAAAATCTGAACCTGTGCTGATAACAGCTTCAGCGCCATAGTGCTCATGCCTGATCGAGGCCTGATGAGCACCGAACCATGAACGATAGAATCCCTTTTACAGAGATCGCCACCGCCCCAGCCAGCATCCAAGCACGACGCACCGACAg                                                                                                         >  2:132077/1‑149 (MQ=255)                                         tAAAATCTGAACCTGTGCTGATAACAGCTTCAGCGCAATAGTGCTCATGCCTGATCGAGGCCTGATGAGCACCGAACCATGAACGATAGAATCCCTTTTACAGAGATTGCCACCGACCCAGCCAGCATACAAGc                                                                                                                        >  1:101648/1‑134 (MQ=255)                                         tAAAATCTGAACCTGTGATGATAACAGCTTCAGCGCCATAGTTCTCATGCCTGATAGAGGCCTGATGAGCACCGAACCATGAACGATAGAATCCCTTTTACAGAGATTGCCACCGCCCCAGCCAGCATCCAAGc                                                                                                                        <  2:101648/134‑1 (MQ=255)                                          aaaaTCTGAACCTGTGCTGATAACAGCTTCAGCGCCATAGTGCTCATGCCTGATCGAGGCCTGATGAGCACCGAACCATGAACGATAGAATCCCTTTTACAGAGATCGCCACCGCCCCAGCCAGCATCCAAGCACGACGCAcc                                                                                                              >  4:242568/1‑143 (MQ=255)                                          aaaaTCTGAACCTGTGCTGATAACAGCTTCAGCGCCATAGTGCTCATGCCTGATCGAGGCCTGATGAGCACCGAACCATGAACGATAGAATCCCTTTTACAGAGATCGCCACCGCCCCAGCCAGCATCCAAGCACGACGCAcc                                                                                                              <  3:242568/143‑1 (MQ=255)                                           aaaTCTGAACCTGTGCTGATAACAGCTTCAGCGCCATAGTGCTCATGCCTGATCGAGGCCTGATGAGCACCGAACCATGAACGATAGAATCCCTTTTACAGAGATCGCCACCGCCCCAGCCAGCATCCAAGCACGACGCACCGACAGc                                                                                                        <  2:43940/148‑1 (MQ=255)                                                     cTGTGCTGATAACAGCTTCAGCGCCATAGTGCTCATGCCTGATCGAGGCCTGATGAGCGCCGAACCATGAACGATAGAATCCCTTTTACAGAGATCGCCACCGCCCCAGCCAGCATCCAAGCACGACGCACCGACAGCGGCATCCATAc                                                                                             >  5:372731/1‑149 (MQ=255)                                                     cTGTGCTGATAACAGCTTCAGCGCCATAGTGCTCATGCCTGATCGAGGCCTGATGAGCCCCGAACCATGAACGATAGACTCCCTATTACAGAGATCGCCACCGCCCCAGCCAGCATCCAAGCACGACGCACCGACAGCGGCATCCATAc                                                                                             >  5:269774/1‑149 (MQ=255)                                                           tGATAACAGCTTCAGCGCCATAGTGCTCATGCCTGATCGAGGCCTGATGAGCACCGAACCATGAACGATAGAATCCCTTTTACAGAGATCGCCACCGCCCCAGCCAGCATCCAAGCACGACGCACCGACAGCGGCATCCATACCCGCAg                                                                                       >  1:14098/1‑149 (MQ=255)                                                             aTAACAGCTTCAGCGCCATAGTGCTCATGCCTGATCGAGGCCTGATGAGCACCGAACCATGAACGATAGAATCCCTTTTACAGAGATCGCCACCGCCCCAGCCAGCATCCAAGCACGACGCACCGACAGCGGCATCCATACCCGCAGc                                                                                      <  8:94532/148‑1 (MQ=255)                                                                 cAGCTTCAGCGCCATAGTGCTCATGCCTGATCGAGGCCTGATGAGCACCGACCCATGAACGATAGCATCCCTTTTACAGAGATCGCCACCGCCCCAGCCAGCATCCAAGCACGACGCACCGACAGCGGCATCCATACCCGCAGCTTTAc                                                                                 >  5:267035/1‑149 (MQ=255)                                                                      tCAGCGCCATAGTGCTCATGCCTGATCGAGGCCTGATGAGCACCGAACCATGAACGATAGAATCCCTTTTACAGAGATCGCCACCGCCCCAGCCAGCATCCAAGCACGACGCACCGACAGCGGCATCCATACCCGCAGCTTTACCGGGc                                                                            >  4:359097/1‑149 (MQ=255)                                                                        aGCGCCATAGTGCTCATGCCTGATCGAGGCCTGATGAGCACCGAACCATGAACGATAGAATCCCTTTTACAGAGATCGCCACCGCCCCAGCCAGCATCCAAGCACGACGCACCGACAGCGGCATCCATACCCGCAGCTTTACCGGGCtt                                                                          >  2:30416/1‑149 (MQ=255)                                                                         gcgcCATAGTGCTCATGCCTGATCGAGGCCTGATGAGCACCGAACCATGAACGATAGAATCCCTTTTACAGAGATCGCCACCGCCCCAGCCAGCATCCAAGCACGACGCACCGACAGCGGCATCCATACCCGCAGCTTTACCGGGCtt                                                                          >  7:210344/1‑148 (MQ=255)                                                                          cgcCATAGTGCTCAAGCCTGATCGAGGCCTGATGAGCACCGAACCATGAACGATAGAATCCCTTTTACAGAGATCGCCACCGCCCCAGCCAGCATCCAAGCACGACGCACCGACAGCGGCATCCATACCCGCAGCTTTACCGGGCTTTa                                                                        <  1:291921/149‑1 (MQ=255)                                                                           gcCATAGTGCTCATGCCTGATCGAGGCCTGATGAGCACCGAACCATGAACGATAGAATCCCTTTTACAGAGATCGCCACCGCCCCAGCCAGCATCCAAGCACGACGCACCGACAGCGGCATCCATAGCCGCAGCTTTACCGGGCTTTAc                                                                       <  8:210344/149‑1 (MQ=255)                                                                             cATAGTGCTCATGCCTGATCGAGGCCTGATGAGCACCGAACCATGAACGATAGAATCCCTTTTACAGAGATCGCCACCGCCCCAGCCAGCATCCAAGCACGACGCACCGACAGCGGCATCCATACCCGCAGCTTTACCGGGCTTTACCg                                                                     >  3:226180/1‑149 (MQ=255)                                                                                  tGCTCATGCCTGATCGAGGCCTGATGAGCACCGAACCATGAACGATAGAATCCCTTTTACAGAGATCGCCACCGCCCCAGCCAGCATCCAAGCACGACGCACCGACAGCGGCATCCATACCCGCAGCTTTACCGGGCTTTACCGTAAcc                                                                <  1:132077/149‑1 (MQ=255)                                                                                    cTCATGCCTGATCGAGGCCTGATGAGCACCGAACCATGAACGATAGAATCCCTTTTACAGAGATCGCCACCGCCCCAGCCAGCATCCAAGCACGACGCACCGACAGCGGCATCCATACCCGCAGCTTTACCGGGCTTTACCGTAACCTg                                                              >  4:52119/1‑149 (MQ=255)                                                                                       aTGCCTGATCGAGGCCTGATGAGCGCCGAACCATGAACGATAGAATCCCTTTTACAGAGATCGCCACCGCCCCAGCCAGCATCCAAGCACGACGCACCGACAGCGGCATCCATACCCGCAGCTTTACCGGGCTTTACCGTAACCTgc                                                             <  6:372731/147‑1 (MQ=255)                                                                                       aTGCCTGATCGAGGCCTGATGAGCACCGAACCATGAACGATAGAATCCCTTTTACAGAGATCGCCACCGCCCCAGCCAGCATCCAAGCACGACGCACCGACAGCGGCATCCATACCCGCAGCTTTACCGGGCTTTACCg                                                                     >  7:170675/1‑139 (MQ=255)                                                                                       aTGCCTGATCGAGGCCTGATGAGCACCGAACCATGAACGATAGAATCCCTTTTACAGAGATCGCCACCGCCCCAGCCAGCATCCAAGCACGACGCACCGACAGCGGCATCCATACCCGCAGCTTTACCGGGCTTTACCg                                                                     <  8:170675/139‑1 (MQ=255)                                                                                         gCCTGATCGAGGCCTGATGAGCACCGAACCATGAACGATAGAATCCCTTTTACAGAGATCGCCACCGCCCCAGCTAGCATCCAAGCACGACGCACCGACAGCGGCATCCATACCCGCAGCTTTACCGGGCTTTACCGTAACCTGCGCAt                                                         >  1:344187/1‑149 (MQ=255)                                                                                          ccTGATCGAGGCCTGATGAGCACCGAACCATGAACGATAGAATCCCTTTTACAGAGATCGCCACCGCCCCAGCCAGCATCCAAGCACGACGCACCGACAGCGGCATCCATACCCGCAGCTTTACCGGGCTTTACCGTAACCTGCGCATc                                                        <  6:271047/149‑1 (MQ=255)                                                                                           cTGATCGAGGCCTGATGAGCACCGAACCATGAACGATAGAATCCCTTTTACAGAGATCGCCACCGCCCCAGCCAGCATCCAAGCACGACGCACCGACAGCGGCATCCATACCCGCAGCTTTACCgg                                                                              <  7:3856/126‑1 (MQ=255)                                                                                           cTGATCGAGGCCTGATGAGCACCGAACCATGAACGATAGAATCCCTTTTACAGAGATCGCCACCGCCCCAGCCAGCATCCAAGCACGACGCACCGACAGCGGCATCCATACCCGCAGCTTTACCgg                                                                              >  8:3856/1‑126 (MQ=255)                                                                                              aTCGAGGCCTGATGAGCACCGAACCATGAACGATAGAATCCCTTTTACATAGATCGCCACCGCCCCAGCCAGCATCCAAGCACGACGCACCGACACCGGCATCCATACCCGCAGCTTTACCGCGCATTACCGTAACCTGCGCATCGGCt                                                    >  3:238823/1‑149 (MQ=255)                                                                                              aTCGAGGCCTGATGAGCACCGAACCATGAACGATAGAATCCCTTTTACAGAGATCGCCACCGCCCCAGCCAGCATCCAAGCACGACGCACCGACAGCGGCATCCATACCCGCAGCTTTACCGGGCTTTACCGTAACCTGCGCATCGGCt                                                    >  3:340666/1‑149 (MQ=255)                                                                                                cGAGGCCTGATGAGCACCGAACCATGAACGATAGAATCCCTTTTACAGAGATCGCCAcc                                                                                                                                            <  6:153783/59‑1 (MQ=255)                                                                                                cGAGGCCTGATGAGCACCGAACCATGAACGATAGAATCCCTTTTACAGAGATCGCCAcc                                                                                                                                            >  5:153783/1‑59 (MQ=255)                                                                                                   tcaggcatgaGCACCGAACCATGAACGATAGAATCCCTTTTACAGAGATCGCCACCGCCCCAGCCAGCATCCAAGCACGACGCACCGACAGCGGCATCCATACCCGCAGCTTTACCGGGCTTTACCGTAACCTGCGCATCGGCTTTGcc                                               <  1:199982/143‑1 (MQ=255)                                                                                                    gCCTGATGAGCACCGAACCATGAACGATAGAATCCCTTTTACAGAGATCGCCACCGCCCCAGCCAGCATCCAAGCACGACGCACCGACAGCGGCATCCATACCCGCAGCTTTACCGGGCTTTACCGTAAc                                                                 <  6:218597/130‑1 (MQ=255)                                                                                                    gCCTGATGAGCACCGAACCATGAACGATAGAATCCCTTTTACAGAGATCGCCACCGCCCCAGCCAGCATCCAAGCACGACGCACCGACAGCGGCATCCATACCCGCAGCTTTACCGGGCTTTACCGTAAc                                                                 >  5:218597/1‑130 (MQ=255)                                                                                                      cggatgaGCACCGAACCCTGAACGCTAGAATCCCTTTTACAGCGATCGCCACCGCCCCAGCCAGCATCCAAGCACGACGCACCGACAGCGGCATCCATACCCGCAGCTTTACCGGGCTTTACCGTAACCTGCGCATCGGCTTTGCCGgc                                            <  6:267035/147‑1 (MQ=255)                                                                                                      cTGATGAGCAGCGAACCATGAACGATAGAATCCCTTTTACAGAGATCGCCACCGCGCCAGCCAGCATCCAAGCACGACGCACCGACAGCGGCATCCATACCCGCAGCTTTACCGGGCTTTACCGTAACCTGCGCATCGGCTTTGCCGgc                                            <  1:287174/149‑1 (MQ=255)                                                                                                       tgatgaGCACCGAACCATGAACGATAGAATCCCTTTTACAGAGATCGCCACCGCCCCAGCCAGCATCCAAGCACGACGCACCGACAGCGGCATCCATACCCGCAGCTTTACCGGGCTTTACCGTAACCTGCGCATCGGCTTTGCCGgcg                                           >  6:333249/1‑149 (MQ=255)                                                                                                         atgaGCACCGAAACATGAACGATAGAATCCCTTTTACAGAGATCGCCACCGCCCCAGCCAGCATCCAAGCACGACGCACCGACAGCGGCATCCATACCCGCAGCTTTACCGGGCTTTACCGTAACCTGCGCATCGGCTTTGCCGGCGcc                                         <  6:269774/149‑1 (MQ=255)                                                                                                                 cGAACCATGAACGATAGAATCCCTTTTACAGAGATCGCCACCGCCCCAGCCAGCATCCAAGCACGACGCACCGACAGCGGCATCCATACCCGCAGCTTTACCGGGCTTTACCGTAACCTGCGCATCGGCTTTGCCGGCGCCCTGTTCGt                                 >  2:161154/1‑149 (MQ=255)                                                                                                                  gAACCATGAACGATAGAATCCCTTTTACAGAGATCGCCACCGCCCCAGCCAGCATCCAAGCACGACGCACCGACAGCGGCATCCATACCCGCAGCTTTACCGGGCTTTACCGTAACCTGCGCATCGGCTTTGCCGGCGCCCTGTTCGTg                                <  4:340666/149‑1 (MQ=255)                                                                                                                        tGAACGATAGAATCCCTTTTACAGAGATCGCCACCGCCCCAGCCAGCATCCAAGCACGACGCACCGACAGCGGCATCCATACCCGCAGCTTTACCGGGCTTTACCGTAACCTGCGCATCGGCTTTGCCGGCGCCCTGTTCGTGCTGttc                          >  6:317298/1‑149 (MQ=255)                                                                                                                             gATAGAATCCCTTTTACAGAGATCGCCACCGCCCCAGCCAGCATCCAAGCACGACGCACCGACAGCGGCATCCATACCCGCAGCTTTACCGGGCTTTACCGTAACCTGCGCATCGGCTTTGCCGGCGCCCTGTTCGTGCTGTTCTTCgg                     <  3:52119/149‑1 (MQ=255)                                                                                                                             gATAGAATCCCTTTTACAGAGATCGCCACCGCCCCAGCCAGCATCCAAGCACGACGCACCGACAGCGGCATCCATACCCGCAGCTTTACCGGGCTTTACCGTAACCTGCGCATCGGCTATGCCGGCGCCCTGTTCGTGCTGTTCTTCgg                     <  1:30416/149‑1 (MQ=255)                                                                                                                                     cccTTTTACAGAGATCGCCACCGCCCCAGCCAGCATCCAAGCACGACGCACCGACAGCGGCATCCATACCCGCAGCTTTACCGGGCTTTACCGTAACCTGCGCATCGGCTTTGCCGGCGCCCTGTTCGTGCTGTTCTTCGGCACGGCCt             >  7:243162/1‑149 (MQ=255)                                                                                                                                       cTTTTACAGAGATCGCCACCGCCCCAGCCAGCATCCAAGCACGACGCACCGACAGCGGCATCCATACCCGCAGCTTTACCGGGCTTTACCGTAACCTGCGCTTCGGCTTTGCCGGCGCCCTGTTCGTGCTGTTCTTCGGCACGGCCTgg           <  2:14098/149‑1 (MQ=255)                                                                                                                                          ttACAGAGATCGCCACCGCCCCAGCCAGCATCCAAGCACGACGCACCGACAGCGGCATCCATACCCGCAGCTTTACCGGGCTTTACCGTAACCTGCGCATCGGCTTTGCCGGCGCCCTGTTCGTGCTGTTCTTCGGCACGGCCTGGCt         <  1:313411/148‑1 (MQ=255)                                                                                                                                          ttACAGAGATCGCCACCGCCCCAGCCAGCATCCAAGCACGACGCACCGACAGCGGCATCCATACCCGCAGCTTTACCGGGCTTTACCGTAACCTGCGCATCGGCTTTGCCGGCGCCCTGTTCGTGCTGTTCTTCGGCACGGCCTGGCt         >  2:313411/1‑148 (MQ=255)                                                                                                                                              agagATCGCCACCGCCCCAGCCAGCATCCAAGCACGACGCACCGACAGCGGCATCCATACCCGCAGCTTTACCGGGCTTTACCGTAACCTGCGCATCGGCTTTGCCGGCGCCCTGTTCGTGCTGTTCTTCGGCACGGCCTGGCTGAACt    >  3:124782/1‑149 (MQ=255)                                                                                                                                                agaTCGCCACCGCCCCAGCCAGCATCCAAGCACGACGCACCGACAGCGGCATCCATACCCGCAGCTTTACCGGGCTTTACCGTAACCTGCGCATCGGCTTTGCCGGCGCCCTGTTCGTGCTGTTCTTCGGCACGGCCTGGCTGAACTgg  <  8:243162/149‑1 (MQ=255)                                                                                                                                                agaTCGCCACCGCCCCAGCCAGCATCCAAGCACGACGCACCGACAGCGGCATCCATACCCGCAGCTTTACCGGGCTTTACCGTAACCTGCGCATCGGCTTTGCCGGCGCCCTGTTCGTGCTGTTCTTCGGCACGGCCTGGCTGAACTg   >  1:116423/1‑148 (MQ=255)                                                                                                                                                agaTCGCCACCGCCCCAGCCAGCATCCAAGCACGACGCACCGACAGCGGCATCCATACCCGCAGCTTTACCGGGCTTTACCGTAACCTGCGCATCGGCTTTGCCGGCGCCCTGTTCGTGCTGTTCTTCGGCACGGCCTGGCTGAACTg   >  2:37053/1‑148 (MQ=255)                                                                                                                                                    cGCCACCGCCCCAGCCAGCATCCAAGCACGACGCACCGACAGCGGCATCCATACCCGCAGCTTTACCGGGCTTTACCGTAACCTGCGCATCGGCTTTGCCGGCGCCCTGTTCGTGCTGTTCTTCGGCACGGCCTGGCt         <  4:158176/138‑1 (MQ=255)                                                                                                                                                    cGCCACCGCCCCAGCCAGCATCCAAGCACGACGCACCGACAGCGGCATCCATACCCGCAGCTTTACCGGGCTTTACCGTAACCTGCGCATCGGCTTTGCCGGCGCCCTGTTCGTGCTGTTCTTCGGCACGGCCTGGCt         >  3:158176/1‑138 (MQ=255)                                                                                                                                                    |                                                                                                                                                 CTGCCGGGGCTATGCGGCAAATCTGATATGGTTTTTATCGTAAAATCTGAACCTGTGCTGATAACAGCTTCAGCGCCATAGTGCTCATGCCTGATCGAGGCCTGATGAGCACCGAACCATGAACGATAGAATCCCTTTTACAGAGATCGCCACCGCCCCAGCCAGCATCCAAGCACGACGCACCGACAGCGGCATCCATACCCGCAGCTTTACCGGGCTTTACCGTAACCTGCGCATCGGCTTTGCCGGCGCCCTGTTCGTGCTGTTCTTCGGCACGGCCTGGCTGAACTGG  >  NC_002947/4014527‑4014818

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