Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,534,762	T→G	11.5%	intergenic (+49/‑106)	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,762	0	T	G	11.5%	148.5 / 9.3	58	intergenic (+49/‑106)	atoB/PP_3124	3‑oxoacid CoA‑transferase subunit B/short chain fatty acid transporter
Reads supporting (aligned to +/- strand):  ref base T (26/25);  new base G (4/3);  total (30/28)
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.79e-01

TCTGCGCGAGACCGCGCCTGGGGTCAGCGTCGAGGAAATCATCGAGAAAACCGCCGGCAAACTGATCGTGCCGGATGATGTGAAGGAAATG‑ACCTTCTAAGGCATGACGTGATCTCACCCTGTAGGCGCGGGTTTACCCGCGAATGCGTCGGTACATTGAC‑C‑GCCGCATTCGCGGGTAAGCCTGTTCCTACAGGTTTTTTGCGTTGCCTGGCAATACCGCTTCATCCCATAAAACCAATAAAAGGAATTCAACCGTGGCCGCTGAAATTCAAGACAGCCGCTCCGCGCGCTTT  >  NC_002947/3534614‑3534906                                                                                                                                                      |                                                                                                                                                   tCTGCGCGAGACCGCGCCTGGGGTCAGCGTCGAGGAAATCATCGAGAAAACCGCCGGCAAACTGATCGTGCCGGATGATGTGAAGGAAATG‑ACCTTCTAAGGCATGACGTGATCTCACCCTGTAGGCGCGGGTTTACCCGCGAATGCGt                                                                                                                                                    >  6:345024/1‑149 (MQ=255)    gcgcGAGACCGCGGCTGGGGTCATCGTCGAGGAAATCATCGATAAAACCCCCCGCAAACTGATCGTGCCGGATGTTGTGAAGGAAATGCCCCTTCTAAGGCATGACGTGATCTCACCCTGTAGGCGCGGGTTTACCCGCGAATGCGTc                                                                                                                                                   <  1:289037/148‑1 (MQ=255)    gcgcGAGACCGCGCCTGGGGTCAGCGTCGAGGAAATCATCGAGAAAACCGCCGGCAAACTGATCGTGCCGGATGATGTGAAGGAAATG‑ACCTTCTAAGGCATGACGTGATCTCACCCTGTAGGCGCGGGTTTACCCGCGAATGcggcgg                                                                                                                                                 >  1:230706/1‑149 (MQ=255)     cgcgAGACCGCGCCTGGGGTCAGCGTCGAGGAAATCATCGAGAAAACCGCCGGCAAACTGATCGTGCCGGATGATGTGAAGGAAATG‑ACCTTCTAAGGCATGACGTGATCTCACCCTGTAGGCGCGGGTTTACCCGCGAATGCGTCGGt                                                                                                                                                <  8:4779/149‑1 (MQ=255)        gagaCCGCGCCTGGGGTCAGCGTCGAGGAAATCATCGAGAAAACCGCCGGCAAACTGATCGTGCCGGATGATGTGAAGGAAATG‑ACCTTCTAAGGCATGACGTGATCTCACCCTGTAGGCGCGGGTTTACCCCCGAAAGCGTCGGTTCa                                                                                                                                             >  5:62670/1‑149 (MQ=255)           aCCGCGCCTGGGGTCAGCGTCGAGGAAATCATCGAGAAAACCGCCGGCAAACTGATCGTGCCGGATGATGTGAAGGAAATG‑ACCTTCTAAGGCATGACGTGATCTCACCCTGTAGGCGCGGGTTTACCCGCGAATGCGTCGGTACAtt                                                                                                                                           >  7:150257/1‑148 (MQ=255)                       gTCAGCGTCGAGGAAATCATCGAGAAAACCGCCGGCAAACTGATCGTGCCGGATGATGTGAAGGAAATG‑ACCTTCTAAGGCATGACGTGATCTCACCCTGTAGGAACAGGCTTACCCGCGAATGCGGCGGTCAATGTAC‑C‑GACGCATTc                                                                                                                            <  1:327784/149‑1 (MQ=25)                          aGCGTCGAGGAAATCATCGAGAAAACCGCCGGCAAACTGATCGTGCCGGATGATGTGAAGGAAATG‑ACCTTCTAAGGCGTGACGTGATCTCACCCTGTAGGCGCGGGTTTACCCGCGAATGCGTCGGTACATTGAC‑C‑GCCGCATTcgcg                                                                                                                         <  3:354104/149‑1 (MQ=255)                          aGCGTCGAGGAAATCATCGAGAAAACCGCCGGCAAACTGATCGTGCCGGATGATGTGAAGGAAATG‑ACCTTCTAAGGCGTGACGTGATCTCACCCTGTAGGCGCGGGTTTACCCGCGAATGCGTCGGTACATTGAC‑C‑GCCGCATTcgcg                                                                                                                         <  3:354100/149‑1 (MQ=255)                           tcGTCGAGGAAATCATCTAGAACACCGCCGGCAAACTGATCGTGCCGGACGCTCTGAAGGAAATC‑ACCTTCTAATGCATGACGTGATCTCACCCTGTAGGCGCGGGTTTACCCGCGAATGCGTCGGTACATTGAC‑C‑GCCGCATTCGCgg                                                                                                                        <  1:308698/148‑1 (MQ=255)                           gCGTCGAGGAAATCATCGAGAAAACCGCCGGCAAACTGATCGTGCCGGATGATGTGAAGGAAATG‑ACCTTCTAAGGCATGACGTGATCTCACCCTGTAGGCGCGGGTTTACCCGCGAATGCGGCGGTCAATGTAC‑C‑GACGCATTCGCgg                                                                                                                        <  2:230706/149‑1 (MQ=38)                            cGTCGAGGAAATCATCGAGAAAACCGCCGGCAAACTGATCGTGCCGGATGATGTGAAGGAAATG‑ACCTTCTAAGGCATGACGTGATCTCACCCTGTAGGCGCGGGTTTACCCGCGAATGCGTCGGTACATTGAC‑C‑GCCGCATTCGCggg                                                                                                                       <  2:12981/149‑1 (MQ=255)                             gTCGAGGAAATCATCGAGAAAACCGCCGGCAAACTGATCGTGCCGGATGATGTGAAGGAAATG‑ACCTTCTAAGGCATGACGTGATCTCACCCTGTAGGCGCGGGTTTACCCGCGAATGCGGCGGTCAATGTAC‑C‑GACGCATTcgcg                                                                                                                         <  3:75776/146‑1 (MQ=38)                             gTCGAGGAAATCATCGAGAAAACCGCCGGCAAACTGATCGTGCCGGATGATGTGAAGGAAATG‑ACCTTCTAAGGCATGACGTGATCTCACCCTGTAGGCGCGGGTTTACCCGCGAATGCGGCGGTCAATGTAC‑C‑GACGCATTcgcg                                                                                                                         >  4:75776/1‑146 (MQ=38)                                gAGGAAATCATCGAGAAAACCGCCGGCAAACTGATCGTGCCGGATGATGTGAAGGAAATG‑ACCTTCTAAGGCATGACGTGATCTCACCCTGTAGGCGCGGGTTTACCCGCGAATGCGTCGGTACATTGAC‑C‑GCCGCATTCGCGGGTaa                                                                                                                    <  7:276118/148‑1 (MQ=255)                                gAGGAAATCATCGAGAAAACCGCCGGCAAACTGATCGTGCCGGATGATGTGAAGGAAATG‑ACCTTCTAAGGCATGACGTGATCTCACCCTGTAGGCGCGGGTTTACACGCGAATGCGGCGGGACATTGAC‑C‑GCCGCATACGCGGGTacg                                                                                                                   >  2:19523/1‑147 (MQ=255)                                 aGGAAATCATCGAGAAAACCGCCGGCAAACTGATCGTGCCGGATGATGTGAAGGAAATG‑ACCTTCTAAGGCATGACGTGATCTCACCCTGTAGGCGCGGGTTTACCCGCGAATGCGTCGGTACATTGAC‑A‑GCAGCATTCTCGTGTaatc                                                                                                                  >  8:171216/1‑147 (MQ=255)                                      atcatcGAGTAAACCGCCGCCAAACTGATCGTGCCGGATGATGTGAAGGAAATC‑ACCTTCTAAGCCATGTCGTGATCTCACCCTGTAGGCGCGGGATTACCCGCGAATGCGTCGGTACATTGAC‑C‑GCCGCATTCGCGGGTAAGCCTGtt                                                                                                             <  1:19523/149‑1 (MQ=255)                                                    cgccgGCAAACTGCTCGTGCCGGATGATGTGAAGGAAATG‑ACCTTCTAAGGCATGACGTGATCTCACCCTTTAGGCGCTGGTTTACCCGCGAATGCGTCGGTACATTGAC‑C‑GCCGCATTCGCGGGTAAGCCTGTTCCTACAGGtttttt                                                                                               <  7:244853/149‑1 (MQ=255)                                                    cgccgGCAAACTGATCGTGCCGGATGATGTGAAGGAAATG‑AGCTTCTAAGGCATGACGTGATCTGACCCTGTAGGCGCGGGTTTACCCGCGAATGCGTCGGTACATTGAC‑C‑GCCGCTTTCGCGGGTAAGCCTGTTCCTACAGGtttttt                                                                                               <  7:278233/149‑1 (MQ=255)                                                                                  gAAGGAAATG‑ACCTTCTAAGGCATGACGTGATCTCACCCTGTAGGCGCGGGTTTACCCGCGAATGCGTCGGTACATTGAC‑C‑GCCGCATTCGCGGGTAAGCCTGTTCCTACAGGTTTTTTGCGTTGCCTgg                                                                                    <  7:11454/130‑1 (MQ=255)                                                                                  gAAGGAAATG‑ACCTTCTAAGGCATGACGTGATCTCACCCTGTAGGCGCGGGTTTACCCGCGAATGCGTCGGTACATTGAC‑C‑GCCGCATTCGCGGGTAAGCCTGTTCCTACAGGTTTTTTGCGTTGCCTgg                                                                                    >  8:11454/1‑130 (MQ=255)                                                                                   aaGGAAATG‑ACCTTCTAAGGCATGACGTGATCTCACCCTGTAGGCGCGGGTTTACCCGCGAATGCGTCGGTACATTGACTCAGCCGCATTCGCGGGTAAGCCTGTTCCTACAGGTTTTTTGCGTTGCCTGGCAATACCGCTTCATccc                                                                   <  1:298635/148‑1 (MQ=255)                                                                                       aaaTG‑ACCTTCTAAGGCATGACGTGATCTCACCCTGTAGGCGCGGGTTTACCCGCGAATGCGTCGGTACATTGAC‑C‑GCCGCATTCGCGGGTAAGCCTGTTCCTACAGGTTTTTTGCGTTGCCTGGCAATACCGCTTCATCCCATAAAAc                                                            <  5:345024/149‑1 (MQ=255)                                                                                        aaTG‑ACCTTCTAAGGCATGACGTGATCTCACCCTGTAGGCGCGGGTTTACCCGCGAATGCGTCGGTACATTGAC‑C‑GCCGCATTCGCGGGTAAGCCTGTTCCTACAGGTTTTTTGCGTTGCCTGGCAATACCGCTTCATCCCATAAAAcc                                                           <  8:221634/149‑1 (MQ=255)                                                                                          tG‑ACCTTCTAAGGCATGACGTGATCTCACCCTGTAGGGGCGTGTTTACCCGCGAATGCGTCGGTACATTGAC‑C‑GCCGCATTCGCGGGTAAGCCTGTTCCTACAGGTTTTTTGCGTTGCCTGGCAATACCGCTTCATCCCATAAAACCaa                                                         <  8:150257/149‑1 (MQ=255)                                                                                             aCCTTCTAAGGCATGACGTGATCTCACCCTGTAGGCGCGGGTTTACCCGCGAATGCGTCGGTACATTGAC‑C‑GCCGCATTCGCGGGTAAGCCTGTTCCTACAGGTTTTTTGCGTTGCCTGGCAATACCGCTTCATCCCATAAAACCAATa                                                       <  6:36381/149‑1 (MQ=255)                                                                                                ttCTAAGGCATGACGTGATCTCACCCTGTAGGCGCGGGTTTACCCGCGAATGCATCGGTACATTGAC‑C‑GCCGCATTCGCGGGTAAGCCTGTTCCTACAGGTTTTTTGCGTTGCCTGGCAATACCGCTTCATCCCATAAAACCAATaa                                                      >  6:22982/1‑147 (MQ=255)                                                                                                  cTAAGGCATGACGTGATCTCACCCTGTAGGCGCGGGTTTACCCGCGAATGCGTCGGTACATTGAC‑C‑GCCGCATTCGCGGGTAAGCCTGTTCCTACAGGTTTTTTGCGTTGCCTGGCAATACCGCTTCATCCCATAAAACCAATaaaa                                                    <  4:301919/147‑1 (MQ=255)                                                                                                    aaGGCATGACGTGATCTCACCCTGTAGGCGCGGGTTTACCCGCGAATGCGTCGGTACATTGAC‑C‑GCCGCATTCGCGGGTAAGCCTGTTCCTACAGGTTTTTTGCGTTGCCTGGCAATACCGCTTCATCCCATAAAACCAATAAAAGGaa                                                >  5:139684/1‑149 (MQ=255)                                                                                                      ggCATGACGTGATCTCACCCTGTAGGCGCGGGTTTACCCGCGAATGCGTCGGTACATTGAC‑C‑GCCGCATTCGCGGGTAAGCCTGTTCCTACAGGTTTTTTGCGTTGCCTGGCAATACCGCTTCATCCCATAAAACCAATAAAAGGAAt                                               >  3:160221/1‑148 (MQ=255)                                                                                                      ggCATGACGTGATCTCACCCTGTAGGCGCGGGTTTACCCGCGAATGCGTCGGTACATTGAC‑C‑GCCGCATTCGCGGGTAAGCCTGTTCCTACAGGTTTTTTGCGTTGCCTGGCAATACCGCTTCATCCCATAAAAACAATAAAAGGAAt                                               >  5:8658/1‑148 (MQ=255)                                                                                                        cATGACGTGATCTCACCCTGTAGGCGCGGGTTTACCCGCGAATGCGTCGGTACATTGAC‑C‑GCCGCATTCGCGGGTAAGCCTGTTCCTACAGGTTTTTTGCGTTGCCTGGCAATACCGCTTCATCCCATAAAACCAATAAAAGGAATTc                                             >  7:175542/1‑148 (MQ=255)                                                                                                              gTGATCTCACCCTGTAGGCGCGGGTTTACCCGCGAATGCGTCGGTACATTGAC‑C‑GCCGCATTCGCGGGTAAGCCTGTTCCTACAGGTTTTTTGCGTTGCCTGGCAATACCGCTTCATCCCATAAAACCAATAAAAGGAATTCAACCGTg                                      >  5:167350/1‑149 (MQ=255)                                                                                                              gTGATCTCACCCTGTAGGCGCGGGTTTACCCGCGAATGCGTCGGTACATTGAC‑C‑GCCGCATTCGCGGGTAAGCCTGTTCCTACAGGTTTTTTGCGTTGCCTGGCAATACCGCTTCATCCCATAAAACCAATAAAAGGAATTCAACCGTg                                      >  7:265166/1‑149 (MQ=255)                                                                                                                gATCTCACCCTGTAGGCGCGGGTTTACCCGCGAATGCGTCGGTACATTGAC‑C‑GCCGCATTCGCGGGTAAGCCTGTTCCTACAGGTTTTTTGCGTTg                                                                                         >  2:132284/1‑96 (MQ=255)                                                                                                                gATCTCACCCTGTAGGCGCGGGTTTACCCGCGAATGCGTCGGTACATTGAC‑C‑GCCGCATTCGCGGGTAAGCCTGTTCCTACAGGTTTTTTGCGTTg                                                                                         <  1:132284/96‑1 (MQ=255)                                                                                                                 aTCTCACCCTGTAGGCGCGGGTTTACCCGCGAATGCGTCGGTACATTGAC‑C‑GCCGCATTAGCGGGTAAGCCTGTTCCTACAGGTTTTTTGCGTTGCCTGGCAATACCGCTTCATCCCATAAAACCAATAAAAGGAATTCAACCGTGGcc                                   >  3:88098/1‑149 (MQ=255)                                                                                                                      aCCCTGTAGGCGCGGGTTTACCCGCGAATGCGTCGGTACATTGAC‑C‑GCCGCATTCGCGGGTAAGCCTGTTGCTACAGGGTTTTTGCga                                                                                           >  2:165977/1‑87 (MQ=255)                                                                                                                       cccTGTAGGCGCGGGTTTACCCGCGAATGCGTCGGTACATTGAC‑C‑GCCGCATTCGCGGGTAAGCCTGTTCCTACAGGTTTTTTGCGTTGCCTGGCAATACCGCTTCATCCCATAAAACCAATAAAAGGAATTCAACCGTGCCCGCTGa                              >  2:255896/1‑148 (MQ=255)                                                                                                                       cccTGTAGGCGCGGGTTTACCCGCGAATGCATCGGTACATTGAC‑C‑GCCGCATTCGCGGGTAAGCCTGTTCCTACAGGTTTTTTGCGTTGCCTGGCAATACCGCTTCATCCCATAAAACCAATAAAAGGAATTCAACCGTGGCCGCTGaa                             <  5:22982/149‑1 (MQ=255)                                                                                                                        ccTGTAGGCGCGGGTTTACCCGCGAATGCGTCGGTACATTGAC‑C‑GCCGCATTAGCGGGTAAGCCTGTTCCTACAGGTTTTTTGCGTTGCCTGGCAATACCGCTTCATCCCATAAAACCAATAAAAGGAATTCAACCGTGGCCGCTGaa                             >  4:231229/1‑148 (MQ=255)                                                                                                                          tGTAGGCGCGGGTTTACCCGCGAATGCGTCGGTACATTGAC‑C‑GCCGCATTCGCGGGTAAGCCTGTTCCTACAGGTTTTTTGCGTTGCCTGGCAATACCGCTTCATCCCATAAAACCAATAAAAGGAATTCAACCGTGGCCGCTGAAAt                           >  4:66624/1‑148 (MQ=255)                                                                                                                            tAGGCGCGGGTTTACCCGCGAATGCGTCGGTACATTGAC‑C‑GCCGCATTCGCGGGTAAGCCTGTTCCTACAGGTTTTTTGCGTTGCCTGGCAATACCGCTTCATCCCATAAAACCAATAAAAGGAATTCAACCGTGCCCGCTGAAAtt                          <  1:255896/147‑1 (MQ=255)                                                                                                                               gcgcGGGTTTACCCGCGAATGCGTCGGTACATTGAC‑C‑GCCGCATTCGCGGGTAAGCCTGTTCCTACAGGTTTTTTGCGTTGCCTGGCAATACCGCTTCATCCCATAAAACCAATAAAAGGAATTCAACCGTGGCCGCTGAAATTCAAGa                     <  8:175542/149‑1 (MQ=255)                                                                                                                                 acaggcttACCCGCGAATGCGGCGGTCAATGTAC‑C‑GACGCATTCGCGGGTAAGCCTGTTCCTACAGGTTTTTTGCGTTGCCTGGCAATACCGCTTCATCCCATAAAACCAATAAAAGGAATTCAACCGTGGCCGCTGAAATTCAAGa                     >  2:267005/7‑147 (MQ=37)                                                                                                                                   aggcttACCCGCGAATGCGGCGGTCAATGTAC‑C‑GACGCATTCGCGGGTAAGCCTGTTCCTACAGGTTTTTTGCGTTGCCTGGCAATACCGCTTCATCCCATAAAACCAATAAAAGGAATTCAACCGTGGCCGCTGAAATTCAAGACAGc                 >  8:51931/5‑149 (MQ=37)                                                                                                                                     gTTTACCCGCTAATGCGTCGGTACATTGAC‑C‑GCCGCATTCGCGGGTAAGCCTGTTCCTACAGGTTTTTTGCGTTGCCTGGCAATACCGCTTCATCCCATAAAACCAATAAAAGGAATTCAACCGTGGCCGCTGAAATTCAAGACAGc                 <  5:280596/147‑1 (MQ=255)                                                                                                                                     gTTTACCCGCGAATTCGTCGGTACATTGAC‑C‑GCCGCATTGGGGGGTAAGCCTGTTGCTACAGGTTTTTTGAGATGCCTGGCAATACCGCATCATCCCATAAAACCAATAAAAGGAATTAAACCGTGGCCGCTGAAATTAAAGAAAGCCg               >  6:224347/1‑149 (MQ=255)                                                                                                                                      tttACCCGCGAATGCGTCGGTACATTGAC‑C‑GCCGCATTCGCGGGTAAGCCTGTTCCTACAGGTTTTTTGCGTTGCCTGGCAATACCGCTTCATCCCATAAAACCAATAAAAGGAATTCAACCGTGGCCGCTGaaa                            >  2:196963/1‑135 (MQ=255)                                                                                                                                      tttACCCGCGAATGCGTCGGTACATTGAC‑C‑GCCGCATTCGCGGGTAAGCCTGTTCCTACAGGTTTTTTGCGTTGCCTGGCAATACCGCTTCATCCCATAAAACCAATAAAAGGAATTCAACCGTGGCCGCTGaaa                            <  1:196963/135‑1 (MQ=255)                                                                                                                                          cccGCGAATGCGTCGGTACATTGAC‑C‑GCCGCATTCGCGGGTAAGCCTGTTCCTACAGGTTTTTTGCGTTGCCTGGCAATACCGCTTCATACCATAAAACCAATAAAAGGAATTCAACCGTGGCCGCTGAAATTCAAGACAGCCGCTCcg          >  1:113750/1‑149 (MQ=255)                                                                                                                                           ccGCGAATGCGTCGGTACATTGAC‑C‑GCCGCATTCGCGGGTAAGCCTGTTCCTACAGGTTTTTTGCGTTGCCTGGCAATACCGCTTCATCCCATAAAACCAATAAAAGGAATTCAACCGTGGCCGCTGAAATTCAAGACAGCCGCTCcgc         >  4:335367/1‑149 (MQ=255)                                                                                                                                           ccGCGAATGCGTCGGTACATTGAC‑C‑GCCGCATTCGCGGGTAAGCCTCTTCCTACAGGTTTTTTGCGTTGCCTGGCAATACCGCTTCATCCCATAAAACCAATAAAAGGAATTCAACCGTGCCAGCTGAAATTCAAGAAAGCCGCTCcgc         >  6:73465/1‑149 (MQ=255)                                                                                                                                            cgcgAATGCGTCGGTACATTGAC‑C‑GCCGCATTCGCGGGTAAGCCTGTTCCTACAGGTTTTTTGCGTTGCCTGGCAATACCGCTTCATCCCATAAAACCAATAAAAGGAATTCAACCGTGGCCGCTGAAATTCAAGACAGCCGCTCcgcg        >  3:88131/1‑149 (MQ=255)                                                                                                                                            cgcgAATGCGTCGGTACATTGAC‑C‑GCCGCATTCGCGGGTAAGCCTGTTCCTACAGGTTTTTTGCGTTGCCTGGCAATACCGCTTCATCCCATAAAACCAATAAAAGGAATTCAACCGTGGCCGCTGAAATTCAAGACAGCCGCTCcgcg        >  1:131735/1‑149 (MQ=255)                                                                                                                                            cgcaAATGCGTCGGTACATTGAC‑C‑GCCGCATTCGCGGGTAAGCCTGTTCCTACAGGTTTTTTGCGTTGCCTGGCAATACCGCTTCATCCCATAAAACCAATAAAAGGAATTCAACCGTGGCCGCTGAAATTCAAGACAGCCGCTCcgcg        >  6:73258/1‑149 (MQ=255)                                                                                                                                             gcgAATGCGTCGGTACATTGAC‑C‑GCCGCATTCGCGGGTAAGCCTGTTCCTACAGGTTTTTTGCGTTGCCTGGCAATACCGCTTCATCCCATAAATCCAATAAAAGGAATTCAACCGTGGCCGCTGAAATTCAAGACAGCCGCTCcgcgc       >  3:222776/1‑149 (MQ=255)                                                                                                                                                aaTGCGTCGGTACATTGAC‑C‑GCCGCATTCGCGGGTAAGCCTGTTCCTACAGGTTTTTTGCGTTGCCTGGCAATACCGCTTCATCCCATAAAACCAATAAAAGGAATTCAACCGTGGCCGCTGAAATTCAAGACAGCCGCTCCGCGCGCt    <  7:171216/149‑1 (MQ=255)                                                                                                                                                  tGCGTCGGTACATTGAC‑C‑GCCGCATTCGCGGGTAAGCCTGTTCCTACAGGTTTTTTGCGTTGCCTGGCAATACCGCTTCATCCCATAAAACCAATAAAAGGAATTCAACCGTGGCCGCTGAAATTCAAGACAGCCGCTCCGCGCGCttt  <  4:160221/149‑1 (MQ=255)                                                                                                                                                      |                                                                                                                                                   TCTGCGCGAGACCGCGCCTGGGGTCAGCGTCGAGGAAATCATCGAGAAAACCGCCGGCAAACTGATCGTGCCGGATGATGTGAAGGAAATG‑ACCTTCTAAGGCATGACGTGATCTCACCCTGTAGGCGCGGGTTTACCCGCGAATGCGTCGGTACATTGAC‑C‑GCCGCATTCGCGGGTAAGCCTGTTCCTACAGGTTTTTTGCGTTGCCTGGCAATACCGCTTCATCCCATAAAACCAATAAAAGGAATTCAACCGTGGCCGCTGAAATTCAAGACAGCCGCTCCGCGCGCTTT  >  NC_002947/3534614‑3534906

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑