Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,729,636	T→A	6.3%	L42Q (CTG→CAG)	PP_3296 →	putative transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,729,636	0	T	A	6.3%	180.7 / 2.7	63	L42Q (CTG→CAG)	PP_3296	putative transporter
Reads supporting (aligned to +/- strand):  ref base T (27/32);  new base A (2/2);  total (29/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 = 7.34e-01

GCATATACTTTCAGGAGCATGAAATGAAGCGGATCAGTGGCTTTATCGCAGGGTTGCTGTTCGGCATGGGAATATTGCTATCCGGCATGGCAAACCCCGCCAAGGTGCTGGGGTTTCTCGACATTGCCGGGCAGTGGGACCCGTCGCTGGGACTGGT‑‑CATGGTGGGTGCGATAGGTGCCGCACTGTTGCCGATGACTTGGGCGCGACGTCATACCACGGCCTTGCTGGGAAGTAAGATACAGCTGCCTGCCAGGCGCGACATCGACCGTCGACTCGTCGGGGGGAGCCTG  >  NC_002947/3729489‑3729778                                                                                                                                                    |                                                                                                                                                 gCATATACTTTCAGGAGCATGAAATGAAGCGGATCAGTGGCTTTATCGCAGGGTTGCTGTTCGGCATGGGAATATTGCTATCCGGCATGGCAAACCCCGCCAAGGTGCTGGGGTTTCTCGACATTGCCGGGCAGTGGGACCCGTCGCTg                                                                                                                                                 >  4:105672/1‑149 (MQ=255)  cATATACTTTCAGGAGCATGAAATGAAGCGGATCAGTGGCTTTATCGCAGGGTTGCTGTTCGGCATGGGAATATTGCTATCCGGCATGGCAAACCCCGCCAAGGTGCTGGGGTTTCTCGACATTGCCGGGCAGTGGGACCCGTCGCTgg                                                                                                                                                >  8:265595/1‑149 (MQ=255)   atatACTTTCAGGAGCATGAAATGAAGCGGATCAGTGGCTTTATCGCAGGGTTGCTGTTCGGCATGGGAATATTGCTATCCGGCATGGCAAACCCCGCCAAGGTGCTGGGGTTTCTCGACATTGCCGGGCAGTGGGACCCGTCGCTgg                                                                                                                                                >  2:14475/1‑148 (MQ=255)   atatACTTTCAGGAGCATGAAATGAAGCGGATCAGTGGCTTTATCGCAGGGTTGCTGTTCGGCATGGGAATATTGCTATCCGGCATGGCAAACCCCGCCAAGGTGCTGGGGTTTCTCGACATTGCCGGGCAGTGGGACCCGTCGCTgg                                                                                                                                                >  1:95908/1‑148 (MQ=255)   atatAATTTCAGGAGCATGAAATGAAGCGGATAAGTGGATTTATCGAAGTGTTGATGTTCGTCATGGGAATATTGCTATCCGGCATGGCAAACCACTCCAAGGTGATGTGGTTTATCGACATTGCAGGGCAGTGTGACCCGTCGCTgg                                                                                                                                                >  1:91217/1‑148 (MQ=255)      taCTTTCAGGAGCATGAAATGAAGCGGATCAGTGGCTTTATCGCAGGGTTGCTGTTCGGCATGGGAATATTGCTATCCGGCATGGCAAACCCCGCCAAGGTGCTGGGGTTTCTCGACATTGCCGGGCAGTGGGACCCGTCGCTGGGAc                                                                                                                                             <  3:105672/148‑1 (MQ=255)             cttaGCATGAAATGAATCGGATAAGTGGCTTTATCGCAGGGTTGCGTTTATGAATGGTAATATTGCTATCCGGCATGGCAAACCCCTCCAATGTGCTGTGGTTTCTCGACATTGCCGGGCAGTGGGACCCGTCGCTGGGACTGGT‑‑CAt                                                                                                                                    <  2:91217/145‑1 (MQ=255)             aGGAGCATGAAATGAAGCGGATCAGTGGCTTTATCGCAGGGTTGCTGTTCGGCATGGGAATATTGCTATCCGGCATGGCAAACCCCGCCAAGGTGCTGGGGTTTCTCGACATTGCCGGGCAGTGGGACCCGTCGCTGGGACTGGT‑‑CAt                                                                                                                                    <  2:95908/148‑1 (MQ=255)                   atgaaatgaaGCGGATCAGTGGCTTTATCGCAGGGTTGCTGTTCGGCATGGGAATATTGCTATCCGGCATGGCAAACCCCGCCAAGGTGCTGGGGTTTCTCGACATTGCCGGGCAGTGGGACCCGTCGCTGGGACTGGT‑‑CATggtgggt                                                                                                                             >  8:189120/1‑149 (MQ=255)                      aaatgaaGCGGATCAGTGGCTTTATCGCAGGGTTGCTGTTCGGCATGGGAATATTGCTATCCGGCATGACAAACCCCGCCAAGGTGCTGGGGTTTCTCGACATTGCCGGGCAGTGGGACCCGTCGCTGGGACTGGT‑‑CATGGTGGGTGc                                                                                                                           >  7:320758/1‑148 (MQ=255)                       aatgaaGCGGATCAGTGGCTTTATCGCAGGGTTGCTGTTCGGCATGGGAATATTGCTATCCGGCATGGCAAACCCCGCCAAGGTGCTGGGGTTTCTCGACATTGCCGGGCAGTGGGACCCGTCGCTGGGACTGGT‑‑CATGGTGGGTGCGa                                                                                                                         <  8:142587/149‑1 (MQ=255)                        atgaaGCGGATCAGTGGCTTTATCGCAGGGTTGCTGTTCGGCATGGGAATATTGCTATCCGGCATGGCAAACCCCGCCAAGGTGCTGGGGTTTCTCGACATTGCCGGGCAGTGGGACCCGTCGCTGGGACTGGT‑‑CATGGTGGGTGCGa                                                                                                                         >  3:319190/1‑148 (MQ=255)                         tgaaGCGGATCAGTGGCTTTATCGCAGGGTTGCTGTTCGGCATGGGAATATTGCTATCCGGCATGGCAAACCCCGCCAAGGTGCTGGGGTTTCTCGACATTGCCGGGCAGTGGGACCCGTCGCTGGGACTGGT‑‑CATGGTGGGTGCGATa                                                                                                                       <  7:265595/149‑1 (MQ=255)                         tgaaGCGGATCAGTGGCTTTATCGCAGGGTTGCTGTTCGGCATGGGAATATTGCTATCCGGCATGGCAAACCCCGCCAAGGTGCTGGGGTTTCTCGACATTGCCGGGCAGTGGGACCCGTCGCTGGGACTGGT‑‑CATGGTGGGTGCGATa                                                                                                                       <  4:95852/149‑1 (MQ=255)                           aaGCTGATCAGTGGCTTTATCGCAGGGTTGCTGTTCGGCATGGGAATATTGCTATCCGGCATGGCAAACCCCGCCAAGGTGCTGGGGTTTCTCGACATTGCCGGGCAGTGGGACCCGTCGCTGGGACTGGT‑‑CATGGTGGGTGCGATAg                                                                                                                      <  4:233702/148‑1 (MQ=255)                           aaGCGGATCAGTGGCTTTATCGCAGGGTTGCTGTTCGGCATGGGAATATTGCTATCCGGCATGGCAAACCCCGCCAAGGTGCTGGGGTTTCTCGACATTGCCGGGCAGTGGGACCCGTCGCTGGGACTGGT‑‑CATGGTGGGTGCGATAgg                                                                                                                     <  6:327813/149‑1 (MQ=255)                                    aGTGGCTTTATCGCAGGGTTGATGTTCGGCATGTGAATATTGCTATCCTGCATGACAAACCCCGCCAAGGTGCTGGGGTTTCTCGACATTTCCGGGCAGTGGGACCCGTCGCTGGGACTGGT‑‑CATGGTGGGTGCGATAGGTGCCGCACt                                                                                                            <  8:320758/149‑1 (MQ=255)                                     gTGGCTTTATCGCAGGGTTGCTGTTCGGCATGGGAATATTGCTATCCGGCATGGCAAACCCCGCCAAGGTGCTGGGGTTTCTCGACATTGCCGGGCAGTGGGACCCGTCGCTGGGACTGGT‑‑CAtggtgg                                                                                                                               <  4:17001/129‑1 (MQ=255)                                     gTGGCTTTATCGCAGGGTTGCTGTTCGGCATGGGAATATTGCTATCCGGCATGGCAAACCCCGCCAAGGTGCTGGGGTTTCTCGACATTGCCGGGCAGTGGGACCCGTCGCTGGGACTGGT‑‑CAtggtgg                                                                                                                               >  3:17001/1‑129 (MQ=255)                                     gTGGCTTTATCGCAGGGTTGCTGTTCGGCATGGGAATATTGCTATCCGGCATGGCAAACCCCGCCAAGGTGCTGGGGTTTCTCGACATTGCCGGGCAGTGGGACCCGTCGCTGGGACTGGT‑‑CATGGTGGGTGCGATAGGTGCCGCACTg                                                                                                           <  1:297296/149‑1 (MQ=255)                                     gTGGCTTTATCGCAGGGTTGCTGTTCGGCATGGGAATATTGCTATCCGGCATGGCAAACCCCGCCAAGGTGCTGGGGTTTCTCGACATTGCCGGGCAGTGGGACCAGTCCCAGCGACGGGT‑‑cc                                                                                                                                     <  8:115991/123‑2 (MQ=255)                                     gTGGCTTTATCGCAGGGTTGCTGTTCGGCATGGGAATATTGCTATCCGGCATGGCAAACCCCGCCAAGGTGCTGGGGTTTCTCGACATTGCCGGGCAGTGGGACCAGTCCCAGCGACGGGT‑‑cc                                                                                                                                     >  7:115991/1‑122 (MQ=255)                                        gCTTTATCGCAGGGTTGCTGTTCGGCATGGGAATATTGCTATCCGGCATGGCAAACCCCGCCAAGGTGCTGGGGTTTCTCGACATTGCCGGGCAGTGGGACCCGTCGCTGGGACTGGT‑‑CATGGTGGGTGCGATAGGTGCCGCACTGTTg                                                                                                        >  2:322959/1‑149 (MQ=255)                                        gCTTTATCGCAGGGTTGCTGTTCGGCATGGGAATATTGCTATCAGGCATGGCAAACCCCGCCAAGGTGCTGGGGTTTCTCGACATTGCCGGGCAGTGGGACCCGTCGCTGGGACTGGT‑‑CATGGTGGGTGCGATAGGTGCCGCACTGTTg                                                                                                        >  3:28738/1‑149 (MQ=255)                                             aTCGCAGGGTTGCTGTTCGGCATGGGAATATTGCTATCCGGCATGGCAAACCCCGCCAAGGTGCTGGGGTTTCTCGACATTGCCGGGCAGTGGGACCCGTCGCTGGGACTGGT‑‑CATGGTGGGTGCGATAGGTGCCGCACTGTTGCCGa                                                                                                    <  7:110442/148‑1 (MQ=255)                                              tCGCAGGGTTGCTGTTCGGCATGGGAATATTGCTATCCGGCATGGCAAACCCCGCCAAGGTGCTGGGGTTTCTCGACATTGCCGGGCAGTGGGACCCGTCGCTGGGACTGGT‑‑CATGGTGGGTGCGATAGGTGCCGCACTGTTGCCGAt                                                                                                   <  7:226644/148‑1 (MQ=255)                                                 cAGGGTTGCTGTTCGGCATGGGAATATTGCTATCCGGCATGGCAAACCCCGCCAAGGTGCTGGGGTTTCTCGACATTGCCGGGCAGTGGGACCCGTCGCTGGGACTGGT‑‑Ca                                                                                                                                     >  1:69202/1‑111 (MQ=255)                                                 cAGGGTTGCTGTTCGGCATGGGAATATTGCTATCCGGCATGGCAAACCCCGCCAAGGTGCTGGGGTTTCTCGACATTGCCGGGCAGTGGGACCCGTCGCTGGGACTGGT‑‑Ca                                                                                                                                     <  2:69202/111‑1 (MQ=255)                                                           gTTCGGCATGGGAATATTGCTATCCGGCATGGCAAACCCCGCCAAGGTGCTGGGGTTTCTCGACATTGCCGGGCAGTGGGACCCGTCGCTGGGACTGGT‑‑CATGGTGGGTGCGATAGGTGCCGCACTGTTGCCGATGACTTGGGCGCGAc                                                                                     <  1:322959/149‑1 (MQ=255)                                                            ttCGGCATGGGAATATTGCTATCAGGCATGGCAAACCCCGCCAAGGTGCTGGGGTTTCTCGACATTGCCGGGCAGTGGGACCCGTCGCTGGGACTGGT‑‑CATGGTGGGTGCGATAGGTGCCGCACTGTTGCCGATGACTTGGGCGCGAc                                                                                     <  4:28738/148‑1 (MQ=255)                                                               ggCATGGGAATATTGCTATCCGGCATGGCAAACCCCGCCAAGGTGCTGGGGTTTCTCGACATTGCCGGGCAGTGGGACCCGTCGCTGGGACTGGT‑‑CATGGTGGGTGCGATAGGTGCCGCACTGTTGCCGATGACTTGGGCGCGACGTCa                                                                                 >  6:345784/1‑149 (MQ=255)                                                                gCATGGGAATATTGCTATCCGGCATGGCAAACCCCGCCAAGGTGCTGGGGTTTCTCGACATTGCCGGGCAGTGGGACCCGTCGCTGGGACTGGT‑‑CATGGTGGGTGCGATAGGTGCCGCACTGTTGCCGATGACTTGGGCGCGACGTCa                                                                                 >  3:340960/1‑148 (MQ=255)                                                                gCATGGGAATATTGCTATCCGGCATGGCAAACCCCGCCAAGGTGCTGGGGTTTCTCGACATTGCCGGGCAGTGGGACCCGTCGCTGGGACTGGT‑‑CATGGTGGGTGCGATAGGTGCCGCACTGTTGCCGATGACTTGGGCGCGACGTCAt                                                                                <  7:246553/149‑1 (MQ=255)                                                                  aTGGGAATATTGCTATCCGGCATGGCAAACCCCGCCAAGGTGCTGGGGTTTCTCGACATTGCCGGGCAGTGGGACCCGTCGCTGGGACTGGT‑‑CATGGTGGGTGCGATAGGTGCCGCACTGTTGCCGATGACTTGGGCGCGACGTCATAc                                                                              >  4:187045/1‑149 (MQ=255)                                                                      gAATATTGCTATCCGGCATGGCAAACACCGCCAAGGTGCTGGGGTTTCTCGACATTGACGGGCAGTGGGACCAGTCGCTGGGACTGGT‑‑CATGGTGGGTGCGATAGGTGCCGCACTGTTGCAGATGACGTGGGCGCGACGTCATAcgccg                                                                          >  2:134181/1‑145 (MQ=255)                                                                               tATCCGGCATGGCAAACCCCGCCAAGGTGCTGGGGTTTCTCGACATTGCCGGGCAGTGGGACCCGTCGCTGGGACTGGT‑‑CATGGTGGGTGCGATAGGTGCCGCACTGTTGCCGATGACTTGGGCGCGACGTCATACCACGGCCTTGCTg                                                                 >  5:91545/1‑149 (MQ=255)                                                                                aTCCGGCATGGCAAACCCCGCCAAGGTGCTGGGGTTTCTCGACATTGCCGGGCAGTGGGACCCGTCGCTGGGACTGGT‑‑CATGGTGGGTGCGATAGGTGCCGCACTGTTGCCGATGACTTGGGCGCGACGTCATACCACGGCCTTGCTgg                                                                >  1:156503/1‑149 (MQ=255)                                                                                   cGGCATGGCAAACCCCGCCAAGGTGCTGGGGTTTCTCGACATTGCCGGGCAGTGGGACCCGTCGCTGGGACTGGT‑‑CATGGTGGGTGCGATAGGTGCCGCACTGTTACCGATGACTTGGGCGCGACGTCATACCACGGCCTTGCTGGGa                                                              <  4:319190/148‑1 (MQ=255)                                                                                    ggCATGGCAAACCCCGCCAAGGTGCTGGGGTTTCTCGACATTGCCGGGCAGTGGGACACGTCGCTGGGACTGGT‑‑CATGGTGGGTGCGATAGGTGCCGCACTGTTGCCGATGACTTGGGCGCGACGTCATACCACGGCCTTGCTGGGAAg                                                            >  3:327623/1‑149 (MQ=255)                                                                                       aTGGCAAACCCCGCCAAGGTGCTGGGGTTTCTCGACATTGCCGGGCAGTGGGACCCGTCGCTGGGACTGGT‑‑CATGGTGGGTGCGATAGGTGCCGCACTGTTGCCGATGACTTGGGCGCGACGTCATACCACGGCCTTGCTGGGAAGTaa                                                         <  4:340960/149‑1 (MQ=255)                                                                                       aTGGCAAACCCCGCCAAGGTGCTGGGGTTTCTCGACATTGCCGGGCAGTGGGACCCGTCGCTGGGACTGGT‑‑CATGGTGGGTGCGATAGGTGCCGCACTGTTGCCGATGACTTGGGCGCGACGTCATACCACGGCCTTGCTGGGAAGTaa                                                         <  3:187045/149‑1 (MQ=255)                                                                                          gCAAACCCCGCCAAGGTGCTGGGGTTTCTCGACATTGCCGGGCAGTGGGACCCGTCGCTGGGACTGGT‑‑CATGGTGGGTGCGATAGGTGCCGCACTGTTGCCGATGACTTGGGCGCGACGTCATACCACGGCCTTGCTGGGAAGTAAGa                                                       >  8:165401/1‑148 (MQ=255)                                                                                              aCCCCGCCAAGGTGCTGGGGTTTCTCGACATTGCCTGGCAGTGGGACCCGTCGCTGGGACTGGT‑‑CATGGTGGGTGCGATAGGTGCCGCACTGTTGCCGATGACTTGGGCGCGACGTCATACCACGGCCTTGCTGGGAAGTAAGATACAg                                                  <  2:156503/149‑1 (MQ=255)                                                                                              aCCCCGCCAAGGTGCTGGGGTTTCTCGACATTGCCGGGCAGTGGGACCCGTCGCTGGGACTGGTGGCATGGTGGGTGCGATAGGTGCCGCACTGTTGCCGATGACTTGGGCGCGACGTCATACCACGGCCTTGCTGGGAAGTAAGATa                                                     >  1:286169/1‑148 (MQ=255)                                                                                                        ggTGCTGGGGTTTCTCGACATTGCCGGGCAGTGGGACCCGTCGCTGGGACTGGT‑‑CATGGTGGGTGCGATAGGTGCCGCACTGTTGCCGATGACTTGGGCGCGACGTCATACCACGGCCTTGCTGGGAAGTAAGATACAGTTGCCTGCCa                                        <  7:165401/149‑1 (MQ=255)                                                                                                                 gTTTCTCGACATTGCCGGGCAGTGGGACCCGTCCCTCGGACTCGT‑‑CATGGTGGGTGCGATAGGTGCCGCACTGATGCCGATGACTTGGGCTCGACGTCATACCACGTCGTTGCAGGGAAGTATGATACAGCTGCCTGACAGGcgcgcc                                >  5:85045/1‑146 (MQ=255)                                                                                                                  tttCTCGACATTGCAGGGCAGTGGGACACGTCGCTGGGACTGGT‑‑CAGGGTGGGTGCGATAGGTGCCGCACTGTTGCCGATGACTTGGGCGCGACGTCATACCACGGCCTTGCTGTGAAGTAAGATACAGCTGCCTGCCAGGCGCGACAt                              <  1:134181/149‑1 (MQ=255)                                                                                                                     ctcGACATTGCCGGGCAGTGGGACCCGTCGCTGGGACTGGT‑‑CATGGTGGGTGCGATAGGTGCCGCACTGTTGCCGATGACTTGGGCGCGACGTCATACCACGGCCTTGCTGGGAAGTAAGATACAGCTGCCTGCCAGGCGCGACATCg                            <  5:345784/148‑1 (MQ=255)                                                                                                                        gACATTGCCGGGCAGTGGGACCCGTCGCTGGGACTGGT‑‑CATGGTGGGTGCGATAGGTGCCGCACTGTTGCCGATGACTTGGGCGCGACGTCATACCACGGCCTTGCTGGGAAGTAAGATa                                                     >  3:187119/1‑120 (MQ=255)                                                                                                                        gACATTGCCGGGCAGTGGGACCCGTCGCTGGGACTGGT‑‑CATGGTGGGTGCGATAGGTGCCGCACTGTTGCCGATGACTTGGGCGCGACGTCATACCACGGCCTTGCTGGGAAGTAAGATa                                                     <  4:187119/120‑1 (MQ=255)                                                                                                                        gACATTGCCGGGCAGTGGGACCCGTCGCTGGGACTGGT‑‑CATGGTGGGTGCGATAGGTGCCGCACTGTTGCCGATGACTTGGGCGCGACGTCATACCACGGCCTTGCTGGGAAGTAAGATACAGCTGCCTGCCAGGCGCGACATCGa                           <  1:14475/146‑1 (MQ=255)                                                                                                                         aCATTGCCGGGCAGTGGGACCCGTCGCTGGGACTGGT‑‑CATGGTGGGTGCGATAGGTGCCGCACTGTTGCCGATGACTTGGGCGCGACGTCATACCACGGCCTTGCTGGGAAGTAAGATACAGCTGCCTGCCAGGCGCGACATCGACCGt                       <  2:286169/149‑1 (MQ=255)                                                                                                                               ccGGGCAGTGGGACCCGTCGCTGGGACTGGT‑‑CATGGTGGGTGCGATAGGTGCCGCACTGTTGCCGATGACTTGGGCGCGACGTCATACCACGGCCTTGCTGGGAAGTAAGATACAGCTGCCTGCCAGGCGCGACATCGACCGTCGACtc                 <  5:93967/149‑1 (MQ=255)                                                                                                                                cGGGCAGTGGGACCCGTCGCTGGGACTGGT‑‑CATGGTGGGTGCGATAGGTGCCGCACTGTTGCCGATGACTTGGGCGCGACGTCATACCACGGCCTTGCt                                                                  <  6:291182/99‑1 (MQ=255)                                                                                                                                cGGGCAGTGGGACCCGTCGCTGGGACTGGT‑‑CATGGTGGGTGCGATAGGTGCCGCACTGTTGCCGATGACTTGGGCGCGACGTCATACCACGGCCTTGCt                                                                  >  5:291182/1‑99 (MQ=255)                                                                                                                                   gCAGTGGGACCCGTCGCTGGGACTGGT‑‑CATGGTGGGTGCGATAGGTGCCGCACTGTTGCCGATGACTTGGGCGCGACGTCATACCACGGCCTTGCTGGGAAGTAAGATACAGCTGCCTGCCAGGCGCGACATCGACCGTCGACtcgtcg             >  7:13414/1‑149 (MQ=255)                                                                                                                                   gCAGTGGGACCCGTCGCTGGGACTGGT‑‑CATGGTGGGTGCGATAGGTGCCGCACTGTTGCCGATGACTTGGGCGCGACGTCATACCACGGCCTTGCTGGGAAGTAAGATACAGCTGCCTGCCAGGCGCGACATCGACCGTCGACtcgtcg             <  7:189120/149‑1 (MQ=255)                                                                                                                                       catgACCAGTCCCAGCGACGGGT‑‑CATGGTGGGTGCGATAGGTGCCGCACTGTTGCCGATGACTTGGGCGCGACGTCATACCACGGCCTTGCTGGGAAGTAAGATACAGCTGCCTGCCAGGCGCGACATCGACCGTCGACTCGTCgggg          <  4:229171/145‑1 (MQ=255)                                                                                                                                       catgACCAGTCCCAGCGACGGGT‑‑CATGGTGGGTGCGATAGGTGCCGCACTGTTGCCGATGACTTGGGCGCGACGTCATACCACGGCCTTGCTGGGAAGTAAGATACAGCTGCCTGCCAGGCGCGACATCGACCGTCGACTCGTCgggg          >  3:229171/4‑148 (MQ=255)                                                                                                                                           aCCCGTCGCTGGGACTGGT‑‑CATGGTGGGTGCGATAGGTGCCGCACTGTTGCCGATGACTTGGGCGCGACGTCATACCACGGCCTTGCTGGGAAGTAAGATACAGCTGCCTGCCAGGCGCGa                                 >  6:18475/1‑121 (MQ=255)                                                                                                                                           aCCCGTCGCTGGGACTGGT‑‑CATGGTGGGTGCGATAGGTGCCGCACTGTTGCCGATGACTTGGGCGCGACGTCATACCACGGCCTTGCTGGGAAGTAAGATACAGCTGCCTGCCAGGCGCGa                                 <  5:18475/121‑1 (MQ=255)                                                                                                                                             ccGTCGCTGGGACTGGT‑‑CATGGTGGGTGCGATAGGTGCCGCACTGTTGCCGATGACTTGGGCGCGACGTCATACCACGGCCTTGCTGGGAAGTAAGATACAGCTGCCTGCCAGGCGCGACATCGACCGTCGACTCGTCGGGGGGAGCCt   >  1:138740/1‑149 (MQ=255)                                                                                                                                              cGTCGCTGGGACTGGT‑‑CATGGTGGGTGCGATAGGCGCCGCACTGTTGCCGATGACTTGGGCGCGACGTCATACCACGGCCTTGCTGGGAAGTAAGATACAGCTGCCTGCCAGGCGCGACATCGACCGTCGACTCGTCGGGGGGAGCCTg  <  7:141703/149‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                                 GCATATACTTTCAGGAGCATGAAATGAAGCGGATCAGTGGCTTTATCGCAGGGTTGCTGTTCGGCATGGGAATATTGCTATCCGGCATGGCAAACCCCGCCAAGGTGCTGGGGTTTCTCGACATTGCCGGGCAGTGGGACCCGTCGCTGGGACTGGT‑‑CATGGTGGGTGCGATAGGTGCCGCACTGTTGCCGATGACTTGGGCGCGACGTCATACCACGGCCTTGCTGGGAAGTAAGATACAGCTGCCTGCCAGGCGCGACATCGACCGTCGACTCGTCGGGGGGAGCCTG  >  NC_002947/3729489‑3729778

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