Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,729,630	C→A	7.6%	P40Q (CCG→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,630	0	C	A	7.6%	181.3 / 4.4	62	P40Q (CCG→CAG)	PP_3296	putative transporter
Reads supporting (aligned to +/- strand):  ref base C (26/31);  new base A (3/2);  total (29/33)
Fisher's exact test for biased strand distribution p-value = 6.58e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.59e-01

CAGGCATATACTTTCAGGAGCATGAAATGAAGCGGATCAGTGGCTTTATCGCAGGGTTGCTGTTCGGCATGGGAATATTGCTATCCGGCATGGCAAACCCCGCCAAGGTGCTGGGGTTTCTCGACATTGCCGGGCAGTGGGACCCGTCGCTGGGACTGGT‑‑CATGGTGGGTGCGATAGGTGCCGCACTGTTGCCGATGACTTGGGCGCGACGTCATACCACGGCCTTGCTGGGAAGTAAGATACAGCTGCCTGCCAGGCGCGACATCGACCGTCGACTCGTCGGGGGGAGCCTG  >  NC_002947/3729486‑3729778                                                                                                                                                 |                                                                                                                                                       cAGGCATATACTTTCAGGAGCATGAAATGAAGCGGATCAGTGGCTTTATCGCAGGGTTGCTGTTCGGCATGGGAATATTGCTATCCGGCATGGCAAACCCCGCCAAGGTGCTGGGGTTTCTCGACATTGCCGGGCAGTGGGACCCGTCg                                                                                                                                                    >  7:142587/1‑149 (MQ=255)    gCATATACTTTCAGGAGCATGAAATGAAGCGGATCAGTGGCTTTATCGCAGGGTTGCTGTTCGGCATGGGAATATTGCTATCCGGCATGGCAAACCCCGCCAAGGTGCTGGGGTTTCTCGACATTGCCGGGCAGTGGGACCCGTCGCTg                                                                                                                                                 >  4:105672/1‑149 (MQ=255)     cATATACTTTCAGGAGCATGAAATGAAGCGGATCAGTGGCTTTATCGCAGGGTTGCTGTTCGGCATGGGAATATTGCTATCCGGCATGGCAAACCCCGCCAAGGTGCTGGGGTTTCTCGACATTGCCGGGCAGTGGGACCCGTCGCTgg                                                                                                                                                >  8:265595/1‑149 (MQ=255)      atatACTTTCAGGAGCATGAAATGAAGCGGATCAGTGGCTTTATCGCAGGGTTGCTGTTCGGCATGGGAATATTGCTATCCGGCATGGCAAACCCCGCCAAGGTGCTGGGGTTTCTCGACATTGCCGGGCAGTGGGACCCGTCGCTgg                                                                                                                                                >  1:95908/1‑148 (MQ=255)      atatACTTTCAGGAGCATGAAATGAAGCGGATCAGTGGCTTTATCGCAGGGTTGCTGTTCGGCATGGGAATATTGCTATCCGGCATGGCAAACCCCGCCAAGGTGCTGGGGTTTCTCGACATTGCCGGGCAGTGGGACCCGTCGCTgg                                                                                                                                                >  2:14475/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                                                                                                                               >  3:17001/1‑129 (MQ=255)                                        gTGGCTTTATCGCAGGGTTGCTGTTCGGCATGGGAATATTGCTATCCGGCATGGCAAACCCCGCCAAGGTGCTGGGGTTTCTCGACATTGCCGGGCAGTGGGACCCGTCGCTGGGACTGGT‑‑CAtggtgg                                                                                                                               <  4:17001/129‑1 (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                                                                                                                                     <  2:69202/111‑1 (MQ=255)                                                    cAGGGTTGCTGTTCGGCATGGGAATATTGCTATCCGGCATGGCAAACCCCGCCAAGGTGCTGGGGTTTCTCGACATTGCCGGGCAGTGGGACCCGTCGCTGGGACTGGT‑‑Ca                                                                                                                                     >  1:69202/1‑111 (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                                                         <  3:187045/149‑1 (MQ=255)                                                                                          aTGGCAAACCCCGCCAAGGTGCTGGGGTTTCTCGACATTGCCGGGCAGTGGGACCCGTCGCTGGGACTGGT‑‑CATGGTGGGTGCGATAGGTGCCGCACTGTTGCCGATGACTTGGGCGCGACGTCATACCACGGCCTTGCTGGGAAGTaa                                                         <  4:340960/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                                                     <  4:187119/120‑1 (MQ=255)                                                                                                                           gACATTGCCGGGCAGTGGGACCCGTCGCTGGGACTGGT‑‑CATGGTGGGTGCGATAGGTGCCGCACTGTTGCCGATGACTTGGGCGCGACGTCATACCACGGCCTTGCTGGGAAGTAAGATa                                                     >  3:187119/1‑120 (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                                                                  >  5:291182/1‑99 (MQ=255)                                                                                                                                   cGGGCAGTGGGACCCGTCGCTGGGACTGGT‑‑CATGGTGGGTGCGATAGGTGCCGCACTGTTGCCGATGACTTGGGCGCGACGTCATACCACGGCCTTGCt                                                                  <  6:291182/99‑1 (MQ=255)                                                                                                                                      gCAGTGGGACCCGTCGCTGGGACTGGT‑‑CATGGTGGGTGCGATAGGTGCCGCACTGTTGCCGATGACTTGGGCGCGACGTCATACCACGGCCTTGCTGGGAAGTAAGATACAGCTGCCTGCCAGGCGCGACATCGACCGTCGACtcgtcg             <  7:189120/149‑1 (MQ=255)                                                                                                                                      gCAGTGGGACCCGTCGCTGGGACTGGT‑‑CATGGTGGGTGCGATAGGTGCCGCACTGTTGCCGATGACTTGGGCGCGACGTCATACCACGGCCTTGCTGGGAAGTAAGATACAGCTGCCTGCCAGGCGCGACATCGACCGTCGACtcgtcg             >  7:13414/1‑149 (MQ=255)                                                                                                                                          catgACCAGTCCCAGCGACGGGT‑‑CATGGTGGGTGCGATAGGTGCCGCACTGTTGCCGATGACTTGGGCGCGACGTCATACCACGGCCTTGCTGGGAAGTAAGATACAGCTGCCTGCCAGGCGCGACATCGACCGTCGACTCGTCgggg          >  3:229171/4‑148 (MQ=255)                                                                                                                                          catgACCAGTCCCAGCGACGGGT‑‑CATGGTGGGTGCGATAGGTGCCGCACTGTTGCCGATGACTTGGGCGCGACGTCATACCACGGCCTTGCTGGGAAGTAAGATACAGCTGCCTGCCAGGCGCGACATCGACCGTCGACTCGTCgggg          <  4:229171/145‑1 (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)                                                                                                                                                 |                                                                                                                                                       CAGGCATATACTTTCAGGAGCATGAAATGAAGCGGATCAGTGGCTTTATCGCAGGGTTGCTGTTCGGCATGGGAATATTGCTATCCGGCATGGCAAACCCCGCCAAGGTGCTGGGGTTTCTCGACATTGCCGGGCAGTGGGACCCGTCGCTGGGACTGGT‑‑CATGGTGGGTGCGATAGGTGCCGCACTGTTGCCGATGACTTGGGCGCGACGTCATACCACGGCCTTGCTGGGAAGTAAGATACAGCTGCCTGCCAGGCGCGACATCGACCGTCGACTCGTCGGGGGGAGCCTG  >  NC_002947/3729486‑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: ‑