Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	6,004,528	A→C	6.2%	T153P (ACC→CCC)	PP_5259 →	LysR family transcriptional regulator

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	6,004,528	0	A	C	6.2%	203.6 / 4.7	65	T153P (ACC→CCC)	PP_5259	LysR family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base A (31/30);  new base C (2/2);  total (33/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.11e-01

CGAGTTTTGGCGTGCGCTGGTTGATCCCGCACCTCAAGGGGTTTGGCAAACGGCACGGCAATATCCACCTGGA‑CATTCGCAACGAGATGGAGCCGTTCGCCCTGCTGCAAGGCTCGGCTGACGTGGTGTTCTTCTATGGCCAGGGCACCTGGCCGGGGGCCACTTGCGTGGAGTTGTTTCGCGAGGAGGTGGTGCCGGTGTGTGCGCCAGAGCTGCTGGCCGGGCGCAAGCTGGCCGGTGCCGAGGAGCTGGCCGAGCTGGTGCTGCTGCAAAGCACCTCGCGCCCGGAGGCCT  >  NC_002947/6004382‑6004675                                                                                                                                                    |                                                                                                                                                    cGAGTTTTGGCGTGCGCTGGTTGATCCCGCACCTCAAGGGGTTTGGCAAACGGCACGGCAATATCCACCTGGA‑CATTCGCAACGAGATGGAGCCGTTCGCCCTGCTGCAAGGCTCGGCTGACGTGGTGTTCTTCTATGGCCAGGGCAcc                                                                                                                                                   >  6:52549/1‑149 (MQ=255)    gTTATGGCGTGCCCAGGATGATCCCGCACCTCAAGGGGTATGGAAAACGGCACGGCAATATCCACCTGGA‑CATTCGAAACGAGATCGAGCCGTTCGACCTGCTGCAAGGCTCGGCTGACGTGGTGTTCTTCTAAGGCCAGGGCACCTgg                                                                                                                                                >  8:17626/1‑149 (MQ=255)        tGGCGTGCGCTGGTTGATCCCGCACCTCAAGGGGTTTGGCAAACGGCACGGCAATATCCACCTGGA‑CATTCGCAACGAGATGGAGCCGTTCGCCCTGCTGCAAGGCTCGGCTGACGTGGTGTTCTTCTATGGCCAGGTGCCCTGGcca                                                                                                                                             <  5:140144/148‑2 (MQ=255)          gCGTGCGGTTGTTTATCCAGCACCTCAAGGGGGTTGGGAAACGGCAGGGCAATATCAAGCTGGA‑GATTCGCAACGATATGGAGCCGTTCGCCCTGCTGCAAGGCTGGGCTGACGTGTTGTTCTTCTTTGGCCATGGCACCTGGCCggg                                                                                                                                           <  7:17626/148‑1 (MQ=255)          gCGTGCGCTGGTTGATCCCGCACCTCAGGGGGTTTGGCAAACGGCACGGCAATATCCACCTGGA‑CATTCGCAACGAGATGGAGCCGTTCGCCCTGCTGCAAGGCTCGGCTGACGTGGTGTTCTTCTATGGCCAGGGCACCTGGCCgggg                                                                                                                                          <  1:305573/149‑1 (MQ=255)          gCGTGCGCTGGTTGATCCCGCACCTCAAGGGGTTTGGCAAACGGCACGGCAATATCCACCTGGA‑CATTCGCAACGAGATGGAGCCGTTCGCCCTGCTGCAAGGCTCGGCTGACGTGGTGTTCTTCTATGGCCAGGGCACCTGGCCgggg                                                                                                                                          <  5:52549/149‑1 (MQ=255)           cGTGCGCTGGTTGATCCCGCACCTCAAGGGGTTTGGCAAACGGCACGGCAATATCCACCTGGA‑CATTCGCAACGAGATGGAGCCGTTCGCCCTGCTGCAAGGCTCGGCTGACGTGGTGTTCTTCTATGGCCAGGGCACCTGGCCggggg                                                                                                                                         <  3:74385/149‑1 (MQ=255)              gcgcTGGTTGATCCCGCACCTCAAGGGGTTTGGCAAACGCCACGGCAATATCCACCTGGA‑CATTCGCAACGAGATGGAGCCGTTCGCCCTGCTGCAAGGCTCGGCTGACGTGGTGTTCTTCTATGGCCAGGGCACCTGGCCGGGGGCCa                                                                                                                                      <  2:154458/149‑1 (MQ=255)               cgcTGGTTGATCCCGCACCTCAAGGGGTTTGGCAAACGGCACGGCAATATCCACCTGGA‑CATTCGCAACGAGATGGAGCCGTTCGCCCTGCTGCAAGGCTCGGCTGACGTGGTGTTCTTCTATGGCCAGGGCACCTGGCCGGGGGCCAc                                                                                                                                     <  7:106058/149‑1 (MQ=255)                     ttGATCCCGCACCTCAAGGGGTTTGGCAAACGGCACGGCAATATCCACCTGGA‑CATTCGCAACGAGATGGAGCCGTTCGCCCTGCTGCAAGGCTCGGCTGACGTGGTGTTCTTCTATGGCCAGGGCACCTGGCCGGGGGCCACTTGCGt                                                                                                                               <  8:70851/149‑1 (MQ=255)                       gATCCCGCACCTCAAGGGGTTTGGCAAACGGCACGGCAATATCCACCTGGA‑CATTCGCAACGAGATGGAGCCGTTCGCCCTGCTGCAAGGCTCGGCTGACGTGGTGTTCTTCTATGGCCAGGGCACCTGGCCGGGGGCCACTTGCGTg                                                                                                                              >  6:108819/1‑148 (MQ=255)                           ccGCACCTCAAGGGGTTTGGCAAACGGCACGGCAATATCCACCTGGA‑CATTCGCAACGAGATGGAGCCGTTCGCCCTGCTGCAAGGCTCGGCTGACGTGGTGTTCTTCTATGGCCAGGGCACCTGGCCGGGGGCCACTTGCGTGGAgtt                                                                                                                         >  2:22949/1‑149 (MQ=255)                             gCACCTCAAGGGGTTTGGCAAACGGCACGGCAATATCCACCTGGA‑CATTCGCAACGAGATGGAGCCGTTCGCCCTGCTGCAAGGCTCGGCTGACGTGGTGTTCTTCTATGGCCAGGGCACCTGGCCGGGGGCCACTTGCGTGGAgttgt                                                                                                                       <  1:22949/149‑1 (MQ=255)                                ccTCAAGGGGTTTGGCAAACGGCACGGCAATATCCACCTGGA‑CATTCGCAACGAGATGGAGCCGTTCGCCCTGCTGCAAGGCTCGGCTGACGTGGTGTTCTTCTATGGCCAGGGCACCTGGCCGGGGGCCACTTGCGTGGAGTTGTTTc                                                                                                                    <  8:150503/149‑1 (MQ=255)                                ccTCAAGGGGTTTGGCAAACGGCACGGCAATATCCACCTGGA‑CATTCGCAACGAGATGGAGCCGTTCGCCCTGCTGCAAGGCTCGGCTGACGTGGTGTTCTTCTATGGCCAGGGCACCTGGCCGGGGGCCACTTGAGTGGAGTTGTTTc                                                                                                                    >  1:265505/1‑149 (MQ=255)                                    aaGGGGTTTGGCAAACGGCACGGCAATATCCACCTGGA‑CATTCGCAACGAGATGGAGCCGTTCGCCCTGCTGCAAGGCTCGGCTGACGTGGTGTTCTTCTATGGCCAGGGCACCTGGCCGGGGGCCACTTGCGTGGAGTTGTTTcgc                                                                                                                  >  5:326914/1‑147 (MQ=255)                                    aaGGGGTTTGGCAAACGGCACGGCAATATCCACCTGGA‑CATTCGCAACGAGATGGAGCCGTTCGCCCTGCTGCAAGGCTCGGCTGACGTGGTGTTCTTCTATGGCCAGGGCACCTGGCCGGGGGCCACTTGCGTGGAGTTGTTTCGCga                                                                                                                <  2:276102/149‑1 (MQ=255)                                    aaGGGGTTTGGCAAACGGCACGGCAATATCCACCTGGA‑CATTCGCAACGAGATGGAGCCGTTCGCCCTGCTGCAAGGCTCGGCTGACGTGGTGTTCTTCTATGGCCAGGGCACCTGGCCGGGGGCCACTTGCGTGGAGTTGTTTCGCga                                                                                                                <  1:230128/149‑1 (MQ=255)                                            tGGCAAACGGCACGGCAATATCCACCTGGA‑CATTCGCAACGAGATGGAGCCGTTCGCCCTGCTGCAAGGCTCGGCTGACGTGGTGTTCTTCTATGGCCAGGGCACCTGGCCGGGGGCCACTTGCGTGGAGTTGTTTCGCGAGGAggtg                                                                                                         <  6:80797/148‑1 (MQ=255)                                               cAAACGGCACGGCAATATCCACCTGGA‑CATTCGCAACGAGATGGAGCCGTTCGCCCTGCTGCAAGGCGCGGCTGACGTGGTGTTCTTCTATGGCCAGGGCACCTGGCCGGGGGCCACTTGCGTGGAGTTGTTTCGCGAGGAGGTGGTGc                                                                                                     <  4:103905/149‑1 (MQ=255)                                                 aaCGGCACGGCAATATCCACCTGGATAATTCGCAACGAGATGGAGCCGTTCGCCCTGCTGCAAGGCTCGGCTGACGTGGTGTTCTTCTATGGCCAGGGCACCTGGCCGGGGGCCACTTGCGTGGa                                                                                                                            >  1:164185/1‑125 (MQ=255)                                                         ggcaATATCCACCTGGA‑CATTCGCAACGAGATGGAGCCGTTCGCCCTGCTGCAAGGCTCGGCTGACGTGGTGTTCTTCTATGGCCAGGGCACCTGGCCGGGGGCCACTTGCGTGGAGTTGTTTCGCGAGGAGGTGGTGCCGGTGTGTgc                                                                                           <  5:210237/149‑1 (MQ=255)                                                         ggcaATATCCACCTGGA‑CATTCGCAACGAGATGGAGCCGTTCGCCCTGCTGCAAGGCTCGGCTGACGTGGTGTTCTTCTATGGCCAGGGCACCTGGCCGGGGGCCACTTGCGTGGAGTTGTTTCGCGAGGAGGTGGTGCCGGTGTGTgc                                                                                           <  8:144163/149‑1 (MQ=255)                                                                           cATTCGCAACGAGATGGAGCCGTTCGCCCTGCTGCAAGGCTCGGCTGACGTGGTGTTCTTCTATGGCCAGGGCACCTGGCCGGGGGCCACTTGCGTGGAGTTGTTTCGCGAGGAGGTGGTGCCGGTGTGTGCGCCAGAGCTGCTGGCCg                                                                          >  3:61892/1‑149 (MQ=255)                                                                            aTTCGCATCGAGATGGAGCCGTTCGCCCTGCTGCAAGGCGCGGCTGACGTGGTGTTCTTCTATGGCCAGGGCACCTGGCCGGGGGCCACTTGCGTGGAGTTGTTTCGCGAGGAGGTGGTGCCGGTGTGTGCGCCAGAGCTGCTGGCCgg                                                                         <  5:183134/149‑1 (MQ=255)                                                                                gCAACGAGATGGAGCCGTTCGCCCTGCTGCAAGGCTCGGCTGACGTGGTGTTCTTCTATGGCCAGGGCACCTGGCCGGGGGCCACTTGCGTGGAGTTGTTTCGCGAGGAGGTGGTGCCGGTGTGTGCGCCAGAGCTGCTGGCCGGgcg                                                                      >  4:307733/1‑148 (MQ=255)                                                                                      agaTGGAGCCGTTCGCCCTGCTGCAAGGCTCGGCTGACGTGGTGTTCTTCTATGGCCAGGGCACCTGGCCGGGGGCCACTTGCGTGGAGTTGTTTCGCGAGGAGGTGGTGCCGGTGTGTGCGCCAGAGCTGCTGGCCGGGCGCAAGCTg                                                               >  8:119232/1‑149 (MQ=255)                                                                                         tGGAGCCGTTCGCCCTGCTGCAAGGCTCGGCTGACGTGGTGTTCTTCTATGGCCAGGGCACCTGGCCGGGGGCCACTTGCGTGGAGTTGTTTCGCGAGGAGGTGGTGCCGGTGTGTGCGCCAGAGCTGCTGGCCGGGCGCAAGCTGGcc                                                            >  2:57759/1‑149 (MQ=255)                                                                                          ggAGCCGTTCGCCCTGCTGCAAGGCTCGGCTGACGTGGTGTTCTTCTATGGCCAGGGCACCTGGCCGGGGGCCACTTGCGTGGAGTTGTTTCGCGAGGAGGTGGTGCCGgt                                                                                                 >  3:160470/1‑111 (MQ=255)                                                                                          ggAGCCGTTCGCCCTGCTGCAAGGCTCGGCTGACGTGGTGTTCTTCTATGGCCAGGGCACCTGGCCGGGGGCCACTTGCGTGGAGTTGTTTCGCGAGGAGGTGGTGCCGgt                                                                                                 <  4:160470/111‑1 (MQ=255)                                                                                          ggAGCCGTTCGCCCTGCTGCAAGGCTCGGCTGACGTGGTGTTCTTCTATGGCCAGGGCACCTGGCCGGGGGCCACTTGCGTGGAGTTGTTTCGCGAGGAGGTGGTGCCGGTGTGTGCGCCAGAGCTGCTGGCCGGGCGCAAGCTGGCCg                                                           >  8:210938/1‑149 (MQ=255)                                                                                              ccGTTCGCCCTGCTGCAAGGCTCGGCTGACGTGGTGTTCTTCTATGGCCAGGGCACCTGGCCGGGGGCCACTTGCGTGGAGTTGTTTCGCGAGGAGGTGGTGCCGGTGTGTGCGCCAGAGCTGCTGGCCGGGCGCAAGCTGGCCGGTGc                                                       <  6:326914/149‑1 (MQ=255)                                                                                                      ccTGCTGCAAGGCTCGGCTGACGTGGTGTTCTTCTATGGCCAGGGCACCTGGCCGGGGGCCACTTGCGTGGAGTTGTTTCGCGAGGAGGTGGTGCCGGTGTGTGCGCCAGAGCTGCTGGCCGGGCGCAAGCTGGCCGGTGCCGAGGAGc                                               >  7:289906/1‑149 (MQ=255)                                                                                                       ctgctgCAAGGCTCGGCTGACGTGGTGTTCTTCTATGGCCAGGGCACCTGGCCGGGGGCCACTTGCGTGGAGTTGTTTCGCGAGGAGGTGGTGCCGGTGTGTGCGCCAGAGCTGCTGGCCGGGCGCAAGCTGGCCGGTGCCGAGGAGCt                                              <  7:119232/149‑1 (MQ=255)                                                                                                         gctgcAAGGCTCGGCTGACGTGGTGTTCTTCTATGGCCAGGGCACCTGGCCGGGGGCCACTTGCGTGGAGTTGTTTCGCGAGGAGGTGGTGCCGGTGTGTGCGCCAGAGCTGCTGGCCGGGCGCAAGCTGGCCGGTGCCGAGGAGCTgg                                            >  8:96545/1‑149 (MQ=255)                                                                                                                 gCTCGGCTGACGTGGTGTTCTTCTATGGCCAGGGCACCTGGCCGGGGGCCACTTGCGTGGAGTTGTTTCGCGAGGAGGTGGTGCCGGTGTGTGCGCCAGAGCTGCTGGCCGGGCGCAAGCTGGCCGGTGCCGAGGAGCTGGCCGAGccg                                    <  6:241774/149‑3 (MQ=255)                                                                                                                   tCGGCTTACTTGGGGTTCTTCTATGGCCTGGGCACCTGGCCGGGGGCCACTTGCGTGGAGTTGTTTCGCGAGGAGGTGTTGCCGGTGTGTGCGCCAGAGCTGCTGGCCGGGCGCAAGCTGGCCGGTGCCGAGGAGCTGGCCGAGCTGGt                                  <  2:109145/149‑1 (MQ=255)                                                                                                                   tCGGCTGACGTGGTGTTCTTCTATGGCCAGGGCACCTGGCCGGGGGCCACTTGCGTGGAGTTGTTTCGCGAGGAGGTGGTGCCGGTGTGTGCGCCAGAGCTGCTGGcc                                                                           <  1:299307/108‑1 (MQ=255)                                                                                                                   tCGGCTGACGTGGTGTTCTTCTATGGCCAGGGCACCTGGCCGGGGGCCACTTGCGTGGAGTTGTTTCGCGAGGAGGTGGTGCCGGTGTGTGCGCCAGAGCTGCTGGcc                                                                           >  2:299307/1‑108 (MQ=255)                                                                                                                   tCGGCTGACGTGGTGTTCTTCTATGGCCAGGGCACCTGGCCGGGGGCCACTTGCGTGGAGTTGTTTCGCGAGGAGGTGGTGCCGGTGTGTGCGCCAGAGCTGCTGGCCGGGCGCAAGCTGGcc                                                            >  2:244004/1‑123 (MQ=255)                                                                                                                     ggCTGACGTGGTGTTCTTCTATGGCCAGGGCACCTGGCCGGGGGCCACTTGCGTGGAGTTGTTTCGCGAGGAGGTGGTGCCGGTGTGTGCGCCAGAGCTGCTGGCCGGGCGCAAGCTGGCCGGTGCCGAGGAGCTGGCCGAGc                                      >  1:170201/1‑143 (MQ=255)                                                                                                                     ggCTGACGTGGTGTTCTTCTATGGCCAGGGCACCTGGCCGGGGGCCACTTGCGTGGAGTTGTTTCGCGAGGAGGTGGTGCCGGTGTGTGCGCCAGAGCTGCTGGCCGGGCGCAAGCTGGCCGGTGCCGAGGAGCTGGCCGAGc                                      <  2:170201/143‑1 (MQ=255)                                                                                                                     ggCTGACGTGGTGTTCTTCTATGGCCAGGGCACCTGGCCGGGGGCCACTTGCGTGGAGTTGTTTCGCGAGGAGGTGGTGCCGGTGTGTGCGCCAGAGCTGCTGGCCGGGCGCAAGCTGGCCGGTGCCGAGGAGCTGGCCGAGCTGGtgc                                >  4:148742/1‑149 (MQ=255)                                                                                                                       cTGCCGTGGTGTTCTTCTATGGCCAGGGCACCTGGCCGGGGGCCACTTGCGTGGAGTGGTTTCGCGAGGAGGTGGTGCCGGTGTGTGCGCCAGAGCTGCTGGCCGGGCGCAAGCTGGCCGGTGCCGAGGAGCTGGCCGAGCTGGtgctg                              >  6:229290/1‑149 (MQ=255)                                                                                                                          aCGTGGTGTTCTTCTATGGCCAGGGCACCTGGCCGGGGGCCACTTGCGTGGAGTTGTTTCGCGAGGAGGTGGTGCCGGTGTGTGCGCCAGAGCTGCTGGCCGGGCGCAAGCTGGCCGGTGCCgagga                                                 >  5:293657/1‑127 (MQ=255)                                                                                                                          aCGTGGTGTTCTTCTATGGCCAGGGCACCTGGCCGGGGGCCACTTGCGTGGAGTTGTTTCGCGAGGAGGTGGTGCCGGTGTGTGCGCCAGAGCTGCTGGCCGGGCGCAAGCTGGCCGGTGCCgagga                                                 <  6:293657/127‑1 (MQ=255)                                                                                                                             tggtgTTCTTCTATGGCCAGGGCACCTGGCCGGGGGCCACTTGCGTGGAGTTGTTTCGCGAGGAGGTGGTGCCGGTGTGTGCGCCAGAgctg                                                                                 <  1:124344/92‑1 (MQ=255)                                                                                                                             tggtgTTCTTCTATGGCCAGGGCACCTGGCCGGGGGCCACTTGCGTGGAGTTGTTTCGCGAGGAGGTGGTGCCGGTGTGTGCGCCAGAgctg                                                                                 >  2:124344/1‑92 (MQ=255)                                                                                                                               gtgtTCTTCTATGGCCAGGGCACCTGGCCGGGGGCCACTTGCGTGGAGTTGTTTCGCGAGGAGGTGGTGCCGGTGTGTGCGCCAGAGCTGCTGGCCGGGCGCAAGCTGGCCGGTGCCGAGGAGCTGGCCGAGCTGGTGCTGCTGCaaa                       <  3:148742/148‑1 (MQ=255)                                                                                                                                    cttctATGGCCAGGGCACCTGGCCGGGGGCCACTTGCGTGGAGTTGTTTCGCGAGGAGGTGGTGCCGGTGTGTGCGCCAGAGCTGCTGGCCGGGCGCAAGCTGGCCGGTGCCGAGGAGCTGGCCGAGCTGGTGCTGCTGCAAAGCAcc                  <  3:307733/148‑1 (MQ=255)                                                                                                                                    cttctATGGCCAGGGCACCTGGCCGGGGGCCACTTGCGTGGAGTTGTTTCGCGAGGAGGTGGTGCCGGTGTGTGCGCCAGAGCTGCTGGCCGGGCGCAAGCTGGCCGGTGCCGAGGAGCTGGCCGAGCTGGTGCTGCTGCAAAGCACCt                 >  3:208803/1‑149 (MQ=255)                                                                                                                                      tctATGGCCAGGGCACCTGGCCGGGGGCCACTTGCGTGGAGTTGTTTCGCGAGGAGGTGGTGCCGGTGTGTGCGCCAGAGCTGCTGGCCGGGCGCAAGCTGGCCGGTGCCGAGGAGCTGGCCGAGCTGGTGCTGCTGCAAAGCACCTcg               >  5:313701/1‑149 (MQ=255)                                                                                                                                      tctATGGCCAGGGCACCTGGCCGGGGGCCACTTGCGTGGAGTTGTTTCGCGAGGAGGTGGTGCCGGTGTGTGCGCCAGAGCTGCTGGCCGGGCGCAAGCTGGCCGGTGCCGAGGAGCTGGCCGAGCTGGTGCTGCTGCAAAGCACCTcg               >  5:145170/1‑149 (MQ=255)                                                                                                                                       ctATGGCCAGGGCACCTGGCCGGGGGCCACTTGCGTGGAGTTGTTTCGCGAGGAGGTGGTGCCGGTGTGTGCGCCAGAGCTGCTGGCCGGGCGCAAGCTGGCCGGTGCCGAGGAGCTGGCCGAGCTGGTGCTGCTGCAAAGCACCTcgc              >  3:239861/1‑149 (MQ=255)                                                                                                                                       ctATGGCCAGGGCACCTGGCCGGGGGCCACTTGCGTGGAGTTGTTTCGCGAGGAGGTGGTGCCGGTGTGTGCGCCAGAGCTGCTGGCCGGGCGCAAGCTGGCCGGTGCCGAGGAGCTGGCCGAGCTGGTGCTGCTGCAAAGCACCTc                <  7:210938/147‑1 (MQ=255)                                                                                                                                       ccccggCCAGGTGCCCTGGCCGGGGGCCACTTGCGTGGAGTTGTTTCGCGAGGAGGTGGTGCCGGTGTGTGCGCCAGAGCTGCTGGCCGGGCGCAAGCTg                                                               <  2:24900/96‑1 (MQ=255)                                                                                                                                       ccccggCCAGGTGCCCTGGCCGGGGGCCACTTGCGTGGAGTTGTTTCGCGAGGAGGTGGTGCCGGTGTGTGCGCCAGAGCTGCTGGCCGGGCGCAAGCTg                                                               >  1:24900/5‑100 (MQ=255)                                                                                                                                         ccggCCAGGTGCCCTGGCCGGGGGCCACTTGCGTGGAGTTGTTTCGCGAGGAGGTGGTGCCGGTGTGTGCGCCAGAGCTGCTGGCCGGGCGCAAGCTGGCCGGTGCCGAGGAGCTGGCCGAGCTGGTGCTGCTGCAAAGCACCTCgcgc            >  5:101079/3‑149 (MQ=255)                                                                                                                                            gCCAGGGCACCTGGCCGGGGGCCACTTGCGTGGAGTTGTTTCGCGAGGAGGTGGTGCCGGTGTGTGCGCCAGAGCTGCTGGCCGGGCGCAAGCTGGCCGGTGCCGAGGAGCTGGCCGAGCTGGTGCTGCTGCAAAGCACCTCGCGCCCg         >  1:227194/1‑149 (MQ=255)                                                                                                                                              cAGGGCACCTGGCCGGGGGCCACTTGCGTGGAGTTGTTTCGCGAGGAGGTGGTGCCGGTGTGTGCGCCAGAGCTGCTGGCCGGGCGCAAGCTGGCCGGTGCCGAGGAGCTGGCCGAGCTGGTGCTGCTGCAAAGCACCTCGCGCCCg         >  3:198651/1‑147 (MQ=255)                                                                                                                                              cAGGGCACCTGGCCGGGGGCCACTTGCGTGGAGTTGTTTCGCGAGGAGGTGGTGCCGGTGTGTGCGCCAGAGCTGCTGGCCGGGCGCAAGCTGGCCGGTGCCGAGGAGCTGGCCGAGCTGGTGCTGCTGCAAAGCACCTCGCGCCCg         >  7:180031/1‑147 (MQ=255)                                                                                                                                              cAGGGCACCTGGCCGGGGGCCACTTGCGTGGAGTTGTTTCGCGAGGAGGTGGTGCCGGTGTGTGCGCCAGAGCTGCTGGCCGGGCGCAAGCTGGCCGGTGCCGAGGAGCTGGCCGAGCTGGTGCTGCTGCAAAGCACCTCGCGCCCGGa       <  4:61892/149‑1 (MQ=255)                                                                                                                                               aGGGCACCTGGCCGGGGGCCACTTGCGTGGAGTTGTTTCGCGAGGAGGTGGTGCCGGTGTGTGCGCCAGAGCTGCTGGCCGGGCGCAAGCTGGCCGGTGCCGAGGAGCTGGCCGAGCTGGTGCTGCTGCAAAGCACCTCGCGCCCGGAg      >  7:212964/1‑149 (MQ=255)                                                                                                                                                   cACCTGGCCGGGGGCCACTTGCGTGGAGTTGTTTCGCGAGGAGGTGGTGCCGGTGTGTGCGCCAGAGCTGCTGGCCGGGCGCAAGCTGGCCGGTGCCGAGGAGCTGGCCGAGCTGGTGCTGCTGCAAAGCACCTCGCGCCCGGAGGCCt  >  5:2639/1‑149 (MQ=255)                                                                                                                                                   cACCGGGCCGGGGGCCACTTGCGGGGAGTGGTTTCGCGAGGACGTGGTGCCGGGGTGTCCGCCAGAGCTGCTGGCCGGGCGCAAGCTGGCCGGTGCCGAGGAGCTGGCCGAGCTGGTGCTGCTGCAAAGCACCTCGCGCCCGGAGGCCt  <  4:309897/149‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                                    CGAGTTTTGGCGTGCGCTGGTTGATCCCGCACCTCAAGGGGTTTGGCAAACGGCACGGCAATATCCACCTGGA‑CATTCGCAACGAGATGGAGCCGTTCGCCCTGCTGCAAGGCTCGGCTGACGTGGTGTTCTTCTATGGCCAGGGCACCTGGCCGGGGGCCACTTGCGTGGAGTTGTTTCGCGAGGAGGTGGTGCCGGTGTGTGCGCCAGAGCTGCTGGCCGGGCGCAAGCTGGCCGGTGCCGAGGAGCTGGCCGAGCTGGTGCTGCTGCAAAGCACCTCGCGCCCGGAGGCCT  >  NC_002947/6004382‑6004675

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