Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	163,978	Δ1 bp	5.6%	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	.	5.6%	215.6 / 5.9	71	intergenic (‑145/+42)	scpC/pntB	propionyl‑CoA:succinate CoA transferase/pyridine nucleotide transhydrogenase subunit beta
Reads supporting (aligned to +/- strand):  ref base A (34/33);  new base . (2/2);  total (36/35)
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.59e-01

ATGGATTGTTATCTCGGGCTACTGAAACTACTGCAGCGCAGTCTAGAGATTTCGCCCATTGCCGTCCCACGACCATGGTCGTATGAGAAGCCTTGTATTTGAGTCGTTGATCCGGATCAATAAAAGAAAAGCCCCGG‑CAGATTCTGACCGGGGCTTCCTTTATATCAGCTGGTTTGCAGCGGCTGGTGTCACTCCACCGCCTTGACCATGTCCTCGATGACCTTCTTGGCGTCACCGAACACCATCATGGTCTTGTCCAGATAGAACAGTTCGTTGTCCAGGCCGGCATAGCCG  >  NC_002947/163832‑164125                                                                                                                                                    |                                                                                                                                                    aTGGATTGTTATCTCGGGCTACTGAAACTACTGCAGCGCAGTCTAGAGATTTCGCCCATTGCCGTCCCACGACCATGGTCGTATGAGAAGCCTTGTATTTGAGTCGTTGATCCGGATCAATAAAAGAAAAGCCCCGG‑CAGATTCTGACCg                                                                                                                                                  <  1:1900312/150‑1 (MQ=255)    gATTGTTATCTCGGGCTACTGAAACTACTGCAGCGCAGTCTAGAGATTTCGCCCATTGCCGTCCCACGACCATGGTCGTATGAGAAGCCTTGTATTTGAGTCGTTGATCCGGATCAATAAAAGAAAAGCCCCGG‑CAGATTCTGACCgggg                                                                                                                                               <  1:256828/150‑1 (MQ=255)         ttATCTCGGGCTACTGAAACTACTGCAGCGCAGTCTAGAGATTTCGCCCATTGCCGTCCCACGACCATGGTCGTATGAGAAGCCTTGTATTTGAGTCGTTGATCCGGATCAATAAAAGAAAAGCCCCGG‑CAGATTCTGa                                                                                                                                                     <  1:1390901/139‑1 (MQ=255)         ttATCTCGGGCTACTGAAACTACTGCAGCGCAGTCTAGAGATTTCGCCCATTGCCGTCCCACGACCATGGTCGTATGAGAAGCCTTGTATTTGAGTCGTTGATCCGGATCAATAAAAGAAAAGCCCCGG‑CAGATTCTGa                                                                                                                                                     >  2:1390901/1‑139 (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                                                                                                                               >  2:178520/1‑150 (MQ=255)                    tACTGAAACTACTGCAGCGCAGTCTAGAGATTTCGCCCATTGCCGTCCCACGACCATGGTCGTATGAGAAGCCTTGTATTTGAGTCGTTGATCCGGATCAATAAAAGAAAAGCCCCGG‑CAGATTCTGACCGGGGCTTCCTTTATATCAGc                                                                                                                               <  1:468350/150‑1 (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                                                                                                                    >  2:1463221/1‑145 (MQ=255)                                    gcgcAGTCTAGAGATTTCGCCCATTTCCGTCCCACGACCATGGTCGTATGAGAAGCCTTGTATTTGAGTCGTTGATCCGGATCAATAAAAGAAAAGCCCCGG‑CAGATTCTGACCGGGGCTTCCTTTATATCAGCTGGTTTGCAGc                                                                                                                    <  1:1463221/145‑1 (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                                                                                                      >  1:1257271/1‑114 (MQ=255)                                                                                 gTATGAGAAGCCTTGTATTTGAGTCGTTGATCCGGATCAATAAAAGAAAAGCCCCGG‑CAGATTCTGACCGGGGCTTCCTTTATATCAGCTGGTTTGCAGCGGCTGGTGTCACTc                                                                                                      <  2:1257271/114‑1 (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                                        <  1:844201/150‑1 (MQ=255)                                                                                                           ttGATCCGGATCAATAAAAGAAAAGCCCCGG‑CAGATTCTGACCGGGGCTTCCTTTATATCAGCTGGTTTGCAGCGGCTGGTGTCACTCCACCGCCTTGACCATGTCCTCGATGACCTTCTTGGCGTCACCGAACACCATCATGGTCTTGt                                        <  2:483656/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:141378/1‑150 (MQ=255)                                                                                                                         tAAAAGAAAAGCCCCGG‑CAGATTCTGACCGGGGCTTCCTTTATATCAGCTGGTTTGCAGCGGCTGGTGTCACTCCACCGCCTTGACCATGTCCTCGATGACCTTCTTGGCGTCACCGAACACCATCATGGTCTTGTCCAGATAGAACAGt                          >  2:1830298/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)                                                                                                                                                  tGACCGGGGCTTCCTTTATATCAGCTGGTTTGCAGCGGCTGGTGTCACTCCACCGCCTTGACCATGTCCTCGATGACCTTCTTGGCGTCACCGAACACCATCATGGTCTTGTCCAGATAGAACAGTTCGTTGTCCAGGCCGGCCTAGCCg  >  1:392001/1‑150 (MQ=255)                                                                                                                                                    |                                                                                                                                                    ATGGATTGTTATCTCGGGCTACTGAAACTACTGCAGCGCAGTCTAGAGATTTCGCCCATTGCCGTCCCACGACCATGGTCGTATGAGAAGCCTTGTATTTGAGTCGTTGATCCGGATCAATAAAAGAAAAGCCCCGG‑CAGATTCTGACCGGGGCTTCCTTTATATCAGCTGGTTTGCAGCGGCTGGTGTCACTCCACCGCCTTGACCATGTCCTCGATGACCTTCTTGGCGTCACCGAACACCATCATGGTCTTGTCCAGATAGAACAGTTCGTTGTCCAGGCCGGCATAGCCG  >  NC_002947/163832‑164125

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