Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,758,372	T→A	7.1%	*264C (TGA→TGT)	PP_2412 ←	ParA family protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,758,372	0	T	A	7.1%	210.5 / 3.0	70	*264C (TGA→TGT)	PP_2412	ParA family protein
Reads supporting (aligned to +/- strand):  ref base T (33/32);  new base A (3/2);  total (36/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.91e-01

GCAGTCGGGCCGCGCCATCGTTATTGTGCTTGCCAGTAGATAGCCGAGCGCTGCAAGGCGCTTTGTCGAGGGCGACCTCGGTATACCTGATTCGACCGGTGCAGCAGGTTACAGCCGGAGTGCCCCCAACGGCCTGTCAGGCCGCACGCTCAAGCACATCCAGCAAATCGATGAATTCCAGGGCCAGTTTAT‑GCCGCGGCGCCAAATGCACCAATGGCACCGAGGCATGGTGCGATTCGCGCATCTTGATCGAACTGCTCAGGTACACCGGCAGCACCGGCA  >  NC_002947/2758236‑2758517                                                                                                                                         |                                                                                                                                                   gcAGTCGGGCCGCGCCATCGTTATTGTGCTTGCCAGTAGATAGCCGAGCGCTGCAAGGCGCTTTGTCGAGGGCGACCTCGGTATACCTGATTCGACCGGTGCAGCAGGTTACAGCCGGAGTGCCCCCAACGGCCTGTCAGGCCGCACGc                                                                                                                                        >  5:151678/1‑149 (MQ=255)  cAGTCGGGCCGCGCCATCGTTATTGTGCTTGCCAGTAGATAGCCGAGCGCTGCAAGGCGCTTTGTCGAGGGCGACCTCGGTATACCTGATTCGACCGGTGCAGCAGGTTACAGCCGGAGTGCCCCCAACGGCCTGTCAGGCCGCACGCt                                                                                                                                       <  5:259534/149‑1 (MQ=255)   aGTCGGGCCGCGCCATCGTTATTGTGCTTGCCAGTAGATAGCCGAGCGCTGCAAGGCGCTTTGTCGAGGGCGACCTCGGTATACCTGATTCGACCGGTGCAGCAGGTTACAGCCGGAGTGCCCCCAACGGCCTGTCAGGCCGCACGCTc                                                                                                                                      <  6:151678/149‑1 (MQ=255)           cgcgCCATCGTTATTGTGCTTGCCAGTAGATAGCCGAGCGCTGCAAGGCGCTTTGTCGAGGGCGACCTCGGTATACCTGATTCGACCGGTGCAGCAGGTTACAGCCGGAGTGCCCCCAACGGCCTGTCAGGCCGCACGCTCAAGCACAt                                                                                                                              >  2:120633/1‑149 (MQ=255)                cATCGTTATTGTGCTTGCCAGTAGATAGCCGAGCGCTGCAAGGCGCTTTGTCGAGGGCGACCTCGGTATACCTGATTCGACCGGTGCAGCAGGTTACAGCCGGAGTGCCCCCAACGGCCTGTCAGGCCGCACGCTCAAGCACATCCAg                                                                                                                          >  3:170289/1‑148 (MQ=255)                        ttGTGCTTGCCAGTAGATAGCCGAGCGCTGCAAGGCGCTTTGTCGAGGGCGACCTCGGTATACCTGATTCGACCGGTGCAGCAGGTTACAGCCGGAGTGCCCCCAACGGCCTGTCAGGCCGCACGCTCAAGCACATCCAGCAAATCGAt                                                                                                                 <  3:195696/149‑1 (MQ=255)                        ttGTGCTTGCCAGTAGATAGCCGAGCGCTGCAAGGCGCTTTGTCGAGGGCGACCTCGGTATACCTGATTCGACCGGTGCAGCAGGTTACAGCCGGAGTGCCCCCAACGGCCTGTCAGGCCGCACGCTCAAGCACATCCAGCAAATCGAt                                                                                                                 <  6:47852/149‑1 (MQ=255)                                   aGTAGATAGCCGAGCGCTGCAAGGCGCTTTGTCGAGGGCGACCTCGGTATACCTGATTCGACCGGTGCAGCAGGTTACAGCCGGAGTGCCCCCAACGGCCTGTCAGGCCGCACGCTCAAGCACATCCAGCAAATCGATGAATTCCAggg                                                                                                      <  4:170289/149‑1 (MQ=255)                                            ccGAGCGCTGCAAGGCGCTTTGTCGAGGGCGCCCTCGGTATACCTGATTCGACCGTTGCAGCAGGTTACAGCCGGAGTGCGCCAAACGGCCTGTCAGGCCGCACGCTCAAGCACATCCAGCAAATCCATGAATTCCAGGGCCAGTTTAt                                                                                             >  7:240342/1‑149 (MQ=255)                                            ccGAGCGCTGCAAGGCGCTTTGTCGAGGGCGACCTCGGTATACCTGATTCGACCGGTGCAGCAGGTTACAGCCGGAGTGCCCCCAACGGCCTGTCAGGCCGCACGCTCAAGCACATCCAGCAAATCGATGAATTCCAGGGCCAGTTTAt                                                                                             >  5:125816/1‑149 (MQ=255)                                                gcgcTGCAAGGCGCTTTGTCGAGGGCGACCTCGGTATACCTGATTCGACCGGTGCAGCAGGTTACAGCCGGAGTGCCCCCAACGGCCTGTCAGGCCGCACGCTCAAGcaca                                                                                                                               >  3:73313/1‑111 (MQ=255)                                                gcgcTGCAAGGCGCTTTGTCGAGGGCGACCTCGGTATACCTGATTCGACCGGTGCAGCAGGTTACAGCCGGAGTGCCCCCAACGGCCTGTCAGGCCGCACGCTCAAGcaca                                                                                                                               <  4:73313/111‑1 (MQ=255)                                                gcgcTGCAAGGCGCTTTGTCGAGGGCGACCTCGGTATACCTGATTCGACCGGTGCAGCAGGTTACAGCCGGAGTGCCCCCAACGGCCTGTCAGGCCGCACGCTCAAGCACATCCAGCAAATCGAt                                                                                                                 >  3:191625/1‑125 (MQ=255)                                                gcgcTGCAAGGCGCTTTGTCGAGGGCGACCTCGGTATACCTGATTCGACCGGTGCAGCAGGTTACAGCCGGAGTGCCCCCAACGGCCTGTCAGGCCGCACGCTCAAGCACATCCAGCAAATCGAt                                                                                                                 <  4:191625/125‑1 (MQ=255)                                                    tGCAAGGCGCCTTGTCGAGGGCGACCTCGGTATACCTGATTCGACCGGTGCAGCAGGTTACAGCCGGAGTGCCCCCAACGGCCTGTCAGGCCGCACGCTCAAGCACATCCAGCAAATCGATGAATTCCAGGGCCAGTTTAT‑GCCgcggc                                                                                    >  1:312090/1‑149 (MQ=255)                                                      cAAGGCGCTTTGTCGAGGGCGACCTCGGTATACCTGATTCGACCGGTGCAGCAGGTTACAGCCGGAGTGCCCCCAACGGCCTGTCAGGCCGCACGCTCAAGCACATCCAGCAAATCGATGAATTCCAgg                                                                                                       >  5:182679/1‑129 (MQ=255)                                                      cAAGGCGCTTTGTCGAGGGCGACCTCGGTATACCTGATTCGACCGGTGCAGCAGGTTACAGCCGGAGTGCCCCCAACGGCCTGTCAGGCCGCACGCTCAAGCACATCCAGCAAATCGATGAATTCCAgg                                                                                                       <  6:182679/129‑1 (MQ=255)                                                      cAAGGCGCCTTGTCGAGGGCGACCTCGGTATACCTGATTCGACCGGTGCAGCAGGTTACAGCCGGAGTGCCCCCAACGGCCTGTCAGGCCGCACGCTCAAGCACATCCAGCAAATCGATGAATTCCAGGGCCAGTTTAT‑GCCgcggc                                                                                    <  2:312090/147‑1 (MQ=255)                                                       aaGGCGCTTTGTCGAGGGCGACCTCGGTATACCTGATTCGACCGGTGCAGCAGGTTACAGCCGGAGTGCCCCCAACGGCCTGTCAGGCCGCACGCTCAAGCACATCCTGCAAATCGATGAATTCCAGGGCCAGTTTAT‑GCCTCGgtgc                                                                                  >  5:308395/1‑148 (MQ=255)                                                        aGGCGCTTTGTCGAGGGCGACCTCGGTATACCTGATTCGACCGGTGCAGCAGGTTACAGCCGGAGTGCCCCCAACGGCCTGTCAGGCCGCACGCTCAAGCACATCCAGCAAATCGATGAATTCCAGGGCCAGTTTAT‑GCCGCGGCGcc                                                                                 >  7:134169/1‑148 (MQ=255)                                                                    gAGGGCGACCTCGGTATACCTGATTCGACCGGTGCAGCAGGTTACAGCCGGAGTGCCCCCAACGGCCTGTCAGGCCGCACGCTCAAGCACATcc                                                                                                                            <  7:173358/94‑1 (MQ=255)                                                                    gAGGGCGACCTCGGTATACCTGATTCGACCGGTGCAGCAGGTTACAGCCGGAGTGCCCCCAACGGCCTGTCAGGCCGCACGCTCAAGCACATcc                                                                                                                            >  8:173358/1‑94 (MQ=255)                                                                    gAGGGCGACCTCGGTATACCTGATTCGACCGGTGCAGCAGGTTACAGCCGGAGTGCCCCCAACGGCCTGTCAGGCCGCACGCTCAAGCACATCCAGCAAATCGATGAATTCCAGGGCCAGTTTAT‑GCCGCGGCGCCAAATGCACCAAt                                                                     <  1:174624/148‑1 (MQ=255)                                                                        gCGACCTCGGTATACCTGATTCGACCGGTGCAGCAGGTTACAGCCGGAGTGCCCCCAACGGCCTGTCAGGCCGCACGCTCAAGCACATCCAGCAAATCGATGAATTCCAGGGCCAGTTTAT‑GCCGCGGCGCCAAATGCACCAATGGCAc                                                                >  2:278742/1‑149 (MQ=255)                                                                        gCGACCTCGGTATACCTGATTCGACCGGTGCAGCAGGTTACAGCCGGAGTGCCCCCAACGGCCTGTCAGGCCGCACGCTCAAGCACATCCAGCAAATCGATGAATTCCAGGGCCAGTTTAT‑GCCGCGGCGCCAAATGCACCAATGGCAc                                                                <  2:70631/149‑1 (MQ=255)                                                                        gCGACCTCGGTATACCTGATTCGACCGGTGCAGCAGGTTACAGCCGGAGTGCCCCCAACGGCCTGTCAGGCCGCACGCTCAAGCACATCCAGCAAATCGATGAATTCCAGGGCCAGTTTAT‑GCCGCGGCGCCAAATGCACCAATGGCAc                                                                <  8:134169/149‑1 (MQ=255)                                                                           aCCTCGGTATACCTGATTCGACCGGTGCAGCAGGTTACAGCCGGAGTGCCCCCAACGGCCTGTCAGGCCGCACGCTCAAGCACATCCAGCAAATCGATGAATTCCAGGGCCAGTTTAT‑GCCGCGGCGCCAAATGCACCAATGGCACCGa                                                             <  1:120633/149‑1 (MQ=255)                                                                             cTCGGTATACCTGATTCGACCGGTGCAGCAGGTTACAGCCGGAGTGCCCCCAACGGCCTGTCAGGCCGCACGCTCAAGCACATCCAGCAAATCGATGAATTCCAGGGCCAGTTTAT‑GCCGCGGCGCCAAATGCACCAATGGCACCGAgg                                                           >  3:40671/1‑149 (MQ=255)                                                                                     aCCTGATTCGACCGGTGCAGCAGGTTACAGCCGGAGTGCCCCCAACGGCCTGTCAGGCCGCACGCTCAAGCACATCCAGCAAATCGATGAATTCCAGGGCCAGTTTAT‑GCCGCGGCGCCAAATGCACCAATGGCACCGAGGCATGGTGc                                                   >  2:14553/1‑149 (MQ=255)                                                                                      ccTGATTCGACCGGTGCAGCAGGTTACAGCCGGAGTGCCCCCAACGGCCTGTCAGGCCGCACGCTCAAGCACATCCAGCAAATCGATGAATTCCAGGGCCAGTTTAT‑GCcg                                                                                        <  4:341556/111‑1 (MQ=255)                                                                                       cTGATTCGACCGGTGCAGCAGGTTACAGCCGGAGTGCCCCCAACGGCCTGTCAGGCCGCACGCTCAAGCACATCCAGCAAATCGATGAATTCCAGGGCCAGTTTAT‑GCCGCGGCGCCAAATGCACCAATGGCACCGAGGCATGGTGCGa                                                 >  4:145488/1‑149 (MQ=255)                                                                                       cTGATTCGACCGGTGCAGCAGGTTACAGCCGGAGTGCCCCCAACGGCCTGTCAGGCCGCACGCTCAAGCACATCCAGCAAATCGATGAATTCCAGGGCCAGTTTAT‑GCCGCGGCGCCAAATGCACCAATGGCACCGAGGCATGGTGCGa                                                 >  4:284848/1‑149 (MQ=255)                                                                                       cTGATTCGACCGGTGCAGCAGGTTACAGCCGGAGTGCCCCCAACGGCCTGTCAGGCCGCACGCTCAAGCACATCCAGCAAATCGATGAATCCCAGGGCCAGTTTAT‑GCCGCGGCGCCAAATGCACCAATGGCACCGAGGc                                                          <  7:249938/140‑1 (MQ=255)                                                                                       cTGATTCGACCGGTGCAGCAGGTTACAGCCGGAGTGCCCCCAACGGCCTGTCAGGCCGCACGCTCAAGCACATCCAGCAAATCGATGAATCCCAGGGCCAGTTTAT‑GCCGCGGCGCCAAATGCACCAATGGCACCGAGGc                                                          >  8:249938/1‑140 (MQ=255)                                                                                        tGATTCGACCGGTGCAGCAGGTTACAGCCGGAGTGCCCCCAACGGCCTGTCAGGCCGCACGCTCAAGCACATCCAGCAAATCGATGAATTCCAGGGCCAGTTTAT‑GCCGCGGCGCCAAATTCACAAATGGCACCGAGGCATGGTGCGAt                                                >  2:301809/1‑149 (MQ=255)                                                                                                  ggTGCAGCAGGTTACAGCCGGAGTGCCCCCAACGGCCTGTCAGGCCGCACGCTCAAGCACATCCAGCAAATCGATGAATTCCAGGGCCAGTTTAT‑GCCGCGGCGCCAAATGCACCAATGGCACCGAGGCATGGTGCGATTCGCGCATCt                                      <  6:125816/149‑1 (MQ=255)                                                                                                       agcagGTTACAGCCGGAGTGCCCCCAACGGCCTGTCAGGCCGCACGCTCAAGCACATCCAGCAAATCGATGAATTCCAGGGCCAGTTTAT‑GCCGCGGCGCCAAATGCAc                                                                         <  2:71142/109‑1 (MQ=255)                                                                                                       agcagGTTACAGCCGGAGTGCCCCCAACGGCCTGTCAGGCCGCACGCTCAAGCACATCCAGCAAATCGATGAATTCCAGGGCCAGTTTAT‑GCCGCGGCGCCAAATGCACCAATGGCACCGAGGCATGGTGCGATTCGCGCATCTTGATc                                 <  4:40671/149‑1 (MQ=255)                                                                                                          agGTTACAGCCGGAGTGCCCCCAACGGCCTGTCAGGCCGCACGCTCAAGCACATCCAGCAAATCGATGAATTCCAGGGCCAGTTTAT‑GCCGCGGCGCCAAATGCACCAATGGCACCGAGGCATg                                                       >  1:237559/1‑124 (MQ=255)                                                                                                          agGTTACAGCCGGAGTGCCCCCAACGGCCTGTCAGGCCGCACGCTCAAGCACATCCAGCAAATCGATGAATTCCAGGGCCAGTTTAT‑GCCGCGGCGCCAAATGCACCAATGGCACCGAGGCATg                                                       <  2:237559/124‑1 (MQ=255)                                                                                                               actgCCGGAGTGCCCCCAACGGCCTGTCAGGCCGCACGCTCAAGCACATCCAGCAAATCGATGAATTCCAGGGCCAGTTTATCTCCGCGGCGCCAAATGCACCAATGGCACCGAGGCATGGTGCGATTCGCGCATCTTGATCGAACTg                           <  1:278742/145‑1 (MQ=255)                                                                                                                 aGCCGGAGTGCCCCCAACGGCCTGTCAGGCCGCACGCTCAAGCACATCCAGCAAATCGATGAATTCCAGGGCCAGTTTAT‑GCCGCGGTGCCAAATGCACCAATGGCACCGAGGCATGGTGCGATTCGCGCATCTTGATCGAACTGCt                         >  3:190356/1‑147 (MQ=255)                                                                                                                 aGCCGGAGTGCCCCCAACGGCCTGTCAGGCCGCACGCTCAAGCACATCCAGCAAATCGATGAATTCCAGGGCCAGTTTAT‑GCCGCGGCGCCAAATGCACCAATGGCACCGAGGCATGGTGCGATTCGCGCATCTTGATCGaa                              >  8:125908/1‑142 (MQ=255)                                                                                                                 aGCCGGAGTGCCCCCAACGGCCTGTCAGGCCGCACGCTCAAGCACATCCAGCAAATCGATGAATTCCAGGGCCAGTTTAT‑GCCGCGGCGCCAAATGCACCAATGGCACCGAGGCATGGTGCGATTCGCGCATCTTGATCGaa                              <  7:125908/142‑1 (MQ=255)                                                                                                                  gCCGGAGTGCCCCCAACGGCCTGTCAGGCCGCACGCTCAAGCACATCCAGCAAATCGATGAATTCCAGGGCCAGTTTAT‑GCCGCGGCGCCAAATGCACCAATGGCACCGAGGCATGGTGCGATTCGCGCATCTTGATCGAACTGCTCAg                      <  3:145488/149‑1 (MQ=255)                                                                                                                      gAGTGCCCCCAACGGCCTGTCAGGCCGCACGCTCAAGCACATCCAGCAAATCGATGAATTCCAGGGCCAGTTTAT‑GCCGCGGCGCCAAATGCACCAATGGCACCGAGGCATGGTGCGATTCGCGCATCTTGATCGAACTGCTCAgg                     >  8:56297/1‑146 (MQ=255)                                                                                                                      gAGTGCCCCCAACGGCCTGTCAGGCCGCACGCTCAAGCACATCCAGCAAATCGATGAATTCCAGGGCCAGTTTAT‑GCCGCGGCGCCAAATGCACCAATGGCACCGAGGCATGGTGCGATTCGCGCATCTTGATCGAACTGCTCAgg                     <  7:56297/146‑1 (MQ=255)                                                                                                                           cccccAACGGCCTGTCAGGCCGCACGCTCAAGCACATCCAGCAAATCGATGAATTCCAGGGCCAGTTTAT‑GCCGCGGCGCCAAATGCACCAATGGCACCGAGGCATGGTGCGATTCGCGCATCt                                      >  8:311712/1‑124 (MQ=255)                                                                                                                           cccccAACGGCCTGTCAGGCCGCACGCTCAAGCACATCCAGCAAATCGATGAATTCCAGGGCCAGTTTAT‑GCCGCGGCGCCAAATGCACCAATGGCACCGAGGCATGGTGCGATTCGCGCATCt                                      <  7:311712/124‑1 (MQ=255)                                                                                                                           cccccAACGGCCTGTCAGGCCGCACGCTCAAGCACATCCAGCAAATCGATGAATTCCAGGGCCAGTTTAT‑GCCGCGGCGCCAAATGCACCAATGGCACCGAGGCATGGTGCGATTCGCGCATCTTGATCGAACTGCTCa                       <  8:246966/139‑1 (MQ=255)                                                                                                                           cccccAACGGCCTGTCAGGCCGCACGCTCAAGCACATCCAGCAAATCGATGAATTCCAGGGCCAGTTTAT‑GCCGCGGCGCCAAATGCACCAATGGCACCGAGGCATGGTGCGATTCGCGCATCTTGATCGAACTGCTCa                       >  7:246966/1‑139 (MQ=255)                                                                                                                            ccccAACGGCCTGTCAGGCCGCACGCTCAAGCACATCCAGCAAATCGATGAATTCCAGGGCCAGTTTAT‑GCCGCGGCGCCAAATGCACCAATGGCACCGAGGCATGGTGCGATTCGCGCATCTTGATCGAACTGCTCAGGTACACCgg             >  7:141125/1‑148 (MQ=255)                                                                                                                             gcgtgcGGCCTGACAGGCCGCACGCTCAAGCACATCCAGCAAATCGATGAATTCCAGGGCCAGTTTAT‑GCCGCGGCGCCAAATGCACCAATGGCACCGAGGCATGGTGCGATTCGCGCATCTTGATCGAACTGCTCAGGTACACCGgc            >  1:184219/6‑148 (MQ=255)                                                                                                                             gcgtgcGGCCTGACAGGCCGCACGCTCAAGCACATCCAGCAAATCGATGAATTCCAGGGCCAGTTTAT‑GCCGCGGCGCCAAATGCACCAATGGCACCGAGGCATGGTGCGATTCGCGCATCTTGATCGAACTGCTCAGGTACACCGgc            <  2:184219/143‑1 (MQ=255)                                                                                                                             cccAACGGCCTGTCAGGCCGCACGCTCAAGCACATCCAGCAAATCGATGAATTCCAGGGCCAGTTTAT‑GCCGCGGCGCCAAATGCACCAATGGCACCGAGGCATGGTGCGATTCGCGCATCt                                      >  8:179891/1‑122 (MQ=255)                                                                                                                             cccAACGGCCTGTCAGGCCGCACGCTCAAGCACATCCAGCAAATCGATGAATTCCAGGGCCAGTTTAT‑GCCGCGGCGCCAAATGCACCAATGGCACCGAGGCATGGTGCGATTCGCGCATCTTGATCGaa                              <  3:221055/130‑1 (MQ=255)                                                                                                                             cccAACGGCCTGTCAGGCCGCACGCTCAAGCACATCCAGCAAATCGATGAATTCCAGGGCCAGTTTAT‑GCCGCGGCGCCAAATGCACCAATGGCACCGAGGCATGGTGCGATTCGCGCATCTTGATCGaa                              >  1:288094/1‑130 (MQ=255)                                                                                                                             cccAACGGCCTGTCAGGCCGCACGCTCAAGCACATCCAGCAAATCGATGAATTCCAGGGCCAGTTTAT‑GCCGCGGCGCCAAATGCACCAATGGCACCGAGGCATGGTGCGATTCGCGCATCTTGATCGaa                              <  2:288094/130‑1 (MQ=255)                                                                                                                             cccAACGGCCTGTCAGGCCGCACGCTCAAGCACATCCAGCAAATCGATGAATTCCAGGGCCAGTTTAT‑GCCGCGGCGCCAAATGCACCAATGGCACCGAGGCATGGTGCGATTCGCGCATCTTGATCGaa                              >  4:221055/1‑130 (MQ=255)                                                                                                                             cccAACGGCCTGTCAGGCCGCACGCTCAAGCACATCCAGCAAATCGATGAATTCCAGGGCCAGTTTAT‑GCCGCGGCGCCAAATGCACCAATGGCACCGAGGCATGGTGCGATTCGCGCATCTTGATCGAACTGCTCAGGTACACCGgca           <  1:301809/149‑1 (MQ=255)                                                                                                                              ccAACGGCCTGTCAGGCCGCACGCTCAAGCACATCCAGCAAATCGATGAATTCCAGGGCCAGTTTAT‑GCCGCGGCGCCAAATGCACCAATGGCACCGAGGCATGGTGCGATTCGCGCATCTTGATCGaa                              <  4:5715/129‑1 (MQ=255)                                                                                                                              ccAACGGCCTGTCAGGCCGCACGCTCAAGCACATCCAGCAAATCGATGAATTCCAGGGCCAGTTTAT‑GCCGCGGCGCCAAATGCACCAATGGCACCGAGGCATGGTGCGATTCGCGCATCTTGATCGaa                              >  3:5715/1‑129 (MQ=255)                                                                                                                              ccAACGGCCTGTCAGGCCGCACGCTCAAGCACATCCAGCAAATCGATGAATTCCAGGGCCAGTTTAT‑GCCGCGGCGCCAAATGCACCAATGGCACCGAGGCATGGTGCGATTCGCGCATCTTGATCGAACTGCTCAGGTACACCGgcag          >  7:156920/1‑149 (MQ=255)                                                                                                                               gtgcGGCCTGACAGGCCGCACGCTCAAGCACATCCAGCAAATCGATGAATTCCAGGGCCAGTTTAT‑GCCGCGGCGCCAAATGCACCAATGGCACCGAGGCATGGTGCGATTCGCGCATCTTGATCGAACTGCTCAGGTACACCGgcagc         >  3:178400/4‑149 (MQ=255)                                                                                                                                tgcGGCCTGACAGGCCGCACGCTCAAGCACATCCAGCAAATCGATGAATTCCAGGGCCAGTTTAT‑GCCGCGGCGCCAAATGCACCAATGGCACCGAGGCATGGTGCGATTCGCGCATc                                       >  2:107566/3‑118 (MQ=255)                                                                                                                                tgcGGCCTGACAGGCCGCACGCTCAAGCACATCCAGCAAATCGATGAATTCCAGGGCCAGTTTAT‑GCCGCGGCGCCAAATGCACCAATGGCACCGAGGCATGGTGCGATTCGCGCATc                                       <  1:107566/116‑1 (MQ=255)                                                                                                                                aaCGGCCTGTCAGGCCGCACGCTCAAGCACATCCAGCAAATCGATGAATTCCAGGGCCAGTTTAT‑GCCGCGGCGCCAAATGCACCAATGGCACCGAGGCATGGTGCGATTCGCGCATCTTGATCGAACTGCTCa                       <  4:185212/134‑1 (MQ=255)                                                                                                                                aaCGGCCTGTCAGGCCGCACGCTCAAGCACATCCAGCAAATCGATGAATTCCAGGGCCAGTTTAT‑GCCGCGGCGCCAAATGCACCAATGGCACCGAGGCATGGTGCGATTCGCGCATCTTGATCGAACTGCTCa                       >  3:185212/1‑134 (MQ=255)                                                                                                                                     ccTGTCAGGCCGCACGCTCAAGCACATCCAGCAAATCGATGAATTCCAGGGCCAGTTTAT‑GCCGCGGCGCCAAATGCACCAATGGCACCGAGGCATGGTGCGATTCGCGCATCTTGATCGAACTGCTCAGGTACACCGGCAGCACCGGc   <  8:141125/149‑1 (MQ=255)                                                                                                                                      cTGTCAGGCCGCACGCTCAAGCACATCCAGCAAATCGATGAATTCCAGGGCCAGTTTAT‑GCCGCGGCGCCAAATGCACCAATGGCACCGAGGCATGGTGCGATTCGCGCATCTTGATCGAACTGCTCAGGTACACCGGCAGCACCGGCa  >  1:244323/1‑149 (MQ=255)                                                                                                                                         |                                                                                                                                                   GCAGTCGGGCCGCGCCATCGTTATTGTGCTTGCCAGTAGATAGCCGAGCGCTGCAAGGCGCTTTGTCGAGGGCGACCTCGGTATACCTGATTCGACCGGTGCAGCAGGTTACAGCCGGAGTGCCCCCAACGGCCTGTCAGGCCGCACGCTCAAGCACATCCAGCAAATCGATGAATTCCAGGGCCAGTTTAT‑GCCGCGGCGCCAAATGCACCAATGGCACCGAGGCATGGTGCGATTCGCGCATCTTGATCGAACTGCTCAGGTACACCGGCAGCACCGGCA  >  NC_002947/2758236‑2758517

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