Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,534,777	C→A	8.8%	intergenic (+64/‑91)	atoB → / → PP_3124	3‑oxoacid CoA‑transferase subunit B/short chain fatty acid transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,534,777	0	C	A	8.8%	175.1 / 5.4	57	intergenic (+64/‑91)	atoB/PP_3124	3‑oxoacid CoA‑transferase subunit B/short chain fatty acid transporter
Reads supporting (aligned to +/- strand):  ref base C (30/22);  new base A (3/2);  total (33/24)
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.19e-01

CTGGGGTCAGCGTCGAGGAAATCATCGAGAAAACCGCCGGCAAACTGATCGTGCCGGATGATGTGAAGGAAATGACCTTCTAAGGCATGACGTGATCTCACCCTGTAGGCGCGGGTTTACCCGCGAATGCGTCGGTACATTGACCGCCGCATTCGCGGGT‑AAGCCTGTTCCTACAGG‑TTTTTTGCGTTGCCTGGCAATACCGCTTCATCCCATAAAACCAATAAAAGGAATTCAACCGTGGCCGCTGAAATTCAAGACAGCCGCTCCGCGCGCTTTGCCTTGCGCTGCTCCAA  >  NC_002947/3534631‑3534923                                                                                                                                                   |                                                                                                                                                     cTGGGGTCAGCGTCGAGGAAATCATCGAGAAAACCGCCGGCAAACTGATCGTGCCGGATGATGTGAAGGAAATGACCTTCTAAGGCATGACGTGATCTCACCCTGTAGGCGCGGGTTTACCCGCGAATGCGTCGGTACATTGAccgccg                                                                                                                                                    <  8:306785/149‑1 (MQ=255)    gggTCAGCGTCGAGGAAATCATCGAGAAAACCGCCGGCAAACTGATCGTGCCGGATGATGTGAAGGAAATGACCTTCTAAGGCATGACGTGATCTCACCCTGTAGGCGCGGGTTTACCCGCGAATGCGTCGGTACATTGACCGCCGCAt                                                                                                                                                 >  5:289646/1‑149 (MQ=255)         aGCGTCGAGGAAATCATCGAGAACACCGCCGGCAAACTGATCGTGCCGGATGATGTGAAGGAAATGACCTTCTAAGGCATGACGTGATCTCACCCTGTAGGCGCGGGTTTACCCGCGAATGCGTCGGTACATTGACCGCCGCATTcgcg                                                                                                                                            <  5:118788/149‑1 (MQ=255)          gCGTCGAGGAAATCATCGAGAAAACCGCCGGCAAACTGATCGTGCCGGATGATGTGAAGGAAATGACCTTCTAAGGCATGACGTGATCTCACCCTGTAGGCGCGGGTTTACCCGCGAATGCGTCGGTACATTGACCGCCGCATTcgc                                                                                                                                             >  6:361564/1‑147 (MQ=255)          gCGTCGAGGAAATCATCGAGAAAACCGCCGGCAAACTGATCGTGCCGGATGATGTGAAGGAAATGACCTTCTAAGGCATGACGTGATCTCACCCTGTAGGCGCGGGTTTACCCGCGAATGCGTCGGTACATTGACCGCCGCATTcgc                                                                                                                                             <  5:361564/147‑1 (MQ=255)               gAGGAAATCATCGAGAAAACCGCCGGCAAACTGATCGTGCCGGATGATGTGAAGGAAATGACCTTCTAAGGCATGACGTGATCTCACCCTGTAGGCGCGGGTTTACCCGCGAATGCGTCGGTACATTGACCGCCGCATTCGCGGGT‑AAg                                                                                                                                     >  2:330395/1‑149 (MQ=255)               gAGGAAATCATCGAGAAAACCGCCGGCAAACTGATCGTGCCGGATGATGTGAAGGAAATGACCTTCTAAGGCATGACGTGATCTCACCCTGTAGGCGCGGGTTTACCCGCGAATGCGTCGGTACATTGACCGCCGCATTCGCGGGT‑AAg                                                                                                                                     >  2:231946/1‑149 (MQ=255)                      tcatcGAGAAAACCGCCGGCAAACTGATCGTGCCGGATGATGTGAAGGAAATGACCTTCTAAGGCATGACGTGATCTCACCCTGTAGGCGCGGGTTTACCCGCGAATGCGTCGGTACATTGACCGCCGCATTCGCGGGT‑AAGCCTGTTc                                                                                                                              <  2:148582/149‑1 (MQ=255)                        atcGAGAAAACCGCCGGCAAACTGATCGTGCCGGATGATGTGAAGGAAATGACCTTCTAAGGCATGACGTGATCTCACCCTGTAGGCGCGGGTTTACCCGCGAATGCGTCGGTACATTGACCGCCGCATTCGCGGGT‑AAGCCTGTTcc                                                                                                                             >  6:111620/1‑148 (MQ=255)                         tcGAGAAAACCGCCGGCAAACTGATCGTGCCGGATGATGTGAAGGAAATGACCTTCTAAGGCATGACGTGATCTCACCCTGTAGGCGCGGGTTTACCCGCGAATGCGTCGGTACATTGACCGCCGCATTCGCGGGT‑AAGCCTGTTCCTa                                                                                                                           >  1:363432/1‑149 (MQ=255)                         tcGAGAAAACCGCCGGCAAACTGATCGTGCCGGATGATGTGAAGGAAATGACCTTCTAAGGCATGACGTGATCTCACCCTGTAGGCGCGGGTTTACCCGCGAATGCGTCGGTACATTGACCGCCGCATTCGCGGGT‑AAGCCTGTTCCTa                                                                                                                           >  1:265005/1‑149 (MQ=255)                          cGAGAAAACCGCCGGCAAACTGATCGTGCCGGATGATGTGAAGGAAATGACCTTCTAAGGCATGACGTGATCTCACCCTGTAGGCGCGGGTTTACCCGCGAATGCGTCGGTACATTGACCGCCGCATTCGCGGGT‑AAGCCTGTTCCTAc                                                                                                                          <  5:111620/149‑1 (MQ=255)                           gagaAAACCGCCGGCAAACTGATCGTGCCGGATGATGTGAAGGAAATGACCTTCTAAGGCATGACGTGATCTCACCCTGTAGGCGCGGGTTTACCCGCGAATGCGTCGGTACATTGACCGCCGCATTCGCGGGT‑AAGCCTGTTCCTAc                                                                                                                          <  6:353981/148‑1 (MQ=255)                           gagaAAACCGCCGGCAAACTGATCGTGCCGGATGATGTGAAGGAAATGACCTTCTAAGGCATGACGTGATCTCACCCTGTAGGCGCGGGTTTACCCGCGAATGCGTCGGTACATTGACCGCCGCATTCGCGGGT‑AAGCCTGTTCCTACa                                                                                                                         <  2:265005/149‑1 (MQ=255)                                aaCCGCCGGCAAACTGATCGTGCCGGATGATGTGAAGGAAATGACCTTCTAAGGCATGACGTGATCTCACCCTGTAGGCGCGGGTTTACCCGCGAATGCGTCGGTACATTGACCGCCGCATTCGCGGGT‑AAGCCTGTTCCTACAGG‑ttt                                                                                                                   >  8:155101/1‑149 (MQ=255)                                    gccgGCAAACTGATCGTGCCGGATGATGTGAAGGAAATGACCTTCTAAGGCATGACGTGATCTCACCCTGTAGGCGCGGGTTTACCCGCGAATGCGTCGGTACATTGACCGCCGCATTCGCGGGT‑AAGCCTGTTCCTACAGG‑TTTTTTg                                                                                                               <  1:152015/149‑1 (MQ=255)                                       ggCAAACTGATCGTGCCGGATGATGTTAAGGAAATGACCTTCTAAGGCATGACGTGATCTAACCCTGTAGGCGCGGGTTTACCCGCGAATGCGTCGGTACATAGAACCCCGCATTCGCGGGT‑AACCCTGTTCCTACAGG‑TTTTTTGCGt                                                                                                            >  5:312631/1‑149 (MQ=38)                                                aTCGTGCCGGATGATGTGAAGGAAATGACCTTCTAAGGCATGACGTGATCTCACCCTGTAGGCGCGGGTTTACCCGCGAATGCGTCGGTACATTGACCGCCGCATTCGCGGGTAAAGCCTGTTCCTACAGG‑TTTTTTGCGTTGCCTgg                                                                                                     <  1:231946/148‑1 (MQ=255)                                                 tCGTGCCGGATGATGTGAAGGAAATGACCTTCTAAGGCATGACGTGATCTCACCCTGTAGGCGCGGGTTTACCCGCGAATGCGTCGGTACATTGACCGCCGCATTCGCGGGT‑AAGCCTGTTACTACAGG‑TTTTTTGCGTTGCCTGGCaa                                                                                                  >  6:94428/1‑149 (MQ=255)                                                          atgatgTGAAGGAAATGACCTTCTAAGGCATGACGTGATCTCACCCTGTAGGCGCGGGTTTACCCGCGAATGCGTCGGTACATTGACCGCCGCATTCGCGGGT‑AAGCCTGTTCCTACAGG‑TTTTTTGCGTTGCCTGGCAATACCGCTTc                                                                                         <  5:174408/149‑1 (MQ=255)                                                                  aaGGAAATGACCTTCTAAGGCATGACGTGATCTCACCCTGTAGGCGCGAGTTTACCCCCGAATGCGTCGGTACATTGACCGCCGCATTCGCGGGT‑AAGCCTGTTCCTACAGG‑TTTTTTGCGTTGCCTGGCAATACCGCTTCATCCCAt                                                                                  >  7:317574/1‑148 (MQ=255)                                                                      aaaTGACCTTCTAAGGCATGACGTGATCTCACCCTGTAGGCGCGGGTTTACCCGCGAATGCGTCGGTACATTGACCGCCGCATTCGCGGGT‑AAGCCTGTTCCTACAGG‑TTTTTTGCGTTGCCTGGCAATACCGCTTCATCCCATAAAAc                                                                             <  5:94428/149‑1 (MQ=255)                                                                       aaTGACCTTCTAAGGCATGACGTGATCTCACCCTGTAGGCGCGGGTTTACCCGCGAATTAGTCGGTAAAATTACCCCCGCAATAGCGGGT‑AAGCCTGTTAGAAAAAG‑TTTTTTTAGTTGACAGGCCATAAACCTTAATACCaaaaaaa                                                                             >  2:349792/1‑147 (MQ=255)                                                                        aTGACCTTCTAAGGCATGACGTGATCTCACCCTGTAGGCGCGGGTTTACCCGCGAATGCGTCGGTACATTGACCGCCGCATTCGCGGGT‑AAGCCTGTTCCTACAGG‑TTTTTTGCGTTGCCTGGCAATACCGCTTCATCCCATAAAACCa                                                                           <  6:289646/149‑1 (MQ=255)                                                                                cTAAGGCATGACGTGATCTCACCCTGTAGGCGCGGGTTTACCCGCGAATGCGTCGGTACATTGACCGCCGCATTCGCGGGT‑AAGCCTGTTCCTACAGG‑TTTTTTGCGTTGCCTGGCAATACCGCTTCATCCCATAAACCCAATAAAAg                                                                    >  7:273019/1‑148 (MQ=255)                                                                                cTAAGGCATGACGTGATCTCACCCTGTAGGCGCGGGTTTACCCGCGAATGCGTCGGTACATTGACCGCCGCATTCGCGGGT‑AAGCCTGTTCCTACAGG‑TTTTTTGCGTTGCCTGGCAATACCGCTTCATCCCATAAAACCAATAAAAg                                                                    >  1:232748/1‑148 (MQ=255)                                                                                cTAAGGCATGACGTGATCTCACCCTGTAGGCGCGGGTTTACCCGCGAATGCGTCGGTACATTGACCGCCGCATTCGCGGGT‑AAGCCTGTTCCTACAGG‑TTTTTTGCGTTGCCTGGCAATACCGCTTCATCCCATAAAACCAATAAAAg                                                                    >  2:301765/1‑148 (MQ=255)                                                                                   aGGCATGACGTGATCTCACCCTGTAGGCGCGGGTTTACCCGCGAATGCGTCGGTACATTGACCGCCGCATTCGCGGGT‑AAGCCTGTTCCTACAGG‑TTTTTTGCGTTGCCTGGCAATACCGCTTCATc                                                                                      >  1:259745/1‑127 (MQ=255)                                                                                   aGGCATGACGTGATCTCACCCTGTAGGCGCGGGTTTACCCGCGAATGCGTCGGTACATTGACCGCCGCATTCGCGGGT‑AAGCCTGTTCCTACAGG‑TTTTTTGCGTTGCCTGGCAATACCGCTTCATc                                                                                      <  2:259745/127‑1 (MQ=255)                                                                                       aTGACGTGATCTCACCCTGTAGGCGCGGGTTTACCCGCGAATGCGTCGGTACATTGACCGCCGCATTCGCGGGT‑AAGCCTGTTCCTACAGG‑TTTTTTGCGTTGCCTGGCAATACCGCTTCATCCCATAAAACCaa                                                                          <  6:175643/135‑1 (MQ=255)                                                                                       aTGACGTGATCTCACCCTGTAGGCGCGGGTTTACCCGCGAATGCGTCGGTACATTGACCGCCGCATTCGCGGGT‑AAGCCTGTTCCTACAGG‑TTTTTTGCGTTGCCTGGCAATACCGCTTCATCCCATAAAACCaa                                                                          >  5:175643/1‑135 (MQ=255)                                                                                            gTGATCTCACCCTGTAGGCGCGGGTTTACCCGCGAATGCGTCGGTACATTGACCGCCGCATTCGCGGGT‑AAGCCTGTTCCTACAGG‑TTTTTTGCGTTGCCTGGCAATACCGCTTCATCCCATAAAACCAATAAAAGGAATTCAACCGTg                                                       >  4:256419/1‑149 (MQ=255)                                                                                            gTGATCTCACCCTGTAGGCGCGGGTTTACCCGCGAATGCGTCGGTACATTGACCGCCGCATTCGCGGGT‑AAGCCTGTTCCTACAGG‑TTTTTTGCGTTGCCTGGCAATACCGCTTCATCCCATAAAACCAATAAAAGGAATTCAACCGTg                                                       <  1:255453/149‑1 (MQ=255)                                                                                              gATCTCACCCTGTAGGCGCGGGTTTACCCGCGAATGCGTCGGTACATTGACCGCCGCATTCGCGGGT‑AAGCCTGTTCCTACAGG‑TTTTTTGCGTTGCCTGGCAATACCGCTTCATCCCATAAAACCAATAAAAGGAATTCAACCGTgg                                                      <  2:363432/148‑1 (MQ=255)                                                                                               aTCTCACCCTGTAGGCGCGGGTTTACCCGCGAATGCGTCGGTACATTGACCGCCGCATTCGCGGGT‑AAGCCTGTTCCTACAGG‑TTTTTTGCGTTGCCTGGCAATACCGCTTCATCCCATAAAACCAATAAAAGGAATTCAACCGTGGc                                                     <  3:256419/148‑1 (MQ=255)                                                                                                  tcACCCTGTAGGCGCGGGTTTACCCGCGAATGCGTCGGTACATTGACCGCCGCATTCGCGGGT‑AAGCCTGTTCCTACAGG‑TTTTTTGCGTTGCCTGGCAATACCGCTTCATCCCATAAAACCAATAAAAGGAATTCAACCGTGGCCGCt                                                 >  6:166472/1‑149 (MQ=255)                                                                                                  tcACCCTGTAGGCGCGGGTTTACCCGCGAATGCGTCGGTACATTGACCGCCGCATTCGCGGGT‑AAGCCTGTTCCTACAGG‑TTTTTTGCGTTGCCTGGCAATACCGCTTCATCCCATAAAACCAATAAAAGGAATTCAACCGTGGCCGCt                                                 >  3:164517/1‑149 (MQ=255)                                                                                                    aCCCTGTAGGCGCGGGTTTACCCGCGAATGTGTCGGTACATTGACCGCCGCATTCGCGGGT‑AAGCCTGTTCCTACAGG‑TTTTTTGCGTTGCCTGGCAATACCGCTTCATCCCATAAAACCAATAAAAGGAATTCAACCGTGGCCGCTg                                                >  3:358649/1‑148 (MQ=255)                                                                                                    aCCCTGTAGGCGCGGGTTTACCCGCGAATGCGTCGGTACATTGACCGCCGCATTCGCGGGT‑AAGCCTGTTCCTACAGG‑TTTTTTGCGTTGCCTGGCAATACCGCTTCATCCCATAAAACCAATAAAAGGAATTCAACCGTGGCCGCTg                                                >  4:313838/1‑148 (MQ=255)                                                                                                      ccTGTAGGCGCGGGTTTACCCGCGAATGCGTCGGTACATTGACCGCCGCATTCGCGGGT‑AAGCCTGTTCCTACAGG‑TTTTTTGCGTTGCCTGGCAATACCGCTTCATCCCATAAAAc                                                                             <  6:289596/117‑1 (MQ=255)                                                                                                      ccTGTAGGCGCGGGTTTACCCGCGAATGCGTCGGTACATTGACCGCCGCATTCGCGGGT‑AAGCCTGTTCCTACAGG‑TTTTTTGCGTTGCCTGGCAATACCGCTTCATCCCATAAAAc                                                                             >  5:289596/1‑117 (MQ=255)                                                                                                          tAGGCGCGGGTTTACCCGCGAATGTGTCGGTACATTGACCGCCGCATTCGCGGGT‑AAGCCTGTTCCTACAGG‑TTTTTTGCGTTGCCTGGCAATACCGCTTCATCCCATAAAACCAATAAAAGGAATTCAACCGTGGCCGCTGAAATTCa                                         <  4:358649/149‑1 (MQ=255)                                                                                                           attcgcgGGTAAACCCGCGAATGCGTCGGTACATTGACCGCCGCATTCGCGGGT‑AAGCCTGTTCCTACAGG‑TTTTTTGCGTTGCCTGGCAATACCGCTTCATCACATAAAACCAATAAAAGGAATTCAACCGTGGCCGCTGAAATTCa                                         >  8:155766/4‑148 (MQ=255)                                                                                                                        cccGCGAATGCGTCGGTACATTGACCGCCGCATTCGCGGGT‑AAGCCTGTTCCTACAGG‑TTTTTTGCGTTGCCTGGCAATACCGCTTCATCCCATAAAACCAATAAAAGGAATTCAACCGTGGCCGCTGAAATTCAAGACAGCCGCTCcg                           >  3:81175/1‑149 (MQ=255)                                                                                                                         ccGCGAATGCGTCGGTACATTGACCGCCGCATTCGCGGGT‑AAGCCTGTTCCTACAGG‑TTTTTTGCGTTGCCTGGCAATACCGCTTCATCCCATAAAACCAATAAAAGGAATTCAACCGTGGCCGCTGAAATTCAAGACAGCCGCTCcgc                          >  8:249900/1‑149 (MQ=255)                                                                                                                         ccGCGAATGCGTCGGTACATTGACCGCCGCATTCGCGGGT‑AAGCCTGTTCCTACAGG‑TTTTTTGCGTTGCCTGGCAATACCGCTTCATCCCATAAAACCAATAAAAGGAATTCAACCGTGGCCGCTGAAATTCAAGACAGCCGCTCagc                          >  2:44774/1‑149 (MQ=255)                                                                                                                         ccGCGAATGCGGCGGTCAATGTACCGACGCATTCGCGGGT‑AAGCCTGTTCCTACAGG‑TTTTTTGCGTTGCCTGGCAATACCGCTTCATCCCATAAAACCAATaaaa                                                                     <  4:231867/106‑1 (MQ=25)                                                                                                                         ccGCGAATGCGGCGGTCAATGTACCGACGCATTCGCGGGT‑AAGCCTGTTCCTACAGG‑TTTTTTGCGTTGCCTGGCAATACCGCTTCATCCCATAAAACCAATaaaa                                                                     >  3:231867/1‑106 (MQ=25)                                                                                                                         ccGCGAATGCGGCGGTCAATGTACCGACGCATTCGCGGGT‑AAACCCGCGCCTACAGG‑TTTTTTGCGTTGCCTGGCAATACCGCTTCATCCCATAAAACCAATAAAAGGAATTCAACCGTGGCCGCTGAAATTCaa                                        <  2:291953/135‑1 (MQ=21)                                                                                                                         ccGCGAATGCGGCGGTCAATGTACAGACGCATTCGCGGGT‑AAACCTGTTCCTACAGG‑TTTTTTGCGTTGCCTGGCAATACCGCTTCATCCCATAAAACCAATAAAAGGAATTAAAACGTGGCAGGTGAAATTCaa                                        >  7:54125/1‑135 (MQ=25)                                                                                                                               aTGCGGCGGTCAATGTACCGACGCATTCGCGGGT‑AAGCCTGTTCCTACAGG‑TTTTTTGCGTTGCCTGGCAATACCGCTGCATCCCATAAAACCAATAAAAGGAATTCAACCGTGGCCGCTGAAATTCAAGACAGCCGCTCCGCGCGCtt                    >  5:344348/1‑149 (MQ=38)                                                                                                                                tGCGTCGGTACATTGACCGCCGCATTCGCGGGT‑AAGCCTGTTCCTACAGG‑TTTTTTGCGTTGCCTGGCAATACCGCTTCATCCCATAAAACCAATAAAAGGAATTCAACCGTGGCCGCTGAAATTCAAGACAGCCGCTCCGCGCGCttt                   <  4:81175/149‑1 (MQ=255)                                                                                                                                tGCGTCGGTACATTGACCGCCGCATTCGCGGGT‑AAGCCTGTTCCTACAGG‑TTTTTTGCGTTGCCTGGCAATACCGCTTCATCCCATAAAACCAATAAAAGGAATTCAACCGTGGCCGCTGAAATTCAAGACAGCCGCTCCGCGCGCttt                   <  3:313838/149‑1 (MQ=255)                                                                                                                                  cGTCGGTACATTGACCGCCGCATTCGCGGGT‑AAGCCTGTTCCTACAGG‑TTTTTTGCGTTGCCTGGCAATACCGCTTCATCCCATAAAACCAATAAAAGGAATTCAACCGTGGCCGCTGAAATTCAAGACAGCCGCTCCGCGCGCTTTGc                 <  5:83283/149‑1 (MQ=255)                                                                                                                                   ttCGGTACATTGACCTCCGCATTCGCGGGT‑AAGCCTGTTCCTACAGG‑TTTTTTGCGTTGCCTGGCAATACCGCTTCATCCCATAAAACCAATAAAAGGAATTCAACCGTGGCCGCTGAAATTCAAGACTGCCGCTCCGCGCGCTTTGcc                <  6:312631/148‑1 (MQ=255)                                                                                                                                       gTACATTGACCGCCGCATTCGCGGGT‑AAGCCTGTTCCTACAGG‑TTTTTTGCGTTGCCTGGCAATACCGCTTCATCCCATAAAACCAATAAAAGGAATTCAACCGTGGCCGCTGAAATTCAAGACAGCCGCTCCGCGCGCTTTGCCTTgc            >  1:156395/1‑149 (MQ=255)                                                                                                                                            ttGACCGCCGCATTCGCGGGT‑AAGCCTGTTCCTACAGG‑TTTTTTGCGTTGCCTGGCAATACCGCTTCATCCCATAAAACCAATAAAAGGAATTCAACCGTGGCCGCTGAAATTCAAGACAGCCGCTCCGCGCGCTTTGCCTTGCgctgc       >  3:192268/1‑149 (MQ=255)                                                                                                                                              gcccgccgCACGCGCGGGT‑AAGCCGGTCCCTCCCGG‑TATTTTGCGTAGCCTCGCACTACCCCTTCATCCCATAACACCACTACAAGGAATTCAACCGCGGCCGCTGAAAATCAAGACAGCCGCTCCGCGCGCTTTGCCTTGCgctgcgc     <  8:273019/147‑3 (MQ=255)                                                                                                                                               tccgccgCATTCGGGGGT‑AATCCCGTTACTACAGGTTTTTTTTCCCAGCCTGGCAATACCGCTTCATCCCATAAAACCAATAAAAGGAATTCAACCGTGGCCGCTGAAATTCAAGACAGCCGCTCCGCGCGCTTTGCCTTGCGCTGCTc     <  1:301765/148‑1 (MQ=255)                                                                                                                                               aCCGCCGCATTCGCGGGT‑AAGCCTGTTCCTACAGG‑TTTTTTGCGTTGCCTGGCAATACCGCTTCATCCCATAAAACCAATAAAAGGAATTCAACCGTGGCCGCTGAAATTCAAGACAGCCGCTCCGCGCGCTTTGCCTTGCGCTGCTc     >  1:100191/1‑148 (MQ=255)                                                                                                                                                 cgccgcATTCGCGGGT‑AAGCCTGTTCCTACAGG‑TTTTTTGCGTTGCCTGGCAATACCGCTTCATCCCATAAAACCAATAAAAGGAATTCAACCGTGGCCGCTGAAATTCAAGACAGCCGCTCCGCGCGCTTTGCCTTGCGCTGCTCCaa  <  8:317574/149‑1 (MQ=255)                                                                                                                                                   |                                                                                                                                                     CTGGGGTCAGCGTCGAGGAAATCATCGAGAAAACCGCCGGCAAACTGATCGTGCCGGATGATGTGAAGGAAATGACCTTCTAAGGCATGACGTGATCTCACCCTGTAGGCGCGGGTTTACCCGCGAATGCGTCGGTACATTGACCGCCGCATTCGCGGGT‑AAGCCTGTTCCTACAGG‑TTTTTTGCGTTGCCTGGCAATACCGCTTCATCCCATAAAACCAATAAAAGGAATTCAACCGTGGCCGCTGAAATTCAAGACAGCCGCTCCGCGCGCTTTGCCTTGCGCTGCTCCAA  >  NC_002947/3534631‑3534923

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