Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	163,973	T→A	5.6%	intergenic (‑140/+47)	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,973	0	T	A	5.6%	220.5 / 2.4	72	intergenic (‑140/+47)	scpC/pntB	propionyl‑CoA:succinate CoA transferase/pyridine nucleotide transhydrogenase subunit beta
Reads supporting (aligned to +/- strand):  ref base T (34/33);  new base A (2/2);  total (36/36)
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.93e-01

ACGGTACATGGATTGTTATCTCGGGCTACTGAAACTACTGCAGCGCAGTCTAGAGATTTCGCCCATTGCCGTCCCACGACCATGGTCGTATGAGAAGCCTTGTATTTGAGTCGTTGATCCGGATCAATAAAAGAAAAGCCCCGG‑CAGATTCTGACCGGGGCTTCCTTTATATCAGCTGGTTTGCAGCGGCTGGTGTCACTCCACCGCCTTGACCATGTCCTCGATGACCTTCTTGGCGTCACCGAACACCATCATGGTCTTGTCCAGATAGAACAGTTCGTTGTCCAGGCCGGCATA  >  NC_002947/163825‑164121                                                                                                                                                      |                                                                                                                                                     aCGGTACATGGATTGTTATCTCGGGCTACTGAAACTACTGCAGCGCAGTCAAGAGATTTCGCCGATTGCCGGCCCACGACCATGGTCGTATGAGAAGCCTTGTATTTGAGTCGTTGATCCGGATCAATAAAAGAAAAACCCCGG‑CAGAtt                                                                                                                                                     >  2:1056352/1‑150 (MQ=255)        aTGGATTGTTATCTCGGGCTACTGAAACTACTGCAGCGCAGTCTAGAGATTTCGCCCATTGCCGTCCCACGACCATGGTCGTATGAGAAGCCTTGTATTTGAGTCGTTGATCCGGATCAATAAAAGAAAAGCCCCGG‑CAGATTCTGACCg                                                                                                                                              <  1:1900312/150‑1 (MQ=255)           gATTGTTATCTCGGGCTACTGAAACTACTGCAGCGCAGTCTAGAGATTTCGCCCATTGCCGTCCCACGACCATGGTCGTATGAGAAGCCTTGTATTTGAGTCGTTGATCCGGATCAATAAAAGAAAAGCCCCGG‑CAGATTCTGACCgggg                                                                                                                                           <  1:256828/150‑1 (MQ=255)                ttATCTCGGGCTACTGAAACTACTGCAGCGCAGTCTAGAGATTTCGCCCATTGCCGTCCCACGACCATGGTCGTATGAGAAGCCTTGTATTTGAGTCGTTGATCCGGATCAATAAAAGAAAAGCCCCGG‑CAGATTCTGa                                                                                                                                                 >  2:1390901/1‑139 (MQ=255)                ttATCTCGGGCTACTGAAACTACTGCAGCGCAGTCTAGAGATTTCGCCCATTGCCGTCCCACGACCATGGTCGTATGAGAAGCCTTGTATTTGAGTCGTTGATCCGGATCAATAAAAGAAAAGCCCCGG‑CAGATTCTGa                                                                                                                                                 <  1:1390901/139‑1 (MQ=255)                   tctcGGGCTACTGAAACTACTGCATCGCAGTCTAGAGATTTCGCCCATTGCCGTCCCACGACCATGGTCGTATGAGAAGCCTTGTATTTGAGTCGTTGATCCGGATCAATAAAAGAAAAGCCCCGG‑CAGATTCTGACCGGGGCTTCCttt                                                                                                                                   >  2:1785988/1‑150 (MQ=255)                    ctcGGGCTACTGAAACTACTGCAGCGCAGTCTAGAGATTTCGCCCATTGCCGTCCCACGACCATGGTCGTATGAGAAGCCTTGTATTTGAGTCGTTGATCCGGATCAATAAAAGAAAAGCCCCGGTCAGAATCTG‑CCGGGGCttttcttt                                                                                                                                  <  2:1876403/150‑7 (MQ=255)                     tcGGGCTACTGAAACTACTGCAGCGCAGTCTAGAGATTTCGCCCATTGCCGTCCCACGACCATGGTCGTATGAGAAGCCTTGTATTTGAGTCGTTGATCCGGATCAATAAAAGAAAAGCCCCGG‑CAGATTCTGACCGGGGCTTCCTTtat                                                                                                                                 <  1:406686/150‑1 (MQ=255)                           tACTGAAACTACTGCAGCGCAGTCTAGAGATTTCGCCCATTGCCGTCCCACGACCATGGTCGTATGAGAAGCCTTGTATTTGAGTCGTTGATCCGGATCAATAAAAGAAAAGCCCCGG‑CAGATTCTGACCGGGGCTTCCTTTATATCAGc                                                                                                                           <  1:468350/150‑1 (MQ=255)                           tACTGAAACTACTGCAGCGCAGTCTAGAGATTTCGCCCATTGCCGTCCCACGACCATGGTCGTATGAGAAGCCTTGTATTTGAGTCGTTGATCCGGATCAATAAAAGAAAAGCCCCGG‑CAGATTCTGACCGGGGCTTCCTTTATATCAGc                                                                                                                           >  2:178520/1‑150 (MQ=255)                             cTGAAACTACTGCAGCGCAGTCTAGAGATTTCGCCCATTGCCGTCCCACGACCATGGTCGTATGAGAAGCCTTGTATTTGAGTCGTTGCTCCGGATCAATAAAAGAAAAGCCCCGG‑CAGATTCTGACCGGGGCTTCCTTTATATCAGCTg                                                                                                                         <  2:1324106/150‑1 (MQ=255)                             cTGAAACTACTGCAGCGCAGTCTAGAGATTTCGCCCATTGCCGTCCCACGACCATGGTCGTATGAGAAGCCTTGTATTTGAGTCGTTGATCCGGATCAATAAAAGAAAAGCCCCGG‑CAGATTCTGACCGGGGCTTCCTTTATATCAGCTg                                                                                                                         >  2:1693532/1‑150 (MQ=255)                               gAAACTACTGCAGCGCAGTCTAGAGATTTCGCCCATTGCCGTCCCACGACCATGGTCGTATGAGAAGCCTTGTATTTGAGTCGTTGATCCGGATCAATAAAAGAAAAGCCCCGG‑CAGATTCTGACCGGGGCTTCCTTTATATCAGCTGGt                                                                                                                       >  2:844201/1‑150 (MQ=255)                                     actGCAGCGCAGTCTAGAGATTTCGCCCATTGCCGTCCCACGACCATGGTCGTATGAGAAGCCTTGTATTTGAGTCGTTGATCCGGATCAATAAAAGAAAAGCCCCGG‑CAGATTCTGACCGGGGCTTCCTTTATATCAGCTGGTTTGCAg                                                                                                                 <  1:1700845/150‑1 (MQ=255)                                        gCAGCGCAGTCTAGAGATTTCGCCCATTGCCGTCCCACGACCATGGTCGTATGAGAAGCCTTGTATTTGAGTCGTTGATCCGGATCAATAAAAGAAAAGCCCCCC‑CAGACTCTGACTTGTTCTTCCTTTATATCAGCTGGTTTGCAGCgg                                                                                                              <  1:832890/150‑1 (MQ=255)                                           gcgcAGTCTAGAGATTTCGCCCATTTCCGTCCCACGACCATGGTCGTATGAGAAGCCTTGTATTTGAGTCGTTGATCCGGATCAATAAAAGAAAAGCCCCGG‑CAGATTCTGACCGGGGCTTCCTTTATATCAGCTGGTTTGCAGc                                                                                                                <  1:1463221/145‑1 (MQ=255)                                           gcgcAGTCTAGAGATTTCGCCCATTTCCGTCCCACGACCATGGTCGTATGAGAAGCCTTGTATTTGAGTCGTTGATCCGGATCAATAAAAGAAAAGCCCCGG‑CAGATTCTGACCGGGGCTTCCTTTATATCAGCTGGTTTGCAGc                                                                                                                >  2:1463221/1‑145 (MQ=255)                                             gcAGTCTAGAGATTTCGCCCATTGCCGTCCCACGACCATGGTCGTATGAGAAGCCTTGTATTTGAGTCGTTGATCCGGATCAATAAAAGAAAAGCCCCGG‑CAGATTCTGACCGGGGCTTCCTTTATATCAGCTGGTTTGCAGCGGCTGgt                                                                                                         <  1:92057/150‑1 (MQ=255)                                             gcAGTCTAGAGATTTCGCCCATTGCCGTCCCACGACCATGGTCGTATGAGAAGCCTTGTATTTGAGTCGTTGATCCGGATCAATAAAAGAAAAGCCCCGG‑CAGATTCTGACCGGGGCTTCCTTTATATCAGCTGGTTTGCAGCGGCTGgt                                                                                                         <  2:63806/150‑1 (MQ=255)                                                 tCTAGAGATTTCGCCCATTGCCGTCCCACGACCATGGTCGTATGAGAAGCCTTGTATTTGAGTCGTTGATCCGGATCAATAAAAGAAAAGCCCCGG‑CAGATTCTGACCGGGGCTTCCTTTATATCAGCTGGTTTGCAGCGGCTGGTGTCa                                                                                                     >  1:1509623/1‑150 (MQ=255)                                                 tCTAGAGATTTCGCCCATTGCCGTCCCACGACCATGGTCGTATGAGAAGCCTTGTATTTGAGTCGTTGATCCGGATCAATAAAAGAAAAGCCCCGG‑CAGATTCTGACCGGGGCTTCCTTTATATCAGCTGGTTTGCAGCGGCTGGTGTCa                                                                                                     >  1:1022722/1‑150 (MQ=255)                                                   tAGAGATTTCGCCCATTGCCGTCCCACGACCATGGTCGTATGAGAAGCCTTGTATTTGAGTCGTTGATCCGGATCAATAAAAGAAAAGCCCCGGTCAGAATCTG‑CCGGGGCttttctt                                                                                                                                   <  2:1188056/118‑6 (MQ=38)                                                   tAGAGATTTCGCCCATTGCCGTCCCACGACCATGGTCGTATGAGAAGCCTTGTATTTGAGTCGTTGATCCGGATCAATAAAAGAAAAGCCCCGGTCAGAATCTG‑CCGGGGCttttctt                                                                                                                                   >  1:1188056/1‑113 (MQ=38)                                                   tAGAGATTTCGCCCATTGCCGTCCCACGACCATGGTCGTATGAGAAGCCTTGTATTTGAGTCGTTGATCCGGATCAATAAAAGAAAAGCCCCGG‑CAGATTCTGACCGGGGCTTCCTTTATATCAGCTGGTTTGCAGCGGCTGGTGTCACt                                                                                                   >  1:1320389/1‑150 (MQ=255)                                                    agagATTTCGCCCATTGCCGTCCCACGACCATGGTAGTATGAGAAGCCTTGTATTTGAGTCGTTGATCCGGATCAATAAAAGAAAAGCCCCGG‑CAGATTCTGACCGGGGCTTCCTTTATATCAGCTGGTTTGCAGCGGCTGGTGTCACTc                                                                                                  >  1:410333/1‑150 (MQ=255)                                                     gagaTTTCGCCCATTGCCGTCCCACGACCATGGTCGTATGAGAAGCCTTGTATTTGAGTCGTTGATCCGGATCAATAAAAGAAAAGCCCCGG‑CAGATTCTGACCGGGGCTTCCTTTATATCAGCTGGTTTGCAGCGGCTGGTGTCACTcc                                                                                                 >  2:153770/1‑150 (MQ=255)                                                     cccaagtCGCCCATTGCCGTCCCACGACCATGGTCGTATGAGAAGCCTTGTATTTGAGTCGTTGATCCGGATCAATAAAAGAAAAGCCCCGG‑CAGATTCTGACCGGGGCTTCCTTTATATCAGCTGGTTTGCAGCGGCTGGTGTCACTcc                                                                                                 >  2:1472644/7‑150 (MQ=255)                                                             gCCCATTGCCGTCCCACGACCATGGTCGTATGAGAAGCCTTGTATTTGAGTCGTTGATCCGGATCAATAAAAGAAAAGCCCCGG‑CAGATTCTGACCGGGGCTTCCTTTATATCAGc                                                                                                                           <  2:922913/116‑1 (MQ=255)                                                             gCCCATTGCCGTCCCACGACCATGGTCGTATGAGAAGCCTTGTATTTGAGTCGTTGATCCGGATCAATAAAAGAAAAGCCCCGG‑CAGATTCTGACCGGGGCTTCCTTTATATCAGc                                                                                                                           >  1:922913/1‑116 (MQ=255)                                                             gCCCATTGCCGTCCCACGACCATGGTAGTGTGAGAAGCCTTGTATTTGAGTCGTTGATCCGGATCAATAAAAGAAAAGCCCCGG‑CAGATTCTGACCGGGGCTTCCTTTATATCAGCTGGTTTGCAGCGGCTGGTGTCACTCCACCGCCtt                                                                                         >  1:529294/1‑150 (MQ=255)                                                               ccATTGCCGTCCCACGACCATGGTCGTATGAGAAGCCTTGTATTTGAGTCGTTGATCCGGATCAATAAAAGAAAAGCCCCGG‑CAGATTCTGACCGGGGCTTCCTTTATATCAGCTGGTTTGCAGCGGCTGGTGTCACTCCACCGCCTTGa                                                                                       >  2:314869/1‑150 (MQ=255)                                                                 aTTGCCGTCCCACGACCATGGTCGTATGAGAAGCCTTGTATTTGAGTCGTTGATCCGGATCAATAAAAGAAAAGCCCCGG‑CAGATTCTGACCGGGGCTTCCTTTATATCAGCTGGTTTGCAGCGGCTGGTGTCACTCCACCGCCTTGAcc                                                                                     <  1:1693532/150‑1 (MQ=255)                                                                     ccGTCCCACGACCATGGTCGTATGAGAAGCCTTGTATTTGAGTCGTTGATCCGGATCAATAAAAGAAAAGCCCCGG‑CAGATTCTGACCGGGGCTTCAATTATACGAGCTGGTTTTCAGCGGCTGGTGTCACTCCACCGCCTTGACCATGt                                                                                 >  1:1563271/1‑150 (MQ=255)                                                                            aCGACCATGGTCGTATGAGAAGCCTTGTATTTGAGTCGTTGATCCGGATCAATAAAAGAAAAGCCCCGG‑CAGATTCTGACCGGGGCTTCCTTTATATCAGCTGGTTTGCAGCGGCTGGTGTCACTCCACCGCCTTGACCATGTCCTCGAt                                                                          >  1:825561/1‑150 (MQ=255)                                                                              gACCATGGTCGTATGAGAAGCCTTGTATTTGAGTCGTTGATCCGGATCAATAAAAGAAAAGCCCCGG‑CAGATTCTGACCGGGGCTTCCTTTATATCAGCTGGTTTGCAGCGGCTGGTGTCACTCCACCGCCTTGACCATGTCCTCGATGa                                                                        >  1:992255/1‑150 (MQ=255)                                                                               aCCATGGTCGTATGAGAAGCCTTGTATTTGAGTCGTTGATCCGGATCAATAAAAGAAAAGCCCCGG‑CAGATTCTGACCGGGGCTTCCTTTATATCAGCTGGTTTGCAGCGGCTGGTGTCACTCCACCGCCTTGACCATGTCCTCGATGAc                                                                       <  1:314869/150‑1 (MQ=255)                                                                                ccATGGTCGTATGAGAAGCCTTGTATTTGAGTCGTTGATCCGGATCAATCCGAGAAAAGCCCCGG‑CAGATTCTGACCGGGGCTTCCTTTATATCAGCTGGTTTGCAGCGGCTGGTGTCACTCCACCGCCTTGACCATGTCCTCGATGAcc                                                                      <  1:117036/150‑1 (MQ=255)                                                                                ccATGGTCGTATGAGAAGCCTTGTATTTGAGTCGTTGATCCGGATCAATAAAAGAAAAGCCCCGG‑CAGATTCTGACCGGGGCTTCCTTTATATCAGCTGGTTTGCAGCGGCTGGTGTCACTCCACCGCCTTGACCATGTCCTCGATGAcc                                                                      >  2:428027/1‑150 (MQ=255)                                                                                     gTCGTATGAGAAGCCTTGTATTTGAGTCGTTGATCCGGATCAATAAAAGAAAAGCCCCGG‑CAGATTCTGACCGGGGCTTCCTTTATATCAGCTGGTTTGCAGCGGCTGGTGTCACTCCACCGCCTTGACCATGTCCTCGATGACcttctt                                                                 <  1:1785988/150‑1 (MQ=255)                                                                                      tCGTATGAGAAGCCTTGTATTTGAGTCGTTGATCCGGATCAATAAAAGAAAAGCCCCGG‑CAGATTCTGACCGGGGCTTCCTTTATATCAGCTGGTTTGCAGCGGCTGGTGTCACTCCACCGCCTTGACCATGTCCTCGATGACCTTCTTg                                                                <  1:428027/150‑1 (MQ=255)                                                                                      tCGTATGAGAAGCCTTGTATTTGAGTCGTTGATCCGGATCAATAAAAGAAAAGCCCCGG‑CAGATTCTGACCGGGGCTTCCTTTATATCAGCTGGTTTGCAGCGGCTGGTGTCACTCCACCGCCTTGACCATGTCCTCGATGACCTTCTTg                                                                >  1:1487814/1‑150 (MQ=255)                                                                                        gTATGAGAAGCCTTGTATTTGAGTCGTTGATCCGGATCAATAAAAGAAAAGCCCCGG‑CAGATTCTGACCGGGGCTTCCTTTATATCAGCTGGTTTGCAGCGGCTGGTGTCACTc                                                                                                  <  2:1257271/114‑1 (MQ=255)                                                                                        gTATGAGAAGCCTTGTATTTGAGTCGTTGATCCGGATCAATAAAAGAAAAGCCCCGG‑CAGATTCTGACCGGGGCTTCCTTTATATCAGCTGGTTTGCAGCGGCTGGTGTCACTc                                                                                                  >  1:1257271/1‑114 (MQ=255)                                                                                            gagaAGCCTTGTATTTGAGTCGTTGATCCGGATCAATAAAAGAAAAGCCCCAG‑CAGATTCTGACCGGGGCTTCCTTTATATCAGCTGGTTTGCAGCGGCTGGTGTCACTCCACCGCCTTGACCATGTCCTCGATGACCTTCTTGGCGTCa                                                          <  2:529294/150‑1 (MQ=255)                                                                                              gaAGCCTTGTATTTGAGTCGTTGATCCGGATCAATCACTAAAAAGCCCCGG‑CAGATTCTGACCGGGGCTTCCTTTATATCAGCTGGTTTGCAGCGGCTGGTGTCACTCCACCGCCTTGACCATGTCCTCGATGACCTTCTTGGCGTCAcc                                                        <  1:99352/150‑1 (MQ=255)                                                                                              gaAGCCTTGTATTTGAGTCGTTGATCCGGATCAATAAAAGAAAAGCCCCGG‑CAGATTCTGACCGGGGCTTCCTTTATATCAGCTGGTTTGCAGCGGCTGGTGTCACTCCACCGCCTTGACCATGTCCt                                                                              <  2:1334458/128‑1 (MQ=255)                                                                                              gaAGCCTTGTATTTGAGTCGTTGATCCGGATCAATAAAAGAAAAGCCCCGG‑CAGATTCTGACCGGGGCTTCCTTTATATCAGCTGGTTTGCAGCGGCTGGTGTCACTCCACCGCCTTGACCATGTCCt                                                                              >  1:1334458/1‑128 (MQ=255)                                                                                               aaGCCTTGTATTTGAGTCGTTGATCCGGATCAATAAAAGAAAAGCCCCGG‑CAGATTCTGACCGGGGCTTCCTTTATATCAGCTGGTTTGCAGCGGCTGGTGTCACTCCACCGCCTTGACCATGTCCTCGATGACCTTCTTGGCGTCACCg                                                       <  1:153770/150‑1 (MQ=255)                                                                                                  ccTTGTATTTGAGTCGTTGATCCGGATCAATAAAAGAAAAGCCCCGG‑CAGATTCTGACCGGGGCTTCCTTTATATCAGCTGGTTTGCAGCGGCTGGTGTCACTCCACCGCCTTGACCATGTCCTCGATGACCTTCTTGGCGTCACCGAac                                                    <  1:1518812/150‑1 (MQ=255)                                                                                                         tttGAGTCGTTGATCCGGATCAATAAAAGAAAAGCCCCGG‑CAGATTCTGACCGGGGCTTCCTTTATATCAGCTGGTTTGCAGCGGCTGGTGTCACTCCACCGCCTTGACCATGTCCTCGATGACCTTCTTg                                                                >  2:1338505/1‑131 (MQ=255)                                                                                                         tttGAGTCGTTGATCCGGATCAATAAAAGAAAAGCCCCGG‑CAGATTCTGACCGGGGCTTCCTTTATATCAGCTGGTTTGCAGCGGCTGGTGTCACTCCACCGCCTTGACCATGTCCTCGATGACCTTCTTg                                                                <  1:1338505/131‑1 (MQ=255)                                                                                                         tttGAGTCGTTGATCCGGATCAATAAAAGAAAAGCCCCGG‑CAGATTCTGACCGGGGCTTCCTTTATATCAGCTGGTTTGCAGCGGCTGGTGTCACTCCACCGCCTTGACCATGTCCTCGATGACCTTCTTGGCGTCACCGAACACcatca                                             <  1:407917/150‑1 (MQ=255)                                                                                                         tttGAGTCGTTGATCCGGATCAATAAAAGAAAAGCCCCGG‑CAGATTCTGACCGGGGCTTCCTTTATATCAGCTGGTTTGCAGCGGCTGGTGTCACTCCACCGCCTTGACCATGTCCTCGATGACCTTCTTGGCGTCACCGAACACcatca                                             >  1:1121197/1‑150 (MQ=255)                                                                                                               tCGTTGATCCGGATCAATAAAAGACAAGCCCCGG‑CAGATTCTGACCGGGGCTTCCTTTATATCAGCTGGTTTGCAGCGGCTGGTGTCACTCCACCGCCTTGACCATGTCCTCGATGACCTTCTTGGCGTCACCGAACACCATCATGGTCt                                       >  2:1284346/1‑150 (MQ=255)                                                                                                                 gTTGATCCGGATCAATAAAAGAAAAGCCCCGG‑CAGATTCTGACCGGGGCTTCCTTTATATCAGCTGGTTTGCAGCGGCTGGTGTCACTCCACCGCCTTGACCATGTCCTCGATGACCTTCTTGGCGTCACCGAACACCATCATGGTCTTg                                     <  1:46831/150‑1 (MQ=255)                                                                                                                  ttGATCCGGATCAATAAAAGAAAAGCCCCGG‑CAGATTCTGACCGGGGCTTCCTTTATATCAGCTGGTTTGCAGCGGCTGGTGTCACTCCACCGCCTTGACCATGTCCTCGATGACCTTCTTGGCGTCACCGAACACCATCATGGTCTTGt                                    <  2:483656/150‑1 (MQ=255)                                                                                                                  ttGATCCGGATCAATAAAAGAAAAGCCCCGG‑CAGATTCTGACCGGGGCTTCCTTTATATCAGCTGGTTTGCAGCGGCTGGTGTCACTCCACCGCCTTGACCATGTCCTCGATGACCTTCTTGGCGTCACCGAACACCATCATGGTCTTGt                                    <  1:844201/150‑1 (MQ=255)                                                                                                                     aTCCGGATCAATAAAAGAAAAGCCCCGG‑CAGATTCTGACCGGGGCTTCCTTTATATCAGCTGGTTTGCAGCGGCTGGTGTCACTCCACCGCCTTGACCATGTCCTCGATGACCTTCTTGGCGTCACCGAACACCATCATGGTCTTGTCCa                                 >  1:267238/1‑150 (MQ=255)                                                                                                                        cGGATCAATAAAAGAAAAGCCCCGG‑CAGATTCTGACCGGGGCTTCCTTTATATCAGCTGGTTTGCAGCGGCTGGTGTCACTCCACCGCCTTGACCATGTCCTCGATGACCTTCTTGGCGTCACCGAACACCATCATGGTCTTGTCCAGAt                              <  2:1509623/150‑1 (MQ=255)                                                                                                                        cGGATCAATAAAAGAAAAGCCCCGG‑CAGATTCTGACCGGGGCTTCCTTTATATCAGCTGGTTTGCAGCGGCTGGTGTCACTCCACCGCCTTGACCATGTCCTCGATGACCTTCTTGGCGTCACCGAACACCATCATGGTCTTGTCCAGAt                              <  2:410333/150‑1 (MQ=255)                                                                                                                         ggATCAATAGAAGAAAAGCCCCGG‑CAGATTCTGACCGGGGCTTCCTTTATATCAGCTGGTTTGCAGCGGCTGGTGTCACTCCACCGCCTTGACCATGTCCTCGATGACCTTCTTGGCGTCACCGAACACCATCATGGTCTTGTCCAGATa                             <  2:1661388/150‑1 (MQ=255)                                                                                                                             cAATAAAAGAAAAGCCCCGG‑CAGATACTGACCGGGGCTTCCTTTATATCAGCTGGTTTGCAGCGGCTGGTGTCACTCCACCGCCTTGACCATGTCCTCGATGACCTTCTTGGCGTCACCGAACACCATCATGGTCTTGTCCAGATAGAAc                         >  2:643983/1‑150 (MQ=255)                                                                                                                               aTAAAAGAAAAGCCCCGG‑CAGATTCTGACCGGGGCTTCCTTTATATCAGCTGGTTTGCAGCGGCTGGTGTCACTCCACCGCCTTGACCATGTCCTCGATGACCTTCTTGGCGTCACCGAACACCATCATGGTCTTGTCCAGATCGAACAg                       >  2:306596/1‑150 (MQ=255)                                                                                                                                tAAAAGAAAAGCCCCGG‑CAGATTCTGACCGGGGCTTCCTTTATATCAGCTGGTTTGCAGCGGCTGGTGTCACTCCACCGCCTTGACCATGTCCTCGATGACCTTCTTGGCGTCACCGAACACCATCATGGTCTTGTCCAGATAGAACAGt                      >  2:1830298/1‑150 (MQ=255)                                                                                                                                tAAAAGAAAAGCCCCGG‑CAGATTCTGACCGGGGCTTCCTTTATATCAGCTGGTTTGCAGCGGCTGGTGTCACTCCACCGCCTTGACCATGTCCTCGATGACCTTCTTGGCGTCACCGAACACCATCATGGTCTTGTCCAGATAGAACAGt                      >  2:141378/1‑150 (MQ=255)                                                                                                                                tAAAAGAAAAGCCCCGG‑CAGATTCTGACCGGGGCGTCCTTTATATCAGCTGGTTTGCAGCGGCTGGTGTCACTCCACCGCCTTGACCATGTCCTCGATGACCTTCTTGGCGTCACCGAACACCATCATGGTCTTGTCCAGATAGAACAGt                      >  2:55136/1‑150 (MQ=255)                                                                                                                                 aaaagaaaagCCCCGG‑CAGATTCTGACCGGGGCTTCCTTTATATCCGCTGGTGTGCAGCGGCTGGTGTCACTCCACCGCCTTGACCATGTCCTCGATGACCTTCTTGGCGTCACCGAACACCATCATGGTCTTGTCCAGATAGAACAGtt                     <  1:1056352/150‑1 (MQ=255)                                                                                                                                 aaaagaaaagCCCCGG‑CAGATTCTGACCGGGGCTTCCTGTATATCAGCTGGTTTGCAGCGGCTGGTGTCACTCCACCGCCTTGACCATGTCCTCGATGACCTTCTTGGCGTCACCGAACACCATCATGGTCTTGTCCAGATAGAACAGtt                     >  2:1914659/1‑150 (MQ=255)                                                                                                                                     gaaaagCCCCGG‑CAGATTCTGACCGGGGCTTCCTTTATATCAGCTGGTTTGCAGCGGCTGGTGTCACTCCACCGCCTTGACCATGTCCTCGATGACCTTCTTGGCGTCACCGAACACCATCATGGTCTTGTCCAGATAGAACAGTTCGtt                 <  1:75476/150‑1 (MQ=255)                                                                                                                                           ccccGGTCAGAATCTG‑CCGGGGCTTCCTTTATATCAGCAGGTTTGCTGCTACTGGTGTCACTCCACCGCCTTGACCATGTCCTCGATGACCTTCTTGGCGTCACCGAACACCATCATGGTCTTGTCCAGATAGAACAGTTCGTTGTCCAg           >  1:1281049/1‑150 (MQ=255)                                                                                                                                              cGG‑CAGATTCTGACCGGGGCTTCCTTTATATCAGCTGGTTTGCAGCGGCTGGTGTCACTCCACCGCCTTGACCATGTCCTCGATGACCTTCTTGGCGTCACCGAACACCATCATGGTCTTGTCCAGATAGAACAGTTCGTTGTCCAGGcc        >  2:402113/1‑150 (MQ=255)                                                                                                                                                g‑CAGATTCTGACCGGGGCTTCCTTTATATCAGCTGGTTTGCAGCGGCTGGTGTCACTCCACCGCCTTGACCATGTCCTCGATGACCTTCTTGGCGTCACCGAACACCATCATGGTCTTGTCCAGATAGAACAGTTCGTTGTCCAGGCCgg      <  2:1487814/150‑1 (MQ=255)                                                                                                                                                     aTTCTGACCGGGGCTTCCTTTATATCAGCTGGTTTGCAGCGGCTGGTGTCACTCCACCGCCTTGACCATGTCCTCGATGACCTTCTTGGCGTCACCGAACACCATCATGGTCTTGTCCAGATAGAACAGTTCGTTGTCCAGGCCGGCATa  >  1:1086524/1‑150 (MQ=255)                                                                                                                                                      |                                                                                                                                                     ACGGTACATGGATTGTTATCTCGGGCTACTGAAACTACTGCAGCGCAGTCTAGAGATTTCGCCCATTGCCGTCCCACGACCATGGTCGTATGAGAAGCCTTGTATTTGAGTCGTTGATCCGGATCAATAAAAGAAAAGCCCCGG‑CAGATTCTGACCGGGGCTTCCTTTATATCAGCTGGTTTGCAGCGGCTGGTGTCACTCCACCGCCTTGACCATGTCCTCGATGACCTTCTTGGCGTCACCGAACACCATCATGGTCTTGTCCAGATAGAACAGTTCGTTGTCCAGGCCGGCATA  >  NC_002947/163825‑164121

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