Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,587,764	T→C	5.5%	S462G (AGC→GGC)	parE ←	DNA topoisomerase IV subunit B

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,587,764	0	T	C	5.5%	265.5 / 4.3	91	S462G (AGC→GGC)	parE	DNA topoisomerase IV subunit B
Reads supporting (aligned to +/- strand):  ref base T (42/44);  new base C (3/2);  total (45/46)
Fisher's exact test for biased strand distribution p-value = 6.77e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.30e-01

GAACAGCGCACACAGCAATGTGGCGATGTGCAGACCATCGGAGTCAGCGTCGGCGAGGATGCAGATCTTGCCATAGCGCAGCTGCGCAAGGTCGGTTGCACCCG‑GA‑TCGACGCCAATGGCCA‑CGGCAATGTTGTGCACTTCCTGGCTGGCCAGGACTTCGCCCCCGTCCA‑‑CTTCC‑‑CAGGTGTTGAGGATCTTGCCGCGCAGCGGCAGGATGGCCTGGAATTCCTTGTCCCGCGCCTGCTTGGCCGAACCACCGGCCGAGTCACCCTCGACGAGGAACAACTCGGCAC  >  NC_002947/5587618‑5587904                                                                                                                                                      |                                                                                                                                                 gaacAGCGCACACAGCAATGTGGCGATGTGCAGACCATCGGAGTCAGCGTCGGCGAGGATGCAGATCTTGCCATAGCGCAGCTGCGCAAGGTCGGTTGCACCCG‑GA‑TCGACGCCAATGGCCA‑CGGCAATGTTGTGCACTTCctggctgg                                                                                                                                                <  3:371066/149‑1 (MQ=255)   acAGCGCACACAGCAATGTGGCGATGTGCAGACCATCGGAGTCAGCGTCGGCGAGGATGCAGATCTTGCCATAGCGCAGCTGCGCAAGGTCGGTTGCACCCG‑GA‑TCGACGCCAATGGCCA‑CGGCAATGTTGTGCACTTCCtggctggc                                                                                                                                               >  2:278117/1‑148 (MQ=255)    cAGCGCACACAGCAATGTGGCGATGTGCAGACCATCGGAGTCAGCGTCGGCGAGGATGCAGATCTTGCCATAGCGCAGCTGCGCAAGGTCGGTTGCACCCG‑GA‑TCGACGCCAATGGCCA‑CGGCAATGTTGTGCACTTCCtggctggc                                                                                                                                               >  3:263793/1‑147 (MQ=255)    cAGCGCACACAGCAATGTGGCGATGTGCAGACCATCGGAGTCAGCGTCGGCGAGGATGCAGATCTTGCCATAGCGCAGCTGCGCAAGGTCGGTTGCACCCG‑GA‑TCGACGCCAATGGCCA‑CGGCAATGTTGTGCACTTCCTGGCTGGCCa                                                                                                                                             <  7:60934/149‑1 (MQ=255)     aGCGCACACAGCAATGTGGCGATGTGCAGACCATCGGAGTCAGCGTCGGCGAGGATGCAGATCTTGCCATAGCGCAGCTGCGCAAGGTCGGTTGCACCCG‑GA‑TCGACGCCAATGGCCA‑CGGCAATGTTGTGCACTTCCTGGCTGGCCAg                                                                                                                                            >  3:361162/1‑149 (MQ=255)     aGCGCACACAGCAATGTGGCGATGTGCAGACCATCGGAGTCAGCGTCGGCGAGGATGCAGATCTTGCCATAGCGCAGCTGCGCAAGGTCGGTTGCACCCG‑GA‑TCGACGCCAATGGCCA‑CGGCAATGTTGTGCACTTCCTGGCTGGCCAg                                                                                                                                            >  3:367225/1‑149 (MQ=255)       cgcACACAGCAATGTGGCGATGTGCAGACCATCGGAGTCAGCGTCGGCGAGGATGCAGATCTTGCCATAGCGCAGCTGCGCAAGGTCGGTTGCACCCG‑GA‑TCGACGCCAATGGCCA‑CGGCAATGTTGTGCACTTCCTGGCTGGCCAGGa                                                                                                                                          <  6:181118/149‑1 (MQ=255)       cgcACACAGCAATGTGGCGATGTGCAGACCATCGGAGTCAGCGTCGGCGAGGATGCAGATCTTGCCATAGCGCAGCTGCGCAAGGTCGGTTGCACCCG‑GA‑TCGACGCCAATGGCCA‑CGGCAATGTTGTGCACTTCCTGGCTGGCCAGGa                                                                                                                                          >  6:420541/1‑149 (MQ=255)        gcACACAGCAATGTGGCGATGTGCAGACCATCGGAGTCAGCGTCGGCGAGGATGCAGATCTTGCCATAGCGCAGCTGCGCAAGGTCGGTTGCACCCG‑GA‑TCGACGCCAATGGCCA‑CGGCAATGTTGTGCACTTCCTGGCCAGCCAGGaa                                                                                                                                         <  7:81485/149‑2 (MQ=255)        gcACACAGCAATGTGGCGATGTGCAGACCATCGGAGTCAGCGTCGGCGAGGATGCAGATCTTGCCATAGCGCAGCTGCGCAAGGTCGGTTGCACCCG‑GA‑TCGACGCCAATGGCCA‑CGGCAATGTTGTGCACGTCCTGGCCAGCCAGGaa                                                                                                                                         <  3:133559/149‑2 (MQ=255)          acacaGCAATGTGGCGATGTGCAGACCATCGGAGTCAGCGTCGGCGAGGATGCAGATCTTGCCATAGCGCAGCTGCGCAAGGTCGGTTGCACCCG‑GA‑TCGACGCCAATGGCCA‑CGGCAATGTTGTGCACTTCCTGGCTGGCCAGGACtt                                                                                                                                       >  1:155267/1‑149 (MQ=255)            acaGCAATGTGGCGATGTGCAGACCATCGGAGTCAGCGTCGGCGAGGATGCAGATCTTGCCATAGCGCAGCTGCGCAAGGTCGGTTGCACCCG‑GA‑TCGACGCCAATGGCCA‑CGGCAATGTTGTGCACTTCCTGGCTGGCCAGGACTTCg                                                                                                                                     <  7:225772/149‑1 (MQ=255)                cAATGTGGCGATGTGCAGACCATCGGAGTCAGCGTCGGCGAGGATGCAGATCTTGCCATAGCGCAGCTGCGCAAGGTCGGTTGCACCCG‑GA‑TCGACGCCAATGGCCA‑CGGCAATGTTGTGCACTTCCTGGCTGGCCAGGACTTCGcccc                                                                                                                                 >  2:152688/1‑149 (MQ=255)                cAATGTGGCGATGTGCAGACCATCGGAGTCAGCGTCAGCGAGGATGCAGATCTTACCATAGCGCAGCTGCGCAAGGTCGGTTGCACCCG‑GA‑TCGACGCCAATGGCCA‑CGGCAATGTTGTGCACTTCCTGGCTGGCCAGGACTTCGcccc                                                                                                                                 >  4:145839/1‑149 (MQ=255)                            gtgCAGACCATCGGAGTCAGCGTCGGCGAGGATGCAGATCTTGCCATAGCGCAGCTGCGCAAGTTCGGTTGCACCCG‑GA‑TCTACGCCAATGGCCA‑CGGCAATGTTGTGCACTTCCTGGCTGGCCAGGACTTCGCCCCCTTCCA‑‑CTTcc                                                                                                                    <  6:421584/148‑1 (MQ=255)                             tgCAGACCATCGGAGTCAGCGTCGGCGAGGATGCAGATCTTGCCATAGCGCAGCTGCGCAAGGTCGGTTGCACCCG‑GA‑TCGACGCCAATGGCCA‑CGGCAATGTTGTGCACTTCCTGGCTGGCCAGGACTTCGCCCCCGTCCA‑‑CTTcc‑‑c                                                                                                                 <  4:145047/148‑1 (MQ=255)                               cAGACCATCGGAGTCAGCGTCGGCGAGGATGCAGATCTTGCCATAGCGCAGCTGCGCAAGGTCGGTTGCACCCG‑GA‑TCGACGCCAATGGCCA‑CGGCAATGTTGTGCACTTCCTGGCTGGCCAGGACTTCGCCCCCGTCCA‑‑CTTCC‑‑CAgg                                                                                                              >  7:97332/1‑149 (MQ=255)                               cAGACCATCGGAGTCAGCGTCGGCGAGGATGCAGATCTTGCCATAGCGCAGCTGCGCAAGGTCGGTTGCACCCG‑GA‑TCGACGCCAATGGCCA‑CGGCAATGTTGTGCACTTCCTGGCTGGCCAGGACTTCGCCCCCGTCCA‑‑CTTCC‑‑CAg                                                                                                               <  1:152688/148‑1 (MQ=255)                                   ccATCGGAGTCAGCGTCGGCGAGGATGCAGATCTTGCCATAGCGCAGCTGCGCAAGGTCGGTTGCACCCGATAGAGGACGCCAATGGCCAGCGGCAATGTTGTGCACTTACTGGCTGTCCAGGACTTCGCCCCCGTCCA‑‑CTTCC‑‑CAGgt                                                                                                             >  1:353818/1‑149 (MQ=255)                                   ccATCGGAGTCAGCGTCGGCGAGGATGCAGATCTTGCCATAGCGCAGCTGCGCAAGGTCGGTTGCACCCG‑GA‑TCGACGCCAATGGCCA‑CGGCAATGTTGTGCACTTCCTGGCTGGCCAGGACTTCGCCCCCGTCCA‑‑CTTCC‑‑CAGgtgt                                                                                                           <  1:269587/148‑1 (MQ=255)                                   ccATCGGAGTCAGCGTCGGCGAGGATGCAGATCTTGCCATAGCGCAGCTGCGCAAGGTCGGTTGCACCCG‑GA‑TCGACGCCAATGGCCA‑CGGCAATGTTGTGCACTTCCTGGCTGGCCAGGACTTCGCCCCCGTCCA‑‑CTTCC‑‑CAGGTGtt                                                                                                          <  1:410984/149‑1 (MQ=255)                                    cATCGGAGTCAGCGTCGGCGAGGATGCAGATCTTGCCATAGCGCAGCTGCGCAAGGTCGGTTGCACCCG‑GA‑TCGACGCCAATGGCCA‑CGGCAATGTTGTGCACTTCCTGGCTGGCCAGGACTTCGCCCCCGTCCa                                                                                                                           <  5:31682/135‑1 (MQ=255)                                    cATCGGAGTCAGCGTCGGCGAGGATGCAGATCTTGCCATAGCGCAGCTGCGCAAGGTCGGTTGCACCCG‑GA‑TCGACGCCAATGGCCA‑CGGCAATGTTGTGCACTTCCTGGCTGGCCAGGACTTCGCCCCCGTCCa                                                                                                                           >  6:31682/1‑135 (MQ=255)                                     aTCGGAGTCAGCGTCGGCGAGGATGCAGATCTTGCCATAGCGCAGCTGCGCAAGGTCGGTTGCACCCG‑GA‑TCGACGCCAATGGCCA‑CGGCAATGTTGTGCACTTCCTGGCTGGCCAGGACTTCGCCCCCGTCCA‑‑CTTCC‑‑CAGGTGtt                                                                                                          >  2:318600/1‑147 (MQ=255)                                     aTCGGAGTCAGCGTCGGCGAGGATGCAGATCTTGCCATAGCGCAGCTGCGCAAGGTCGGTTGCACCCG‑GA‑TCGACGCCAATGGCCA‑CGGCAATGTTGTGCACTTCCTGGCTGGCCAGGACTTCGCCCCCGTCCA‑‑CTTCC‑‑CAGGTGtt                                                                                                          >  8:251526/1‑147 (MQ=255)                                     aTCGGAGTCAGCGTCGGCGAGGATGCAGATCTTGCCATAGCGCAGCTGCGCAAGGTCGGTTGCACCCG‑GA‑TCGACGCCAATGGCCA‑CGGCAATGTTGTGCACTTCCTGGCTGGCCAGGACTTCGCCCCCGTCCA‑‑CTTCC‑‑CAGGTGTTGa                                                                                                        <  5:420541/149‑1 (MQ=255)                                     aTCGGAGTCAGCGTCGGCGAGGAGGCAGATCTTGCCATAGCGCAGCTGCGCAAGGTCGGTTGCACCCG‑GA‑TCGACGCCAATGGCCA‑CGGCAATGTTGTGCACTTCCTGGCTGGCCAGGACTTCGCCCCCGTCCA‑‑CTTCC‑‑CAGGTGTTGa                                                                                                        <  8:355321/149‑1 (MQ=255)                                      tCGGAGTCAGCGTCGGCGAGGATGCAGATCTTGCCATAGCGCAGCTGCGCAAGGTCGGTTGCACCCG‑GA‑TCGACGCCAATGGCCA‑CGGCAATGTTGTGCACTTCCTGGCTGGCCAGGACTTCGCCCCCGTCCA‑‑CTTCC‑‑CAGGTGTTGAg                                                                                                       >  6:27816/1‑149 (MQ=255)                                          aGTCAGCGTCGGCGAGGATGCAGATCTTGCCATAGCGCAGCTGCGCAAGGTCGGTTGCACCCG‑GA‑TCGACGCCAATGGCCA‑CGGCAATGTTGTGCACTTCCTGGCTGGCCAGGACTTCGCCCCCGTCCA‑‑CTTCC‑‑CAGGTGTTGAGGATc                                                                                                   <  6:278123/149‑1 (MQ=255)                                            tCAGCGTCGGCGAGGATGCAGATCTTGCCATAGCGCAGCTGCGCAAGGTCGGTTGCACCCG‑GA‑TCGACGCCAATGGCCA‑CGGCAATGTTGTGCACTTCCTGGCTGGCCAGGACTTCGCCCCCGTCCA‑‑CTTCC‑‑CAGGTGTTGAGGATCtt                                                                                                 >  5:168783/1‑149 (MQ=255)                                            tCAGCGTCGGCGAGGATGCAGATCTTGCCATAGCGCAGCTGCGCAAGGTCGGTTGCACCCG‑GA‑TCGACGCCAATGGCCA‑CGGCAATGTTGTGCACTTCCTGGCTGGCCAGGACTTCGCCCCCGTCCA‑‑CTTCC‑‑CAGGTGTTGAGGATCtt                                                                                                 >  4:101753/1‑149 (MQ=255)                                               gCGTCGGCGAGGATGCAGATCTTGCCATAGCGCAGCTGCGCAAGGTCGGTTGCACCCG‑GA‑TCGACGCCAATGGCCA‑CGGCAATGTTGTGCACTTCCTGGCTGGCCAGGACTTCGCCCCCGTCCA‑‑CTTCC‑‑CAGGTGTTGAGGATCTTGcc                                                                                              >  2:141532/1‑149 (MQ=255)                                                   cGGCGAGGATGCAGATCTTGCCATAGCGCAGCTGCGCAAGGTCGGTTGCACCCG‑GA‑TCGACGCCAATGGCCA‑CGGCAATGTTGTGCACTTCCTGGCTGGCCAGGACTTCGCCCCCGTCCA‑‑CTTCC‑‑CAGGTGTTGAGGATCTTGCcgcg                                                                                           >  6:182536/1‑148 (MQ=255)                                                     gCGAGGATGCAGATCTTGCCATAGCGCAGCTGCGCAAGGTCGGTTGCACCCG‑GA‑TCGACGCCAATGGCCA‑CGGCAATGTTGTGCACTTCCTGGCTGGCCAGGACTTCGCCCCCGTCCA‑‑CTTCC‑‑CAGGTGTTGAGGATCTTGCCGCGCAg                                                                                        >  5:217462/1‑149 (MQ=255)                                                       gAGGATGCAGATCTTGCCATAGCGCAGCTGCGCAAGGTCGGTTGCACCCG‑GA‑TCGACGCCAATGGCCA‑CGGCAATGTTGTGCACTTCCTGGCTGGCCAGGACTTCGCCCCCGTCCA‑‑CTTCC‑‑CAGGTGTTGAGGATCTTGCCGCGCAGCg                                                                                      >  3:211701/1‑149 (MQ=255)                                                           aTGCAGATCTTGCCATAGCGCAGCTGCGCAAGGTCGGTTGCACCCG‑GA‑TCGACGCCAATGGCCA‑CGGCAATGTTGTGCACTTCCTGGCTGGCCAGGACTTCGCCCCCGTCCA‑‑CTTCC‑‑CAGGTGTTGAGGATCTTGCCGCGCAGCGGCAg                                                                                  <  7:219374/149‑1 (MQ=255)                                                            tGCAGATCTTGCCATAGCGCAGCTGCGCAAGGTCGGTTGCACCCG‑GA‑TCGACGCCAATGGCCA‑CGGCAATGTTGTGCACTTCCTGGCTGGCCAGGACTTCGCCCCCGTCCA‑‑CTTCC‑‑CAGGTGTTGAGGATCTTGCCGCGCAGCGGCAgg                                                                                 <  6:217462/149‑1 (MQ=255)                                                             gCAGATCTTGCCATAGCGCAGCTGCGCAAGGTCGGTTGCACCCG‑GA‑TCGACGCCAATGGCCA‑CGGCAATGTTGTGCACTTCCTGGCTGGCCAGGACTTCGCCCCCGTCCA‑‑CTTCC‑‑CAGGTGTTGAGGATCTTGCCGCGCAGCGGCAGGa                                                                                >  5:77339/1‑149 (MQ=255)                                                              cAGATCTTACCATAGCGCAGCTGCGCAAGGTCGGTTGCACCCG‑GA‑TCGACGCCAATGGCCA‑CGGCAATGTTGTGCACTTCCTGGCTGGCCAGGACTTCGCCCCCGTCCA‑‑CTTCC‑‑CAGGTGTTGAGGATCTTGCCGCGCAGCGGCAGGa                                                                                <  3:145839/148‑1 (MQ=255)                                                                 aTCTTGCCATAGCGCAGCTGCGCAAGGTCGGTTGCACCCG‑GA‑TCGACGCCAATGGCCA‑CGGCAATGTTGTGCACTTCCTGGCTGGCCAGGACTTCGCCCCCGTCCA‑‑CTTCC‑‑CAGGTGTTGAGGATCTTGCCGCGCAGCGGCAGGATGGc                                                                            >  6:430153/1‑149 (MQ=255)                                                                   cTTGCCATAGCGCAGCTGCGCAAGGTCGGTCGCACCCG‑GA‑TCGACGCCAATGGCCA‑CGGCAATGTTGTGCACTTCCTGGCTGGCCAGGACTTCGCCCCCGTCCA‑‑CTTCC‑‑CAGGTGTTGAGGATCTTGCCGCGCAGCGGCAGGATGGCCt                                                                          >  7:35248/1‑149 (MQ=255)                                                                      gCCATAGCGCAGCTGCGCAAGGTCGGTTGCACCCG‑GA‑TCGACGCCAATGGCCA‑CGGCAATGTTGTGCACTTCCTGGCTGGCCAGGACTTCGCCCCCGTCCA‑‑CTTCC‑‑CAGGTGTTGAGGATCTTGCCGCGCAGCGGCAGGATGGCCt                                                                          <  7:251526/146‑1 (MQ=255)                                                                      gCCATAGCGCAGCTGCGCAAGGTCGGTTGCACCCG‑GA‑TCGACGCCAATGGCCA‑CGGCAATGTTGTGCACTTCCTGGCTGGCCAGGACTTCGCCCCCGTCCA‑‑CTTCC‑‑CAGGTGTTGAGGATCTTGCCGCGCAGCGGCAGGATGGCCTgg                                                                        >  7:348423/1‑148 (MQ=255)                                                                       ccATAGCGCAGCTGCGCAAGGTCGGTTGCACCCG‑GA‑TCGACGCCAATGGCCA‑CGGCAATGTTGTGCACTTCCTGGCTGGCCAGGACTTCGCCCCCGTCCA‑‑CTTCC‑‑CAGGTGTTGAGGATCTTGCCGCGCAGCGGCAGGATGGCCTGGa                                                                       <  6:168783/148‑1 (MQ=255)                                                                        cATAGCGCAGCTGCGCAAGGTCGGTTGCACCCG‑GA‑TCGACGCCAATGGCCA‑CGGCAATGTTGTGCACTTCCTGGCTGGCCAGGACTTCGCCCCCGTCCA‑‑CTTCC‑‑CAGGTGTTGAGGATCTTGCCGCGCAGCGGCAGGATGGCCt                                                                          <  2:369631/144‑1 (MQ=255)                                                                        cATAGCGCAGCTGCGCAAGGTCGGTTGCACCCG‑GA‑TCGACGCCAATGGCCA‑CGGCAATGTTGTGCACTTCCTGGCTGGCCAGGACTTCGCCCCCGTCCA‑‑CTTCC‑‑CAGGTGTTGAGGATCTTGCCGCGCAGCGGCAGGATGGCCt                                                                          >  1:369631/1‑144 (MQ=255)                                                                        cATAGCGCAGCTGCGCAAGGTCGGTTGCACCCG‑GA‑TCGACGCCAATGGCCA‑CGGCAATGTTGTGCACTTCCTGGCTGGCCAGGACTTCGCCCCCGTCCA‑‑CTTCC‑‑CAGGTGTTGAGGATCTTGCCGCGCAGCGGCAGGATGGCCTGGAAt                                                                     <  5:430153/149‑1 (MQ=255)                                                                        cATAGCGCAGCTGCGCAAGGTCGGTTGCACCCG‑GA‑TCGACGCCAATGGCCA‑CGGCAATGTTGTGCACTTCCTGGCTGGCCAGGACTTCGCCCCCGTCCA‑‑CTTCC‑‑CAGGTGTTGAGGATCTTGCCGCGCAGCGGCAGGATGGCCTGGAAt                                                                     >  6:327473/1‑149 (MQ=255)                                                                         aTAGCGCAGCTGCGCAAGGTCGGTTGCACCCG‑GA‑TCGACGCCAATGGCCA‑CGGCAATGTTGTGCACTTCCTGGCTGGCCAGGACTTCGCCCCCGTCCA‑‑CTTCC‑‑CAGGTGTTGAGGATCTTGCCGCGCAGCGGCAGGATGGCCTGGAAtt                                                                    >  3:50339/1‑149 (MQ=255)                                                                         aTAGCGCAGCTGCGCAAGGTCGGTCGCACCCG‑GA‑TCGACGCCAATGGCCA‑CGGCAATGTTGTGCACTTCCTGGCTGGCCAGGACTTCGCCCCCGTCCA‑‑CTTCC‑‑CAGGTGTTGAGGATCTTGCCGCGCAGCGGCAGGATGGCCTGGAAtt                                                                    <  8:35248/149‑1 (MQ=255)                                                                            gcgcAGCTGCGCAAGGTCGGTTGCACCCG‑GA‑TCGACGCCAATGGCCA‑CGGCAATGTTGTGCACTTCCTGGCTGGCCAGGACTTCGCCCCCGTCCA‑‑CTTCC‑‑CAGGTGTTGAGGATCTTGCCGCGCAGCGGCAGGATGGCCTGGAATTcc                                                                  <  4:361162/148‑1 (MQ=255)                                                                            gcgcAGCTGCGCAAGGTCGGTTGCACCCG‑GA‑TCGACGCCAATGGCCA‑CGGCAATGTTGTGCACTTCCTGGCTGGCCAGGACTTCGCCCCCGTCCA‑‑CTTCC‑‑CAGGTGTTGAGGATCTTGCCGCGCAGCGGCAGGATGGCCTGGAATTcc                                                                  <  4:367225/148‑1 (MQ=255)                                                                             cgcAGCTGCGCAAGGTCGGTTGCACCCG‑GA‑TCGACGCCAATGGCCA‑CGGCAATGTTGTGCACTTCCTGGCTGGCCAGGACTTCGCCCCCGTCCA‑‑CTTCC‑‑CAGGTGTTGAGGATCTTGCCGCGCAGCGGCAGGATGGCCTGGAATTCCtt                                                                <  5:75031/149‑1 (MQ=255)                                                                              gcAGCTGCGCAAGGTCGGTTGCACCCG‑GA‑TCGACGCCAATGGCCA‑CGGCAATGTTGTGCACTTCCTGGCTGGCCAGGACTTCGCCCCCGTCCA‑‑Ctt                                                                                                                      <  8:158122/96‑1 (MQ=255)                                                                              gcAGCTGCGCAAGGTCGGTTGCACCCG‑GA‑TCGACGCCAATGGCCA‑CGGCAATGTTGTGCACTTCCTGGCTGGCCAGGACTTCGCCCCCGTCCA‑‑Ctt                                                                                                                      >  7:158122/1‑96 (MQ=255)                                                                               cAGCTGCGCAAGGTCGGTTGCACCCG‑GA‑TCGCCGCCAATGGCCA‑CGGCAATGTTGTGCACTTCCTGGCTGGCCAGGACTTCGCCCCCGTCCA‑‑CTTCC‑‑CAGGTGTTGAGGATCTTGCCGCGCACCGGCAGGATGGCCTGGAATTCCTTGt                                                              >  3:284801/1‑149 (MQ=255)                                                                                        aaGGTCGGTTGCACCCG‑GA‑TCGACGCCAATGGCCA‑CGGCAATGTTGTGCACTTCCTGGCTGGCCAGGACTTCGCCCCCGTCCA‑‑CTTCC‑‑CAGGTGTTGAGGATCTTGCCGCGCAGCGGCAGGATGGCCTGGAATTCCTTGTCCCGCGCCt                                                     >  2:18604/1‑149 (MQ=255)                                                                                        aaGGTCGGTTGCACCCG‑GA‑TCGACGCCAATGGCCA‑CGGCAATGTTGTGCACTTCCTGGCTGGCCAGGACTTCGCCCCCGTCCA‑‑CTTCC‑‑CAGGTGTTGAGGATCTTGCCGCGCAGCGGCAGGATGGCCTGGAATTCCTTGTCCCGCGCCt                                                     >  2:274058/1‑149 (MQ=255)                                                                                          ggTCGGTTGCACCCG‑GA‑TCGACGCCAATGGCCA‑CGGCAATGTTGTGCACTTCCTGGCTGGCCAGGACTTCGCCCCCGTCCA‑‑CTTCC‑‑CAGGTGTTGAGGATCTTGCCGCGCAGCGGCAGGATGGCCTGGAATTCCTTGTCCCGCGCCTGc                                                   <  4:50339/149‑1 (MQ=255)                                                                                           gTCGGTTGCACCCG‑GA‑TCGACGCCAATGGCCA‑CGGCAATGTTGTGCACTTCCTGGCTGGCCAGGACTTCGCCCCCGTCCA‑‑CTTCC‑‑CAGGTGTTGAGGATCTTGCCGCGCAGCGGCAGGATGGCCTGGAATTCCTTGTCCCGCGCCTGCt                                                  <  4:211701/149‑1 (MQ=255)                                                                                               gTTGCACCCG‑GA‑TCGACGCCAATGGCCA‑CGGCAATGTTGTGCACTTCCTGGCTGGCCAGGACTTCGCCCCCGTCCA‑‑CTTCC‑‑CAGGTGTTGAGGATCTTGCCGCGCAGCGGCAGGATGGCCTGGAATTCCTTGTCCCGCGCCTGCTTGGc                                              >  7:255809/1‑149 (MQ=255)                                                                                               gTTGCACCCG‑GA‑TCGACGCCAATGGCCA‑CGGCAATGTTGTGCACTTCCTGGCTGGCCAGGACTTCGCCCCCGTCCA‑‑CTTCC‑‑CAGGTGTTGAGGATCTTGCCGCGCAGCGGCAGGATGGCCTGGAATTCCTTGTCCCGCGCCTGCTTGGc                                              <  2:353818/149‑1 (MQ=255)                                                                                                ttGCACCCG‑GA‑TCGACGCCAATGGCCA‑CGGCAATGTTGTGCACTTCCTGGCTGGCCAGGACTTCGCCCCCGTCCA‑‑CTTCC‑‑CAGGTGTTGAGGATCTTGCCGCGCAGCGGCAGGATGGCCTGGAATTCCTTGTCCCGCGCCTGCTTGGcc                                             <  3:27280/149‑1 (MQ=255)                                                                                                  gCACCCG‑GA‑TCGACGCCAATGGCCA‑CGGCAATGTTGTGCACTTCCTGGCTGGCCAGGACTTCGCCCCCGTCCA‑‑CTTCC‑‑CAGGTGTTGAGGATCTTGCCGCGCAGCGGCAGGATGGCCTgg                                                                        >  1:233290/1‑120 (MQ=255)                                                                                                  gCACCCG‑GA‑TCGACGCCAATGGCCA‑CGGCAATGTTGTGCACTTCCTGGCTGGCCAGGACTTCGCCCCCGTCCA‑‑CTTCC‑‑CAGGTGTTGAGGATCTTGCCGCGCAGCGGCAGGATGGCCTgg                                                                        <  2:233290/120‑1 (MQ=255)                                                                                                   cACCCG‑GA‑TCGACGCCAATGGCCA‑CGGCAATGTTGTGCACTTCCTGGCTGGCCAGGACTTCGCCCCCGTCCA‑‑CTTCC‑‑CAGGTGTTGAGGATCTTGCCGCGCAGCGGCAGGATGGCCt                                                                          <  5:146728/117‑1 (MQ=255)                                                                                                   cACCCG‑GA‑TCGACGCCAATGGCCA‑CGGCAATGTTGTGCACTTCCTGGCTGGCCAGGACTTCGCCCCCGTCCA‑‑CTTCC‑‑CAGGTGTTGAGGATCTTGCCGCGCAGCGGCAGGATGGCCt                                                                          >  6:146728/1‑117 (MQ=255)                                                                                                     cccG‑GA‑TCGACGCCAATGGCCA‑CGGCAATGTTGTGCACTTCCTGGCTGGCCAGGACTTCGCCCCCGTCCA‑‑CTTCC‑‑CAGGTGTTGAGGATCTTGCCGCGCAGCGGCAGGATGGCCTGGAATTCCTTGTCCCGCGCCTGCTTGGCCGAac                                         >  4:57213/1‑148 (MQ=255)                                                                                                      ccG‑GA‑TCGACGCCAATGGCCA‑CGGCAATGTTGTGCACTTCCTGGCTGGCCAGGACTTCGCCCCCGTCCA‑‑CTTCC‑‑CAGGTGTTGAGGATCTTGCCGCGCAGCGGCAGGATGGCCTGGAATTCCTTGTCCCGCGCCTGCTTGGCCGAacca                                       >  4:111207/1‑149 (MQ=255)                                                                                                       cG‑GA‑TCGACGCCAATGGCCA‑CGGCAATGTTGTGCACTTCCTGGCTGGCCAGGACTTCGCCCCCGTCCA‑‑CTTCC‑‑CAGGTGTTGAGGATCTTGCCGCGCAGCGGCAGGATGGCCTGGAATTCCTTGTCCCGCGCCTGCTTGGCCGAaccac                                      <  4:152497/149‑1 (MQ=255)                                                                                                        g‑gA‑TCGACGCCAATGGCCA‑CGGCAATGTTGTGCACTTCCTGGCTGGCCAGGACTTCGCCCCCGTCCA‑‑CTTCC‑‑CAGGTGTTGAGGATCTTGCCGCGCAGCGGCAGGATGGCCTGGAATTCCTTGTCCCGCGCCTGCTTGGCCGAaccacc                                     <  1:17097/149‑1 (MQ=255)                                                                                                           a‑TCGACGCCAATGGCCA‑CGGCAATGTTGTGCACTTCCTGGCTGGCCAGGACTTCGCCCCCGTCCA‑‑CTTCC‑‑CAGGTGTTGAGGATCTTGCCGCGCAGCGGCAGGATGGCCTGGAATTCCTTGTCCCGCGCCTGCTTGGCCGAACCACCgg                                   <  3:55327/149‑1 (MQ=255)                                                                                                             tCGACGCCAATGGCCA‑CGGCAATGTTGTGCACTTCCTGGCTGGCCAGGACTTCGCCCCCGTCCA‑‑CTTCC‑‑CAGGTGTTGAGGATCTTGCCGCGCAGCGGCAGGATGGCCTGGAATTCCTTGTCCCGCGCCTGCTTGGCCGAACCACCGGc                                  <  8:255809/149‑1 (MQ=255)                                                                                                             tCGACGCCAATGGCCA‑CGGCAATGTTGTGCACTTCCTGGCTGGCCAGGACTTCGCCCCCGTCCA‑‑CTTCC‑‑CAGGTGTTGAGGATCTTGCCGCGCAGCGGCAGGATGGCCTGGAATTCCTTGTCCCGCGCCTGCTTGGCCGAACCACCGGc                                  <  1:278117/149‑1 (MQ=255)                                                                                                               gACGCCAATGGCCA‑CGGCAATGTTGTGCACTTCCTGGCTGGCCAGGACTTCGCCCCCGTCCA‑‑CTTCC‑‑CAGGTGTTGAGGATCTTGCCGCGCAGCGGCAGGATGGCCTGGAATTCCTTGTCCCGCGCCTGCTTGGCCGAACCACCGGcc                                 <  4:263793/148‑1 (MQ=255)                                                                                                               gACGCCAATGGCCA‑CGGCAATGTTGTGCACTTCCTGGCTGGCCAGGACTTCGCCCCCGTCCA‑‑CTTCC‑‑CAGGTGTTGAGGATCTTGCCGCGCAGCGGCAGGATGGCCTGGAATTCCTTGTCCCGCGCCTGCTTGGCCGAACCACCGGcc                                 <  3:111207/148‑1 (MQ=255)                                                                                                                aCGCCAATGGCCA‑CGGCAATGTTGTGCACTTCCTGGCTGGCCAGGACTTCGCCCCCGTCCA‑‑CTTCC‑‑CAGGTGTTGAGGATCTTGCCGCGCAGCGGCAGGATGGCCTGGAATTCCTTGTCCCGCGCCTGCTTGGCCGAACCACCGGcc                                 >  4:35186/1‑147 (MQ=255)                                                                                                                  gCCAATGGCCA‑CGGCAATGTTGTGCACTTCCTGGCTGGCCAGGACTTCGCCCCCGTCCA‑‑CTTCCGGCAGGTGTTGAGGATCTTGCCGCGCAGCGGCAGGATGGCCTGGAATTCCTTGTCCCGCGCCTGCTTGGCCGAACCACCGGcc                                 >  1:348781/1‑147 (MQ=255)                                                                                                                    cAATGGCCA‑CGGCAATGTTGTGCACTTCCTGGCTGGCCAGGACTTCGCCCCCGTCCA‑‑CTTCC‑‑CAGGTGTTGAGGATCTTGCCGCGCAGCGGCAGGATGGCCTGGAATTCCTTGTCCCGCGCCTGCTTGGCCGAACCACCGGCCGAGTCa                           >  6:422151/1‑149 (MQ=255)                                                                                                                    cAATGGCCA‑CGGCAATGTTGTGCACTTCCTGGCTGGCCAGGACTTCGCCCCCGTCCA‑‑CTTCC‑‑CAGGTGTTGAGGATCTTGCCGCGCAGCGGCAGGATGGCCTGGAATTCCTTGTCCCGCGCCTGCTTGGCCGAACCACCGGCCGAGTCa                           >  8:252663/1‑149 (MQ=255)                                                                                                                       tGGCCA‑GGTCAATGTTGTGCACTTCCTGGCTGGCCAGGACTTCGCCCCCGTCCA‑‑CTTCC‑‑CAGGTGTTGAGGATCTTGCCGCGCAGCGGCAGGATGGCCTGGAATTCCTTGTCCCGCGCCTGCTTGGCCGAACCACCGGCCGAGTCAccc                        <  3:56889/149‑1 (MQ=255)                                                                                                                          ccA‑CGGCAATGTTGTGCACTTCCTGGCTGGCCAGGACTTCGCCCCCGTCCA‑‑CTTCC‑‑CAGGTGTTGAGGATCTTGCCGCGCAGCGGCAGGATGGCCTGGAATTCCTTGTCCCGCGCCTGCTTGGCCGAACCACCGGCCGAGTCACCCTcg                     <  3:35186/149‑1 (MQ=255)                                                                                                                           cA‑CGGCAATTTTTTGCAATTCCTGGCTGGCCAGGAATTCGCCCCCGTCCACTCTTCC‑‑CCGGACTTCTCCCATTTGCCGCGCAGCGGCAGGATGGCCTGGAATTCCTCGTCCCGCGCCTGCTTGGCCGAACCACCGGCCGAGTCACCCTc                      <  1:141532/149‑1 (MQ=255)                                                                                                                               ggCAATGTTGTGCACTTCCTGGCTGGCCAGGACTTCGCCCCCGTCCA‑‑CTTCC‑‑CAGGTGTTGAGGATCTTGCCGCGCAGCGGCAGGATGGCCTGGAATTCCTTGTCCCGCGCCTGCTTGGCCGAACCACCGGCCGAGTCACCCTcgacga                 <  2:155267/149‑1 (MQ=255)                                                                                                                                 cAATGTTGTGCACTTCCTGGCTGGCCAGGACTTCGCCCCCGTCCA‑‑CTTCC‑‑CAGGTGTTGAGGATCTTGCCGCGCAGCGGCAGGATGGCCTGGAATTCCTTGTCCCGCGCCTGCTTGGCCGAACCACCGGCCGAGTCACCCTCGACGAgg               <  5:422151/149‑1 (MQ=255)                                                                                                                                  aaTGTTGTGCACTTCCTGGCTGGCCAGGACTTCGCCCCCGTCCA‑‑CTTCC‑‑CAGGTGTTGAGGATCTTGCCGCGCAGCGGCAGGATGGCCTGGAATTCCTTGTCCCGCGCCTGCTTGGCCGAACCACCGGCCGAGTCACCCTCGACGAgg               >  4:94219/1‑148 (MQ=255)                                                                                                                                    tgttgtGCACTTCCTGGCTGGCCAGGACTTCGCCCCCGTCCA‑‑CTTCC‑‑CAGGTGTTGAGGATCTTGCCGCGCAGCGGCAGGATGGCCTGGAATTCCTTGTCCCGCGCCTGCTTGGCCGAACCACCGGCCGAGTCACCCTCGACGAGGaac            <  2:348781/149‑1 (MQ=255)                                                                                                                                    tgttgtGCACTTCCTGGCTGGCCAGGACTTCGCCCCCGTCCA‑‑CTTCC‑‑CAGGTGTTGAGGATCTTGCCGCGCAGCGGCAGGATGGCCTGGAATTCCTTGTCCCGCGCCTGCTTGGCCGAACCACCGGCCGAGTCACCCTCGACGAGGaa             >  1:174688/1‑148 (MQ=255)                                                                                                                                          gCACTTCCTGGCTGGCCAGGACTTCGCCCCCGTCCA‑‑CTTCC‑‑CAGGTGTTGAGGATCTTGCCGCGCAGCGGCAGGATGGCCTGGAATTCCTTGTCCCGCGCCTGCTTGGCCGAACCACCGGCCGAGTCACCCTCGACGAGGAACAACTCg      >  8:264801/1‑149 (MQ=255)                                                                                                                                           cACTTCCTGGCTGGCCAGGACTTCGCCCCCGTCCA‑‑CTTCC‑‑CAGGTGTTGAGGATCTTGCCGCGCAGCGGCAGGATGGCCTGGAATTCCTTGTCCCGCGCCTGCTTGGCCGAACCACCGGCCGAGTCACCCTCGACGAGGAACAACTCg      <  7:264801/148‑1 (MQ=255)                                                                                                                                             agtCCTGGCCAGCCAGGACTTCGCCCCCGTCCA‑‑CTTCC‑‑CAGGTGTTGAGGATCTTGCCGCGCAGCGGCAGGATGGCCTGGAATTCCTTGTCCCGCGCCTGCTTGGCCGAACCACCGGCCGAGTCACCCTCGACGAGGAACAACTCGGCa   >  1:218697/3‑149 (MQ=255)                                                                                                                                             agtCCTGGCCAGCCAGGACTTCGCCCCCGTCCA‑‑CTTCC‑‑CAGGTGTTGAGGATCTTGCCGCGCAGCGGCAGGATGGCCTGGAATTCCTTGTCCCGCGCCTGCTTGGCCGAACCACCGGCCGAGTCACCCTCGACGAGGAACAACTCGGCa   >  1:153772/3‑149 (MQ=255)                                                                                                                                              ttCCTGGCTGGCCAGGACTTCGCCCCCGTCCA‑‑CTTCC‑‑CAGGTGTTGAGGATCTTGCCGCGCAGCGGCAGGATGGCCTGGAATTCCTTGTCCCGCGCCTGCTTGGCCGAACCACCGGCCGAGTCACCCTCGACGAGGAACAACTCGGCAc  >  8:55211/1‑149 (MQ=255)                                                                                                                                              gtCCTGGCCAGCCAGGACTTCGCCCCCGTCCA‑‑CTTCC‑‑CAGGTGTTGAGGATCTTGCCGCGCAGTGGCAGGATGGCCTGGAATTCCTTGTCCCGCGCCTGCTTGGCCGAACCACCGGCCGAGTCACCCTCGACGAGGAACAACTCGGCAc  >  6:203607/2‑149 (MQ=255)                                                                                                                                                      |                                                                                                                                                 GAACAGCGCACACAGCAATGTGGCGATGTGCAGACCATCGGAGTCAGCGTCGGCGAGGATGCAGATCTTGCCATAGCGCAGCTGCGCAAGGTCGGTTGCACCCG‑GA‑TCGACGCCAATGGCCA‑CGGCAATGTTGTGCACTTCCTGGCTGGCCAGGACTTCGCCCCCGTCCA‑‑CTTCC‑‑CAGGTGTTGAGGATCTTGCCGCGCAGCGGCAGGATGGCCTGGAATTCCTTGTCCCGCGCCTGCTTGGCCGAACCACCGGCCGAGTCACCCTCGACGAGGAACAACTCGGCAC  >  NC_002947/5587618‑5587904

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