Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	163,978	Δ1 bp	6.7%	intergenic (‑145/+42)	scpC ← / ← pntB	propionyl‑CoA:succinate CoA transferase/pyridine nucleotide transhydrogenase subunit beta

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	163,978	0	A	.	6.7%	190.1 / 7.0	60	intergenic (‑145/+42)	scpC/pntB	propionyl‑CoA:succinate CoA transferase/pyridine nucleotide transhydrogenase subunit beta
Reads supporting (aligned to +/- strand):  ref base A (25/31);  new base . (2/2);  total (27/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.59e-01

ATGGATTGTTATCTCGGGCTACTGAAACTACTGCAGCGCAGTCTAGAGATTTCGCCCATTGCCGTCCCACGACCATGGTCGTATGAGAAGCCTTGTATTTGAGTCGTTGATCCGGATCAATAAAAGAAAAGCCCCGG‑CAGATTCTGACCGGGGC‑TTCCTTTATATCAGCTGGTTTGCAGCGGCTGGTGTCACTCCACCGCCTTGACCATGTCCTCGATGACCTTCTTGGCGTCACCGAACACCATCATGGTCTTGTCCAGATAGAACAGTTCGTTGTCCAGGCCGGCATAG  >  NC_002947/163832‑164122                                                                                                                                                    |                                                                                                                                                  aTGGATTGTTATCTCGGGCTACTGAAACTACTGCAGCGCAGTCTAGAGATTTCGCCCATTGCCGTCCCACGACCATGGTCGTATGAGAAGCCTTGTATTTGAGTCGTTGATCCGGATCAATAAAAGAAAAGCCCCGG‑CAGATTCTGAcc                                                                                                                                                 >  7:251539/1‑149 (MQ=255)     aTTGTTATCTCGGGCTACTGAAACTACTGCAGCGCAGTCTAGAGATTTCGCCCATTGCCGTCCCACGACCATGGTAGTATGAGAAGCCTTGTATTTGAGTCGTTGATCCGGATCAATAAAAGAAAAGCCCCGG‑CAGATTCTGACCgggg                                                                                                                                             >  5:297837/1‑149 (MQ=255)        gTTATCTCGGGCTACTGAAACTACTGCAGCGCAGTCTAGAGATTTCGCCCATTGCCGTCCCACGACCATGGTCGTATGAGAAGCCTTGTATTTGAGTCGTTGATCCGGATCAATAAAAGAAAAGCCCCGGTCAGAATCTG‑CCGGGGC‑tt                                                                                                                                         >  5:129122/1‑149 (MQ=255)         ttATCTCGGGCTACTGAAACTACTGCAGCGCAGTCTAGAGATTTCGCCCATTGCCGTCCCACGACCATGGTCGTATGAGAAGCCTTGTATTTGAGTCGTTGATCCGGATCAATAAAAGAAAAGCCCCGGTCAGAATCTG‑CCGGGGC‑ttt                                                                                                                                        <  3:302507/149‑2 (MQ=255)             ctcGGGCTACTGAAACTACTGCAGCGCAGTCTAGAGATTTCGCCCATTGCCGTCCCACGACCATGGTCGTATGAGAAGCCTTGTATTTGAGTCGTTGATCCGGATCAATAAAAGAAAAGCCCCGG‑CAGATTCTGACCGGGGC‑TTCCtt                                                                                                                                     >  2:90352/1‑148 (MQ=255)                gggCTACTGAAACTACTGCAGCGCAGTCTAGAGATTTCGCCCATTGCCGTCCCACGACCATGGTCGTATGAGAAGCCTTGTATTTGAGTCGTTGATCCGGATCAATAAAAGAAAAGCCCCGGTCAGAATCTG‑CCGGGGC‑ttttcttt                                                                                                                                   <  6:129122/147‑7 (MQ=255)                gggCTACTGAAACTACTGCAGCGCAGTCTAGAGATTTCGCCCATTGCCGTCCCACGACCATGGTCGTATGAGAAGCCTTGTATTTGAGTCGTTGATCCGGATCAATAAAAGAAAAGCCCCGG‑CAGATTCTGACCGGGGC‑TTCCtt                                                                                                                                     <  2:249890/145‑1 (MQ=255)                 ggCTACTGAAACTACTGCAGCGCAGTCTAGAGATTTCGCCCATTGCCGTCCCACGACCATGGTCGTATGAGAAGCCTTGTATTTGAGTCGTTGATCAGGATCAATAAAAGAAAAGCCCCGG‑CAGATTCTGACCGGGGC‑TTCCTTAatat                                                                                                                                >  6:221203/1‑149 (MQ=255)                      cTGAAACTACTGCAGCGCAGTCTAGAGATTTCGCCCATTGCCGTCCCACGACCATGGTCGTATGAGAAGCCTTGTATTTGAGTCGTTGATCCGGATCAATAAAAGAAAAGCCCCGG‑CAGATTCTGACCGGGGC‑TTCCTTTATATCAg                                                                                                                             <  7:214797/147‑1 (MQ=255)                      cTGAAACTACTGCAGCGCAGTCTAGAGATTTCGCCCATTGCCGTCCCACGACCATGGTCGTATGAGAAGCCTTGTATTTGAGTCGTTGATCCGGATCAATAAAAGAAAAGCCCCGG‑CAGATTCTGACCGGGGC‑TTCCTTTATATCAGCt                                                                                                                           <  5:240222/149‑1 (MQ=255)                      cTGAAACTACTGCAGCGCAGTCTAGAGATTTCGCCCATTGCCGTCCCACGACCATGGTCGTATGAGAAGCCTTGTATTTGAGTCGTTGATCCGGATCAATAAAAGAAAAGCCCCGG‑CAGATTCTGACCGGGGC‑TTCCTTTATATCAGCt                                                                                                                           >  6:49216/1‑149 (MQ=255)                      cTGAAACTACTGCAGCGCAGTCTAGAGATTTCGCCCATTGCCGTCCCACGACCATGGTCGTATGAGAAGCCTTGTATTTGAGTCGTTGATCCGGATCAATAAAAGAAAAGCCCCGG‑CAGATTCTGACCGGGGC‑TTCCTTTATATCAGCt                                                                                                                           >  1:181841/1‑149 (MQ=255)                        gAAACTACTGCAGCGCAGTCTAGAGATTTCGCCCATTGCCGTCCCACGACCATGGTCGTATGAGAAGCCTTGTATTTGAGTCGTTGATCCGGATCAATAAAAGAAAAGCCCCGG‑CAGATTCTGACCGGGGC‑TTCCTTTATATCAGCTgg                                                                                                                         >  6:240636/1‑149 (MQ=255)                           actactGCAGCGCAGTCTAGAGATTTCGCCCATTGCCGTCCCACGACTATGGTCGTATGAGAAGCCTTGTATTTGAGTCGTTGATCCGGATCAATAAAAGAAAAGCCCCGG‑CAGATTCTGACCGGGGC‑TTCCTTTATATCAGCTGGttt                                                                                                                      >  2:215832/1‑149 (MQ=255)                           actactGCAGCGCAGTCTAGAGATTTCGCCCATTGCCGTCCCACGACTATGGTCGTATGAGAAGCCTTGTATTTGAGTCGTTGATCCGGATCAATAAAAGAAAAGCCCCGG‑CAGATTCTGACCGGGGC‑TTCCTTTATATAAGCTGGttt                                                                                                                      >  2:211510/1‑149 (MQ=255)                           actactGCAGCGCAGTCTAGAGATTTCGCCCATTGCCGTCCCACGACCATGGTCGTATGAGAAGCCTTGTATTTGAGTCGTTGATCCGGATCAATAAAAGAAAAGCCCCGG‑CAGATTCTGACCGGGGC‑TTCCTTTATATCAGCTGGttt                                                                                                                      <  4:186562/149‑1 (MQ=255)                           actactGCAGCGCAGTCTAGAGATTTCGCCCATTGCCGTCCCACGACCATGGTCGTATGAGAAGCCTTGTATTTGAGTCGTTGATCCGGATCAATAAAAGAAAAGCCCCGG‑CAGATTCTGACCGGGGC‑TTCCTTTATATCAGCTGGttt                                                                                                                      >  3:277271/1‑149 (MQ=255)                           actactGCAGCGCAGTCTAGAGATTTCGCCCATTGCCGTCCCACGACCATGGTCGTATGAGAAGCCTTGTATTTGAGTCGTTGATCCGGATCAATAAAAGAAAAGCCCCGG‑CAGATTCTGACCGGGGC‑TTCCTTTATATCAGCTGGtt                                                                                                                       >  5:55034/1‑148 (MQ=255)                                  cAGCGCAGTCTAGAGATTTCGCCCATTGCCGTCCCACGACCATGGTCGTATGAGAAGCCTTGTATTTGAGTCGTTGATCCGGATCAATAAAAGAAAAGCCCCGG‑CAGATTCTGACCGGGGC‑TTCCTTTATATCAGCTGGTTTGCAGCg                                                                                                                <  7:280784/148‑1 (MQ=255)                                       cAGTCTAGAGATTTCGCCCATTGCCGTCCCACGACCATGGTCGTATGAGAAGCCTTGTATTTGAGTCGTTGATCCGGATCAATAAAAGAAAAGCCCCGG‑CAGATTCTGACCGGGGC‑TTCCTTTATATCAGCTGGTTTGCAGCGGCTGgt                                                                                                          <  6:184921/149‑1 (MQ=255)                                       cAGTCTAGAGATTTCGCCCATTGCCGTCCCACGACCATGGTCGTATGAGAAGCCTTGTATTTGAGTCGTTGATCCGGATCAATAAAAGAAAAGCCCCGG‑CAGATTCTGACCGGGGC‑TTCCTTTATATCAGCTGGTTTGCAGCGGCTGgt                                                                                                          <  5:49216/149‑1 (MQ=255)                                            tAGAGATTTCGCCCATTGCCGTCCCACGACCATGGTCGTATGAGAAGCCTTGTATTTGAGTCGTTGATCCGGATCAATAAAAGAAAAGCCCCGG‑CAGATTCTGACCGGGGC‑TTCCTTTATATCAGCTGGTTTGCAGCGGCTGGTGTCAc                                                                                                     >  5:2649/1‑149 (MQ=255)                                               agaTTTCGCCCATTGCCGTCCCACGACCATGGTCGTATGAGAAGCCTTGTATTTGAGTCGTTGATCCGGATCAATAAAAGAAAAGCCCCGG‑CAGATTCTGACCGGGGC‑TTCCTTTATATCAGCTGGTTTGCAGCGGCTGGTGTCACTc                                                                                                   >  7:281912/1‑148 (MQ=255)                                               agaTTTCGCCCATTGCCGTCCCACGACCATGGTCGTATGAGAAGCCTTGTATTTGAGTCGTTGATCCGGATCAATAAAAGAAAAGCCCCGG‑CAGATTCTGACCGGGGC‑TTCCTTTATATCAGCTGGTTTGCAGCGGCTGGTGTCACTc                                                                                                   <  8:281912/148‑1 (MQ=255)                                                gaTTTCGCCCATTGCCGTCCCACGACCATGGTCGTATGAGAAGCCTTGTATTTGAGTCGTTGATCCGGATCAATAAAAGAAAAGCCCCGG‑CAGATTCTGACCGGGGC‑TTCCTTTATATCAGCTGGTTTGCAGCGGCTGGTGTCACTCCa                                                                                                 <  4:53946/149‑1 (MQ=255)                                                  tttCGCCCATTGCCGTCCCACGACCATGGTCGTATGAGAAGCCTTGTATTTGAGTCGTTGATCCGGATCAATAAAAGAAAAGCCCCGG‑CAGATTCTGACCGGGGC‑TTCCTTTATATCAGCTGGTTTGCAGCGGCTTGTGTCACTCCAcc                                                                                               <  6:297837/149‑1 (MQ=255)                                                           ttGCCGTCCCACGACCATGGTCGTATGAGAAGCCTTGTATTTGAGTCGTTGATCCGGATCAATAAAAGAAAAGCCCCGG‑CAGATTCTGACCGGGGC‑TTCCTTTATATCAGCTGGTTTTCAGCGGCTGGTGTCACTCCACCGCCTTGAcc                                                                                      <  6:55034/149‑1 (MQ=255)                                                              ccGTCCCACGACCATGGTCGTATGAGAAGCCTTGTATTTGAGTCGTTGATCCGGATCAATAAAAGAAAAGCCCCGG‑CAGATTCTGACCGGGGC‑TTCCTTTATATCAGCTGGTTTGCAGCGGCTGGTGTCACTCCACCGCCTTGACCATg                                                                                   >  7:15660/1‑149 (MQ=255)                                                                                 gTATGAGAAGCCTTGTATTTGAGTCGTTGATCCGGATCAATAAAAGAAAAGCCCCGG‑CAGATTCTGACCGGGGC‑TTCCTTTATATCAGCTGGTTTGCAGCGGCTGGTGTCACTCCACCGCCTTGACCATGTCCTCGATGACCTTCTTg                                                                 <  2:168803/148‑1 (MQ=255)                                                                                  tATGAGAAGCCTTGTATTTGAGTCGTTGATCCGGATCAATAAAAGAAAAGCCCCGG‑CAGATTCTGACCGGGGC‑TTCCTTTatat                                                                                                                                <  4:154650/84‑1 (MQ=255)                                                                                  tATGAGAAGCCTTGTATTTGAGTCGTTGATCCGGATCAATAAAAGAAAAGCCCCGG‑CAGATTCTGACCGGGGC‑TTCCTTTatat                                                                                                                                >  3:154650/1‑84 (MQ=255)                                                                                       gaAGCCTTGTATTTGAGTCGTTGATCCGGATCAATAAAAGAAAAGCCCCGG‑CAGATTCTGACCGGGGC‑TTCCTTTATATCAGCTGGTTTGCAGCGGCTGGTGTCACTCCACCGCCTTGACCATGTCCTCGATGACCTTCTTGGCGTCAc                                                          >  7:150774/1‑149 (MQ=255)                                                                                       gaAGCCTTGTATTTGAGTCGTTGATCCGGATCAATAAAAGAAAAGCCCCGG‑CAGATTCTGACCGGGGC‑TTCCTTTATATCAGCTGGTTTGCAGCGGCTGGTGTCACTCCACCGCCTTGACCATGTCCTCGATGACCTTCTTGGCGTCAc                                                          <  2:181841/149‑1 (MQ=255)                                                                                            cTTGTATTTGAGTCGTTGATCCGTATCAATAAAAGAAAAGCCCCGG‑CAGATTCTGACCGGGGC‑TTCCTTTATATCAGCTGTTTTGCAGCGGCTGGTGTCACTCCACCGCCTTGACCATGTCCTCGATGACCTTCTTGGCGTCACCGAac                                                     <  8:251539/149‑1 (MQ=255)                                                                                             ttGTATTTGAGTCGTTGATCCGGATCAATAAAAGAAAAGCCCCGG‑CAGATTCTGACCGGGGC‑TTCCTTTATATCAGCTGGTTTGCAGCGGCTGGTGTCACTCCACCGCCTTGACCATGTCCTCGATGACCTTCTTGGCGTCACCGAaca                                                    >  8:186145/1‑149 (MQ=255)                                                                                             ttGTATTTGAGTCGTTGATCCGGATCAATAAAAGAAAAGCCCCGG‑CAGATTCTGACCGGGGC‑TTCCTTTATATCAGCTGGTTTGCAGCGGCTGGTGTCACTCCACCGCCTTGACCATGTCCTCGATGACCTTCTTGGCGTCACCGAaca                                                    >  2:250279/1‑149 (MQ=255)                                                                                             ttGTATTTGAGTCGTTGATCCGGATCAATAAAAGAAAAGCCCCGG‑CAGATTCTGACCGGGGC‑TTCCTTTATATCAGCTGGTTTGCAGCGGCTGGTGTCACTCCACCGCCTTGACCATGTCCTCGATGACCTTCTTGGCGTCACCGAaca                                                    >  2:265306/1‑149 (MQ=255)                                                                                             ttGTATTTGAGTCGTTGATCCGGATCAATAAAAGAAAAGCCCCGG‑CAGATTCTGACAGGGGC‑TTCCTTTATATAAGCTGGTTTGCAGCGGCTGGTGTCACTCCACCGCCTTGACCAt                                                                                    >  3:162412/1‑117 (MQ=255)                                                                                                  tttGAGTCGTTGATCCGGATCAATAAAAGAAAAGCCCCGG‑CAGATTCTGACCGGGGC‑TTCCTTTATATGAGCTGGTTTGCAGCGGCTGGTGTCACTCCACCGCCTTGACCATGTCCTCGATGACCTTCTTGGCGTCACCGAACACcat                                                >  3:287162/1‑148 (MQ=255)                                                                                                       gTCGTTGATCCGGATCAATAAAAGAAAAGCCCCGG‑CAGATTCTGACCGGGGC‑TTCCTTTATATCAGCTGGTTTGCAGCGGCTGGTGTCACTCCACCGCCTTGACCATGTCCTCGATGACCTTCTTGGCGTCACCGAACACCATCATGGt                                          <  3:302889/149‑1 (MQ=255)                                                                                                        tCGTTGATCCGGATCAATAAAAGAAAAGCCCCGG‑CAGATTCTGACCGGGGC‑TTCCTTTATATCAGCTGGTTTGCAGCGGCTGGTGTCACTCCACCGCCTTGACCATGTCCTCGATGACCTTCTTGGCGTCACCGAACACCATCATGGTc                                         <  5:240636/149‑1 (MQ=255)                                                                                                        tCGTTGATCCGGATCAATAAAAGAAAAGCCCCGG‑CAGATTCTGACCGGGGC‑TTCCTTTATATCAGCTGGTTTGCAGCGGCTGGTGTCACTCCACCGCCTTGACCATGTCCTCGATGACCTTCTTGGCGTCACCGAACACCATCATGGTc                                         <  2:154219/149‑1 (MQ=255)                                                                                                         cGTTGATCCGGATCAATAAAAGAAAAGCCCCGG‑CAGATTCTGACCGGGGC‑TTCCTTTATATCAGCTGGTTTGCAGCGGCTGGTGTCACTCCACCGCCTTGACCATGTCCTCGATGACCTTCTTGGCGTCACCGAACACCATCATGGTc                                         <  8:150774/148‑1 (MQ=255)                                                                                                          gTTTATCCTGATCAATAAAAGAATAGCCCGGG‑CAGATTCTGACCGGGGC‑TTCCTTTATATCAGCTGGTTTGCAGCGGCTGGTGTCACTCCTCGGCCTTGACCATGTCCTCGATGACCTTCTTGGCGTCACCGAACACCATCATGGTCtt                                       <  5:221203/149‑1 (MQ=255)                                                                                                          gTTGATCCGGATCAATAAAAGAAAAGCCCCGG‑CAGATTCTGACCGGGGC‑TTCCTTTATATCAGCTGGTTTGCAGCGGCTGGTGTCACTCCACCGCCTTGACCATGTCCTCGATGACCTTCTTGGCGTCACCGAACACCATCATGGTCtt                                       <  6:2649/149‑1 (MQ=255)                                                                                                                 cGGATCAATAAAAGAAAAGCCCCGG‑CAGATTCTGACCGGGGC‑TTCCTTTATATCAGCTGGTTTGCAGCGGCTGGTGTCACTCCACCGCCTTGACCATGTCCTCGATGACCTTCTTGGCGTCACCGAACACCATCATGGTCTTGTCCAGa                                <  1:285028/149‑1 (MQ=255)                                                                                                                  ggATCAATAAAAGAAAAGCCCCGG‑CAGATTCTGACCGGGGC‑TTCCTTTATATCAGCTGGTTTGCAGCGGCTGGTGTCACTCCACCGCCTTGACCATGTCCTCGATGACCTTCTTGGCGTCACCGAACACCATCATGGTCTTGTCCTGAt                               <  2:126173/149‑1 (MQ=255)                                                                                                                  ggATCAATAAAAGAAAAGCCCCGG‑CAGATTCTGACCGGGGC‑TTCCTTTATATCAGCTGGTTTGCAGCGGCTGGTGTCACTCCACCGCCTTGACCATGTCCTCGATGACCTTCTTGGCGTCACCGAACACCATCATGGTCTTGTCCAGAt                               <  4:277271/149‑1 (MQ=255)                                                                                                                  ggATCAATAAAAGAAAAGCCCCGG‑CAGATTCTGACCGGGGC‑TTCCTTTATATCAGCTGGTTTGCAGCGGCTGGTGTCACTCCACCGCCTTGACCATGTCCTCGATGACCTTCTTGGCGTCACCGAACACCATCATGGTCTTGTCCAGAt                               <  7:186145/149‑1 (MQ=255)                                                                                                                   gATCAATAAAAGAAAAGCCCCGG‑CAGATTCTGACCGGGGC‑TTCCTTTATATCAGCTGGTTTGCAGCGGCTGGTGTCACTCCACCGCCTTGACCATGTCCTCGATGACCTTCTTGGCGTCACCGAACACCATCATGGTCTTGTCCAGATa                              <  1:211510/149‑1 (MQ=255)                                                                                                                   gATCAATAAAAGAAAAGCCCCGG‑CAGATTCTGACCGGGGC‑TTCCTTTATATCAGCTGGTTTGCAGCGGCTGGTGTCACTCCACCGCCTTGACCATGTCCTCGAAGACCTTCTTGGCGTCACCGAACACCATCATGGTCTTGTCCAGATa                              <  1:215832/149‑1 (MQ=255)                                                                                                                        aTAAAAGAAAAGCCCCGG‑CAGATTCTGACCGGGGC‑TTCCTTTATATCAGCTGGTTTGCAGCGGCTGGTGTCACTCCACCGCCTTGACCATGTCCTCGATGACCTTCTTGGCGTCACCGAACACCATCATGGTCTTGTCCAGATAGaa                           >  5:11587/1‑147 (MQ=255)                                                                                                                          aaaagaaaagCCCCGG‑CAGATTCTGACCGGGGC‑TTCCTTTATATCAGCTGGTTTGCAGCGGCTGGTGTCACTCCACCGCCTTGACCATGTCCTCGATGACCTTCTTGGCGTCACCGAACACCATCATGGTCTTGTCCAGATAGAACAGt                       >  1:145024/1‑149 (MQ=255)                                                                                                                           aaagaaaagCCCCGG‑CAGATTCTGACCGGGGC‑TTCCTTTATATCAGCTGGTTTGCAGCGGCTGGTGTCACTCCACCGCCTTGACCATGTCCTCGATGACCTTCTTGGCGTCACCGAACACcatc                                               <  2:106209/124‑1 (MQ=255)                                                                                                                             agaaaagCCCCGG‑CAGATTCTGACCGGGGC‑TTCCTTTATATCAGCTGGTTTGCAGCGGCTGGTGTCACTCCACCGCCTTGACCATGTCCTCGATGACCTTCTTGGCGTCACCGAACACCATCATGGTCTTGTCCAGATAGAACAGTTCg                    <  2:145024/149‑1 (MQ=255)                                                                                                                               aaaagCCCCGG‑CAGATTCTGACCGGGGC‑TTCCTTTATATCAGCTGGTTTGCAGCGGCTGGTGTCACTCCACCGCCTTGACCATGTCCTCGATGACCTTCTTGGCGTCACCGAACACCATCATGGTCTTGTCCAGATAGAACAGTTCGtt                  <  2:98836/149‑1 (MQ=255)                                                                                                                                  agCCCCGG‑CAGATTCTGACCGGGGC‑TTCCTTTATATCAGCTGGTTTGCAGCGGCTGGTGTCACTCCACCGCCTTGACCATGTCCTCGATGACCTTCTTGGCGTCACCGAACACCATCATGGTCTTGTCCAGATAGAACAGTTCGTTGTc               <  1:90352/149‑1 (MQ=255)                                                                                                                                    ccccGGTCAGAATCTG‑CCGGGGCTTTTCTTTATATCAGCTGGTTTGCAGCGGCTGGTGTCACTCCACCGCCTTGACCATGTCCTCGATGACCTTCTTGGCGTCACCGAACACCATCATGGTCTTGTCCAGATAGAACAGTTCGTTGTc               >  8:40468/1‑148 (MQ=255)                                                                                                                                        gg‑CAGATTCTGACCGGGGC‑TTCCTTTATATCAGCTGGTTTGCAGCGGCTGGTGTCACTCCACCGCCTTGACCATGTCCTCGATGACCTTCTTGGCGTCACCGAACACCATCATGGTCTTGTCCAGATAGAACAGTTCGTTGTCCAGGcc         >  1:179367/1‑149 (MQ=255)                                                                                                                                                tCTGACCGGGGC‑TTCCTTTATATCAGCTGGTTTGCAGCGGCTGGTGTCACTCCACCGCCTTGACCATGTCCTCGATGACCTTCTTGGCGTCACCGAACACCATCATGGTCTTGTCCAGATAGAACAGTTCGTTGTCCAGGCCGGCATAg  >  6:128896/1‑149 (MQ=255)                                                                                                                                                    aCCGGGGC‑TTCCTTTATATCAGCTGGTTTGCAGCGGCTGGTGTCACTCCACCGCCTTGACCATGTCCTCGATGACCTTCTTGGCGTCACCGAACACcatc                                               >  6:128555/1‑100 (MQ=255)                                                                                                                                                    aCCGGGGC‑TTCCTTTATATCAGCTGGTTTGCAGCGGCTGGTGTCACTCCACCGCCTTGACCATGTCCTCGATGACCTTCTTGGCGTCACCGAACACcatc                                               <  5:128555/100‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                                  ATGGATTGTTATCTCGGGCTACTGAAACTACTGCAGCGCAGTCTAGAGATTTCGCCCATTGCCGTCCCACGACCATGGTCGTATGAGAAGCCTTGTATTTGAGTCGTTGATCCGGATCAATAAAAGAAAAGCCCCGG‑CAGATTCTGACCGGGGC‑TTCCTTTATATCAGCTGGTTTGCAGCGGCTGGTGTCACTCCACCGCCTTGACCATGTCCTCGATGACCTTCTTGGCGTCACCGAACACCATCATGGTCTTGTCCAGATAGAACAGTTCGTTGTCCAGGCCGGCATAG  >  NC_002947/163832‑164122

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